Search

Abstract Several large-scale whole-exome sequencing studies in patients with schizophrenia (SCZ) and autism spectrum disorder (ASD) have identified rare variants with modest or strong effect size as genetic risk factors. Dysregulation of cellular calcium homeostasis might be involved in SCZ/ASD pathogenesis, and genes encoding L-type voltage-gated calcium channel (VGCC) subunits Cav1.1 (CACNA1S), Cav1.2 (CACNA1C), Cav1.3 (CACNA1D), and T-type VGCC subunit Cav3.3 (CACNA1I) recently were identified as risk loci for psychiatric disorders. We performed a screening study, using the Ion Torrent Personal Genome Machine (PGM), of exon regions of these four candidate genes (CACNA1C, CACNA1D, CACNA1S, CACNA1I) in 370 Japanese patients with SCZ and 192 with ASD. Variant filtering was applied to identify biologically relevant mutations that were not registered in the dbSNP database or that have a minor allele frequency of less than 1% in East-Asian samples from databases; and are potentially disruptive, including nonsense, frameshift, canonical splicing site single nucleotide variants (SNVs), and non-synonymous SNVs predicted as damaging by five different in silico analyses. Each of these filtered mutations were confirmed by Sanger sequencing. If parental samples were available, segregation analysis was employed for measuring the inheritance pattern. Using our filter, we discovered one nonsense SNV (p.C1451* in CACNA1D), one de novo SNV (p.A36V in CACNA1C), one rare short deletion (p.E1675del in CACNA1D), and 14 NSstrict SNVs (non-synonymous SNV predicted as damaging by all of five in silico analyses). Neither p.A36V in CACNA1C nor p.C1451* in CACNA1D were found in 1871 SCZ cases, 380 ASD cases, or 1916 healthy controls in the independent sample set, suggesting that these SNVs might be ultra-rare SNVs in the Japanese population. The neuronal splicing isoform of Cav1.2 with the p.A36V mutation, discovered in the present study, showed reduced Ca2+-dependent inhibition, resulting in excessive Ca2+ entry through the mutant channel. These results suggested that this de novo SNV in CACNA1C might predispose to SCZ by affecting Ca2+ homeostasis. Thus, our analysis successfully identified several ultra-rare and potentially disruptive gene variants, lending partial support to the hypothesis that VGCC-encoding genes may contribute to the risk of SCZ/ASD.

Abstract A preceding experiment suggested that a compound, which inhibits binding of the REST/NRSF segment to the cleft of a receptor protein mSin3B, can be a potential drug candidate to ameliorate many neuropathies. We have recently developed an enhanced conformational sampling method, genetic-algorithm-guided multi-dimensional virtual-system-coupled canonical molecular dynamics, and in the present study, applied it to three systems consisting of mSin3B and one of three compounds, sertraline, YN3, and acitretin. Other preceding experiments showed that only sertraline inhibits the binding of REST/NRSF to mSin3B. The current simulation study produced the spatial distribution of the compounds around mSin3B, and showed that sertraline and YN3 bound to the cleft of mSin3B with a high propensity, although acitretin did not. Further analyses of the simulation data indicated that only the sertraline–mSin3B complex produced a hydrophobic core similar to that observed in the molecular interface of the REST/NRSF-mSin3B complex: An aromatic ring of sertraline sunk deeply in the mSin3B’s cleft forming a hydrophobic core contacting to hydrophobic amino-acid residues located at the bottom of the cleft. The present study proposes a step to design a compound that inhibits competitively the binding of a ligand to its receptor.

Abstract Background Rare genetic variants contribute to the etiology of both autism spectrum disorder (ASD) and schizophrenia (SCZ). Most genetic studies limit their focus to likely gene-disrupting mutations because they are relatively easier to interpret their effects on the gene product. Interpretation of missense variants is also informative to some pathophysiological mechanisms of these neurodevelopmental disorders; however, their contribution has not been elucidated because of relatively small effects. Therefore, we characterized missense variants detected in NRXN1, a well-known neurodevelopmental disease-causing gene, from individuals with ASD and SCZ. Methods To discover rare variants with large effect size and to evaluate their role in the shared etiopathophysiology of ASD and SCZ, we sequenced NRXN1 coding exons with a sample comprising 562 Japanese ASD and SCZ patients, followed by a genetic association analysis in 4273 unrelated individuals. Impact of each missense variant detected here on cell surface expression, interaction with NLGN1, and synaptogenic activity was analyzed using an in vitro functional assay and in silico three-dimensional (3D) structural modeling. Results Through mutation screening, we regarded three ultra-rare missense variants (T737M, D772G, and R856W), all of which affected the LNS4 domain of NRXN1α isoform, as disease-associated variants. Diagnosis of individuals with T737M, D772G, and R856W was 1ASD and 1SCZ, 1ASD, and 1SCZ, respectively. We observed the following phenotypic and functional burden caused by each variant. (i) D772G and R856W carriers had more serious social disabilities than T737M carriers. (ii) In vitro assay showed reduced cell surface expression of NRXN1α by D772G and R856W mutations. In vitro functional analysis showed decreased NRXN1α-NLGN1 interaction of T737M and D772G mutants. (iii) In silico 3D structural modeling indicated that T737M and D772G mutations could destabilize the rod-shaped structure of LNS2-LNS5 domains, and D772G and R856W could disturb N-glycan conformations for the transport signal. Conclusions The combined data suggest that missense variants in NRXN1 could be associated with phenotypes of neurodevelopmental disorders beyond the diagnosis of ASD and/or SCZ.

Geometric features of macromolecular shapes are important for binding with other molecules. Kawabata, T. and Go, N. (2007) defined a pocket as a space into which a small probe can enter, but a large probe cannot. In 2010, mathematical morphology (MM) was introduced to provide a more rigorous definition, and the program GHECOM was developed using the grid-based representation of molecules. This method was simple, but effective in finding the binding sites of small compounds on protein surfaces. Recently, many 3D structures of large macromolecules have been determined to contain large internal hollow spaces. Identification and size estimation of these spaces is important for characterizing their function and stability. Therefore, we employ the MM definition of pocket proposed by Manak, M. (2019)—a space into which an internal probe can enter, but an external probe cannot enter from outside of the macromolecules. This type of space is called a “cave pocket”, and is identified through molecular grid-representation. We define a “cavity” as a space into which a probe can enter, but cannot escape to the outside. Three types of spaces: cavity, pocket, and cave pocket were compared both theoretically and numerically. We proved that a cave pocket includes a pocket, and it is equal to a pocket if no cavity is found. We compared the three types of spaces for a variety of molecules with different-sized spherical probes; cave pockets were more sensitive than pockets for finding almost closed internal holes, allowing for more detailed representations of internal surfaces than cavities provide.

Ethoxy (EO) chain nonylphenol dehydrogenase (NPEO-DH) from Ensifer sp. AS08 and EO chain octylphenol dehydrogenase from Pseudomonas putida share common molecular characteristics with polyethylene glycol (PEG) dehydrogenases (PEG-DH) and comprise a PEG-DH subgroup in the family of glucose-methanol-choline (GMC) oxidoreductases that includes glucose/alcohol oxidase and glucose/choline dehydrogenase. Three-dimensional (3D) molecular modeling suggested that differences in the size, secondary structure and hydropathy in the active site caused differences in their substrate specificities toward EO chain alkylphenols and free PEGs. Based on 3D molecular modeling, site-directed mutagenesis was utilized to introduce mutations into potential catalytic residues of NPEO-DH. From steady state and rapid kinetic characterization of wild type and mutant NPEO-DHs, we can conclude that His465 and Asn507 are directly involved in the catalysis. Asn507 mediates the transfer of proton from a substrate to FAD and His465 transfers the same proton from the reduced flavin to an electron acceptor.