SynGap Missense Server

Table of SynGAP1 Isoform α2 (UniProt Q96PV0-1) Missense Variants.

c.dna Variant SGM Consensus Domain ClinVar gnomAD ESM1b AlphaMissense REVEL FoldX Rosetta Foldetta PremPS PROVEAN PolyPhen-2 HumDiv PolyPhen-2 HumVar FATHMM SIFT PAM Physical SASA Normalized B-factor backbone Normalized B-factor sidechain SynGAP Structural Annotation DOI
Clinical Status Review Subm. ID Allele count Allele freq. LLR score Prediction Pathogenicity Class Optimized Score Prediction Average ΔΔG Prediction StdDev ΔΔG Prediction ΔΔG Prediction ΔΔG Prediction Score Prediction pph2_prob Prediction pph2_prob Prediction Nervous System Score Prediction Prediction Status Conservation Sequences PAM250 PAM120 Hydropathy Δ MW Δ Average Δ Δ StdDev Δ StdDev Secondary Tertiary bonds Inside out GAP-Ras interface At membrane No effect MD Alert Verdict Description
c.1202G>AR401Q
(3D Viewer)
Likely PathogenicC2Uncertain 16-33438107-G-A-11.213Likely Pathogenic0.969Likely PathogenicLikely Pathogenic0.780Likely Pathogenic0.96Ambiguous0.11.50Ambiguous1.23Ambiguous1.20Destabilizing-3.69Deleterious0.999Probably Damaging0.978Probably Damaging5.47Benign0.04Affected3.3827111.0-28.06
c.600G>CL200F
(3D Viewer)
PHUncertain 16-33435242-G-C21.24e-6-7.606In-Between0.592Likely PathogenicLikely Benign0.094Likely Benign1.00Ambiguous0.51.45Ambiguous1.23Ambiguous0.43Likely Benign-1.97Neutral0.997Probably Damaging0.916Probably Damaging4.02Benign0.17Tolerated3.46920-1.034.02250.4-15.10.60.20.50.0XUncertainLeu200, a hydrophobic residue located in the N-terminal loop before the first anti-parallel β sheet strand (res. Ile205-Pro208), is replaced by another hydrophobic residue, phenylalanine. Both the phenyl group of Phe200 and the branched iso-butyl hydrocarbon sidechain of Leu200 occupy an inward hydrophobic niche (e.g., Leu246, Val222, Phe231) during the simulations. However, since the model ends abruptly at the N-terminus, no definite conclusions can be drawn from the simulations.
c.775C>TR259W
(3D Viewer)
Likely PathogenicC2Uncertain 1-12.186Likely Pathogenic0.985Likely PathogenicLikely Pathogenic0.691Likely Pathogenic1.95Ambiguous0.80.51Ambiguous1.23Ambiguous0.51Ambiguous-7.35Deleterious1.000Probably Damaging0.993Probably Damaging5.76Benign0.00Affected3.39152-33.630.03254.040.00.20.20.20.4XXXPotentially PathogenicThe guanidinium group of Arg259, located at the beginning of an anti-parallel β sheet strand (res. Arg259-Arg272), forms salt bridges with the carboxylate groups of Asp684 at the end of an α helix (res. Ile683-Gln702, GAP domain) and Asp261 on the same β strand. The Arg259 side chain also frequently forms hydrogen bonds with the backbone carbonyl groups of Ser257, Asn256, and Asp255. In the variant simulations, the indole ring of the Trp259 side chain cannot form salt bridges or maintain hydrogen bonding with the carboxylate group of Asp684 or other nearby residues. Notably, the amino group of the Lys254 side chain maintains a salt bridge with Asp684 and Glu244 throughout the variant simulations, while it forms a cation-π bond with the indole ring of Trp259 in the variant. This salt bridge is not maintained in the WT simulations. Additionally, the partially or loosely α helical conformation of a lysine-containing loop (res. Lys251-Ser257), which extends to a nearby α helix (res. Met414-Asn426), could be stabilized due to the residue swap. Moreover, the bulky size of the Trp259 side chain requires nearby residues to adjust their positioning to accommodate the introduced residue, weakening the tertiary structure assembly between the C2, PH, and GAP domains. The residue swap potentially causes more severe effects during protein folding or for the SynGAP-membrane interaction than the solvent-only simulations imply.
c.895C>TR299C
(3D Viewer)
Likely PathogenicC2Conflicting 26-33437800-C-T31.86e-6-6.326Likely Benign0.572Likely PathogenicLikely Benign0.344Likely Benign1.85Ambiguous0.40.61Ambiguous1.23Ambiguous0.76Ambiguous-3.54Deleterious1.000Probably Damaging0.998Probably Damaging1.65Pathogenic0.06Tolerated3.3919-4-37.0-53.05210.791.30.10.00.00.2XXPotentially PathogenicThe guanidinium group of Arg299, located in a β hairpin loop linking two anti-parallel β sheet strands (res. Met289-Pro298, res. Thr305-Asn315), forms hydrogen bonds that stabilize the tight turn. In the WT simulations, the Arg299 side chain hydrogen bonds with the loop backbone carbonyl groups (e.g., Ser302, Thr305, Leu274, Gly303), the hydroxyl group of Ser300, and even forms a salt bridge with the carboxylate group of Asp304.In the variant simulations, the thiol group of the Cys299 side chain is unable to form any of these well-coordinated or strong interactions, which could affect the initial formation of the secondary hairpin loop during folding. β hairpins are potential nucleation sites during the initial stages of protein folding, so even minor changes in them could be significant. Moreover, the positively charged Arg299 side chain faces the polar head group region of the inner leaflet membrane and could directly anchor the C2 domain to the membrane. In short, the residue swap could negatively affect both protein folding and the stability of the SynGAP-membrane association.
c.1606T>GL536V
(3D Viewer)
Likely PathogenicGAPUncertain 1-9.014Likely Pathogenic0.269Likely BenignLikely Benign0.586Likely Pathogenic1.25Ambiguous0.31.22Ambiguous1.24Ambiguous1.20Destabilizing-2.81Deleterious0.998Probably Damaging0.992Probably Damaging-1.34Pathogenic0.09Tolerated3.3734210.4-14.03204.726.40.20.0-0.20.2XPotentially BenignLeu536 is located on an α-helix (res. Ala533-Val560) at the membrane interface. The iso-butyl group of Leu536 interacts with nearby hydrophobic residues in the preceding loop (e.g., Val526, Pro528, Cys531). In the variant simulations, the iso-propyl side chain of Val536 forms similar hydrophobic interactions as Leu536 in the WT, causing no negative structural effects.
c.762G>CK254N
(3D Viewer)
Likely PathogenicPHUncertain 1-13.306Likely Pathogenic0.999Likely PathogenicLikely Pathogenic0.757Likely Pathogenic0.73Ambiguous0.21.87Ambiguous1.30Ambiguous1.19Destabilizing-4.23Deleterious0.384Benign0.070Benign5.93Benign0.01Affected3.3915100.4-14.07215.3-21.0-1.01.70.20.3XPotentially PathogenicThe amino group of Lys254, located in an α-β loop connecting the PH and C2 domains (res. Lys251-Arg258), forms salt bridges with the carboxylate groups of Glu244 and Asp684. Since the neutral carboxamide group of the Asn254 side chain cannot form salt bridges with acidic residues, the residue swap potentially weakens the tertiary structure assembly and/or influences the loop positioning. Regardless, in both the variant and WT simulations, all hydrogen bonds formed by the residue’s side chain were broken, and the residue rotated outwards. The partially α helical conformation of the loop, which extends to a nearby α helix (res. Met414-Asn426), is dynamic, making it unclear if the mutation affects it.
c.773G>AR258H
(3D Viewer)
C2Benign/Likely benign 36-33437678-G-A106.20e-6-10.533Likely Pathogenic0.525AmbiguousLikely Benign0.830Likely Pathogenic1.60Ambiguous0.61.00Ambiguous1.30Ambiguous1.47Destabilizing-4.06Deleterious1.000Probably Damaging0.991Probably Damaging5.77Benign0.01Affected3.3915201.3-19.05212.581.80.10.0-0.50.2XPotentially PathogenicThe guanidinium group of Arg258, located at the end of an α-β loop connecting the PH domain to the C2 domain (res. Lys251-Arg258), forms hydrogen bonds with the carboxamide groups of Asn727 and Asn729 side chains, as well as with the backbone carbonyl groups of Ala724, Leu725, and Asn727 in the WT simulations. Although the imidazole group of His258 can act as a hydrogen bond donor/acceptor, the swapped residue is unable to maintain an equally well-coordinated hydrogen bond network for linking the C2 and GAP domains in the variant simulations.
c.1771G>AA591T
(3D Viewer)
Likely PathogenicGAPConflicting 36-33440823-G-A181.12e-5-9.572Likely Pathogenic0.704Likely PathogenicLikely Benign0.270Likely Benign1.61Ambiguous0.21.00Ambiguous1.31Ambiguous1.19Destabilizing-3.40Deleterious0.955Possibly Damaging0.209Benign3.48Benign0.01Affected3.373510-2.530.03202.9-43.40.20.00.70.1XPotentially BenignThe methyl group of the Ala591 side chain, located in the middle of an α helix (res. Glu582-Met603), packs against hydrophobic residues (e.g., Ile483, Phe484) of an opposing partially helical loop (res. Phe476-Asn487).In the variant simulations, the hydroxyl group of Thr591 can form hydrogen bonds with the backbone carbonyl of Ile843 in the opposing loop or the backbone carbonyl group of Arg587. These interactions could either reinforce the tertiary assembly or weaken the α helix unity. Additionally, the Thr591 side chain can hydrogen bond with the guanidinium group of the Arg587 side chain, potentially strengthening the α helix unity.Overall, the residue swap does not seem to cause any major negative effects on the protein structure.
c.1126G>TG376CC2Uncertain 1-7.686In-Between0.125Likely BenignLikely Benign0.560Likely Pathogenic2.56Destabilizing0.50.22Likely Benign1.39Ambiguous0.16Likely Benign-1.15Neutral1.000Probably Damaging1.000Probably Damaging1.32Pathogenic0.01Affected-3-32.946.09
c.844T>AC282S
(3D Viewer)
Likely PathogenicC2Uncertain 1-11.846Likely Pathogenic0.958Likely PathogenicLikely Pathogenic0.460Likely Benign1.55Ambiguous0.11.23Ambiguous1.39Ambiguous1.62Destabilizing-9.19Deleterious0.997Probably Damaging0.994Probably Damaging1.64Pathogenic0.03Affected3.39180-1-3.3-16.06233.214.8-0.10.0-0.20.3XPotentially BenignThe thiol-containing side chain of Cys282, located at the beginning of an anti-parallel β sheet strand (res. Arg279-Leu286), packs against multiple hydrophobic residues (e.g., Ile268, Leu284, Trp308, Leu327). In the variant simulations, the hydroxyl-containing side chain of Ser282 is more hydrophilic and, hence, not as favorable as Cys282 for this hydrophobic niche. Due to this polarity difference, the residue swap could potentially weaken the hydrophobic packing of the C2 domain during the folding process.Moreover, because the C2 domain interacts with the membrane, there could also be a negative effect on the stability of the SynGAP-membrane association. However, no large-scale structural changes were observed during the variant simulations. The hydroxyl group of Ser282 forms a hydrogen bond with the backbone carbonyl group of His326 in another β strand (res. Ala322-Arg329), which competes directly with the backbone amide group of Glu283 within the secondary structure element.
c.1304T>GL435W
(3D Viewer)
Likely PathogenicGAPUncertain 1-14.889Likely Pathogenic0.992Likely PathogenicLikely Pathogenic0.572Likely Pathogenic2.11Destabilizing0.10.69Ambiguous1.40Ambiguous1.66Destabilizing-5.63Deleterious1.000Probably Damaging0.998Probably Damaging3.15Benign0.00Affected3.3729-2-2-4.773.05242.2-25.20.00.00.30.1XPotentially PathogenicThe iso-butyl side chain of Leu435, located in an α helix (res. Met414-Glu436), packs against other hydrophobic residues in an interhelix space (e.g., Val699, Val447, Leu489, Leu439) in the WT simulations. In the variant simulations, the indole ring of Trp435 fits into the same niche despite its considerably bulkier size. Additionally, the side chain forms an H-bond with the backbone carbonyl of Leu696 in an α helix (res. Asp684-Gln702). Although no apparent negative changes are observed during the variant simulation, the size difference between the swapped residues could affect the protein folding process.
c.892C>TP298S
(3D Viewer)
Likely BenignC2Benign 16-33437797-C-T53.10e-6-6.342Likely Benign0.144Likely BenignLikely Benign0.189Likely Benign1.38Ambiguous0.21.41Ambiguous1.40Ambiguous0.58Ambiguous-1.20Neutral0.991Probably Damaging0.898Possibly Damaging2.03Pathogenic0.85Tolerated3.3920-110.8-10.04
c.1465C>TL489F
(3D Viewer)
Likely PathogenicGAPUncertain 26-33438497-C-T16.20e-7-12.066Likely Pathogenic0.965Likely PathogenicLikely Pathogenic0.724Likely Pathogenic1.72Ambiguous0.51.14Ambiguous1.43Ambiguous0.56Ambiguous-3.76Deleterious1.000Probably Damaging0.997Probably Damaging-1.51Pathogenic0.01Affected3.373520-1.034.02246.4-17.80.00.00.60.1XPotentially BenignThe iso-butyl side chain of Leu489, located in the α-helix (res. Leu489-Glu519) within an inter-helix space of four helices (res. Ala461-Phe476, res. Val441-Ser457, and res. Met414-Glu436), packs with hydrophobic residues (e.g., Cys432, Ala448, Lys444, Ala493, Val447, Met468) in the inter-helix space. In the variant simulations, the phenyl ring of the Phe489 side chain can also pack favorably in the hydrophobic region. However, due to the size difference, the aromatic side chain of Phe489 tends to reposition to escape the tight region to accommodate the larger side chain, stacking with Lys444. Although no apparent negative changes are observed during the variant simulation, the size difference between the swapped residues could affect the protein folding process.
c.1789T>CF597L
(3D Viewer)
Likely PathogenicGAPUncertain 1-10.173Likely Pathogenic0.998Likely PathogenicLikely Pathogenic0.929Likely Pathogenic0.74Ambiguous0.12.12Destabilizing1.43Ambiguous1.20Destabilizing-5.97Deleterious0.999Probably Damaging0.994Probably Damaging-2.06Pathogenic0.13Tolerated201.0-34.02
c.913A>GT305A
(3D Viewer)
Likely BenignC2Conflicting 26-33437818-A-G138.05e-6-4.307Likely Benign0.078Likely BenignLikely Benign0.144Likely Benign1.30Ambiguous0.61.55Ambiguous1.43Ambiguous0.77Ambiguous-2.10Neutral0.939Possibly Damaging0.645Possibly Damaging1.76Pathogenic0.12Tolerated3.4020102.5-30.03177.943.5-0.20.10.40.0UncertainThe hydroxyl group of Thr305, located at the beginning of an anti-parallel β strand (res. Thr305-Asn315), hydrogen bonds with the carboxylate groups of Glu270 and Asp304 in the anti-parallel β strand and the adjacent β hairpin loop, respectively. In the variant simulations, the methyl group of the Ala305 side chain cannot hydrogen bond with either of the acidic residues, which could weaken the integrity of the tertiary structure and the β hairpin loop. Indeed, the guanidinium group of Arg299 does not acquire its central hairpin loop position due to the residue swap.β hairpins are potential nucleation sites during the initial stages of protein folding, so even minor changes in them could be significant. Due to its location near the membrane surface, the residue swap could also affect the C2 loop dynamics and SynGAP-membrane association. However, this is beyond the scope of the solvent-only simulations to unravel.
c.1835A>CQ612P
(3D Viewer)
Likely PathogenicGAPUncertain 1-9.684Likely Pathogenic0.673Likely PathogenicLikely Benign0.671Likely Pathogenic-0.19Likely Benign0.33.06Destabilizing1.44Ambiguous0.56Ambiguous-5.84Deleterious1.000Probably Damaging1.000Probably Damaging-1.31Pathogenic0.19Tolerated0-11.9-31.01
c.745G>AA249T
(3D Viewer)
Likely BenignPHUncertain 1-3.564Likely Benign0.805Likely PathogenicAmbiguous0.487Likely Benign1.50Ambiguous0.61.39Ambiguous1.45Ambiguous0.30Likely Benign-0.96Neutral0.990Probably Damaging0.815Possibly Damaging5.65Benign0.40Tolerated3.391510-2.530.03214.5-43.30.00.00.50.2XPotentially BenignThe methyl group of Ala249, located on the surface of an α helix (res. Ala236-Val250) facing an anti-parallel β sheet strand (res. Ile205-Val209), packs against nearby hydrophobic residues such as Leu200, Leu246, and Val250. In the variant simulations, the hydroxyl group of Thr249, which is not suitable for hydrophobic packing, forms a stable hydrogen bond with the backbone carbonyl of Asn245 in the same helix. Although this interaction could theoretically weaken the structural integrity of the α helix, this destabilizing effect is not observed in the variant simulations.
c.1832T>CM611T
(3D Viewer)
Likely BenignGAPUncertain 16-33440884-T-C16.19e-7-5.696Likely Benign0.101Likely BenignLikely Benign0.240Likely Benign1.98Ambiguous0.20.94Ambiguous1.46Ambiguous0.87Ambiguous-2.40Neutral0.034Benign0.038Benign-1.19Pathogenic0.29Tolerated3.3735-1-1-2.6-30.09
c.1717C>TR573W
(3D Viewer)
Likely PathogenicGAPConflicting 8-14.078Likely Pathogenic0.995Likely PathogenicLikely Pathogenic0.758Likely Pathogenic2.37Destabilizing0.70.57Ambiguous1.47Ambiguous0.88Ambiguous-6.94Deleterious1.000Probably Damaging0.997Probably Damaging-1.48Pathogenic0.00Affected3.37352-33.630.03257.639.00.10.00.20.0XXPotentially PathogenicThe guanidinium group of Arg573, located in an α-helix (res. Arg563-Glu578), forms a salt bridge with the carboxylate groups of Glu582 and/or Asp586 from a nearby α-helix (res. Glu582-Met603) in the WT simulations. Additionally, the Arg573 side chain stacks planarly with the aromatic phenol ring of Tyr665 and hydrogen bonds with the hydroxyl group of Ser668 from another α-helix (res. Ser641-Ser668). In the variant simulations, the indole ring of the Trp573 side chain is unable to maintain the same level of coordination as the positively charged Arg573 side chain. Indeed, Trp573 is seen hydrogen bonding only briefly with the carboxylate group of Glu582. Consequently, the integrity of the opposing α-helix end (res. Glu582-Met603) is weakened. Overall, the residue swap has the potential to substantially affect the tertiary structure assembly during the protein folding process.
c.772C>TR258C
(3D Viewer)
Likely PathogenicC2Uncertain 16-33437677-C-T16.20e-7-10.285Likely Pathogenic0.790Likely PathogenicAmbiguous0.771Likely Pathogenic1.17Ambiguous0.41.76Ambiguous1.47Ambiguous0.87Ambiguous-6.79Deleterious1.000Probably Damaging0.993Probably Damaging5.77Benign0.00Affected3.3915-3-47.0-53.05
c.1997A>GE666G
(3D Viewer)
Likely PathogenicGAPLikely Benign 16-33441256-A-G106.20e-6-12.261Likely Pathogenic0.911Likely PathogenicAmbiguous0.522Likely Pathogenic1.57Ambiguous0.11.46Ambiguous1.52Ambiguous0.93Ambiguous-6.25Deleterious1.000Probably Damaging0.970Probably Damaging3.37Benign0.02Affected3.38280-23.1-72.06173.998.50.00.0-0.70.0XPotentially PathogenicIn the WT simulations, the carboxylate group of Glu666, located on the α-helix (res. Ser641-Glu666), is involved in a highly coordinated hydrogen-bonding network between residues from two α-helices (res. Ser641-Glu666 and res. Arg563-Glu578) and from the α-α loop connecting the two α-helices (res. Ser641-Glu666 and res. Leu685-Val699), such as Lys566, Thr672, and Asn669. In the variant simulations, the carbonyl group of Gly666 occasionally forms hydrogen bonds with Lys566 and Asn669. However, Gly666 lacks a side chain and thus cannot maintain as well-coordinated a hydrogen-bond network as Glu666 in the WT, which may affect the tertiary structure assembly.
c.662A>GE221G
(3D Viewer)
Likely PathogenicPHUncertain 1-12.221Likely Pathogenic0.992Likely PathogenicLikely Pathogenic0.863Likely Pathogenic1.40Ambiguous0.11.74Ambiguous1.57Ambiguous0.71Ambiguous-5.56Deleterious0.596Possibly Damaging0.201Benign5.79Benign0.00Affected0-23.1-72.06
c.1260T>GF420L
(3D Viewer)
Likely PathogenicGAPUncertain 1-8.432Likely Pathogenic0.998Likely PathogenicLikely Pathogenic0.146Likely Benign1.76Ambiguous0.01.41Ambiguous1.59Ambiguous1.04Destabilizing-5.39Deleterious0.009Benign0.005Benign4.22Benign0.39Tolerated3.3729201.0-34.02231.113.20.00.0-0.10.0XPotentially BenignIn the WT, the phenyl ring of the Phe420 side chain, located on an α helix (res. Met414-Glu436), packs against hydrophobic residues in the interhelix area of the GAP domain (e.g., Leu689, Leu714, Leu717, Leu718). In the variant simulations, the iso-butyl side chain of Leu420 also packs into the hydrophobic inter-helix niche, but due to its smaller size, the resulting steric interactions are not as favorable as with phenylalanine. In short, the residue swap does not cause severe effects on the protein structure based on the variant simulations.
c.1552T>CY518H
(3D Viewer)
Likely PathogenicGAPUncertain 1-9.797Likely Pathogenic0.943Likely PathogenicAmbiguous0.496Likely Benign2.39Destabilizing0.40.82Ambiguous1.61Ambiguous1.31Destabilizing-4.74Deleterious1.000Probably Damaging1.000Probably Damaging3.40Benign0.08Tolerated02-1.9-26.03
c.2050G>AD684N
(3D Viewer)
Likely PathogenicGAPUncertain 1-13.155Likely Pathogenic0.985Likely PathogenicLikely Pathogenic0.382Likely Benign1.47Ambiguous0.81.76Ambiguous1.62Ambiguous0.37Likely Benign-4.99Deleterious0.999Probably Damaging0.746Possibly Damaging3.39Benign0.01Affected210.0-0.98
c.986G>AR329H
(3D Viewer)
Likely PathogenicC2Uncertain 16-33437891-G-A21.24e-6-10.154Likely Pathogenic0.769Likely PathogenicLikely Benign0.155Likely Benign2.53Destabilizing0.70.71Ambiguous1.62Ambiguous0.82Ambiguous-3.17Deleterious0.995Probably Damaging0.778Possibly Damaging4.04Benign0.05Affected3.4115201.3-19.05220.481.40.10.10.20.3UncertainThe guanidinium group of Arg329, located at the end of an anti-parallel β sheet strand (res. Ala322-Asp330), faces the negatively charged lipid bilayer surface. While the residue swap does not cause any apparent negative effects on the protein structure in the variant simulations, it could adversely affect the SynGAP-membrane association in reality. The positively charged Arg329 side chain forms hydrogen bonds with other loop residues (e.g., Ser371, Asp338) that are expected to dynamically interact with the membrane head group region. However, this phenomenon is beyond the scope of the solvent-only simulations to unravel. Notably, histidine can also be double protonated and positively charged, but this alternative protonation state was not considered in the variant simulations.
c.1531G>AG511R
(3D Viewer)
Likely PathogenicGAPLikely Pathogenic 1-11.327Likely Pathogenic0.991Likely PathogenicLikely Pathogenic0.416Likely Benign1.94Ambiguous0.31.32Ambiguous1.63Ambiguous0.94Ambiguous-7.72Deleterious1.000Probably Damaging1.000Probably Damaging3.26Benign0.06Tolerated3.3735-3-2-4.199.14279.4-159.90.00.00.70.1XXPotentially PathogenicGly511 is located in an α-helix (res. Gly502-Tyr518), facing hydrophobic residues in an inter-helix space (e.g., Leu610, Ile514) in the WT simulations. In contrast, in the variant simulations, the bulkier and positively charged guanidinium side chain of Arg511 forms a salt bridge with the carboxylate group of Glu217 or hydrogen bonds with the backbone carbonyl group of Leu610. Although the residue swap introduces a third positively charged residue in close vicinity (Arg511, Lys507, Arg515), the protein structure seems to remain stable in the variant simulations. Importantly, according to ClinVar, the residue swap alters the last nucleotide of an exon and is predicted to destroy the splice donor site, resulting in aberrant splicing and pathogenic status.10.1016/j.ajhg.2020.11.011
c.1531G>CG511R
(3D Viewer)
Likely PathogenicGAPPathogenic 1-11.327Likely Pathogenic0.991Likely PathogenicLikely Pathogenic0.415Likely Benign1.94Ambiguous0.31.32Ambiguous1.63Ambiguous0.94Ambiguous-7.72Deleterious1.000Probably Damaging1.000Probably Damaging3.26Benign0.06Tolerated3.3735-3-2-4.199.14279.4-159.90.00.00.70.1XXPotentially PathogenicGly511 is located in an α-helix (res. Gly502-Tyr518), facing hydrophobic residues in an inter-helix space (e.g., Leu610, Ile514) in the WT simulations. In contrast, in the variant simulations, the bulkier and positively charged guanidinium side chain of Arg511 forms a salt bridge with the carboxylate group of Glu217 or hydrogen bonds with the backbone carbonyl group of Leu610. Although the residue swap introduces a third positively charged residue in close vicinity (Arg511, Lys507, Arg515), the protein structure seems to remain stable in the variant simulations. Importantly, according to ClinVar, the residue swap alters the last nucleotide of an exon and is predicted to destroy the splice donor site, resulting in aberrant splicing and pathogenic status.10.1016/j.ajhg.2020.11.011
c.1484A>GE495G
(3D Viewer)
Likely PathogenicGAPUncertain 16-33438516-A-G16.20e-7-9.400Likely Pathogenic0.923Likely PathogenicAmbiguous0.867Likely Pathogenic1.21Ambiguous0.02.06Destabilizing1.64Ambiguous0.78Ambiguous-6.70Deleterious1.000Probably Damaging0.999Probably Damaging-1.46Pathogenic0.02Affected3.3735-203.1-72.06
c.1256A>GE419G
(3D Viewer)
Likely PathogenicGAPUncertain 1-10.589Likely Pathogenic0.956Likely PathogenicLikely Pathogenic0.469Likely Benign1.41Ambiguous0.01.94Ambiguous1.68Ambiguous0.83Ambiguous-6.42Deleterious1.000Probably Damaging0.997Probably Damaging3.31Benign0.02Affected3.37290-23.1-72.06165.3110.80.00.0-0.10.0XPotentially PathogenicThe carboxylate group of Glu419, located on an α helix (res. Met414-Glu436), forms a salt bridge with the side chain of either Arg716 or Lys418 from an opposing helix (res. Pro713-Arg726). The backbone amide group of Glu419 does not form H-bonds, resulting in a slight bend in the α helix. Thus, although glycine is known as an “α helix breaker,” the residue swap does not disrupt the continuity or integrity of the α helix. However, because Gly419 cannot form a salt bridge with the guanidinium group of the Arg716 side chain, the C2-GAP domain tertiary structure could be compromised during folding.
c.1306G>AE436K
(3D Viewer)
Likely PathogenicGAPUncertain 1-13.869Likely Pathogenic0.997Likely PathogenicLikely Pathogenic0.829Likely Pathogenic0.56Ambiguous0.12.86Destabilizing1.71Ambiguous0.82Ambiguous-3.77Deleterious0.994Probably Damaging0.951Probably Damaging4.71Benign0.02Affected3.372901-0.4-0.94186.839.80.00.0-0.20.0XXXPotentially PathogenicThe carboxylate group of Glu436, located on the α helix (res. Met414-Glu436), forms a salt bridge with the amino group of the Lys444 side chain on an opposing α helix (res. Val441-Ser457). The backbone carbonyl of Glu436 also H-bonds with the Lys444 side chain, which helps keep the ends of the two α helices tightly connected. In contrast, in the variant simulations, the salt bridge formation with Lys444 is not possible. Instead, the repelled Lys436 side chain rotates outward, causing a change in the α helix backbone H-bonding: the amide group of Lys444 H-bonds with the carbonyl of Ala433 instead of the carbonyl of Cys432.
c.1474A>GK492E
(3D Viewer)
Likely PathogenicGAPConflicting 2-16.175Likely Pathogenic0.998Likely PathogenicLikely Pathogenic0.510Likely Pathogenic1.53Ambiguous0.11.90Ambiguous1.72Ambiguous1.42Destabilizing-3.98Deleterious1.000Probably Damaging0.998Probably Damaging2.99Benign0.01Affected3.3735100.40.94
c.1738G>AG580S
(3D Viewer)
Likely PathogenicGAPUncertain 16-33440790-G-A16.20e-7-10.788Likely Pathogenic0.861Likely PathogenicAmbiguous0.644Likely Pathogenic2.84Destabilizing0.20.59Ambiguous1.72Ambiguous0.87Ambiguous-5.73Deleterious1.000Probably Damaging0.999Probably Damaging-1.23Pathogenic0.07Tolerated3.373410-0.430.03233.9-49.30.80.00.60.1XPotentially BenignGly580 is located on the outer surface in a short α-α loop turn connecting two α-helices (res. Arg563-Glu578, res. Glu582-Phe608) in the WT simulations. In the variant simulations, the side chain of Ser580 faces outward, and its hydroxyl group does not make any new or additional interactions compared to Gly580 in the WT simulations that could affect the protein structure.
c.1787G>AR596H
(3D Viewer)
Likely PathogenicGAPLikely Benign 16-33440839-G-A159.29e-6-11.128Likely Pathogenic0.950Likely PathogenicAmbiguous0.717Likely Pathogenic3.00Destabilizing0.90.43Likely Benign1.72Ambiguous1.35Destabilizing-4.97Deleterious1.000Probably Damaging0.999Probably Damaging2.43Pathogenic0.00Affected3.3735201.3-19.05223.580.5-0.10.0-0.10.3XXPotentially PathogenicThe guanidinium group of Arg596, located in an α helix (res. Glu582-Met603), forms a salt bridge with the carboxylate group of Glu495 from another α helix (res. Leu489-Glu519). In the WT simulations, the side chain of Arg596 hydrogen bonds with the backbone carbonyl groups of Asn487, Glu486, Arg485, and Phe484. Additionally, Arg596 can hydrogen bond with the carboxamide group of the Asn487 side chain on an opposing loop that links two α helices (res. Ala461-Arg475, res. Leu489-Glu519).In the variant simulations, the imidazole ring of His596 can form hydrogen bonds with the same residues as arginine; however, these interactions are not as coordinated or strong in comparison. Thus, the residue swap could affect the tertiary structure assembly more profoundly than observed in the simulations. Notably, Arg596 plays a key role in positioning the aforementioned loop, which is crucial for the placement of the “arginine finger” or the Arg485 side chain during RasGTPase activation.
c.1792C>GL598V
(3D Viewer)
Likely PathogenicGAPUncertain 1-10.002Likely Pathogenic0.578Likely PathogenicLikely Benign0.221Likely Benign1.89Ambiguous0.11.58Ambiguous1.74Ambiguous1.01Destabilizing-2.92Deleterious0.944Possibly Damaging0.786Possibly Damaging3.21Benign0.02Affected3.3735210.4-14.03218.429.60.00.00.80.0XPotentially BenignThe iso-butyl side chain of Leu598, located on an α helix (res. Glu582-Met603), packs hydrophobically with other hydrophobic residues in the inter-helix space (e.g., Ile602, Phe594, Ile510).In the variant simulations, Val598, which has similar size and physicochemical properties to leucine, resides in the inter-helix hydrophobic space in a similar manner to Leu598 in the WT. This causes no negative effects on the protein structure.
c.835C>TR279W
(3D Viewer)
Likely PathogenicC2Uncertain 1-11.417Likely Pathogenic0.942Likely PathogenicAmbiguous0.485Likely Benign2.00Destabilizing0.81.47Ambiguous1.74Ambiguous0.80Ambiguous-6.29Deleterious1.000Probably Damaging0.998Probably Damaging1.88Pathogenic0.00Affected3.39182-33.630.03270.038.30.10.00.30.0UncertainThe guanidinium group of Arg279, located at the beginning of an anti-parallel β sheet strand (res. Arg279-Leu286), can form hydrogen bond with the backbone carbonyl groups of nearby loop residues (e.g., Ser296, Ser331, and As332) and form salt bridges with the carboxylate groups of Asp330 and Asp332. In the WT simulations, Arg279 sporadically forms a salt bridge even with the carboxylate group of Glu613, loosely connecting the C2 domain and GAP domain. Meanwhile, the indole ring of the Trp279 side chain is unable to hydrogen bond with the loop residues in the variant simulations. The lack of hydrogen bond or salt bridge formation with the loop residues could be significant, as Arg279 and the loops face the polar head group region of the membrane. Thus, although Trp279 could interact with the membrane surface as a “lipid anchor,” any changes to the wider loop dynamics could still adversely affect the formation of a stable SynGAP-membrane association. However, no definite conclusions on the effect of the residue swap on the SynGAP-membrane association can be drawn from solvent-only simulations.
c.1873C>GL625VLikely PathogenicGAPUncertain 1-11.319Likely Pathogenic0.833Likely PathogenicAmbiguous0.480Likely Benign1.80Ambiguous0.71.69Ambiguous1.75Ambiguous1.42Destabilizing-2.96Deleterious0.998Probably Damaging0.992Probably Damaging3.07Benign0.01Affected210.4-14.03
c.1240A>GM414VGAPUncertain 1-8.003Likely Pathogenic0.541AmbiguousLikely Benign0.261Likely Benign1.81Ambiguous0.41.73Ambiguous1.77Ambiguous0.95Ambiguous-2.95Deleterious0.999Probably Damaging0.987Probably Damaging3.43Benign0.24Tolerated212.3-32.06
c.1150G>AG384S
(3D Viewer)
Likely BenignC2Uncertain 16-33438055-G-A16.22e-7-5.243Likely Benign0.090Likely BenignLikely Benign0.315Likely Benign1.92Ambiguous0.21.66Ambiguous1.79Ambiguous0.19Likely Benign-0.67Neutral0.980Probably Damaging0.968Probably Damaging1.33Pathogenic0.04Affected4.32210-0.430.03202.4-49.80.51.0-0.20.0UncertainGly384 is located in the Gly-rich Ω loop (res. Pro364-Pro398) between two anti-parallel β sheet strands (res. Thr359-Pro364, res. Ala399-Ile411). Because the Ω loop is assumed to directly interact with the membrane, it moves arbitrarily throughout the WT solvent simulations. The Ω loop potentially plays a crucial role in the SynGAP-membrane complex association, stability, and dynamics. However, this aspect cannot be fully addressed through solvent simulations alone.Ω loops are known to play major roles in protein functions that require flexibility, and so they are rich in glycines, prolines, and, to a lesser extent, small hydrophilic residues to ensure maximum flexibility. Thus, the variant’s Ser384 is potentially tolerated in the Ω loop, although the hydroxyl group of Ser384 forms various hydrogen bonds with several other loop residues in the variant simulations. However, since the effects on Gly-rich Ω loop dynamics can only be studied through the SynGAP-membrane complex, no definite conclusions can be drawn.
c.1487A>GE496G
(3D Viewer)
Likely PathogenicGAPUncertain 1-13.529Likely Pathogenic0.850Likely PathogenicAmbiguous0.825Likely Pathogenic1.83Ambiguous0.11.76Ambiguous1.80Ambiguous0.92Ambiguous-6.16Deleterious1.000Probably Damaging0.999Probably Damaging-1.45Pathogenic0.02Affected3.37350-23.1-72.06173.9103.10.00.0-0.70.0XXPotentially PathogenicGlu496 is located in the α-helix (res. Leu489-Glu519), and its carboxylate group forms salt bridges with the neighbouring residues Lys492 and Arg499 in the WT simulations. Glu496 also forms a hydrogen bond with Ser449 on an opposing helix (res. Val441-Ser457). In the variant simulations, Gly496 cannot form these salt bridges, which could weaken the secondary structure. Additionally, the loss of the hydrogen bond with Ser449 on the opposite helix can weaken the tertiary structure assembly. Moreover, glycine is an α-helix breaker, and it is seen to weaken the integrity of the helix as the hydrogen bonding between the backbone atoms of Gly496 and Ala493 breaks down. Also, due to its location at the GAP-Ras interface, the interaction of Glu496 with Arg499 and Lys492 might play a role in complex association and stability, which cannot be fully addressed using the SynGAP solvent-only simulations.
c.1685C>TP562L
(3D Viewer)
Likely PathogenicGAPPathogenic/Likely path. 106-33440737-C-T-13.438Likely Pathogenic0.996Likely PathogenicLikely Pathogenic0.829Likely Pathogenic3.54Destabilizing0.80.17Likely Benign1.86Ambiguous-0.14Likely Benign-9.95Deleterious1.000Probably Damaging1.000Probably Damaging0.58Pathogenic0.00Affected3.3735-3-35.416.04228.8-68.5-0.10.00.10.2XPotentially PathogenicPro562 is located on an α-α loop between two α-helices (res. Ala533-Val560 and res. Arg563-Glu578). The cyclic pyrrolidine side chain of Pro562 hydrophobically packs with other residues in the inter-helix space, such as Leu565, Ile501, and Phe561. In the variant simulations, Leu562 packs more favorably with the nearby hydrophobic residues, and the backbone amide group of Leu562 (absent in proline) does not form any intra-protein hydrogen bonds. However, prolines are well-suited for unstructured regions like loops, and thus, Pro562 in the WT is necessary at the end of the helix to induce a tight turn during folding. Although no negative structural effects are observed during the simulations, the residue swap could potentially cause extensive damage to the protein structure during folding.10.1016/j.ajhg.2020.11.011
c.1136C>TS379L
(3D Viewer)
Likely BenignC2Benign 16-33438041-C-T84.05e-5-5.641Likely Benign0.173Likely BenignLikely Benign0.469Likely Benign0.39Likely Benign0.23.38Destabilizing1.89Ambiguous-0.52Ambiguous-0.85Neutral0.015Benign0.002Benign3.83Benign0.04Affected4.3211-3-24.626.08251.9-48.10.61.10.00.5UncertainSer379 is located in the Gly-rich Ω loop (res. Pro364-Pro398) between two anti-parallel β sheet strands (res. Thr359-Pro364, res. Ala399-Ile411). Because the Ω loop is assumed to directly interact with the membrane, it moves arbitrarily throughout the WT solvent simulations. The Ω loop potentially plays a crucial role in the SynGAP-membrane complex association, stability, and dynamics. However, this aspect cannot be fully addressed through solvent simulations alone.Ω loops are known to play major roles in protein functions that require flexibility, and thus hydrophobic residues like leucine are rarely tolerated. Although no negative structural effects are observed in the variant simulations, Leu379 may exert drastic effects on the SynGAP-membrane complex dynamics and stability. However, since the effect on Gly-rich Ω loop dynamics can only be studied through the SynGAP-membrane complex, no definite conclusions can be drawn.
c.1108G>AG370S
(3D Viewer)
Likely BenignC2Uncertain 16-33438013-G-A159.31e-6-3.533Likely Benign0.081Likely BenignLikely Benign0.282Likely Benign2.83Destabilizing2.01.05Ambiguous1.94Ambiguous-0.02Likely Benign0.47Neutral0.000Benign0.000Benign1.33Pathogenic0.77Tolerated3.421910-0.430.03196.6-49.60.92.2-0.10.4UncertainGly370 is located in the Gly-rich Ω loop (res. Pro364- Pro398) between two anti-parallel β sheet strands (res. Thr359-Pro364, res. Ala399-Ile411). Because, the Ω loop is assumed to be directly interacting with the membrane, it is only seen to move arbitrarily throughout the WT solvent simulations. The Ω loop is potentially playing a crucial loop in the SynGAP-membrane complex association, stability and dynamics, regardless, this aspect cannot be addressed through the solvent simulations only. The Ω-loops are known to have a major role in protein functions that requires flexibility and thus, they are rich in glycines, prolines and to a lesser extent, hydrophilic residues to ensure maximum flexibility. Thus, Ser370 in the variant is potentially tolerated in the Ω loop. However, since the effect on the Gly-rich Ω loop dynamics can only be well-studied through the SynGAP-membrane complex, no definite conclusions can be withdrawn.
c.1405G>AA469T
(3D Viewer)
Likely PathogenicGAPUncertain 1-9.540Likely Pathogenic0.723Likely PathogenicLikely Benign0.527Likely Pathogenic2.26Destabilizing0.11.90Ambiguous2.08Destabilizing0.34Likely Benign-1.46Neutral0.994Probably Damaging0.986Probably Damaging-1.21Pathogenic0.42Tolerated10-2.530.03
c.1045C>TP349S
(3D Viewer)
C2Uncertain 1-7.654In-Between0.217Likely BenignLikely Benign0.277Likely Benign1.92Ambiguous0.12.28Destabilizing2.10Destabilizing0.87Ambiguous-6.13Deleterious1.000Probably Damaging0.996Probably Damaging1.66Pathogenic0.06Tolerated3.37251-10.8-10.04194.9-18.1-0.10.00.20.1XXPotentially PathogenicThe cyclic pyrrolidine side chain of Pro349, located at the end of an anti-parallel β sheet strand (res. Gly341-Pro349), allows the strand to end and make a tight turn before a short α helical section within a loop connecting to another β strand (res. Thr359-Pro364). In the variant simulations, the hydroxyl group of Ser349 forms a hydrogen bond with the backbone amide group of Ala351 in the short helical section. Conversely, the backbone amide group of Ser349 (absent in proline) does not form any intra-protein hydrogen bonds. However, the β strand end connects to the α helical section in a more stable and consistent manner compared to the WT. Although the residue swap does not cause major adverse effects on the protein structure in the simulations, it is possible that the tight turn at the β strand end could not be created during folding without the presence of proline.
c.2086C>GL696V
(3D Viewer)
Likely PathogenicGAPUncertain 1-11.909Likely Pathogenic0.745Likely PathogenicLikely Benign0.351Likely Benign2.35Destabilizing0.11.85Ambiguous2.10Destabilizing1.46Destabilizing-2.79Deleterious0.992Probably Damaging0.970Probably Damaging3.16Benign0.00Affected3.4613120.4-14.03
c.1502T>CI501T
(3D Viewer)
Likely BenignGAPUncertain 1-5.996Likely Benign0.252Likely BenignLikely Benign0.362Likely Benign2.40Destabilizing0.11.81Ambiguous2.11Destabilizing1.57Destabilizing-3.48Deleterious1.000Probably Damaging1.000Probably Damaging3.44Benign0.16Tolerated3.37350-1-5.2-12.05214.526.90.00.00.50.0XPotentially PathogenicIle501 is located near a hinge in the middle of an α-helix (res. Leu489-Glu519). The sec-butyl side chain of Ile501 is hydrophobically packed with other residues in the inter-helix space (e.g., Leu500, Tyr497, Phe679) in the WT simulations. In the variant simulations, the hydroxyl group of Thr501 forms a hydrogen bond with the backbone atoms of Tyr497 on the same α-helix, which may weaken the α-helix integrity. Additionally, the polar hydroxyl group of Thr501 is not suitable for the hydrophobic inter-helix space, and thus, the residue swap could affect protein folding. However, Ile501 is followed by Gly502, which facilitates a hinge in the middle of the α-helix, making further weakening caused by Thr501 unlikely to be harmful to the α-helix integrity.
c.1718G>AR573Q
(3D Viewer)
Likely PathogenicGAPLikely Pathogenic 1-9.900Likely Pathogenic0.923Likely PathogenicAmbiguous0.733Likely Pathogenic2.28Destabilizing0.81.94Ambiguous2.11Destabilizing1.08Destabilizing-3.16Deleterious1.000Probably Damaging0.995Probably Damaging-1.31Pathogenic0.12Tolerated3.3735111.0-28.06230.149.90.00.0-0.60.0XXPotentially PathogenicThe guanidinium group of Arg573, located in an α-helix (res. Arg563-Glu578), forms a salt bridge with the carboxylate groups of Glu582 and/or Asp586 from a nearby α-helix (res. Glu582-Met603) in the WT simulations. Additionally, the Arg573 side chain stacks planarly with the aromatic phenol ring of Tyr665 and hydrogen bonds with the hydroxyl group of Ser668 from another α-helix (res. Ser641-Ser668). In the variant simulations, although the carboxamide group of the Gln573 side chain can hydrogen bond with the carboxylate group of Glu582 or the hydroxyl group of Ser668, these interactions are not as coordinated, stable, or strong as those of the positively charged Arg573. Consequently, the integrity of the opposing α-helix end (res. Glu582-Met603) is weakened. Overall, the residue swap has the potential to substantially affect the tertiary structure assembly during the protein folding process.
c.930G>CE310D
(3D Viewer)
Likely PathogenicC2Likely Pathogenic1-11.218Likely Pathogenic0.994Likely PathogenicLikely Pathogenic0.666Likely Pathogenic1.87Ambiguous0.52.39Destabilizing2.13Destabilizing1.04Destabilizing-2.76Deleterious0.997Probably Damaging0.992Probably Damaging1.19Pathogenic0.02Affected3.3819320.0-14.03232.627.20.10.00.10.1XPotentially BenignThe carboxylate group of Glu310, located in an anti-parallel β sheet strand (res. Thr305-Asn315), is ideally positioned to interact with the hydroxyl and backbone amide groups of Thr295 on a twisted anti-parallel β strand. Because the carboxylate group can also interact with the GAP domain residues (e.g., Gln612, Tyr614), Glu310 potentially plays a key role in maintaining the tertiary assembly between the C2 and GAP domains. In the variant simulations, the carboxylate group of Asp310 can form the same interactions as glutamate; however, due to its one hydrocarbon shorter length, the connections are less stable or less optimal.
c.821T>AL274Q
(3D Viewer)
Likely PathogenicC2Uncertain 1-15.518Likely Pathogenic0.995Likely PathogenicLikely Pathogenic0.774Likely Pathogenic2.54Destabilizing0.31.74Ambiguous2.14Destabilizing1.97Destabilizing-5.42Deleterious1.000Probably Damaging0.999Probably Damaging0.00Pathogenic0.00Affected3.3819-2-2-7.314.97245.91.80.00.00.10.2XXXPotentially PathogenicThe aliphatic side chain of Leu274, located in a β hairpin loop (res. Glu273-Lys278) connecting two anti-parallel β sheet strands, packs against multiple hydrophobic residues facing the β sheet (e.g., Ala271, Leu327, Tyr280, Val306). The hydrophilic carboxamide group of the Gln274 side chain is not suitable for this hydrophobic niche, causing nearby residues to adjust to make room for the hydrophilic glutamine. Additionally, a new hydrogen bond forms with the backbone carboxyl group of Arg272 in another β strand (res. Glu273-Arg259).As a result, the backbone amide group of Ala399 and the carbonyl group of Arg272, which connect two β strands at the β sheet end, form fewer hydrogen bonds in the variant than in the WT simulations. Although no major secondary structure disruption is observed in the variant simulations, the residue swap could profoundly affect the C2 domain folding, as the hydrophobic packing of Leu274 is crucial for maintaining the loop's contact with the rest of the C2 domain. Lastly, because the Leu274-containing loop faces the membrane surface, the residue swap could also negatively impact the SynGAP-membrane association.
c.2143C>TP715S
(3D Viewer)
GAPLikely Pathogenic 16-33441608-C-T16.20e-7-7.635In-Between0.787Likely PathogenicAmbiguous0.277Likely Benign3.54Destabilizing0.00.81Ambiguous2.18Destabilizing0.94Ambiguous-7.17Deleterious1.000Probably Damaging0.998Probably Damaging3.43Benign0.01Affected3.5091-10.8-10.04231.8-14.0-0.10.0-0.80.1XUncertainPro715, along with Gly712 and Pro713, are located in a hinge region of an α-helix making a ~90-degree turn (res. Lys705-Leu725). In the WT simulations, the pyrrolidine side chain of Pro715, lacking the backbone amide groups altogether, forces the tight helix turn to take place while also hydrophobically packing with nearby residues (e.g., Leu700, Leu708, Leu714, and Leu718). Leu715, with a normal amide backbone, could potentially affect protein folding and turn formation, although this was not observed in the variant simulations. Additionally, the hydroxyl group of the Ser715 side chain can form hydrogen bonds with the backbone carbonyl group of Gly712 and disrupt the hydrophobic packing arrangement of the leucine residues from the neighboring α-helices, impacting the GAP domain tertiary assembly.
c.1947G>CM649I
(3D Viewer)
Likely PathogenicGAPUncertain 1-9.361Likely Pathogenic0.995Likely PathogenicLikely Pathogenic0.449Likely Benign2.42Destabilizing0.21.96Ambiguous2.19Destabilizing1.01Destabilizing-3.99Deleterious0.672Possibly Damaging0.093Benign3.40Benign0.02Affected3.3827212.6-18.03243.721.50.00.10.00.1XPotentially BenignThe thioether side chain of Met649, located on an α helix (res. Ser641-Glu666), bridges Phe652, Phe648, and Phe639 in an inter-helix hydrophobic cavity in the WT simulations. In the variant simulations, the sec-butyl side chain of Ile649 maintains hydrophobic interactions with nearby residues, with no significant effects on the protein structure.However, methionine is known as a bridging motif for aromatic residues, and these Met-aromatic interactions are lost in the variant. Indeed, in the second variant simulation,the bridging of Phe652, Phe648 and Phe639 is completely lost. In reality, the effect could be more severe on the structure during the protein folding.
c.1404G>AM468I
(3D Viewer)
Likely PathogenicGAPUncertain 16-33438436-G-A16.20e-7-8.583Likely Pathogenic0.907Likely PathogenicAmbiguous0.508Likely Pathogenic2.53Destabilizing0.21.89Ambiguous2.21Destabilizing0.37Likely Benign-1.06Neutral0.748Possibly Damaging0.886Possibly Damaging-1.10Pathogenic0.07Tolerated3.3731122.6-18.03
c.1971G>CW657CLikely PathogenicGAPUncertain 1-12.035Likely Pathogenic0.997Likely PathogenicLikely Pathogenic0.463Likely Benign2.74Destabilizing0.31.69Ambiguous2.22Destabilizing1.30Destabilizing-11.06Deleterious1.000Probably Damaging0.982Probably Damaging3.43Benign0.03Affected-8-23.4-83.07
c.1786C>TR596C
(3D Viewer)
Likely PathogenicGAPConflicting 26-33440838-C-T63.72e-6-10.805Likely Pathogenic0.972Likely PathogenicLikely Pathogenic0.633Likely Pathogenic2.94Destabilizing0.01.49Ambiguous2.22Destabilizing-0.03Likely Benign-7.96Deleterious1.000Probably Damaging1.000Probably Damaging2.41Pathogenic0.00Affected3.3735-4-37.0-53.05230.797.9-0.10.0-0.30.4XXPotentially PathogenicThe guanidinium group of Arg596, located in an α helix (res. Glu582-Met603), forms a salt bridge with the carboxylate group of Glu495 from another α helix (res. Leu489-Glu519). In the WT simulations, the side chain of Arg596 hydrogen bonds with the backbone carbonyl groups of Asn487, Glu486, Arg485, and Phe484. Additionally, Arg596 can hydrogen bond with the carboxamide group of the Asn487 side chain on an opposing loop that links two α helices (res. Ala461-Arg475, res. Leu489-Glu519).In the variant simulations, the thiol group of the Cys596 side chain is unable to form salt bridges or any of the hydrogen bonds that the Arg596 side chain can. Thus, the residue swap could affect the tertiary structure assembly more profoundly than observed in the simulations. Notably, Arg596 plays a key role in positioning the aforementioned loop, which is crucial for the placement of the “arginine finger” or the Arg485 side chain during RasGTPase activation.
c.1157G>AG386E
(3D Viewer)
C2Uncertain 16-33438062-G-A-9.286Likely Pathogenic0.686Likely PathogenicLikely Benign0.447Likely Benign3.69Destabilizing2.90.79Ambiguous2.24Destabilizing0.54Ambiguous-0.83Neutral0.860Possibly Damaging0.354Benign3.93Benign0.01Affected4.323-20-3.172.06
c.719A>GD240GLikely PathogenicPHUncertain 1-12.825Likely Pathogenic0.951Likely PathogenicAmbiguous0.912Likely Pathogenic1.85Ambiguous0.12.72Destabilizing2.29Destabilizing0.24Likely Benign-6.19Deleterious0.993Probably Damaging0.984Probably Damaging5.79Benign0.01Affected1-13.1-58.04
c.1408A>GM470V
(3D Viewer)
Likely PathogenicGAPUncertain 1-8.856Likely Pathogenic0.478AmbiguousLikely Benign0.770Likely Pathogenic2.73Destabilizing0.11.88Ambiguous2.31Destabilizing1.31Destabilizing-3.58Deleterious0.999Probably Damaging0.993Probably Damaging-1.20Pathogenic0.15Tolerated3.3734122.3-32.06
c.1214G>AR405H
(3D Viewer)
Likely PathogenicC2Conflicting 26-33438119-G-A42.48e-6-9.081Likely Pathogenic0.706Likely PathogenicLikely Benign0.371Likely Benign2.79Destabilizing0.61.85Ambiguous2.32Destabilizing1.26Destabilizing-4.54Deleterious1.000Probably Damaging0.991Probably Damaging3.65Benign0.01Affected3.3828201.3-19.05214.0102.2-0.10.0-0.70.1XPotentially PathogenicThe guanidinium group of Arg405, located in an anti-parallel β sheet strand of the C2 domain (res. Pro398-Ile411), forms a salt bridge with the carboxylate group of the Glu446 side chain from an opposing α helix (res. Val441-Ser457) in the GAP domain. The positively charged Arg405 side chain also stacks with the aromatic ring of the Phe358 side chain from a loop preceding the β strand (res. Thr359-Thr366), which could assist in maintaining the anti-parallel strand arrangement.In the variant simulations, the imidazole ring of His405 does not stack with the aromatic ring of Phe358 nor form any lasting H-bonds with the loop residues. The imidazole ring of His405 (neutral and epsilon protonated in the simulations) is unable to form a salt bridge with Glu446, which could affect the tertiary structure assembly, although this is not apparent based on the variant simulations.
c.1118G>AG373E
(3D Viewer)
C2Uncertain 1-7.281In-Between0.569Likely PathogenicLikely Benign0.420Likely Benign4.13Destabilizing3.20.52Ambiguous2.33Destabilizing-0.02Likely Benign-0.69Neutral0.001Benign0.000Benign3.90Benign0.01Affected0-2-3.172.06
c.844T>CC282R
(3D Viewer)
Likely PathogenicC2Pathogenic 2-16.378Likely Pathogenic0.999Likely PathogenicLikely Pathogenic0.466Likely Benign3.13Destabilizing0.61.58Ambiguous2.36Destabilizing1.70Destabilizing-11.03Deleterious0.999Probably Damaging0.998Probably Damaging1.63Pathogenic0.00Affected3.3918-4-3-7.053.05297.4-98.2-0.10.10.50.0XXXPotentially PathogenicThe thiol-containing side chain of Cys282, located at the beginning of an anti-parallel β sheet strand (res. Arg279-Leu286), is packed against multiple hydrophobic residues (e.g., Ile268, Leu284, Trp308, Leu327). In the variant simulations, the bulky side chain of Arg282 with its positively charged guanidinium group is not suitable for this hydrophobic niche. Consequently, the hydrophobic residues must either make room to accommodate Arg282 or it must escape the hydrophobic C2 domain core.As a result, new hydrogen bonds are formed with the backbone carbonyl groups of the surrounding β sheet residues Ala399, Leu325, and His326, which decreases the unity of the secondary structure elements. Notably, it is likely that the residue swap causes major problems during the C2 domain folding that are not visible in the variant simulations. In fact, even increased lability in the C2 domain could adversely affect the establishment of a stable SynGAP-membrane association.
c.1942T>CF648LLikely PathogenicGAPUncertain 1-9.296Likely Pathogenic0.999Likely PathogenicLikely Pathogenic0.468Likely Benign2.71Destabilizing0.82.08Destabilizing2.40Destabilizing1.04Destabilizing-5.98Deleterious0.999Probably Damaging0.976Probably Damaging3.45Benign0.08Tolerated201.0-34.02
c.1778T>AL593H
(3D Viewer)
Likely PathogenicGAPUncertain 1-16.504Likely Pathogenic0.998Likely PathogenicLikely Pathogenic0.812Likely Pathogenic2.52Destabilizing0.22.32Destabilizing2.42Destabilizing2.75Destabilizing-6.77Deleterious1.000Probably Damaging1.000Probably Damaging2.77Benign0.00Affected3.3735-2-3-7.023.98222.020.70.00.00.20.0XXPotentially PathogenicThe iso-propyl side chain of Leu593, located in an α helix (res. Glu582-Met603), packs favourably with multiple hydrophobic residues in the inter-helix space (e.g., Leu598, Ile589, Phe594, Phe561).In the variant simulations, His593 retains a similar packing arrangement via its aromatic imidazole ring. However, the polar nitrogen atoms introduce hydrogen bond donors and acceptors into the previously hydrophobic space. The epsilon protonated nitrogen of His593 forms a stable hydrogen bond with the phenol group of the Tyr505 side chain in an α helix (res. Gln503-Tyr518).While the residue swap could affect the tertiary assembly and the underlying protein folding process, it is difficult to determine if the mutation would be tolerated based solely on the variant simulations.
c.1402A>GM468V
(3D Viewer)
GAPUncertain 1-9.461Likely Pathogenic0.361AmbiguousLikely Benign0.570Likely Pathogenic2.69Destabilizing0.12.20Destabilizing2.45Destabilizing0.89Ambiguous-1.66Neutral0.998Probably Damaging0.993Probably Damaging-1.21Pathogenic0.08Tolerated3.3731122.3-32.06
c.896G>AR299H
(3D Viewer)
C2Conflicting 26-33437801-G-A106.20e-6-7.731In-Between0.388AmbiguousLikely Benign0.238Likely Benign3.97Destabilizing1.00.94Ambiguous2.46Destabilizing1.41Destabilizing-3.35Deleterious1.000Probably Damaging0.998Probably Damaging1.69Pathogenic0.02Affected3.3919201.3-19.05211.272.5-0.10.2-0.20.3XPotentially PathogenicThe guanidinium group of Arg299, located in a β hairpin loop linking two anti-parallel β sheet strands (res. Met289-Pro298, res. Thr305-Asn315), forms hydrogen bonds that stabilize the tight turn. In the WT simulations, the Arg299 side chain hydrogen bonds with the loop backbone carbonyl groups (e.g., Ser302, Thr305, Leu274, Gly303), the hydroxyl group of Ser300, and even forms a salt bridge with the carboxylate group of Asp304.In the variant simulations, the imidazole ring of His299 (epsilon protonated state) hydrogen bonds with the carbonyl group of Asp304 and the hydroxyl group of Ser300. However, it does not form as many or as strong interactions as arginine, which could affect the initial formation of the secondary hairpin loop during folding. β hairpins are potential nucleation sites during the initial stages of protein folding, so even minor changes in them could be significant.Additionally, His299 prefers to hydrophobically interact with other hydrophobic residues inside the C2 domain core (e.g., Val306, Leu274), which destabilizes the C2 domain. Indeed, the β strand partially unfolds during the second simulation. Moreover, the positively charged Arg299 side chain faces the polar head group region of the inner leaflet membrane and could directly anchor the C2 domain to the membrane. In short, the residue swap could negatively affect both protein folding and the stability of the SynGAP-membrane association.
c.1393C>GL465V
(3D Viewer)
Likely PathogenicGAPUncertain 1-9.893Likely Pathogenic0.838Likely PathogenicAmbiguous0.276Likely Benign2.46Destabilizing0.12.66Destabilizing2.56Destabilizing1.21Destabilizing-2.98Deleterious0.996Probably Damaging0.992Probably Damaging2.44Pathogenic0.10Tolerated3.3734210.4-14.03204.330.90.00.0-0.40.6XPotentially BenignThe iso-butyl side chain of Leu465, located in the middle of an α helix (res. Ala461–Phe476), packs with hydrophobic residues (e.g., Phe464, Met468, Tyr497, Ile494) in an inter-helix space formed with two other α helices (res. Ala461–Phe476 and res. Thr488-Gly502). In the variant simulations, the iso-propyl side chain of Val465 is equally sized and similarly hydrophobic as the original side chain of Leu465. Hence, the mutation does not exert any negative effects on the protein structure based on the variant simulations.
c.1025A>GY342C
(3D Viewer)
Likely PathogenicC2Benign/Likely benign 26-33437930-A-G211.30e-5-7.596In-Between0.682Likely PathogenicLikely Benign0.404Likely Benign2.48Destabilizing0.12.73Destabilizing2.61Destabilizing0.92Ambiguous-6.67Deleterious1.000Probably Damaging0.999Probably Damaging1.72Pathogenic0.02Affected3.37250-23.8-60.04242.462.80.10.0-0.10.2Potentially PathogenicThe phenol ring of Tyr342, located at the end of an anti-parallel β sheet strand (res. Gly341-Pro349), faces outward in the C2 domain. This phenol ring contributes to a triple tyrosine stack (Tyr342, Tyr328, and Tyr281) that links together three anti-parallel β sheet strands. Additionally, it shields Gly344 from the solvent, reducing its exposure and providing stability for the β-sandwich. This motif also contributes to a twist formation in the β sheet.In the variant simulations, the Cys342 side chain cannot participate in the stack formation. Instead, its thiol group forms a hydrogen bond with the backbone carbonyl group of Leu327. Although these changes in surface interactions could weaken the characteristic twist that strengthens the β sheet fold, no major structural effects are observed in the variant simulations. The residue swap could also affect the SynGAP-membrane association; however, this phenomenon cannot be addressed using solvent-only simulations. Notably, the thiol group of cysteine is not a particularly strong hydrogen-bonding partner, which could mitigate the negative effects of the residue swap.
c.1957C>GL653VLikely BenignGAPUncertain 1-7.050In-Between0.301Likely BenignLikely Benign0.146Likely Benign3.28Destabilizing0.32.18Destabilizing2.73Destabilizing1.32Destabilizing-2.25Neutral0.227Benign0.039Benign3.28Benign0.08Tolerated210.4-14.03
c.1877T>CI626TLikely PathogenicGAPUncertain 1-10.420Likely Pathogenic0.946Likely PathogenicAmbiguous0.640Likely Pathogenic2.94Destabilizing0.12.70Destabilizing2.82Destabilizing2.23Destabilizing-4.18Deleterious1.000Probably Damaging1.000Probably Damaging3.04Benign0.00Affected0-1-5.2-12.05
c.928G>AE310K
(3D Viewer)
Likely PathogenicC2Conflicting 4-14.601Likely Pathogenic0.997Likely PathogenicLikely Pathogenic0.764Likely Pathogenic1.97Ambiguous1.23.66Destabilizing2.82Destabilizing1.02Destabilizing-3.68Deleterious1.000Probably Damaging0.995Probably Damaging1.19Pathogenic0.01Affected3.381901-0.4-0.94213.458.00.10.00.20.1XPotentially PathogenicThe carboxylate group of Glu310, located in an anti-parallel β sheet strand (res. Thr305-Asn315), is ideally positioned to interact with the side chain hydroxyl and backbone amide groups of Thr295 on a twisted anti-parallel β strand (res. Met289-Arg299). Because the carboxylate group can also interact with the GAP domain residues (e.g., Gln612, Tyr614), Glu310 plays a key role in maintaining the tertiary assembly between the C2 and GAP domains. In the variant simulations, the amino group of the Lys310 side chain hydrogen bonds with the GAP domain residues and forms a salt bridge with Glu613. Although no apparent negative effects are seen due to the residue swap, it is possible that the loss of hydrogen bonding with the hydroxyl group of the Thr295 side chain causes problems during folding, potentially compromising the twisting of the β sheet.
c.2131C>GL711V
(3D Viewer)
Likely PathogenicGAPUncertain16-33441596-C-G16.20e-7-10.045Likely Pathogenic0.709Likely PathogenicLikely Benign0.170Likely Benign3.48Destabilizing0.12.22Destabilizing2.85Destabilizing1.40Destabilizing-2.59Deleterious0.992Probably Damaging0.970Probably Damaging3.34Benign0.00Affected3.509120.4-14.03
c.1513T>CY505H
(3D Viewer)
Likely PathogenicGAPLikely Pathogenic 1-11.383Likely Pathogenic0.982Likely PathogenicLikely Pathogenic0.646Likely Pathogenic2.91Destabilizing0.12.88Destabilizing2.90Destabilizing1.60Destabilizing-4.97Deleterious1.000Probably Damaging1.000Probably Damaging2.64Benign0.00Affected3.373520-1.9-26.03
c.791T>AL264Q
(3D Viewer)
Likely PathogenicC2Uncertain 1-15.729Likely Pathogenic0.999Likely PathogenicLikely Pathogenic0.678Likely Pathogenic3.43Destabilizing0.12.41Destabilizing2.92Destabilizing2.48Destabilizing-5.52Deleterious1.000Probably Damaging0.999Probably Damaging0.49Pathogenic0.00Affected3.3818-2-2-7.314.97254.7-7.60.00.00.00.3XXXPotentially PathogenicThe iso-butyl branched hydrocarbon side chain of Leu264, located at the end of an anti-parallel β sheet strand (res. Arg259-Arg272), packs against multiple hydrophobic residues such as Leu266, Phe314, Leu317, and Leu323 in the WT simulations. In the variant simulations, the hydrophilic carboxamide group of the Gln264 side chain is not suitable for the hydrophobic niche, causing the hydrophobic residues to make room for the swapped residue. Additionally, the carboxamide group of Gln264 forms hydrogen bonds with the backbone amide groups of Arg405 and Lys256 in the β sheet and the carbonyl group of Val350 in an α helical section of a nearby loop (res. Pro359-Phe358). The residue swap disrupts the packing of the C2 domain, which could adversely affect the C2 domain structure during folding. This disruption could potentially weaken the stability of the SynGAP-membrane association.
c.1409T>CM470T
(3D Viewer)
Likely PathogenicGAPUncertain 1-8.104Likely Pathogenic0.976Likely PathogenicLikely Pathogenic0.763Likely Pathogenic3.19Destabilizing0.12.68Destabilizing2.94Destabilizing1.49Destabilizing-5.30Deleterious0.996Probably Damaging0.985Probably Damaging-1.08Pathogenic0.24Tolerated3.3734-1-1-2.6-30.09213.846.50.00.0-0.20.2XXPotentially PathogenicThe thioether group of Met470, located in the middle of an α helix (res. Ala461–Phe476), interacts with hydrophobic residues in the inter-helix space (e.g., Val473, Leu558, Cys576, Trp572) formed by two other α helices (res. Ser604–Arg581, res. Pro562–Arg579). In the WT simulations, the Met470 side chain also packs against the positively charged guanidinium groups of Arg575, Arg429, and Arg579, which form salt bridges with the negatively charged carboxylate groups of the Asp474 and Asp467 side chains at the protein surface. In the variant simulations, the hydroxyl group of the Thr470 side chain forms an H-bond with the backbone carbonyl group of Ser466 in the α helix, potentially lowering its structural integrity. Importantly, the hydroxyl group of Thr470 also forms an H-bond with the guanidinium group of Arg575, which helps it form a more permanent salt bridge with Asp467.
c.1025A>CY342S
(3D Viewer)
Likely PathogenicC2Uncertain 2-7.996In-Between0.925Likely PathogenicAmbiguous0.407Likely Benign3.03Destabilizing0.12.87Destabilizing2.95Destabilizing0.93Ambiguous-6.60Deleterious1.000Probably Damaging0.998Probably Damaging1.75Pathogenic0.04Affected3.3725-3-20.5-76.10200.177.80.00.0-0.20.1Potentially PathogenicThe phenol ring of Tyr342, located at the end of an anti-parallel β sheet strand (res. Gly341-Pro349), faces outward in the C2 domain. In the WT simulations, the phenol ring of Tyr342 contributes to a triple tyrosine stack (Tyr342, Tyr328, and Tyr281) that links together three anti-parallel β sheet strands. Additionally, it shields Gly344 from the solvent, reducing its exposure and providing stability for the β-sandwich. This motif also contributes to a twist formation in the β sheet.In the variant simulations, the Ser342 side chain cannot participate in the stack formation. Instead, the hydroxyl group of the Ser342 side chain forms a hydrogen bond with the imidazole ring of His326 in a neighboring β strand (res. Ala322-Asp330). This disrupts the formation of a hydrogen bond between His326 and the carboxylate group of the Glu283 side chain from another β strand (res. Arg279-Cys285). Although these changes in surface interactions could weaken the characteristic twist that strengthens the β sheet fold, no major structural effects are observed in the variant simulations. The residue swap could also affect the SynGAP-membrane association, as the hydroxyl group of Ser342 could form hydrogen bonds with membrane-facing loop residues. However, this phenomenon cannot be addressed using solvent-only simulations.
c.929A>GE310G
(3D Viewer)
Likely PathogenicC2Pathogenic 1-14.132Likely Pathogenic0.995Likely PathogenicLikely Pathogenic0.848Likely Pathogenic2.38Destabilizing0.73.56Destabilizing2.97Destabilizing0.36Likely Benign-6.43Deleterious1.000Probably Damaging0.996Probably Damaging1.12Pathogenic0.00Affected3.3819-203.1-72.06
c.1540A>TI514F
(3D Viewer)
Likely PathogenicGAPUncertain 1-13.383Likely Pathogenic0.962Likely PathogenicLikely Pathogenic0.601Likely Pathogenic2.35Destabilizing0.33.74Destabilizing3.05Destabilizing0.93Ambiguous-3.98Deleterious0.997Probably Damaging0.993Probably Damaging2.89Benign0.00Affected3.373501-1.734.02
c.1118G>TG373V
(3D Viewer)
Likely BenignC2Uncertain 16-33438023-G-T65.03e-6-6.062Likely Benign0.112Likely BenignLikely Benign0.428Likely Benign5.32Destabilizing3.20.82Ambiguous3.07Destabilizing0.09Likely Benign-0.98Neutral0.007Benign0.001Benign3.90Benign0.00Affected3.5316-1-34.642.08207.6-68.11.91.1-0.60.1UncertainGly373 is located in the Gly-rich Ω loop (res. Pro364-Pro398) between two anti-parallel β sheet strands (res. Thr359-Pro364, res. Ala399-Ile411). Because the Ω loop is assumed to directly interact with the membrane, it moves arbitrarily throughout the WT solvent simulations. The Ω loop potentially plays a crucial role in the SynGAP-membrane complex association, stability, and dynamics. However, this aspect cannot be fully addressed through solvent simulations alone.Ω loops are known to play major roles in protein functions that require flexibility, and thus hydrophobic residues like valine are rarely tolerated. Although no negative structural effects are observed in the variant simulations, Val373 may exert drastic effects on the SynGAP-membrane complex dynamics and stability. However, since the effect on the Gly-rich Ω loop dynamics can only be studied through the SynGAP-membrane complex, no definite conclusions can be drawn.
c.1199T>AV400E
(3D Viewer)
Likely PathogenicC2Uncertain 1-13.686Likely Pathogenic0.998Likely PathogenicLikely Pathogenic0.810Likely Pathogenic3.70Destabilizing0.22.46Destabilizing3.08Destabilizing2.29Destabilizing-4.88Deleterious0.920Possibly Damaging0.335Benign5.31Benign0.00Affected3.3827-2-2-7.729.98249.1-38.8-0.10.11.00.0XXXPotentially PathogenicThe iso-propyl side chain of Val400, located in an anti-parallel β sheet strand (res. Ala399-Ile411), hydrophobically packs against hydrophobic residues within the anti-parallel β sheet of the C2 domain (e.g., Ile268, Ala404, Leu325, Leu402). In the variant simulations, the negatively charged carboxylate group of the Glu400 side chain is not suitable for occupying the hydrophobic niche. Consequently, the side chain escapes the center of the C2 domain and interacts with the backbone amide groups of Leu402 in the same β strand and/or Ile269 and Glu270 in a neighboring β strand (res. Arg259-Arg272). This residue swap disrupts the hydrophobic packing and generally has extensive negative effects on the C2 domain structure. At a minimum, the residue swap could affect the C2 domain stability and membrane association.
c.1493T>GM498R
(3D Viewer)
Likely PathogenicGAPLikely Pathogenic 1-8.812Likely Pathogenic0.988Likely PathogenicLikely Pathogenic0.869Likely Pathogenic3.85Destabilizing0.22.35Destabilizing3.10Destabilizing1.76Destabilizing-4.53Deleterious0.464Possibly Damaging0.120Benign-1.36Pathogenic0.00Affected0-1-6.424.99
c.2050G>CD684H
(3D Viewer)
Likely PathogenicGAPUncertain 1-14.194Likely Pathogenic0.998Likely PathogenicLikely Pathogenic0.613Likely Pathogenic3.36Destabilizing1.02.95Destabilizing3.16Destabilizing0.55Ambiguous-6.98Deleterious1.000Probably Damaging0.972Probably Damaging3.36Benign0.00Affected3.4217-110.322.05
c.872A>GY291C
(3D Viewer)
Likely PathogenicC2Uncertain 1-8.997Likely Pathogenic0.967Likely PathogenicLikely Pathogenic0.505Likely Pathogenic2.90Destabilizing0.43.51Destabilizing3.21Destabilizing1.35Destabilizing-7.37Deleterious1.000Probably Damaging0.999Probably Damaging1.76Pathogenic0.01Affected3.38230-23.8-60.04205.266.10.10.0-0.40.4XXPotentially PathogenicThe phenol group of the Tyr291 side chain, located in an anti-parallel β sheet strand (res. Met289-Pro298), packs against hydrophobic residues of the C2 and PH domains (e.g., Leu317, Leu286, Leu284, Pro208, Val209). The phenol ring of Tyr291 also forms favorable Met-aromatic stacking with the methyl group of Met289. In the variant simulation, the thiol group of the Cys291 side chain is not as suitable for the hydrophobic inter-domain space as the phenol ring of Tyr291. Consequently, the structural unity of the PH domain is weakened and ultimately unfolds in the second simulation. Moreover, the residue swap might result in severe detrimental effects on the C2 domain structure and the C2-PH domain tertiary structure assembly during folding.
c.1390T>GF464V
(3D Viewer)
Likely PathogenicGAPUncertain 1-12.254Likely Pathogenic0.994Likely PathogenicLikely Pathogenic0.592Likely Pathogenic3.61Destabilizing0.12.89Destabilizing3.25Destabilizing1.40Destabilizing-6.96Deleterious0.998Probably Damaging0.996Probably Damaging3.36Benign0.04Affected3.3734-1-11.4-48.04210.140.5-0.10.0-0.90.3XPotentially PathogenicThe phenyl ring of Phe464, located in the middle of an α helix (res. Ala461–Phe476), packs against hydrophobic residues (e.g., Met468, Leu451, Leu455, and Tyr428) in the inter-helix space formed with two other α helices (res. Asn440-Lys460 and res. Pro413-Glu436). The iso-propyl side chain of Val464 is similarly hydrophobic but considerably smaller than the original phenyl ring of Phe464. To compensate for the size difference, neighboring residues need to fill in the gap in the variant simulations.The phenolic side chain of Tyr428, located at the middle bend of an α helix (res. Glu436-Pro413), assumes a new position in the inter-helix space or rotates inward next to the third α helix (res. Asn440-Lys460) when the stable H-bond between Tyr428 and Asp467 seen in the WT simulations breaks. The residue swap also leads to the loss of the methionine-aromatic interaction between the Met468 and Phe464 side chains, which could weaken the integrity of the parent α helix (res. Ala461-Phe476). Although the simulations likely underestimate the full adverse effect of the introduced mutation during folding, the two opposing α helices (res. Ala461–Phe476 and res. Glu436-Pro413) move substantially closer to each other in the variant simulations.
c.1403T>AM468K
(3D Viewer)
Likely PathogenicGAPLikely Pathogenic 1-16.982Likely Pathogenic0.978Likely PathogenicLikely Pathogenic0.828Likely Pathogenic3.21Destabilizing0.13.30Destabilizing3.26Destabilizing2.57Destabilizing-4.61Deleterious0.878Possibly Damaging0.922Probably Damaging-1.34Pathogenic0.04Affected3.37310-1-5.8-3.02188.769.30.00.0-0.10.2XXPotentially PathogenicThe thioether group of Met468, located in the middle of an α helix (res. Ala461–Phe476), interacts with hydrophobic residues (e.g., Phe464, Leu465, Leu489) in an inter-helix space formed by two other α helices (res. Ala461–Phe476, res. Thr488–Gly502). In the variant simulations, the positively charged side chain of Lys468 rotates outward to escape the hydrophobic niche, forming an H-bond with the hydroxyl group of the Ser471 side chain and a salt bridge with the carboxylate group of the Glu472 side chain. This residue swap also disrupts the methionine-aromatic stacking with the phenyl ring of the Phe464 side chain. Although no large-scale structural changes are observed during the variant simulations, the importance of hydrophobic packing suggests that the effects could be more pronounced during protein folding.
c.2075T>AL692Q
(3D Viewer)
Likely PathogenicGAPPathogenic 1-13.873Likely Pathogenic0.998Likely PathogenicLikely Pathogenic0.596Likely Pathogenic3.24Destabilizing0.13.27Destabilizing3.26Destabilizing2.76Destabilizing-5.98Deleterious1.000Probably Damaging0.998Probably Damaging3.06Benign0.00Affected3.4217-2-2-7.314.97
c.1084T>CW362R
(3D Viewer)
Likely PathogenicC2Pathogenic 2-14.004Likely Pathogenic1.000Likely PathogenicLikely Pathogenic0.706Likely Pathogenic2.64Destabilizing0.33.90Destabilizing3.27Destabilizing1.10Destabilizing-12.87Deleterious0.999Probably Damaging0.996Probably Damaging1.28Pathogenic0.00Affected3.39242-3-3.6-30.03287.5-34.1-0.20.1-0.60.2XXXPotentially PathogenicThe indole ring of Trp362, located on the surface of an anti-parallel β sheet (res. Thr359-Pro364) in the C2 domain, stacks with nearby residues (e.g., Arg401, Arg272). In the variant simulations, the guanidinium group of the introduced residue Arg362 forms a salt bridge with the carboxylate group of Glu273 and, like Trp362, stacks with other arginine residues (e.g., Arg401, Arg272). This residue is at both the C2-membrane interface and the C2-RasGTPase interface, so the residue swap could potentially affect both interactions. However, these phenomena cannot be addressed using solvent-only simulations. Notably, Arg272, which stacks with both the non-mutated Trp362 and the mutated Arg362, forms a salt bridge directly with Asp105 of Ras in the WT simulations. Therefore, the residue swap could affect the C2 domain stability, the SynGAP-membrane association, and the SynGAP-Ras association.10.1016/j.ajhg.2020.11.011
c.1082A>CQ361P
(3D Viewer)
Likely PathogenicC2Likely Pathogenic 1-13.280Likely Pathogenic0.956Likely PathogenicLikely Pathogenic0.482Likely Benign3.12Destabilizing0.03.45Destabilizing3.29Destabilizing0.38Likely Benign-3.03Deleterious0.996Probably Damaging0.979Probably Damaging1.63Pathogenic0.05Affected3.3725-101.9-31.01
c.1403T>CM468T
(3D Viewer)
Likely PathogenicGAPUncertain 26-33438435-T-C16.20e-7-12.399Likely Pathogenic0.862Likely PathogenicAmbiguous0.801Likely Pathogenic3.47Destabilizing0.13.10Destabilizing3.29Destabilizing1.84Destabilizing-3.85Deleterious0.994Probably Damaging0.985Probably Damaging-1.31Pathogenic0.01Affected3.3731-1-1-2.6-30.09214.647.10.00.00.10.0XPotentially PathogenicThe thioether group of Met468, located in the middle of an α helix (res. Ala461–Phe476), interacts with hydrophobic residues (e.g., Phe464, Leu465, Leu489) in an inter-helix space formed by two other α helices (res. Ala461–Phe476, res. Thr488–Gly502). In the variant simulations, the hydrophilic side chain of Thr468 does not pack favorably in the hydrophobic niche, and the methionine-aromatic stacking is lost. Although the hydroxyl group of Thr468 forms an H-bond with the backbone carbonyl group of Phe464, the integrity of the α helix is not affected in the simulations. No large-scale structural changes are observed during the variant simulations; however, due to the importance of hydrophobic packing, the effects could be more pronounced during protein folding.
c.878G>AR293HLikely PathogenicC2Uncertain 1-13.009Likely Pathogenic0.973Likely PathogenicLikely Pathogenic0.438Likely Benign4.45Destabilizing2.32.12Destabilizing3.29Destabilizing0.32Likely Benign-4.60Deleterious1.000Probably Damaging0.998Probably Damaging1.45Pathogenic0.04Affected201.3-19.05
c.1436G>CR479P
(3D Viewer)
Likely PathogenicGAPUncertain 1-11.795Likely Pathogenic0.938Likely PathogenicAmbiguous0.277Likely Benign2.86Destabilizing0.23.88Destabilizing3.37Destabilizing0.81Ambiguous-3.52Deleterious1.000Probably Damaging1.000Probably Damaging3.41Benign0.18Tolerated0-22.9-59.07
c.1813C>TP605S
(3D Viewer)
Likely PathogenicGAPUncertain 1-10.830Likely Pathogenic0.987Likely PathogenicLikely Pathogenic0.718Likely Pathogenic3.40Destabilizing0.13.34Destabilizing3.37Destabilizing1.00Destabilizing-7.96Deleterious1.000Probably Damaging1.000Probably Damaging0.70Pathogenic0.00Affected3.37351-10.8-10.04213.8-15.4-0.30.20.20.1XXPotentially PathogenicPro605 is located in a short turn between an α helix (res. Glu582-Met603) and a short α helical section (res. Ser606-Phe608). The pyrrolidine side chain of Pro605 packs hydrophobically with nearby hydrophobic residues (e.g., Ile514, Leu623, Leu610) in the inter-helix space. Additionally, proline lacks a free backbone amide group, which breaks the α helix and facilitates the turn in the WT structure.In the variant simulations, the hydroxyl side chain of Ser605 forms hydrogen bonds with the backbone carbonyl groups of Ala601 and Ile602. Importantly, the helix end is more stable than with Pro605 in the WT. Indeed, proline is a more effective secondary structure breaker compared to serine.Thus, the residue swap could have a more profound effect on the actual folding process, for example, by preventing the bending at the α helix end, than what the simulations suggest. Moreover, due to its location at the GAP-Ras interface, the residue swap could affect the GAP-Ras association.
c.1169G>AG390E
(3D Viewer)
C2Uncertain 1-7.913In-Between0.646Likely PathogenicLikely Benign0.575Likely Pathogenic2.61Destabilizing0.94.28Destabilizing3.45Destabilizing0.47Likely Benign-0.87Neutral0.276Benign0.045Benign1.32Pathogenic0.05Affected4.3280-2-3.172.06241.5-108.40.60.5-0.10.1UncertainGly390 is located in the Gly-rich Ω loop (res. Pro364-Pro398) between two anti-parallel β sheet strands (res. Thr359-Pro364 and res. Ala399-Ile411). The Ω loop is assumed to directly interact with the membrane, and it is observed to move arbitrarily throughout the WT solvent simulations. This loop potentially plays a crucial role in the SynGAP-membrane complex association, stability, and dynamics. However, this aspect cannot be fully addressed through solvent simulations alone.Ω loops are known to play significant roles in protein functions that require flexibility, and so they are rich in glycine residues, prolines, and to a lesser extent, small hydrophilic residues to ensure maximum flexibility. Thus, the variant’s Glu390 may not be as well tolerated in the Ω loop. Additionally, the carboxylate group of Glu390 occasionally forms H-bonds with other loop residues in the variant simulations. The interaction between the acidic carboxylate side chain and the acidic membrane lipids may further influence the SynGAP-membrane complex. However, since the effects on the Gly-rich Ω loop dynamics can only be well studied through the SynGAP-membrane complex, no definite conclusions can be drawn.
c.1466T>CL489P
(3D Viewer)
Likely PathogenicGAPConflicting 2-13.520Likely Pathogenic0.997Likely PathogenicLikely Pathogenic0.939Likely Pathogenic2.50Destabilizing0.14.69Destabilizing3.60Destabilizing1.73Destabilizing-6.74Deleterious1.000Probably Damaging1.000Probably Damaging-1.56Pathogenic0.00Affected3.3735-3-3-5.4-16.04209.961.90.10.00.60.1XPotentially PathogenicThe iso-butyl side chain of Leu489, located in the α-helix (res. Leu489-Glu519) within an inter-helix space of four helices (res. Ala461-Phe476, res. Val441-Ser457, and res. Met414-Glu436), packs with hydrophobic residues (e.g., Cys432, Ala448, Lys444, Ala493, Val447, Met468). In the variant simulations, Pro489 is located near the beginning of the α-helix, so the residue swap with Leu489 does not affect the continuity of the secondary structure element. However, the side chain of proline is not as optimal as that of leucine for maintaining hydrophobic packing with nearby residues (e.g., Ala448, Lys444). Additionally, the consistently maintained hydrogen bond interaction between the backbone amide group of Leu489 and the carbonyl of Glu436 is lost due to the residue swap, potentially affecting the tertiary structure integrity.
c.1991T>CL664S
(3D Viewer)
Likely PathogenicGAPLikely Benign 16-33441250-T-C16.20e-7-16.498Likely Pathogenic0.997Likely PathogenicLikely Pathogenic0.543Likely Pathogenic3.75Destabilizing0.23.63Destabilizing3.69Destabilizing2.77Destabilizing-5.99Deleterious1.000Probably Damaging0.996Probably Damaging2.85Benign0.00Affected3.3828-3-2-4.6-26.08215.550.10.00.0-0.20.2XPotentially BenignThe iso-butyl side chain of L664, located on an α-helix (res. Ser641-Glu666), hydrophobically interacts with residues in the inter-helix space between three helices (res. Glu617-Asn635, res. Glu582-Met603, and res. Ser641-Glu666), such as Ile589, Phe663, and Met660. In the variant simulations, the hydroxyl group of Ser664 forms hydrogen bonds with the backbone carbonyl oxygen of another helix residue, such as Met660 or Gln661. This interaction is known to destabilize hydrogen bonding in the α-helix, but this effect was not observed in the simulations. Additionally, Ser664 occasionally forms hydrogen bonds with the carboxylate group of Asp586 on another α-helix (res. Glu582-Met603), which could minimally influence the tertiary structure assembly. Despite these interactions, no major negative effects on the protein structure were observed during the simulations.
c.812C>AA271D
(3D Viewer)
Likely PathogenicC2Pathogenic 1-18.590Likely Pathogenic0.999Likely PathogenicLikely Pathogenic0.706Likely Pathogenic4.71Destabilizing0.42.67Destabilizing3.69Destabilizing1.59Destabilizing-5.52Deleterious1.000Probably Damaging0.999Probably Damaging0.62Pathogenic0.00Affected3.38190-2-5.344.01226.2-63.40.00.00.90.1XXXXPotentially PathogenicThe methyl group of Ala271, located near the end of an anti-parallel β sheet strand (res. Arg259-Arg272), packs against multiple hydrophobic residues such as Val400, Val306, and Leu274 in the WT simulations. In the variant simulations, the carboxylate group of Asp271 is not suitable for the hydrophobic niche, causing the hydrophobic residues to make room for the swapped residue. Additionally, the carboxylate group of the Asp271 side chain forms hydrogen bonds with the backbone amide groups of Arg272 and Ala399 in the β sheet, or even forms a salt bridge with the amino group of the Lys394 side chain. This directly affects the integrity of the anti-parallel β sheet at the end. In short, the residue swap disrupts the C2 domain packing during folding, which could weaken the stability of the SynGAP-membrane association.
c.1529T>GI510S
(3D Viewer)
Likely PathogenicGAPLikely Pathogenic 1-11.661Likely Pathogenic0.955Likely PathogenicAmbiguous0.926Likely Pathogenic4.00Destabilizing0.13.78Destabilizing3.89Destabilizing2.34Destabilizing-4.63Deleterious1.000Probably Damaging0.999Probably Damaging-1.44Pathogenic0.00Affected3.3735-1-2-5.3-26.08201.445.9-0.40.20.00.3XPotentially PathogenicIle510 is located in the middle of an α-helix (res. Gly502-Tyr518) within the inter-helix space of three helices (res. Gly502-Tyr518, Ala533-Val560, and res. Glu582-Met603). In the WT simulations, the sec-butyl side chain of Ile510 hydrophobically packs with other residues in the inter-helix space (e.g., Leu506, Leu610, Ile514, Ile602, Leu598). In the variant simulations, the hydroxyl group of Ser510 forms a hydrogen bond with the backbone atoms of Leu506 and Gly511 in the same α-helix, which could further weaken the α-helix integrity. This α-helix already shows weakness in the WT simulations due to Gly511. Although the simulations do not show large-scale effects, the residue swap could have a substantial impact due to the fundamental role of hydrophobic packing during protein folding.
c.1136C>GS379W
(3D Viewer)
C2Uncertain 16-33438041-C-G-8.898Likely Pathogenic0.388AmbiguousLikely Benign0.520Likely Pathogenic4.32Destabilizing3.43.56Destabilizing3.94Destabilizing0.16Likely Benign-1.02Neutral0.998Probably Damaging0.844Possibly Damaging3.82Benign0.01Affected4.3211-2-3-0.199.14271.3-75.71.41.00.60.5UncertainSer379 is located in the Gly-rich Ω loop (res. Pro364-Pro398) between two anti-parallel β sheet strands (res. Thr359-Pro364, res. Ala399-Ile411). Because the Ω loop is assumed to directly interact with the membrane, it moves arbitrarily throughout the WT solvent simulations. The Ω loop potentially plays a crucial role in the SynGAP-membrane complex association, stability, and dynamics. However, this aspect cannot be fully addressed through solvent simulations alone.Ω loops are known to play major roles in protein functions that require flexibility, and thus hydrophobic residues like tryptophan are rarely tolerated. Although no major negative structural effects are observed in the variant simulations, Trp379 may exert drastic effects on the SynGAP-membrane complex dynamics and stability. However, since the effect on Gly-rich Ω loop dynamics can only be studied through the SynGAP-membrane complex, no definite conclusions can be drawn
c.2162T>GI721S
(3D Viewer)
Likely PathogenicGAPUncertain 1-14.032Likely Pathogenic0.996Likely PathogenicLikely Pathogenic0.466Likely Benign3.91Destabilizing0.13.96Destabilizing3.94Destabilizing2.28Destabilizing-5.26Deleterious1.000Probably Damaging1.000Probably Damaging2.21Pathogenic0.00Affected3.509-1-2-5.3-26.08203.349.3-0.10.0-1.10.0XUncertainThe sec-butyl side chain of Ile721, located on an α-helix (res. Leu714-Arg726), engages in hydrophobic packing with other residues in the hydrophobic inter-helix space, such as Phe420, Tyr417, His693, and Leu717. In the variant simulations, the hydroxyl side chain of Ser721 forms hydrogen bonds with nearby residues, such as Leu717 and His693. Although no major structural changes are observed during the variant simulations, the hydrophilic residue Ser721 could disrupt the hydrophobic packing during folding. However, because the model ends abruptly at the C-terminus, no definite conclusions can be drawn based on the simulations.
c.1763T>AL588H
(3D Viewer)
Likely PathogenicGAPLikely Pathogenic 1-16.947Likely Pathogenic0.999Likely PathogenicLikely Pathogenic0.939Likely Pathogenic4.20Destabilizing0.23.69Destabilizing3.95Destabilizing2.26Destabilizing-6.97Deleterious1.000Probably Damaging1.000Probably Damaging-1.42Pathogenic0.00Affected3.3834-2-3-7.023.98214.320.90.00.00.00.2XXXPotentially PathogenicThe isobutyl group of the Leu588 side chain, located in an α helix (res. Glu582-Met603), packs against hydrophobic residues in the inter-helix hydrophobic space (e.g., Ile584, Trp572, Phe484, Met470, Val473, Ile483).In the variant simulations, the imidazole ring of His588 is aromatic but contains polar delta and epsilon nitrogen atoms that are not suited for the hydrophobic niche. The protonated epsilon nitrogen forms a hydrogen bond with the backbone carbonyl group of Ala469, which can disrupt the continuity of the opposing α helix (res. Phe476-Lys460).While the residue swap could affect the tertiary assembly and the underlying protein folding process, it is difficult to determine if the mutation would be tolerated based solely on the variant simulations.
c.1205T>GL402R
(3D Viewer)
Likely PathogenicC2Likely Pathogenic1-13.800Likely Pathogenic0.997Likely PathogenicLikely Pathogenic0.522Likely Pathogenic4.10Destabilizing0.23.82Destabilizing3.96Destabilizing2.24Destabilizing-4.69Deleterious0.967Probably Damaging0.459Possibly Damaging3.69Benign0.00Affected3.3828-3-2-8.343.03259.5-55.40.00.01.40.0XXXPotentially PathogenicThe iso-butyl side chain of Leu402, located in an anti-parallel β sheet strand (res. Ala399-Ile411), packs with residues inside the hydrophobic core of the C2 domain (e.g., Ile268, Ala404, Leu266, Val400). In the variant simulations, the positively charged guanidinium group of the Arg402 side chain is not suitable for the hydrophobic niche. Consequently, the side chain moves outward from the hydrophobic C2 domain core and stacks with the phenol ring of Tyr363 or forms H-bonds with the carboxamide group of the Gln361 side chain in the β sheet strand (res. Thr359-Tyr364). This movement induces extensive negative effects on the C2 domain structure.
c.680G>AG227E
(3D Viewer)
Likely PathogenicPHConflicting 26-33435531-G-A31.86e-6-9.186Likely Pathogenic0.996Likely PathogenicLikely Pathogenic0.792Likely Pathogenic2.56Destabilizing0.45.36Destabilizing3.96Destabilizing0.94Ambiguous-6.49Deleterious0.906Possibly Damaging0.360Benign5.72Benign0.01Affected3.43120-2-3.172.06237.7-112.10.10.30.00.3XXUncertainThe introduced residue Glu227 is located in a β hairpin loop connecting two anti-parallel β sheet strands (res. Cys219-Thr224 and Thr228-Ala232). In the variant simulations, the carboxylate group of Glu227 frequently forms a salt bridge with the amino group of the neighboring residue Lys229. Despite this interaction, the integrity of the secondary structure element is not compromised. However, the β hairpins are potential nucleation sites during the initial stages of protein folding. Additionally, since the model ends abruptly at the N-terminus, no definite conclusions can be drawn from the simulations.
c.1784T>CL595P
(3D Viewer)
Likely PathogenicGAPUncertain 1-11.856Likely Pathogenic0.999Likely PathogenicLikely Pathogenic0.747Likely Pathogenic2.09Destabilizing0.85.88Destabilizing3.99Destabilizing1.78Destabilizing-6.97Deleterious1.000Probably Damaging1.000Probably Damaging2.72Benign0.00Affected3.3735-3-3-5.4-16.04
c.968T>GL323R
(3D Viewer)
Likely PathogenicC2Likely Pathogenic 1-14.568Likely Pathogenic0.997Likely PathogenicLikely Pathogenic0.692Likely Pathogenic3.75Destabilizing0.44.47Destabilizing4.11Destabilizing2.15Destabilizing-4.70Deleterious0.999Probably Damaging0.969Probably Damaging0.59Pathogenic0.01Affected3.3922-3-2-8.343.03261.8-61.6-0.40.20.80.2XXXPotentially PathogenicThe iso-butyl side chain of Leu323, located at the beginning of an anti-parallel β sheet strand (res. Ala322-Asp330), packs against multiple hydrophobic leucine residues (e.g., Leu264, Leu266, Leu284, Leu286). In contrast, in the variant simulations, the positively charged guanidinium group of the Arg323 side chain is unsuitable for the hydrophobic niche. Consequently, the side chain either rotates away from the center of the C2 domain or, if it remains within the C2 domain core, it reorients nearby residues to form hydrogen bonds. Regardless, the residue swap extensively disrupts the C2 domain structure.
c.1214G>CR405P
(3D Viewer)
Likely PathogenicC2Uncertain 1-14.206Likely Pathogenic0.999Likely PathogenicLikely Pathogenic0.572Likely Pathogenic3.11Destabilizing0.35.19Destabilizing4.15Destabilizing1.26Destabilizing-6.32Deleterious1.000Probably Damaging1.000Probably Damaging3.62Benign0.01Affected3.3828-202.9-59.07
c.1490A>GY497C
(3D Viewer)
Likely PathogenicGAPUncertain 1-11.872Likely Pathogenic0.948Likely PathogenicAmbiguous0.806Likely Pathogenic3.88Destabilizing0.14.76Destabilizing4.32Destabilizing1.40Destabilizing-8.82Deleterious1.000Probably Damaging0.995Probably Damaging-1.65Pathogenic0.03Affected3.37350-23.8-60.04209.959.1-0.10.0-0.30.1XXPotentially PathogenicTyr497 is located in the α-helix (res. Leu489-Glu519) within the inter-helix space of four α-helices (res. Leu489-Ile501, res. Val441-Ser457, res. Arg563-Glu578, res. Ala461-Val473). In the WT simulations, the phenol ring of Tyr497 hydrophobically packs with other residues in the inter-helix space (e.g., Leu465, Leu565, Val568). The hydroxyl group of Tyr497 also alternately forms hydrogen bonds with the carboxylate side chain of Gln456 and the backbone carbonyl of Glu564. Thus, Tyr497 plays a role in the folding and maintenance of the tertiary structure assembly between these four helices.In the variant simulations, the comparatively smaller residue, Cys497, cannot maintain any of the interactions seen with Tyr497 in the WT. Although no severe deleterious consequences are observed in the simulations, the structural effects could be more pronounced during actual protein folding. Indeed, the tertiary structure is seen to slightly break apart in the variant simulations.
c.917T>AV306D
(3D Viewer)
Likely PathogenicC2Uncertain 1-18.289Likely Pathogenic0.986Likely PathogenicLikely Pathogenic0.530Likely Pathogenic4.40Destabilizing0.34.29Destabilizing4.35Destabilizing2.44Destabilizing-5.44Deleterious1.000Probably Damaging0.999Probably Damaging1.74Pathogenic0.00Affected3.3819-2-3-7.715.96212.3-18.3-0.20.40.00.2XXXPotentially PathogenicThe isopropyl group of Val396, located at the beginning of an anti-parallel β sheet strand (res. Thr305-Asn315), packs against multiple hydrophobic residues (e.g., Leu274, Trp308, Ala271) in the WT simulations. However, in the variant simulations, the negatively charged carboxylate group of the Asp306 side chain is not suitable for this hydrophobic niche. Consequently, the side chain moves out to interact with Ser300 in the β strand (res. Met289-Arg299) and the guanidinium group of Arg299 in the β hairpin loop.In the third simulation, the residue swap disrupts the C2 domain secondary structure and tertiary assembly to a large degree when the amino group of the Lys297 side chain rotates to form a salt bridge with Asp306. This drastic effect could potentially reflect the challenge presented by the residue swap during the C2 domain folding. Because the residue swap affects the C2 domain structure, the SynGAP-membrane association could also be impacted. However, this is beyond the scope of the solvent-only simulations to unravel.
c.1631G>CR544P
(3D Viewer)
Likely PathogenicGAPUncertain 2-16.905Likely Pathogenic1.000Likely PathogenicLikely Pathogenic0.762Likely Pathogenic4.70Destabilizing0.14.19Destabilizing4.45Destabilizing1.14Destabilizing-4.88Deleterious1.000Probably Damaging1.000Probably Damaging-1.48Pathogenic0.05Affected3.37350-22.9-59.07192.0123.80.10.0-0.30.0XXPotentially PathogenicArg544 is located in the middle of an α-helix (res. Ala533-Val560). In the WT simulations, the guanidinium side chain of Arg544 forms a salt bridge with the carboxylate groups of Glu548 on the same α-helix, and with Glu651 and Glu656 on an opposing α-helix (res. Glu666-Asp644). In the variant simulations, the pyrrolidine side chain of Pro544 cannot form any of the salt bridges that Arg544 does in the WT, potentially weakening the tertiary structure assembly. Additionally, Pro544 lacks the amide group, and thus, unlike Arg544 in the WT, is unable to form a hydrogen bond with the carbonyl of Gln540. This disruption breaks the continuity of the secondary structure element, causing the α-helix to bend slightly in the variant simulations. These negative structural effects could be more pronounced during protein folding and are likely to be undermined in the MD simulations.
c.1172G>TG391V
(3D Viewer)
Likely BenignC2Likely Benign 16-33438077-G-T31.86e-6-6.642Likely Benign0.133Likely BenignLikely Benign0.595Likely Pathogenic4.23Destabilizing1.34.81Destabilizing4.52Destabilizing-0.11Likely Benign-0.98Neutral0.994Probably Damaging0.887Possibly Damaging1.32Pathogenic0.10Tolerated3.698-1-34.642.08228.6-69.00.00.8-0.50.3UncertainGly387 is located in the Gly-rich Ω loop (res. Pro364-Pro398) between two anti-parallel β sheet strands (res. Thr359-Pro364 and res. Ala399-Ile411). The Ω loop is assumed to directly interact with the membrane, and it is observed to move arbitrarily throughout the WT solvent simulations. This loop potentially plays a crucial role in the SynGAP-membrane complex association, stability, and dynamics. However, this aspect cannot be fully addressed through solvent simulations alone.Ω loops are known to play significant roles in protein functions that require flexibility, and thus hydrophobic residues like valine are rarely tolerated. Although no negative structural effects are visualized in the variant’s simulations, Val391 may exert drastic effects on the SynGAP-membrane complex dynamics and stability. Since the effects on the Gly-rich Ω loop dynamics can only be well studied through the SynGAP-membrane complex, no definite conclusions can be drawn.
c.1349C>AA450E
(3D Viewer)
Likely PathogenicGAPUncertain 1-16.578Likely Pathogenic0.989Likely PathogenicLikely Pathogenic0.653Likely Pathogenic3.86Destabilizing0.25.23Destabilizing4.55Destabilizing1.59Destabilizing-4.67Deleterious0.999Probably Damaging0.992Probably Damaging3.38Benign0.07Tolerated3.37320-1-5.358.04240.1-82.60.00.00.70.0XXPotentially PathogenicThe methyl group of Ala450, located in an α helix (res. Asn440-Thr458), packs against hydrophobic residues in the inter-helix space (e.g., Leu692). In the variant simulations, the carboxylate group of the Glu450 side chain rotates outward, away from the hydrophobic niche, where it does not form any lasting salt bridges or H-bonds. Although the residue swap does not negatively affect the protein structure based on the simulations, it is possible that the introduction of the negatively charged residue adversely affects the folding process or tertiary assembly.
c.1715G>CW572S
(3D Viewer)
Likely PathogenicGAPPathogenic 1-17.461Likely Pathogenic0.999Likely PathogenicLikely Pathogenic0.775Likely Pathogenic5.78Destabilizing0.23.37Destabilizing4.58Destabilizing1.79Destabilizing-12.74Deleterious1.000Probably Damaging1.000Probably Damaging-1.24Pathogenic0.01Affected3.3735-2-30.1-99.14235.176.60.00.0-0.40.1XPotentially PathogenicThe introduced residue Ser572, located in an α-helix (res. Arg563-Glu578), is considerably smaller than the tryptophan it replaced. The indole ring of the Trp572 side chain lies in a hydrophobic inter-helix space, where it makes extensive hydrophobic interactions with nearby residues such as Met470, Phe569, Leu588, and Ile589. In the variant simulations, all these favorable packing interactions are completely removed, as the introduced residue Ser572 is too hydrophilic or small to fill the hydrophobic niche occupied by the indole ring. Moreover, the hydroxyl group of Ser572 forms hydrogen bonds with the carbonyl groups of Glu567 and Val568 within the same α-helix, potentially lowering its integrity. Overall, the residue swap is highly likely to cause critical protein folding problems that are underestimated based on the effects seen in the variant simulations.
c.2068T>CS690P
(3D Viewer)
Likely PathogenicGAPUncertain 1-14.568Likely Pathogenic0.999Likely PathogenicLikely Pathogenic0.431Likely Benign4.84Destabilizing0.34.40Destabilizing4.62Destabilizing1.42Destabilizing-4.77Deleterious0.998Probably Damaging0.790Possibly Damaging3.44Benign0.01Affected3.42171-1-0.810.04207.515.10.10.0-0.10.2XXPotentially PathogenicThe hydroxyl side chain of Ser690, located in an α-helix (res. Leu696-Leu685), forms a hydrogen bond with the backbone carbonyl group of Ser410 in an anti-parallel β-sheet of the C2 domain (res. Ile411-Ala399). In the variant simulations, the pyrrolidine side chain of Pro690 cannot form hydrogen bonds with the C2 domain residue, resulting in the loss of this inter-domain connection. Additionally, prolines lack a free amide group necessary for hydrogen bonding with the carbonyl group of Gly686, introducing a slight bend in the α-helix and compromising its integrity.
c.1406C>AA469D
(3D Viewer)
Likely PathogenicGAPUncertain 1-14.643Likely Pathogenic0.999Likely PathogenicLikely Pathogenic0.738Likely Pathogenic5.09Destabilizing0.24.16Destabilizing4.63Destabilizing1.68Destabilizing-3.48Deleterious0.999Probably Damaging0.996Probably Damaging-1.34Pathogenic0.21Tolerated3.37340-2-5.344.01237.0-58.2-0.20.10.80.1XXPotentially PathogenicThe methyl group of Ala469, located in an α helix (res. Ala461–Phe476), interacts with hydrophobic residues (e.g., Trp572, Leu588, Met470) in an inter-helix space formed by two other α helices (res. Glu582–Ser604, res. Arg563–Gly580). In the variant simulations, Asp469 introduces a negatively charged and bulky side chain into the hydrophobic niche. Consequently, the side chain of Asp469 rotates outward, allowing the carboxylate group to form a salt bridge with the guanidinium group of Arg575 on the protein surface. This interaction affects the continuity of the parent α helix (Ala461–Phe476). Due to the importance of hydrophobic packing, the structural effects could be more pronounced during actual protein folding.
c.980T>CL327P
(3D Viewer)
Likely PathogenicC2Pathogenic 3-16.602Likely Pathogenic0.999Likely PathogenicLikely Pathogenic0.658Likely Pathogenic5.38Destabilizing0.14.00Destabilizing4.69Destabilizing2.62Destabilizing-5.97Deleterious1.000Probably Damaging0.999Probably Damaging1.52Pathogenic0.01Affected3.3823-3-3-5.4-16.04221.769.40.10.00.60.1XPotentially PathogenicThe backbone amide group of Leu327, located in the middle of an anti-parallel β sheet strand (res. Ala322-Asp330), forms a hydrogen bond with the carbonyl group of Gly344 on a neighboring β strand (res. Lys336-Pro349) in the WT simulations. In contrast, in the variant simulations, the introduction of Pro327 destabilizes the hydrogen bonding between the two anti-parallel β strands because proline lacks the backbone amide group altogether. Additionally, in the WT simulations, the iso-butyl side chain of Leu327 packs against multiple hydrophobic residues (e.g., Leu274, V400, Val343), whereas the less bulky cyclic five-membered pyrrolidine ring of Pro327 cannot fill the same space as effectively. Thus, although no large-scale unfolding is observed during the variant simulations, the residue swap is likely to cause severe problems for the correct C2 domain folding, which could also affect the SynGAP-membrane association.10.1016/j.ajhg.2020.11.011
c.1505G>AG502D
(3D Viewer)
Likely PathogenicGAPUncertain 1-14.796Likely Pathogenic0.994Likely PathogenicLikely Pathogenic0.915Likely Pathogenic3.79Destabilizing0.95.69Destabilizing4.74Destabilizing1.38Destabilizing-6.80Deleterious0.999Probably Damaging0.977Probably Damaging-1.66Pathogenic0.00Affected3.37351-1-3.158.04224.2-80.0-0.80.70.60.3XXXPotentially PathogenicGly502 is located in a hinge in the middle of an α-helix (res. Leu489-Glu519). In the WT, Gly502 acts as an α-helix breaker due to its lack of a side chain, facilitating a bend in the middle of the α-helix. In the variant simulations, the carboxylate group of Asp502 forms hydrogen bonds with neighboring residues (e.g., Ser677, Lys504), disrupting the hinge. Additionally, Asp502 struggles to fit into the α-helix hinge and cannot generate a similar bend as Gly502, which would drastically affect the secondary structure during folding. Thus, the deleterious effect seen in the simulations is likely an underestimate of the impact of the residue swap on the protein structure during protein folding.
c.1354G>TV452F
(3D Viewer)
Likely PathogenicGAPUncertain 1-14.769Likely Pathogenic0.975Likely PathogenicLikely Pathogenic0.511Likely Pathogenic9.21Destabilizing0.10.37Likely Benign4.79Destabilizing0.61Ambiguous-4.94Deleterious0.999Probably Damaging0.993Probably Damaging3.29Benign0.00Affected3.3734-1-1-1.448.04249.4-35.70.00.00.40.1XPotentially PathogenicThe iso-propyl side chain of Val452, located in the middle of an α helix (res. Val441-Ser457), packs against hydrophobic residues in the inter-helix space at the intersection of three α helices (e.g., Leu500, His453, Leu465). In the variant simulations, the larger side chain of Phe452 cannot pack against the opposing α helix (res. Leu489-Glu519) as efficiently as valine. Due to space restrictions, the phenol ring adjusts to make room by rotating slightly sideways in the inter-helix space. Besides this small and local shift, no large-scale effects on the protein structure are seen based on the simulations. However, the size difference between the swapped residues could affect the protein folding process.
c.922T>CW308R
(3D Viewer)
Likely PathogenicC2Pathogenic 1-12.264Likely Pathogenic1.000Likely PathogenicLikely Pathogenic0.868Likely Pathogenic5.40Destabilizing0.54.27Destabilizing4.84Destabilizing1.88Destabilizing-12.87Deleterious1.000Probably Damaging0.999Probably Damaging0.48Pathogenic0.00Affected3.38192-3-3.6-30.03290.4-26.7-0.10.10.00.2XXXPotentially PathogenicThe indole ring of Trp308, located in an anti-parallel β sheet strand (res. Thr305-Asn315), packs against multiple hydrophobic residues (e.g., Ile268, Val306, Cys282). The indole group of Trp308 also hydrogen bonds with the backbone atoms of the C2 domain residues forming the anti-parallel β sheet (e.g., Tyr280, Thr294). The guanidinium group of Arg308 is comparably sized to the tryptophan it replaced; however, it is also positively charged.In the variant simulations, the charged side chain remains buried deep in the hydrophobic part of the C2 domain, where it forms new hydrogen bonds with the backbone carbonyl atoms of surrounding residues (e.g., Val306, Ile268). However, the residue swap is likely to disrupt the hydrophobic packing during folding. At a minimum, the residue swap could affect the C2 domain stability and membrane association.
c.1513T>GY505D
(3D Viewer)
Likely PathogenicGAPLikely Pathogenic 1-14.078Likely Pathogenic0.993Likely PathogenicLikely Pathogenic0.718Likely Pathogenic4.98Destabilizing0.14.72Destabilizing4.85Destabilizing2.49Destabilizing-9.95Deleterious1.000Probably Damaging1.000Probably Damaging2.60Benign0.00Affected3.3735-3-4-2.2-48.09
c.1621G>CA541P
(3D Viewer)
Likely PathogenicGAPUncertain 1-14.733Likely Pathogenic0.996Likely PathogenicLikely Pathogenic0.594Likely Pathogenic2.47Destabilizing0.37.26Destabilizing4.87Destabilizing0.86Ambiguous-3.16Deleterious1.000Probably Damaging0.998Probably Damaging-1.34Pathogenic0.07Tolerated3.37351-1-3.426.04170.4-11.20.10.00.10.0XPotentially PathogenicAla541 is located on the outer surface of an α-helix (res. Ala533-Val560). The methyl group of Ala541 is on the surface and does not form any interactions. Proline lacks a free backbone amide group, and thus, Pro541 is unable to form a hydrogen bond with the carbonyl group of Ala537 in the variant simulations. Consequently, Pro541 disrupts the continuity of the secondary structure element, causing the α-helix to bend slightly in the variant simulations.
c.924G>CW308C
(3D Viewer)
Likely PathogenicC2Pathogenic/Likely path. 2-12.791Likely Pathogenic1.000Likely PathogenicLikely Pathogenic0.738Likely Pathogenic5.56Destabilizing0.34.38Destabilizing4.97Destabilizing1.26Destabilizing-11.95Deleterious1.000Probably Damaging0.999Probably Damaging0.48Pathogenic0.00Affected3.3819-8-23.4-83.07230.860.5-0.30.1-0.40.4XPotentially PathogenicThe indole ring of Trp308, located in an anti-parallel β sheet strand (res. Thr305-Asn315), packs against multiple hydrophobic residues (e.g., Ile268, Val306, Cys282). The indole group of Trp308 also hydrogen bonds with the backbone atoms of the C2 domain residues forming the anti-parallel β sheet (e.g., Tyr280, Thr294). The introduced Cys308 is smaller than the tryptophan it replaced. The thiol group of the Cys308 side chain is well-suited for the inner hydrophobic part of the C2 domain. Although the negative effects are essentially missing from the simulations, the side chain size difference between the residues is likely to disrupt the hydrophobic packing during folding. At a minimum, the residue swap could affect the C2 domain stability and membrane association.
c.1481T>GI494R
(3D Viewer)
Likely PathogenicGAPLikely Pathogenic 1-15.758Likely Pathogenic0.995Likely PathogenicLikely Pathogenic0.911Likely Pathogenic6.71Destabilizing0.33.40Destabilizing5.06Destabilizing2.19Destabilizing-6.43Deleterious0.999Probably Damaging0.957Probably Damaging-1.41Pathogenic0.00Affected3.3735-2-3-9.043.03273.9-59.80.00.00.00.1XXXXPotentially PathogenicThe sec-butyl side chain of Ile494, located in an α-helix (res. Leu489-Glu519), packs against hydrophobic residues (e.g., Phe484, Leu465, Trp572, Ala493, Met468) in an inter-helix space (res. Leu489-Glu519 and res. Ala461-Phe476). In the variant simulations, the bulkier and positively charged residue, Arg494, weakens the integrity of the opposing helix. Additionally, the bulkier Arg494 stacks with Phe484, causing the α-helices to move farther apart to accommodate it. This mutation could have substantial negative effects due to the fundamental role of hydrophobic packing, which is disrupted by Arg494 during protein folding.
c.1147G>TG383W
(3D Viewer)
C2Uncertain 16-33438052-G-T16.22e-7-10.161Likely Pathogenic0.439AmbiguousLikely Benign0.469Likely Benign5.81Destabilizing3.64.44Destabilizing5.13Destabilizing0.08Likely Benign-1.01Neutral0.959Probably Damaging0.704Possibly Damaging4.09Benign0.00Affected4.327-2-7-0.5129.16
c.1516C>TL506F
(3D Viewer)
Likely PathogenicGAPUncertain 1-11.262Likely Pathogenic0.883Likely PathogenicAmbiguous0.464Likely Benign4.92Destabilizing0.85.76Destabilizing5.34Destabilizing0.91Ambiguous-3.98Deleterious0.999Probably Damaging0.997Probably Damaging1.62Pathogenic0.01Affected3.373502-1.034.02
c.703T>CS235P
(3D Viewer)
Likely PathogenicPHLikely Pathogenic 1-14.857Likely Pathogenic0.998Likely PathogenicLikely Pathogenic0.870Likely Pathogenic4.02Destabilizing0.16.91Destabilizing5.47Destabilizing1.23Destabilizing-4.24Deleterious0.917Possibly Damaging0.446Benign5.47Benign0.01Affected3.40141-1-0.810.04201.517.00.10.0-0.60.0XPotentially PathogenicIn the WT, the hydroxyl group of Ser235, located in a β-α loop between an anti-parallel β sheet strand (res. Gly227-Phe231) and an α helix (residues Ala236-Val250), forms hydrogen bonds with the GAP domain loop residue Glu680 and with the backbone amide groups of Ala237 and Glu238 from the α helix. In the variant simulations, the pyrrolidine ring of Pro235 cannot stabilize the α helix end or maintain tertiary bonding interactions between the PH and GAP domains via hydrogen bonding as effectively as serine.
c.1714T>CW572R
(3D Viewer)
Likely PathogenicGAPNot provided1-17.511Likely Pathogenic1.000Likely PathogenicLikely Pathogenic0.894Likely Pathogenic4.84Destabilizing0.16.19Destabilizing5.52Destabilizing1.79Destabilizing-12.81Deleterious-1.25Pathogenic0.00Affected3.37352-3-3.6-30.03312.6-37.60.00.0-1.00.0XXPotentially PathogenicThe indole ring of Trp572, located in an α-helix (res. Arg563-Glu578), lies in a hydrophobic inter-helix space, where it makes extensive hydrophobic interactions with nearby residues such as Met470, Phe569, Leu588, and Ile589. The guanidinium group of Arg572 is similarly sized to the tryptophan it replaced; however, it is also positively charged. In the variant simulations, Arg572 forms hydrogen bonds with other residues in the inter-helix space, such as Ser592 and the backbone carbonyl atom of Leu465. Additionally, Arg572 hydrophobically packs its carbon chain with surrounding residues such as Phe569 and Ile589.However, the introduced residue arginine is too hydrophilic and charged for the hydrophobic space, disrupting the hydrophobic packing of the inter-helix space. Indeed, in the second simulation, Arg572 successfully escapes the hydrophobic niche completely, causing the whole protein to partially unfold.Overall, the residue swap is highly likely to cause critical protein folding problems, as evidenced by the effects seen in the variant simulations.
c.1259T>CF420S
(3D Viewer)
Likely PathogenicGAPLikely Pathogenic 1-13.231Likely Pathogenic1.000Likely PathogenicLikely Pathogenic0.544Likely Pathogenic5.34Destabilizing0.15.73Destabilizing5.54Destabilizing2.14Destabilizing-7.43Deleterious0.998Probably Damaging0.938Probably Damaging3.09Benign0.00Affected3.3729-3-2-3.6-60.10213.357.80.00.0-0.40.1XPotentially PathogenicIn the WT, the phenyl ring of the Phe420 side chain, located on an α helix (res. Met414-Glu436), packs against hydrophobic residues in the interhelix area of the GAP domain (e.g., Leu689, Leu714, Leu717, Leu718). Although no large-scale adverse effects are seen in the variant simulations, the polar hydroxyl group of Ser420 is not suitable for the hydrophobic inter-helix space. Thus, the residue swap could affect protein folding. In theory, the introduced hydroxyl group could also lower the α helix integrity by H-bonding with the backbone atoms of neighboring residues in the same α helix. However, no such effect is seen in the variant simulations.
c.1706T>CF569S
(3D Viewer)
Likely PathogenicGAPLikely Pathogenic 2-13.384Likely Pathogenic0.999Likely PathogenicLikely Pathogenic0.916Likely Pathogenic5.70Destabilizing0.15.38Destabilizing5.54Destabilizing2.45Destabilizing-7.97Deleterious1.000Probably Damaging1.000Probably Damaging-1.32Pathogenic0.00Affected3.3734-3-2-3.6-60.10213.767.9-0.10.0-1.00.1XPotentially PathogenicPhe569 is located on an α-helix (res. Arg563-Glu578). In the WT simulations, the phenyl side chain of Phe569 packs with hydrophobic residues such as Trp572, Leu565, Ile589, Ile667, and Phe561, originating from three different α-helices (res. Ala533-Val560, res. Arg563-Glu578, and res. Ser641-Glu666). In the variant simulations, the acceptor/donor hydroxyl group of Ser569 forms hydrogen bonds with the carbonyl groups of Glu567 and Lys566 on the same α-helix, which could affect the α-helix integrity, although this is not observed in the simulations. While the simulations do not show large-scale effects, the residue swap could have a substantial impact on the protein structure due to the fundamental role of hydrophobic packing during protein folding.
c.1771G>CA591P
(3D Viewer)
Likely PathogenicGAPUncertain 1-14.479Likely Pathogenic0.991Likely PathogenicLikely Pathogenic0.404Likely Benign3.78Destabilizing0.37.29Destabilizing5.54Destabilizing1.45Destabilizing-4.41Deleterious0.995Probably Damaging0.853Possibly Damaging3.35Benign0.01Affected3.37351-1-3.426.04191.5-10.10.20.10.40.1XPotentially PathogenicThe methyl group of the Ala591 side chain, located in the middle of an α helix (res. Glu582-Met603), packs against hydrophobic residues (e.g., Ile483, Phe484) of an opposing partially helical loop (res. Phe476-Asn487).In the variant simulations, Pro591 lacks a free backbone amide group and, therefore, cannot form a hydrogen bond with the backbone carbonyl of Arg587 as Ala591 does in the WT. This notably weakens the α helix integrity and compromises the continuity of the helix. In reality, the effect on the structure during protein folding could be more severe.
c.2071A>CT691P
(3D Viewer)
Likely PathogenicGAPLikely Pathogenic 1-13.801Likely Pathogenic0.905Likely PathogenicAmbiguous0.214Likely Benign5.04Destabilizing0.46.09Destabilizing5.57Destabilizing1.27Destabilizing-3.43Deleterious1.000Probably Damaging0.952Probably Damaging3.43Benign0.06Tolerated3.43140-1-0.9-3.99188.933.00.10.0-0.60.0XXPotentially PathogenicThe hydroxyl side chain of Thr691, located in an α-helix (res. Leu696-Leu685), can form hydrogen bonds with the backbone carbonyl and the side chain guanidinium group of Arg687. This interaction facilitates the simultaneous formation of salt bridges between Arg687 and Glu688 on the same α-helix. Additionally, Thr691 occasionally interacts with the thioether side chain of Met409 in an anti-parallel β-sheet of the C2 domain (res. Ile411-Ala399), although this interaction is not consistently maintained throughout the WT simulations. In the variant simulations, the pyrrolidine side chain of Pro691 lacks hydrogen bond donors, making a similar setup impossible. Moreover, proline lacks a free amide group necessary for hydrogen bonding with the carbonyl group of Arg687, introducing a slight bend in the α-helix and compromising its integrity.
c.1142G>TG381V
(3D Viewer)
Likely BenignC2Uncertain 16-33438047-G-T21.25e-6-5.967Likely Benign0.146Likely BenignLikely Benign0.618Likely Pathogenic7.16Destabilizing1.04.10Destabilizing5.63Destabilizing-0.32Likely Benign-0.95Neutral0.386Benign0.157Benign1.32Pathogenic0.10Tolerated4.329-1-34.642.08214.6-68.80.30.7-0.50.3UncertainGly381 is located in the Gly-rich Ω loop (res. Pro364-Pro398) between two anti-parallel β sheet strands (res. Thr359-Pro364, res. Ala399-Ile411). Because the Ω loop is assumed to directly interact with the membrane, it moves arbitrarily throughout the WT solvent simulations. The Ω loop potentially plays a crucial role in the SynGAP-membrane complex association, stability, and dynamics. However, this aspect cannot be fully addressed through solvent simulations alone.Ω loops are known to play major roles in protein functions that require flexibility, and thus hydrophobic residues like valine are rarely tolerated. Although no negative structural effects are observed in the variant simulations, Val381 may exert drastic effects on the SynGAP-membrane complex dynamics and stability. However, since the effects on Gly-rich Ω loop dynamics can only be well studied through the SynGAP-membrane complex, no definite conclusions can be drawn.
c.1726T>CC576R
(3D Viewer)
Likely PathogenicGAPConflicting 2-14.886Likely Pathogenic1.000Likely PathogenicLikely Pathogenic0.579Likely Pathogenic7.20Destabilizing1.04.09Destabilizing5.65Destabilizing1.64Destabilizing-10.88Deleterious0.999Probably Damaging0.996Probably Damaging3.38Benign0.00Affected3.3735-3-4-7.053.05
c.1160G>TG387V
(3D Viewer)
Likely BenignC2Uncertain 16-33438065-G-T221.37e-5-6.199Likely Benign0.153Likely BenignLikely Benign0.390Likely Benign5.13Destabilizing1.86.44Destabilizing5.79Destabilizing-0.33Likely Benign-0.54Neutral0.069Benign0.077Benign1.32Pathogenic0.01Affected4.323-1-34.642.08207.7-68.4-0.70.8-0.50.1UncertainGly387 is located in the Gly-rich Ω loop (res. Pro364-Pro398) between two anti-parallel β sheet strands (res. Thr359-Pro364 and res. Ala399-Ile411). The Ω loop is assumed to directly interact with the membrane, and it is observed to move arbitrarily throughout the WT solvent simulations. This loop potentially plays a crucial role in the SynGAP-membrane complex association, stability, and dynamics. However, this aspect cannot be fully addressed through solvent simulations alone.Ω loops are known to play significant roles in protein functions that require flexibility, and thus hydrophobic residues like valine are rarely tolerated. Although no negative structural effects are visualized in the variant’s simulations, Val387 may exert drastic effects on the SynGAP-membrane complex dynamics and stability. Since the effects on the Gly-rich Ω loop dynamics can only be well studied through the SynGAP-membrane complex, no definite conclusions can be drawn.
c.851T>CL284PLikely PathogenicC2Likely Pathogenic1-15.588Likely Pathogenic1.000Likely PathogenicLikely Pathogenic0.794Likely Pathogenic5.83Destabilizing0.25.81Destabilizing5.82Destabilizing1.89Destabilizing-6.17Deleterious1.000Probably Damaging0.999Probably Damaging1.64Pathogenic0.00Affected-3-3-5.4-16.04
c.968T>CL323P
(3D Viewer)
Likely PathogenicC2Uncertain 1-12.507Likely Pathogenic0.998Likely PathogenicLikely Pathogenic0.762Likely Pathogenic3.39Destabilizing0.68.46Destabilizing5.93Destabilizing2.20Destabilizing-4.80Deleterious0.999Probably Damaging0.977Probably Damaging0.59Pathogenic0.01Affected4.29398-3-3-5.4-16.04201.968.20.00.10.60.3XPotentially PathogenicThe iso-butyl side chain of Leu323, located at the beginning of an anti-parallel β sheet strand (res. Ala322-Asp330), packs against multiple hydrophobic leucine residues (e.g., Leu264, Leu266, Leu284, Leu286). In contrast, in the variant simulations, the less bulky cyclic five-membered pyrrolidine ring of Pro323 cannot fill the hydrophobic space as effectively as the branched hydrocarbon side chain of leucine. Notably, the backbone amide group of Leu323 forms a hydrogen bond with the backbone carbonyl group of Cys285. Pro323 cannot form this bond due to the absence of the backbone amide group, resulting in partial unfolding of the anti-parallel β sheet end in the variant simulations.
c.743G>CR248P
(3D Viewer)
Likely PathogenicPHLikely Pathogenic 1-10.751Likely Pathogenic1.000Likely PathogenicLikely Pathogenic0.848Likely Pathogenic3.09Destabilizing0.68.87Destabilizing5.98Destabilizing1.21Destabilizing-5.97Deleterious0.998Probably Damaging0.878Possibly Damaging5.64Benign0.00Affected3.41140-22.9-59.07223.8126.60.00.0-0.20.1XXPotentially PathogenicThe guanidinium group of Arg248, located on an α helix (residues Ala236-Val250), forms two very stable salt bridges with Asp255 (from a short α helical section, res. Lys254-Asn256) and Glu244 (from a nearby loop) in the WT simulations. In the variant simulations, the pyrrolidine side chain of Pro248 cannot form any salt bridges, which could negatively affect the tertiary structure assembly of the PH domain. Additionally, Pro248 lacks a free amide group needed for hydrogen bonding with the backbone carbonyl group of Asn245, disrupting the continuity of the α helix.
c.1367A>CQ456P
(3D Viewer)
Likely PathogenicGAPUncertain 1-15.250Likely Pathogenic0.993Likely PathogenicLikely Pathogenic0.469Likely Benign3.68Destabilizing0.28.43Destabilizing6.06Destabilizing0.82Ambiguous-5.66Deleterious1.000Probably Damaging0.999Probably Damaging3.34Benign0.07Tolerated3.3734-101.9-31.01
c.791T>CL264P
(3D Viewer)
Likely PathogenicC2Uncertain 1-12.285Likely Pathogenic1.000Likely PathogenicLikely Pathogenic0.767Likely Pathogenic5.73Destabilizing0.36.57Destabilizing6.15Destabilizing2.65Destabilizing-6.43Deleterious1.000Probably Damaging0.999Probably Damaging0.49Pathogenic0.00Affected-3-3-5.4-16.04
c.1802C>AA601E
(3D Viewer)
Likely PathogenicGAPConflicting 2-16.752Likely Pathogenic0.992Likely PathogenicLikely Pathogenic0.588Likely Pathogenic6.68Destabilizing0.85.76Destabilizing6.22Destabilizing1.24Destabilizing-4.98Deleterious1.000Probably Damaging0.999Probably Damaging2.54Benign0.00Affected3.37350-1-5.358.04240.0-82.30.00.00.70.1XXXPotentially PathogenicThe methyl side chain of Ala601, located on an α helix (res. Glu582-Met603), packs hydrophobically against other hydrophobic residues in the inter-helix space (e.g., Phe597, Leu598, Leu506, Phe608).In the variant simulations, the carboxylate group of Glu601 faces the inter-helix space and is forced to shift slightly away from the hydrophobic niche. Additionally, in two of the simulations, Glu601 forms a salt bridge with Arg499, causing the otherwise stable salt bridge between Arg499 and Glu496 at the outer surface of an α helix (res. Leu489-Glu519) to break due to the residue swap.These effects suggest that the protein folding process could be seriously affected. Moreover, due to its location at the GAP-Ras interface, it could also impact the complex formation with the GTPase.
c.1558T>CS520P
(3D Viewer)
Likely PathogenicGAPUncertain 1-12.707Likely Pathogenic0.999Likely PathogenicLikely Pathogenic0.855Likely Pathogenic3.72Destabilizing0.88.86Destabilizing6.29Destabilizing0.83Ambiguous-4.57Deleterious0.997Probably Damaging0.986Probably Damaging-1.32Pathogenic0.01Affected1-1-0.810.04
c.878G>CR293P
(3D Viewer)
Likely PathogenicC2Likely Pathogenic 1-16.275Likely Pathogenic0.999Likely PathogenicLikely Pathogenic0.497Likely Benign3.62Destabilizing0.49.06Destabilizing6.34Destabilizing0.47Likely Benign-6.43Deleterious1.000Probably Damaging0.999Probably Damaging1.45Pathogenic0.01Affected3.38230-22.9-59.07202.3132.00.10.00.10.1XXXPotentially PathogenicThe guanidinium group of the Arg293 side chain, located in an anti-parallel β sheet strand (res. Met289-Pro298), packs against the phenol ring of the Tyr281 side chain or forms a salt bridge with the carboxylate group of Glu283 on the outer side of the C2 domain. In the WT simulations, the positively charged side chain of arginine remains outside the hydrophobic C2 domain, resulting in a twist in the β strand. The backbone amide bond of Arg293 potentially maintains this twist by forming a hydrogen bond with the carbonyl group of His210 or the hydroxyl group of Ser211 in the anti-parallel β sheet.Although this twist is also maintained in the variant simulations, replacing the positively charged residue with proline, which lacks the backbone amide group altogether, causes the β strand to unfold. Because Arg293 is positioned at the C2 and PH domain interface, the residue swap could significantly impact the tertiary structure assembly. Notably, Arg293 is located at the SynGAP-Ras interface, and its role in complex formation cannot be fully understood through solvent-only simulations.
c.1652T>CL551P
(3D Viewer)
Likely PathogenicGAPLikely Pathogenic 1-14.620Likely Pathogenic0.999Likely PathogenicLikely Pathogenic0.953Likely Pathogenic6.66Destabilizing0.16.58Destabilizing6.62Destabilizing2.66Destabilizing-4.70Deleterious1.000Probably Damaging1.000Probably Damaging-1.60Pathogenic0.01Affected3.3735-3-3-5.4-16.04208.660.90.10.0-0.30.0XPotentially PathogenicL551 is located on an α-helix (res. Ala533-Val560). The iso-butyl side chain of Leu551 hydrophobically packs with nearby hydrophobic residues such as Cys547, Phe652, Leu633, and Ile630 in the inter-helix space. In the variant simulations, the pyrrolidine side chain of Pro551 is not as optimal as leucine for hydrophobic packing with the nearby residues. Moreover, Pro551 lacks the amide group, and thus, it cannot form a hydrogen bond with the backbone carbonyl group of Cys547, which disrupts the continuity of the secondary structure element.
c.1639T>CC547R
(3D Viewer)
Likely PathogenicGAPUncertain 1-16.967Likely Pathogenic0.999Likely PathogenicLikely Pathogenic0.900Likely Pathogenic7.76Destabilizing0.85.83Destabilizing6.80Destabilizing1.69Destabilizing-11.60Deleterious1.000Probably Damaging0.998Probably Damaging-1.33Pathogenic0.02Affected3.3735-4-3-7.053.05267.4-90.30.00.0-0.10.1XXXXPotentially PathogenicCys547 is located in an α-helix (res. Ala533-Val560). The thiol side chain of Cys is situated in a hydrophobic inter-helix space, where it packs hydrophobically with other residues such as Ile626, Leu551, and Phe652. Additionally, the thiol side chain of Cys547 weakly hydrogen bonds with the carbonyl group of Leu543 in the same α-helix. In the variant simulations, the bulkier, positively charged guanidinium group of Arg547 must rotate out of the hydrophobic space. Consequently, it forms ionic interactions with the carboxylate groups of Glu548 in the same helix and Glu656 in the neighboring α-helix (res. Glu666-Asp644). This causes the two helices to slightly separate, significantly affecting the secondary structure integrity of the latter helix. These negative structural effects could be more pronounced during protein folding and are likely to be undermined in the MD simulations.
c.1714T>GW572G
(3D Viewer)
Likely PathogenicGAPUncertain 1-17.692Likely Pathogenic0.997Likely PathogenicLikely Pathogenic0.900Likely Pathogenic6.57Destabilizing0.27.57Destabilizing7.07Destabilizing1.83Destabilizing-11.98Deleterious1.000Probably Damaging1.000Probably Damaging-1.24Pathogenic0.00Affected3.3735-7-20.5-129.16195.2127.90.00.0-1.00.0XPotentially PathogenicThe introduced residue Gly572, located in an α-helix (res. Arg563-Glu578), is considerably smaller than the tryptophan it replaced. The indole ring of the Trp572 side chain lies in a hydrophobic inter-helix space, where it makes extensive hydrophobic interactions with nearby residues such as Met470, Phe569, Leu588, and Ile589. In the variant simulations, all these favorable packing interactions are completely removed, as the introduced residue Gly572 essentially lacks a side chain altogether. Although not observed in the simulations, the residue swap could also weaken the integrity of the helix (res. Arg563-Glu578), as glycine is known as an “α-helix breaker.” Overall, the residue swap is highly likely to cause critical protein folding problems that are underestimated based on the effects seen in the variant simulations.
c.1640G>AC547Y
(3D Viewer)
Likely PathogenicGAPPathogenic 1-15.871Likely Pathogenic0.999Likely PathogenicLikely Pathogenic0.874Likely Pathogenic8.53Destabilizing1.86.20Destabilizing7.37Destabilizing0.62Ambiguous-10.57Deleterious1.000Probably Damaging0.998Probably Damaging-1.33Pathogenic0.06Tolerated3.37350-2-3.860.04280.1-54.80.00.00.00.0XXXPotentially PathogenicCys547 is located in an α-helix (res. Ala533-Val560). The thiol side chain of Cys547 is situated in a hydrophobic inter-helix space, where it packs hydrophobically with other residues such as Ile626, Leu551, and Phe652. Additionally, the thiol side chain of Cys weakly hydrogen bonds with the carbonyl group of Leu543 in the same α-helix. In the variant simulations, the bulkier phenol ring of Tyr547, with its polar hydroxyl group, is less suited for the hydrophobic space. Consequently, it moves outside and forms a hydrogen bond with the carbonyl group of Phe652 in the neighboring α-helix (res. Glu666-Asp644). This causes the two helices to slightly separate, negatively affecting the secondary structure integrity of the latter helix. These negative structural effects could be more pronounced during protein folding and are likely to be undermined in the MD simulations.
c.1058T>CL353P
(3D Viewer)
Likely PathogenicC2Uncertain 1-7.913In-Between0.936Likely PathogenicAmbiguous0.464Likely Benign4.63Destabilizing0.110.19Destabilizing7.41Destabilizing2.17Destabilizing-3.70Deleterious0.947Possibly Damaging0.454Possibly Damaging1.29Pathogenic0.02Affected3.3725-3-3-5.4-16.04
c.1499T>CL500P
(3D Viewer)
Likely PathogenicGAPPathogenic 1-15.898Likely Pathogenic0.996Likely PathogenicLikely Pathogenic0.894Likely Pathogenic5.91Destabilizing0.38.90Destabilizing7.41Destabilizing1.92Destabilizing-6.96Deleterious1.000Probably Damaging1.000Probably Damaging-1.37Pathogenic0.01Affected3.3735-3-3-5.4-16.04
c.1030G>AG344S
(3D Viewer)
Likely PathogenicC2Pathogenic 5-11.254Likely Pathogenic0.986Likely PathogenicLikely Pathogenic0.790Likely Pathogenic9.02Destabilizing0.76.08Destabilizing7.55Destabilizing0.98Ambiguous-5.28Deleterious1.000Probably Damaging1.000Probably Damaging-0.45Pathogenic0.04Affected3.372510-0.430.03217.3-51.70.00.10.20.1XXPotentially PathogenicBecause Gly344 lacks a proper side chain, it allows the anti-parallel β sheet strand (res. Gly341-Pro349) to have a slight twist. Within a β strand, side chains normally alternate between outward and inward positions, but glycine is an exception as it allows the alternating pattern to skip a residue. Introducing serine or any other residue with a side chain at position 344 prevents this unique skip in the alternating pattern, causing structural strain or likely preventing correct folding altogether. Additionally, Tyr342 shields Gly344 from the solvent, contributing to twist formation in the β sheet and stabilizing the β-strand.In the variant simulations, the side chain of Ser344 assumes the inward position. However, the hydrophobic niche formed by multiple C2 domain residues (e.g., Val365, Val343, Leu327) is not accommodating for its hydroxyl group. The outward position, not seen in the simulations, would be equally disadvantageous due to the presence of hydrophobic residues on that side as well (e.g., Leu345, Tyr342). Serine is also not well-suited for twist formation, as it tends to suppress twisting and bending in β sheets. At this position, the hydroxyl group of Ser344 could also form hydrogen bonds with the backbone atoms of the Gly-rich Ω loop in the C2 domain (e.g., Thr366, Leu367, Gly378; res. Pro364-Pro398), potentially adversely affecting membrane-loop dynamics and ultimately compromising the stability of the SynGAP-membrane association.
c.1814C>GP605R
(3D Viewer)
Likely PathogenicGAPUncertain 1-13.745Likely Pathogenic0.996Likely PathogenicLikely Pathogenic0.845Likely Pathogenic8.71Destabilizing2.56.46Destabilizing7.59Destabilizing0.92Ambiguous-8.95Deleterious1.000Probably Damaging1.000Probably Damaging0.69Pathogenic0.00Affected3.37350-2-2.959.07281.7-118.1-0.20.00.50.1XXXXPotentially PathogenicPro605 is located in a short turn between an α helix (res. Glu582-Met603) and a short α helical section (res. Ser606-Phe608). The pyrrolidine side chain of Pro605 packs hydrophobically with nearby hydrophobic residues (e.g., Ile514, Leu623, Leu610) in the inter-helix space. Additionally, proline lacks a free backbone amide group, which breaks the α helix and facilitates the turn in the WT structure.In the variant simulations, the guanidinium side chain of Arg605 is bulkier than proline, and its positively charged guanidinium group faces mostly hydrophobic residues (e.g., Ile514, Leu623, Leu610). As a result, it needs to rotate away from the hydrophobic niche. The residue swap could have a more profound effect on the actual folding process, for example, by preventing the bending at the α helix end.Moreover, due to its location at the GAP-Ras interface, the residue swap could affect the GAP-Ras association.
c.1352T>CL451P
(3D Viewer)
Likely PathogenicGAPLikely Pathogenic 1-14.549Likely Pathogenic1.000Likely PathogenicLikely Pathogenic0.750Likely Pathogenic6.92Destabilizing0.28.57Destabilizing7.75Destabilizing2.58Destabilizing-6.81Deleterious1.000Probably Damaging1.000Probably Damaging2.43Pathogenic0.00Affected3.3734-3-3-5.4-16.04
c.1517T>CL506P
(3D Viewer)
Likely PathogenicGAPLikely Pathogenic1-12.088Likely Pathogenic0.998Likely PathogenicLikely Pathogenic0.737Likely Pathogenic5.48Destabilizing0.710.19Destabilizing7.84Destabilizing2.50Destabilizing-6.96Deleterious1.000Probably Damaging1.000Probably Damaging1.55Pathogenic0.00Affected3.3735-3-3-5.4-16.04182.664.90.10.00.20.1XPotentially PathogenicLeu506 is located in the middle of an α-helix (res. Gly502-Tyr518) within the inter-helix space of two helices (res. Gly502-Tyr518 and res. Glu582-Met603). In the WT simulations, the iso-butyl side chain of Leu506 hydrophobically packs with residues in the inter-helix space (e.g., Ile510, Phe597, Leu598, Ala601). In the variant simulations, the cyclic five-membered pyrrolidine ring of Pro506 is not as optimal as Leu506 for hydrophobic packing with nearby residues. Additionally, Pro506 cannot maintain the hydrogen bond with the backbone oxygen of Gly502 as Leu506 does in the WT, which disrupts the secondary structure element.
c.1778T>CL593P
(3D Viewer)
Likely PathogenicGAPUncertain 1-13.961Likely Pathogenic0.999Likely PathogenicLikely Pathogenic0.777Likely Pathogenic5.75Destabilizing0.910.77Destabilizing8.26Destabilizing2.43Destabilizing-6.77Deleterious1.000Probably Damaging1.000Probably Damaging2.77Benign0.00Affected-3-3-5.4-16.04
c.1898T>CL633P
(3D Viewer)
Likely PathogenicGAPPathogenic/Likely path. 2-15.669Likely Pathogenic1.000Likely PathogenicLikely Pathogenic0.693Likely Pathogenic6.60Destabilizing0.210.15Destabilizing8.38Destabilizing2.42Destabilizing-6.97Deleterious1.000Probably Damaging1.000Probably Damaging2.70Benign0.00Affected3.3734-3-3-5.4-16.04193.265.10.00.00.10.0XPotentially PathogenicThe iso-butyl side chain of Leu633, located in the middle of an α helix (res. Glu617-Asn635), packs hydrophobically with nearby residues (e.g., Leu653, Val629, Leu551) in the WT simulations.In the variant simulations, the pyrrolidine side chain of Pro633 is not as optimal for hydrophobic packing as Leu633 in the WT. Additionally, proline lacks a free backbone amide group, so Pro633 cannot form a hydrogen bond with the backbone carbonyl group of Val629, which disrupts the continuity of the secondary structure element.
c.2087T>CL696P
(3D Viewer)
Likely PathogenicGAPLikely Pathogenic 1-16.926Likely Pathogenic1.000Likely PathogenicLikely Pathogenic0.678Likely Pathogenic6.66Destabilizing0.210.84Destabilizing8.75Destabilizing2.13Destabilizing-6.58Deleterious1.000Probably Damaging1.000Probably Damaging3.00Benign0.00Affected3.4613-3-3-5.4-16.04180.665.90.10.0-0.60.1XPotentially PathogenicThe isobutyl side chain of Leu696, located in the middle of an α-helix (res. Leu685-Gln702), engages in hydrophobic packing with nearby residues (e.g., Leu441, Leu431, Leu692, Leu714) in the inter-helix space. Prolines lack a free amide group necessary for hydrogen bonding with the carbonyl group of Leu692 in the same manner as Leu696 in the WT. Consequently, the residue swap with proline disrupts the continuity of the secondary structure element in the variant simulations. Additionally, the side chain of Pro696 is not as optimal as Leu696 for hydrophobic packing in the inter-helix space.
c.1394T>CL465P
(3D Viewer)
Likely PathogenicGAPLikely Pathogenic 1-14.824Likely Pathogenic1.000Likely PathogenicLikely Pathogenic0.778Likely Pathogenic7.18Destabilizing0.310.85Destabilizing9.02Destabilizing2.73Destabilizing-6.96Deleterious1.000Probably Damaging1.000Probably Damaging2.29Pathogenic0.00Affected3.3734-3-3-5.4-16.04211.165.90.10.0-0.20.0XPotentially PathogenicThe iso-butyl side chain of Leu465, located in the middle of an α helix (res. Ala461–Phe476), packs with hydrophobic residues (e.g., Phe464, Met468, Tyr497, Ile494) in an inter-helix space formed with two other α helices (res. Ala461–Phe476 and res. Thr488-Gly502). In the variant simulations, the cyclic five-membered pyrrolidine ring of Pro465 is not as optimal as the side chain of Leu465 for filling the three α helix hydrophobic niche. Although the residue swap does not cause a large-scale conformational shift during the simulations, the H-bond between the backbone amide group of Leu465 and the backbone carbonyl group of Ala461 is lost. This, in turn, breaks the continuity of the α helix secondary structure element.
c.1292T>CL431P
(3D Viewer)
Likely PathogenicGAPLikely Pathogenic 1-14.222Likely Pathogenic0.996Likely PathogenicLikely Pathogenic0.659Likely Pathogenic6.78Destabilizing0.311.59Destabilizing9.19Destabilizing2.29Destabilizing-6.39Deleterious1.000Probably Damaging0.998Probably Damaging2.91Benign0.05Affected3.3729-3-3-5.4-16.04222.462.80.10.00.10.0XPotentially PathogenicThe iso-butyl side chain of Leu431, located in an α helix (res. Met414-Glu436), packs against other hydrophobic residues in an interhelix space (e.g., Val434, Leu435, Leu696, Leu711) in the WT simulations. While the backbone amide group of Leu431 forms an H-bond with the carbonyl group of His427, the cyclic five-membered pyrrolidine ring of Pro431, lacking the necessary amide group, cannot do the same. Thus, although the cyclic five-membered pyrrolidine ring of Pro431 packs almost as favorably as the side chain of Leu431 in the hydrophobic niche, the residue swap causes the α helix to partially unfold in the variant simulations.
c.1763T>CL588P
(3D Viewer)
Likely PathogenicGAPUncertain 1-14.771Likely Pathogenic1.000Likely PathogenicLikely Pathogenic0.932Likely Pathogenic5.61Destabilizing0.512.91Destabilizing9.26Destabilizing2.33Destabilizing-6.97Deleterious1.000Probably Damaging1.000Probably Damaging-1.42Pathogenic0.00Affected3.3834-3-3-5.4-16.04
c.2075T>CL692P
(3D Viewer)
Likely PathogenicGAPUncertain 1-16.447Likely Pathogenic1.000Likely PathogenicLikely Pathogenic0.668Likely Pathogenic9.19Destabilizing0.113.20Destabilizing11.20Destabilizing1.69Destabilizing-6.98Deleterious1.000Probably Damaging0.999Probably Damaging3.06Benign0.00Affected3.4217-3-3-5.4-16.04186.262.8-0.20.1-0.70.3XPotentially PathogenicThe isobutyl side chain of Leu692, located in the middle of an α-helix (res. Leu685-Gln702), engages in hydrophobic packing with nearby residues (e.g., Leu441, Leu431, Leu696) in the inter-helix space. Prolines lack a free amide group necessary for hydrogen bonding with the carbonyl group of Glu688 in the same manner as Leu692 in the WT. Consequently, the residue swap with proline disrupts the continuity of the secondary structure element in the variant simulations. Additionally, the side chain of Pro692 is not as optimal as Leu692 for hydrophobic packing in the inter-helix space.
c.2003C>TS668F
(3D Viewer)
Likely PathogenicGAPLikely Pathogenic 1-15.047Likely Pathogenic0.999Likely PathogenicLikely Pathogenic0.643Likely Pathogenic16.72Destabilizing5.011.07Destabilizing13.90Destabilizing0.00Likely Benign-5.98Deleterious0.999Probably Damaging0.935Probably Damaging3.18Benign0.00Affected3.3828-3-23.660.10250.9-59.6-0.10.10.00.1XXXPotentially PathogenicIn the WT simulations, the hydroxyl side chain of Ser668, located on an α-α loop connecting the two α-helices (res. Ser641-Glu666 and res. Leu685-Val699), forms hydrogen bonds with the backbone carbonyl groups of Leu664, Tyr665, and Glu666, as well as the guanidinium group of Arg573 on a nearby α-helix (res. Arg563-Glu578). In the variant simulations, the side chain of Phe668 cannot maintain the same hydrogen-bond network. Due to its larger size, it moves away to avoid steric hindrance. In the WT simulations, a network of hydrogen bonds between several residues (e.g., Asn669, Lys566, and Glu666) keeps both α-helices and the proceeding loop (res. Asn669-Asp684) tightly connected, but this setup is not present in the variant simulations. Additionally, in the variant simulations, the side chain of Arg573 shifts to form a more stable salt bridge with the carboxylate group of Glu582 instead of hydrogen bonding with Ser668 as in the WT simulations.

Found 757 rows. Show 200 rows per page. Page 4/4 |