Molecular dynamics and mutational analysis of a channelopathy mutation in the IIS6 helix of Ca V 1.2

Autism associated mutations in β2 subunit of voltage-gated calcium channels constitutively activate gene expression

Membrane depolarization triggers gene expression through voltage-gated calcium channels (VGCC) in a process called Excitation-transcription (ET) coupling. Mutations in the channel subunits α₁1.2, or β_2d, are associated with neurodevelopmental disorders such as ASD. Here, we found that two mutations S143F and G113S within the rat Cavβ_2a corresponding to autistic related mutations Cavβ_2d^S197F and Cavβ_2d^G167S in the human Cavβ_2d, activate ET-coupling via the RAS/ERK/CREB pathway. Membrane depolarization of HEK293 cells co-expressing α₁1.2 and α_2δ with Cavβ_2a^S143F or Cavβ_2a^G113S triggers constitutive transcriptional activation, which is correlated with facilitated channel activity. Similar to the Timothy-associated autistic mutation α₁1.2^G406R, constitutive gene activation is attributed to a hyperpolarizing shift in the activation kinetics of Cav1.2. Pulldown of RasGRF2 and RhoGEF by wt and the Cavβ_2a autistic mutants is consistent with Cavβ₂/Ras activation in ET coupling and implicates Rho signaling as yet another molecular pathway activated by Cavα₁1.2/Cavβ2 . Facilitated spontaneous channel activity preceding enhanced gene activation via the Ras/ERK/CREB pathway, appears a general molecular mechanism for Ca²⁺ channel mediated ASD and other neurodevelopmental disorders.

The effect of alterations of schizophrenia-associated genes on gamma band oscillations

Abnormalities in the synchronized oscillatory activity of neurons in general and, specifically in the gamma band, might play a crucial role in the pathophysiology of schizophrenia. While these changes in oscillatory activity have traditionally been linked to alterations at the synaptic level, we demonstrate here, using computational modeling, that common genetic variants of ion channels can contribute strongly to this effect. Our model of primary auditory cortex highlights multiple schizophrenia-associated genetic variants that reduce gamma power in an auditory steady-state response task. Furthermore, we show that combinations of several of these schizophrenia-associated variants can produce similar effects as the more traditionally considered synaptic changes. Overall, our study provides a mechanistic link between schizophrenia-associated common genetic variants, as identified by genome-wide association studies, and one of the most robust neurophysiological endophenotypes of schizophrenia.

Computational Modeling of Genetic Contributions to Excitability and Neural Coding in Layer V Pyramidal Cells: Applications to Schizophrenia Pathology

Pyramidal cells in layer V of the neocortex are one of the most widely studied neuron types in the mammalian brain. Due to their role as integrators of feedforward and cortical feedback inputs, they are well-positioned to contribute to the symptoms and pathology in mental disorders-such as schizophrenia-that are characterized by a mismatch between the internal perception and external inputs. In this modeling study, we analyze the input/output properties of layer V pyramidal cells and their sensitivity to modeled genetic variants in schizophrenia-associated genes. We show that the excitability of layer V pyramidal cells and the way they integrate inputs in space and time are altered by many types of variants in ion-channel and Ca²⁺ transporter-encoding genes that have been identified as risk genes by recent genome-wide association studies. We also show that the variability in the output patterns of spiking and Ca²⁺ transients in layer V pyramidal cells is altered by these model variants. Importantly, we show that many of the predicted effects are robust to noise and qualitatively similar across different computational models of layer V pyramidal cells. Our modeling framework reveals several aspects of single-neuron excitability that can be linked to known schizophrenia-related phenotypes and existing hypotheses on disease mechanisms. In particular, our models predict that single-cell steady-state firing rate is positively correlated with the coding capacity of the neuron and negatively correlated with the amplitude of a prepulse-mediated adaptation and sensitivity to coincidence of stimuli in the apical dendrite and the perisomatic region of a layer V pyramidal cell. These results help to uncover the voltage-gated ion-channel and Ca²⁺ transporter-associated genetic underpinnings of schizophrenia phenotypes and biomarkers.

Alterations in Schizophrenia-Associated Genes Can Lead to Increased Power in Delta Oscillations

Genome-wide association studies have implicated many ion channels in schizophrenia pathophysiology. Although the functions of these channels are relatively well characterized by single-cell studies, the contributions of common variation in these channels to neurophysiological biomarkers and symptoms of schizophrenia remain elusive. Here, using computational modeling, we show that a common biomarker of schizophrenia, namely, an increase in delta-oscillation power, may be a direct consequence of altered expression or kinetics of voltage-gated ion channels or calcium transporters. Our model of a circuit of layer V pyramidal cells highlights multiple types of schizophrenia-related variants that contribute to altered dynamics in the delta-frequency band. Moreover, our model predicts that the same membrane mechanisms that increase the layer V pyramidal cell network gain and response to delta-frequency oscillations may also cause a deficit in a single-cell correlate of the prepulse inhibition, which is a behavioral biomarker highly associated with schizophrenia.

Autism throughout genetics: Perusal of the implication of ion channels

BACKGROUND

Autism spectrum disorder (ASD) comprises a group of neurodevelopmental psychiatric disorders characterized by deficits in social interactions, interpersonal communication, repetitive and stereotyped behaviors and may be associated with intellectual disabilities. The description of ASD as a synaptopathology highlights the importance of the synapse and the implication of ion channels in the etiology of these disorders.

METHODS

A narrative and critical review of the relevant papers from 1982 to 2017 known by the authors was conducted.

RESULTS

Genome-wide linkages, association studies, and genetic analyses of patients with ASD have led to the identification of several candidate genes and mutations linked to ASD. Many of the candidate genes encode for proteins involved in neuronal development and regulation of synaptic function including ion channels and actors implicated in synapse formation. The involvement of ion channels in ASD is of great interest as they represent attractive therapeutic targets. In agreement with this view, recent findings have shown that drugs modulating ion channel function are effective for the treatment of certain types of patients with ASD.

CONCLUSION

This review describes the genetic aspects of ASD with a focus on genes encoding ion channels and highlights the therapeutic implications of ion channels in the treatment of ASD.

Pleiotropic effects of schizophrenia-associated genetic variants in neuron firing and cardiac pacemaking revealed by computational modeling

Schizophrenia patients have an increased risk of cardiac dysfunction. A possible factor underlying this comorbidity are the common variants in the large set of genes that have recently been discovered in genome-wide association studies (GWASs) as risk genes of schizophrenia. Many of these genes control the cell electrogenesis and calcium homeostasis. We applied biophysically detailed models of layer V pyramidal cells and sinoatrial node cells to study the contribution of schizophrenia-associated genes on cellular excitability. By including data from functional genomics literature to simulate the effects of common variants of these genes, we showed that variants of voltage-gated Na⁺ channel or hyperpolarization-activated cation channel-encoding genes cause qualitatively similar effects on layer V pyramidal cell and sinoatrial node cell excitability. By contrast, variants of Ca²⁺ channel or transporter-encoding genes mostly have opposite effects on cellular excitability in the two cell types. We also show that the variants may crucially affect the propagation of the cardiac action potential in the sinus node. These results may help explain some of the cardiac comorbidity in schizophrenia, and may facilitate generation of effective antipsychotic medications without cardiac side-effects such as arrhythmia.

Insights into channel dysfunction from modelling and molecular dynamics simulations

Developments in structural biology mean that the number of different ion channel structures has increased significantly in recent years. Structures of ion channels enable us to rationalize how mutations may lead to channelopathies. However, determining the structures of ion channels is still not trivial, especially as they necessarily exist in many distinct functional states. Therefore, the use of computational modelling can provide complementary information that can refine working hypotheses of both wild type and mutant ion channels. The simplest but still powerful tool is homology modelling. Many structures are available now that can provide suitable templates for many different types of ion channels, allowing a full three-dimensional interpretation of mutational effects. These structural models, and indeed the structures themselves obtained by X-ray crystallography, and more recently cryo-electron microscopy, can be subjected to molecular dynamics simulations, either as a tool to help explore the conformational dynamics in detail or simply as a means to refine the models further. Here we review how these approaches have been used to improve our understanding of how diseases might be linked to specific mutations in ion channel proteins. This article is part of the Special Issue entitled 'Channelopathies.'

CACNA1D de novo mutations in autism spectrum disorders activate Cav1.3 L-type calcium channels

BACKGROUND

Cav1.3 voltage-gated L-type calcium channels (LTCCs) are part of postsynaptic neuronal signaling networks. They play a key role in brain function, including fear memory and emotional and drug-taking behaviors. A whole-exome sequencing study identified a de novo mutation, p.A749G, in Cav1.3 α1-subunits (CACNA1D), the second main LTCC in the brain, as 1 of 62 high risk-conferring mutations in a cohort of patients with autism and intellectual disability. We screened all published genetic information available from whole-exome sequencing studies and identified a second de novo CACNA1D mutation, p.G407R. Both mutations are present only in the probands and not in their unaffected parents or siblings.

METHODS

We functionally expressed both mutations in tsA-201 cells to study their functional consequences using whole-cell patch-clamp.

RESULTS

The mutations p.A749G and p.G407R caused dramatic changes in channel gating by shifting (~15 mV) the voltage dependence for steady-state activation and inactivation to more negative voltages (p.A749G) or by pronounced slowing of current inactivation during depolarizing stimuli (p.G407R). In both cases, these changes are compatible with a gain-of-function phenotype.

CONCLUSIONS

Our data, together with the discovery that Cav1.3 gain-of-function causes primary aldosteronism with seizures, neurologic abnormalities, and intellectual disability, suggest that Cav1.3 gain-of-function mutations confer a major part of the risk for autism in the two probands and may even cause the disease. Our findings have immediate clinical relevance because blockers of LTCCs are available for therapeutic attempts in affected individuals. Patients should also be explored for other symptoms likely resulting from Cav1.3 hyperactivity, in particular, primary aldosteronism.

Timothy mutation disrupts the link between activation and inactivation in Ca(V)1.2 protein

The Timothy syndrome mutations G402S and G406R abolish inactivation of Ca_V1.2 and cause multiorgan dysfunction and lethal arrhythmias. To gain insights into the consequences of the G402S mutation on structure and function of the channel, we systematically mutated the corresponding Gly-432 of the rabbit channel and applied homology modeling. All mutations of Gly-432 (G432A/M/N/V/W) diminished channel inactivation. Homology modeling revealed that Gly-432 forms part of a highly conserved structure motif (G/A/G/A) of small residues in homologous positions of all four domains (Gly-432 (IS6), Ala-780 (IIS6), Gly-1193 (IIIS6), Ala-1503 (IVS6)). Corresponding mutations in domains II, III, and IV induced, in contrast, parallel shifts of activation and inactivation curves indicating a preserved coupling between both processes. Disruption between coupling of activation and inactivation was specific for mutations of Gly-432 in domain I. Mutations of Gly-432 removed inactivation irrespective of the changes in activation. In all four domains residues G/A/G/A are in close contact with larger bulky amino acids from neighboring S6 helices. These interactions apparently provide adhesion points, thereby tightly sealing the activation gate of Ca_V1.2 in the closed state. Such a structural hypothesis is supported by changes in activation gating induced by mutations of the G/A/G/A residues. The structural implications for Ca_V1.2 activation and inactivation gating are discussed.

Channelopathies in Cav1.1, Cav1.3, and Cav1.4 voltage-gated L-type Ca2+ channels

Voltage-gated Ca2+ channels couple membrane depolarization to Ca2+-dependent intracellular signaling events. This is achieved by mediating Ca2+ ion influx or by direct conformational coupling to intracellular Ca2+ release channels. The family of Cav1 channels, also termed L-type Ca2+ channels (LTCCs), is uniquely sensitive to organic Ca2+ channel blockers and expressed in many electrically excitable tissues. In this review, we summarize the role of LTCCs for human diseases caused by genetic Ca2+ channel defects (channelopathies). LTCC dysfunction can result from structural aberrations within their pore-forming alpha1 subunits causing hypokalemic periodic paralysis and malignant hyperthermia sensitivity (Cav1.1 alpha1), incomplete congenital stationary night blindness (CSNB2; Cav1.4 alpha1), and Timothy syndrome (Cav1.2 alpha1; reviewed separately in this issue). Cav1.3 alpha1 mutations have not been reported yet in humans, but channel loss of function would likely affect sinoatrial node function and hearing. Studies in mice revealed that LTCCs indirectly also contribute to neurological symptoms in Ca2+ channelopathies affecting non-LTCCs, such as Cav2.1 alpha1 in tottering mice. Ca2+ channelopathies provide exciting disease-related molecular detail that led to important novel insight not only into disease pathophysiology but also to mechanisms of channel function.

Different pathways for activation and deactivation in CaV1.2: a minimal gating model

Point mutations in pore-lining S6 segments of Ca_V1.2 shift the voltage dependence of activation into the hyperpolarizing direction and significantly decelerate current activation and deactivation. Here, we analyze theses changes in channel gating in terms of a circular four-state model accounting for an activation R–A–O and a deactivation O–D–R pathway. Transitions between resting-closed (R) and activated-closed (A) states (rate constants x(V) and y(V)) and open (O) and deactivated-open (D) states (u(V) and w(V)) describe voltage-dependent sensor movements. Voltage-independent pore openings and closures during activation (A–O) and deactivation (D–R) are described by rate constants α and β, and γ and δ, respectively. Rate constants were determined for 16-channel constructs assuming that pore mutations in IIS6 do not affect the activating transition of the voltage-sensing machinery (x(V) and y(V)). Estimated model parameters of 15 Ca_V1.2 constructs well describe the activation and deactivation processes. Voltage dependence of the “pore-releasing” sensor movement ((x(V)) was much weaker than the voltage dependence of “pore-locking” sensor movement (y(V)). Our data suggest that changes in membrane voltage are more efficient in closing than in opening Ca_V1.2. The model failed to reproduce current kinetics of mutation A780P that was, however, accurately fitted with individually adjusted x(V) and y(V). We speculate that structural changes induced by a proline substitution in this position may disturb the voltage-sensing domain.

Coupled and independent contributions of residues in IS6 and IIS6 to activation gating of CaV1.2

Voltage dependence and kinetics of Ca(V)1.2 activation are affected by structural changes in pore-lining S6 segments of the alpha(1)-subunit. Significant effects are induced by either proline or threonine substitutions in the lower third of segment IIS6 ("bundle crossing region"), where S6 segments are likely to seal the channel in the closed conformation (Hohaus, A., Beyl, S., Kudrnac, M., Berjukow, S., Timin, E. N., Marksteiner, R., Maw, M. A., and Hering, S. (2005) J. Biol. Chem. 280, 38471-38477). Here we report that S435P in IS6 results in a large shift of the activation curve (-25.9 +/- 1.2 mV) and slower current kinetics. Threonine substitutions at positions Leu-429 and Leu-434 induced a similar kinetic phenotype with shifted activation curves (L429T by -6.6 +/- 1.2 and L434T by -12.1 +/- 1.7 mV). Inactivation curves of all mutants were shifted to comparable extents as the activation curves. Interdependence of IS6 and IIS6 mutations was analyzed by means of mutant cycle analysis. Double mutations in segments IS6 and IIS6 induce either additive (L429T/I781T, -34.1 +/- 1.4 mV; L434T/I781T, -40.4 +/- 1.3 mV; L429T/L779T, -12.6 +/- 1.3 mV; and L434T/L779T, -22.4 +/- 1.3 mV) or nonadditive shifts of the activation curves along the voltage axis (S435P/I781T, -33.8 +/- 1.4 mV). Mutant cycle analysis revealed energetic coupling between residues Ser-435 and Ile-781, whereas other paired mutations in segments IS6 and IIS6 had independent effects on activation gating.