Frontiers in Molecular Biosciences

Background: Stress-induced childhood-onset neurodegeneration with variable ataxia and seizures (CONDSIAS) is a rare autosomal recessive disorder caused by biallelic variants in ADPRHL2, which encodes ADP-ribosylhydrolase 3 (ARH3), a key enzyme involved in poly (ADP-ribose) (PAR) metabolism. Although Minocycline has been reported to attenuate PAR-mediated neurotoxicity primarily through modulation of PARP-dependent pathways, whether it may also interact with ARH3 or influence the structural behavior of pathogenic ARH3 variants remains unknown. This study was designed to explore this possibility by integrating clinical observation with computational structural analyses. Methods: Comprehensive clinical evaluation, targeted Sanger sequencing, and in silico pathogenicity analyses were performed. Protein modeling, molecular docking, and 100-ns molecular dynamics simulations were conducted to evaluate the predicted structural consequences of the p.Thr79Pro variant and to explore potential interactions between ARH3 and Minocycline. Results: A homozygous ADPRHL2 variant (NM_017825.3:c.235A>C; p.Thr79Pro) was identified in a child with CONDSIAS. Computational analyses predicted reduced structural stability and increased conformational flexibility of the mutant ARH3 protein relative to the wild-type structure. MM-GBSA calculations estimated differences in binding free energies between the wild-type (-34.51 kcal/mol) and mutant (-39.76 kcal/mol) ARH3-Minocycline complexes, suggesting subtle differences in their predicted energetic profiles. Clinically, neurological progression appeared stable, with improved motor function observed during approximately one year of follow-up and no notable treatment-related adverse effects. Conclusions: By integrating clinical observations with computational structural analyses, this study provides preliminary computational support for the hypothesis that Minocycline may influence ARH3 conformational behavior in addition to its proposed effects on PARP-dependent pathways. Although these findings do not demonstrate direct molecular binding or therapeutic efficacy, they provide a biologically plausible framework for future biochemical, cellular, and functional investigations. Keywords: CONDSIAS; ADPRHL2; ARH3; Minocycline; molecular docking; molecular dynamics simulation; structural bioinformatics; translational medicine

Liquid-liquid phase separation is widely invoked in transcriptional regulation, yet prevailing models attribute condensate formation primarily to intrinsically disordered activation domains rather than structured DNA-binding motifs. Here, we overturn this view by demonstrating that the isolated basic leucine zipper (bZIP) domain of the yeast transcription factor Gcn4 undergoes robust DNA-induced phase separation in the complete absence of its activation domain. Using small-angle X-ray scattering in combination with all-atom molecular dynamics simulations and ensemble optimization, we directly resolve the conformational landscape of the Gcn4 bZIP-DNA complex across coexisting dilute and condensed phases. Beyond the canonical uninterrupted helical conformation captured in crystal structures, we identify a previously unrecognized minor population featuring a pronounced helical kink at the basic region-leucine zipper junction. These findings establish DNA binding as a sufficient physical driver of bZIP phase separation and demonstrate that small-angle scattering can quantitatively interrogate protein conformational ensembles within biomolecular condensates, opening new avenues for the structural chemistry of phase-separated systems.

Elucidating the mechanisms of transport kinetics in plants is crucial to develop crops that can use nutrients efficiently. The plant nitrate transporter NRT1.1 rapidly switches between high- and low-affinity transport modes to maintain an optimal uptake amidst fluctuations in nitrate levels. This functional switch is regulated by NRT1.1 phosphorylation, but the precise mechanisms remain poorly understood. Here, using an integrated molecular and systems-level modeling, we identify mechanisms underlying biphasic behaviour of NRT1.1. Phosphorylation of NRT1.1 and its binding to nitrate impacts its overall flexibility and synergistically modulates its global conformation, impacting the nitrate transport rate. Integrating these observations with a regulatory network involving kinases CIPK8/CIPK23 and calcium binding proteins CBL1/9, reveals that in high nitrate conditions, CIPK8-mediated sequestration of CBL1 disrupts the CIPK23-CBL complex required for NRT1.1 phosphorylation, switching NRT1.1 to a low-affinity mode. Together, our findings untangle the molecular complexity enabling NRT1.1 phosphorylation switch with broader implications in nitrate sensing and molecular-level adaption to fluctuating external nutrient levels.

BackgroundProtein raftophilicity refers to the affinity of proteins for cell biomembrane lipid domains, called rafts. Rafts are fluctuating nanoscale platforms that are enriched in cholesterol and sphingolipids, and that are considered relevant for cell signalling, viral function, and biomembrane trafficking. The dynamic partitioning of proteins into rafts depends on the physical and physico-chemical properties of the biomembranes where such proteins are embedded or attached; however it also depends on specific protein "features", such as acylation, glypidation, specific amino acid sequence motifs, transmembrane hydrophobic length, and surface accessible area to solvent. In this paper we present a method, and the resulting "ProRafts" predictor, that can be used to predict if a given mammal protein may be "raftophilic" or "non-raftophilic", without having an a priori knowledge of the physical and physico-chemical properties of the biomembranes where such protein is embedded or attached. ProRafts is based on a machine-learning algorithm, XGBoost, where data regarding the features of known raftophilic human-proteins fed the algorithm. ResultsProRafts enabled to predict correctly more than 80% of human proteins that are a priori known to be raftophilic; this is a promising result considering the limited size of the training dataset that we could build with data retrieved from protein databases. In addition, although we used protein features of known human raftophilic proteins, it was possible to identify accurately raft-proteins from other mammals than humans, such as mouse and rats. This finding suggests that certain protein features are sufficient to predict raftophilicity of proteins from different species. Moreover, our results indicated that phosphorylation may play a more relevant role for protein raftophilicity than indicated by previous studies. ConclusionRaftophilic proteins can be used as biomarkers in medical research, or can serve as targeting sites for therapeutics. In this respect, the machine learning method presented in this paper is a useful tool to guide experimental validations of raftophilicity of proteins in biomembranes, and facilitate the choice of proteins that can be used for experiments on biomimetic membranes.

SARS-CoV-2 ORF3a is believed to form ion channels, which may be involved in the modulation of virus release, and has been implicated in various cellular processes like the up-regulation of fibrinogen expression in lung epithelial cells, downregulation of type 1 interferon receptor, caspase-dependent apoptosis, and increasing IFNAR1 ubiquitination. ORF3a assemblies as homotetramers, which are stabilized by residue C133. A recent cryoEM structure of a homodimeric complex of ORF3a has been released. A lower-resolution cryoEM map of the tetramer suggests two dimers form it, arranged side by side. The dimers cryoEM structure revealed that each protomer contains three transmembrane helices arranged in a clockwise configuration forming a six helices transmembrane domain. This domains potential permeation pathway has six constrictions narrowing to about 1 [A] in radius, suggesting the structure solved is in a closed or inactivated state. At the cytosol end, the permeation pathway encounters a large and polar cavity formed by multiple beta strands from both protomers, which opens to the cytosolic milieu. We modeled the tetramer following the arrangement suggested by the low-resolution tetramer cryoEM map. Molecular dynamics simulations of the tetramer embedded in a membrane and solvated with 0.5 M of KCl were performed. Our simulations show the cytosolic cavity is quickly populated by both K+ and Cl-, yet with different dynamics. K+ ions moved relatively free inside the cavity without forming proper coordination sites. In contrast, Cl- ions enter the cavity, and three of them can become stably coordinated near the intracellular entrance of the potential permeation pathway by an inter-subunit network of positively charged amino acids. Consequently, the central cavitys electrostatic potential changed from being entirely positive at the beginning of the simulation to more electronegative at the end.

AlphaFold2-Multimer was used to generate structures of the heterotrimeric P2X7 receptors composed of wild-type P2X7A subunits and alternatively spliced subunits (P2X7B, P2X7E, P2X7J, and P2X7L) that have been confirmed in humans. The study supports laboratory research by providing insight into the structure and flexibility of the heterotrimeric alternatively spliced receptors in a simulated environment and may thereby aid structure-guided drug design. AbstractP2X7 receptors (P2X7Rs) are membrane-bound ATP-gated ion channels that are composed of three subunits. Different subunit structures may be expressed due to alternative splicing of the P2RX7 gene, altering the receptors function when combined with the wild-type P2X7A subunits. In this study, the application of the deep-learning method, AlphaFold2-Multimer (AF2M), for the generation of trimeric P2X7Rs was first validated by comparing an AF2M-generated rat wild-type P2X7A receptor with a structure determined by cryogenic electron microscopy (Protein Data Bank Identification: 6U9V). The results suggested AF2M could firstly, accurately predict the structures of P2X7Rs and secondly, accurately identify the highest quality model through the ranking system. Subsequently, AF2M was used to generate models of heterotrimeric alternatively spliced P2X7Rs consisting of one or two wild-type P2X7A subunits in combination with one or two P2X7B, P2X7E, P2X7J, and P2X7L splice variant subunits. The top-ranking models were deemed valid based on AF2Ms confidence measures, stability in molecular dynamics simulations, and consistent flexibility of the conserved regions between the models. Visual analysis of the heterotrimeric receptors identified missing residues in the ATP binding sites of the P2X7E, P2X7J, and P2X7L splice variants, likely translating into dysfunctional binding sites. Overall, the models produced in this study (available as supplementary material) unlock the possibility of structure-based studies into the heterotrimeric P2X7Rs.

Secreted apolipoprotein L1 (APOL1) is well known as an innate immune factor, protecting against African trypanosomiasis. The intracellular form has multiple functions, including regulating autophagy, intracellular vesicle trafficking, and ion channel activity. The APOL1 protein (G0) has two common variants (denoted G1 and G2) in the C-terminal region and are associated with a high risk of chronic kidney disease (CKD) and progression to end-stage kidney disease. Our previous studies using molecular modeling suggested that APOL1 G1 and G2 stabilize an autoinhibited state of the C-terminus, leading to impaired intracellular interactions with SNARE proteins. To characterize the structural consequence of kidney disease-associated APOL1 variants further, we assigned the C-terminal region proteins using 1H, 13C, 15N multidimensional nuclear magnetic resonance (NMR) spectra in solution in the presence of membrane mimetic dodecylphosphocholine micelles. We then derived models for the three-dimensional structure of APOL1-G0, and -G1 and -G2 variant C-terminal regions using the chemical shifts of the main chain nuclei followed by NMR relaxation measurements. The data suggest that changes in the three-dimensional structure of APOL1 C-terminal region induced by kidney disease-associated variants, not least the alteration of key sidechains and their interactions, could disrupt membrane association and the yet to be characterized protein-protein interactions including its binding partners, such as SNARE proteins. Such interactions could underlie the intracellular mechanisms that mediate the pathogenesis of CKD. In the future, one may try to reverse such structural and dynamics changes in the protein by designing agents that may bind and then mitigate APOL1 variant-associated CKD.

BackgroundDelirium is common in COVID-19 intensive care unit (ICU) patients. Biomarkers for prediction, detection, and monitoring are missing. Unbiased omics analyses are warranted to gain a systems biology view on pathophysiology. MethodsThis prospective observational satellite study aims to investigate the proteome signatures of COVID-19 ICU patients, comparing those with delirium to those without. This study was conducted in ICUs of a university hospital between March 2020 and September 2021. ICU patients of legal age with a positive SARS-CoV-2 test were screened daily for oversedation and delirium. Blood samples were taken thrice a week. 457 samples were analyzed using data-in-dependent acquisition mass spectrometry to determine protein levels. A mixed-effects regression model was developed to identify proteins significantly influenced by delirium, accounting for sex and age as confounders. This model also aimed to determine proteins that were either up- or downregulated in association with delirium. Additionally, an enrichment analysis was conducted to examine the biological pathways linked to these delirium-associated proteins. ResultsOut of 360 ICU patients, 69 were analyzed for protein profiling. Out of these 69 patients, 42 patients (60.9%) had delirium on ICU admission, and 27 (39.1%) did not. Based on the multivariate model, the analysis of 204 proteins unfolded 125 (61.3%) to be differentially expressed. In total, 80.8% (n=101) of these 125 proteins were associated with delirium. Of these, 10 proteins were uniquely associated with delirium and were not significant in the multivariate model (SERPING1, SERPINA7, HP, TGFBI, CD5L, IGHV3-7, IGHV1-46, IGHV3-15, IGHV3-23, and "IGHV4-34;IGHV4-38-2"). In the univariate model for delirium, six out of 111 significant proteins showed increased expression with a log2FC > 0.5: PIGR, MST1, LBP, CRP, SAA1, and "SAA1;SAA2"; while three showed decreased expression with a log2FC < - 0.5: HP, PPBP, and "HP;HPR". The enrichment analysis of delirium-influenced proteins revealed three significant pathways: "Network map of SARS-CoV-2 signaling" (M42569/WP5115), "Acute inflammatory response" (M10617), and "Regulation of defense response" (M15277). ConclusionWe identified a unique proteomic signature in COVID-19 ICU patients with delirium, including up- and downregulated proteins. These findings may provide potential biomarker candidates for the assessment of delirium risk and its underlying causes. These findings could be a further step towards the development of personalized, causative treatments for delirium and its monitoring in the ICU. Trial registrationThe study was retrospectively registered in the German Clinical Trials Register on May 13, 2020 (DRKS00021688).

NUDT15 or MTH2 is a member of NUDIX protein family that catalyze the hydrolysis of nucleotides and deoxynucleotides, including thioguanine analogues. NUTD15 has been reported as a DNA sanitizer in humans, and more recent studies have proved that genomic variants are related to a poor prognosis in inmoplastic and immunologic diseases with thioguanines. Despite of this, the role of NUTD15 in physiology and molecular biology is quite unclear, as well as the mechanism of action of this enzyme. The existence of clinically relevant variants has prompted the study of this enzymes, whose capacity to bound and hydrolyze thioguanine nucleotides is still poor understood. By using a combination of biomolecular modelling techniques together with molecular dynamics we have studied the monomeric wild type NUTD15, as well as two important variants R139C and R139H. Our findings reveal not only how nucleotide binding stabilizes the enzyme, but also how two loops are responsible for keeping the enzyme in a packed close conformation. Mutations in 2 helix affect a network of hydrophobic and {pi}-interactions that are responsible of active site enclosing.

A variant of the U1A protein containing four substitutions to ionizable residues was generated serendipitously due to a miscommunication. Biophysical measurements show that this variant has at least twice as much helical structure as the wild-type U1A and is trimeric in solution, in contrast to the monomeric wild type. In sharp contrast, structures predicted by deep-learning AI tools (AlphaFold2 and RoseTTAFold2) and transformer-based tools (OmegaFold and ESMFold) are all highly similar to the wild-type U1A (backbone RMSD < 1 [A]). Even more surprising, two of the substituted ionizable residues are predicted to be fully buried in the non-polar core of the protein, an outcome that contradicts well-established physico-chemical principles, as ionizable residues are normally located on the protein surface. To explore this effect further, we generated sequences containing up to all twelve residues that make up the non-polar core of U1A. Across thousands of sequences, and depending on the AI model used, the majority of predicted structures contained fully buried ionizable residues while still maintaining the overall U1A fold. We then examined two additional proteins of comparable size, acylphosphatase and the de novo-designed TOP7 fold, and observed the same phenomenon: AI models frequently predicted structures with buried ionizable residues that nevertheless retained the parent fold. When these AI-predicted structures were subjected to short (50 ns) molecular dynamics simulations using physics-based force fields such as CHARMM or AMBER, the structures rapidly relaxed into ensembles that exposed ionizable residues. We conclude that while AI-based structure prediction tools perform extremely well on naturally occurring sequences, they do not reliably encode the physico-chemical principles governing the placement of ionizable residues. A straightforward remedy is to include a brief molecular dynamics simulation as a final validation step for AI-generated structures.

Voltage-gated sodium (Nav) channels are pivotal for cellular signaling and mutations in Nav channels can lead to excitability disorders in cardiac, muscular, and neural tissues. A major cluster of pathological mutations localizes in the voltage-sensing domains (VSDs), resulting in either gain-of-function (GoF), loss-of-function (LoF) effects, or both. However, the mechanism behind this functional divergence of mutations at equivalent positions remains elusive. Through hotspot analysis, we identified three gating charges (R1, R2, and R3) as major mutational hotspots in VSDs. The same amino-acid substitutions at equivalent gating-charge positions in VSDI and VSDII of the cardiac sodium channel Nav1.5 show differential gating-property impacts in electrophysiology measurements. We conducted 120 {micro}s molecular dynamics (MD) simulations on wild-type and six mutants to elucidate the structural basis of their differential impacts. Our s-scale MD simulations with applied external electric fields captured VSD state transitions and revealed the differential structural dynamics between equivalent R-to-Q mutants. Notably, we observed transient leaky conformations in some mutants during structural transitions, offering a detailed structural explanation for gating-pore currents. Our salt-bridge network analysis uncovered VSD-specific and state-dependent interactions among gating charges, countercharges, and lipids. This detailed analysis elucidated how mutations disrupt critical electrostatic interactions, thereby altering VSD permeability and modulating gating properties. By demonstrating the crucial importance of considering the specific structural context of each mutation, our study represents a significant leap forward in understanding structure- function relationships in Nav channels. Our work establishes a robust framework for future investigations into the molecular basis of ion channel-related disorders.

Brown adipose tissue expresses uncoupling protein 1 (UCP1), a mitochondrial transporter that uncouples respiration from ATP synthesis and dissipates energy as heat, making it a target for treating obesity and related metabolic disorders. Here, we combine molecular dynamics simulations with mitochondrial respiration assays to investigate how purine nucleotides inhibit respiration uncoupling by UCP1. Simulations predict that GDP binds UCP1 in the common substrate binding site in an upright orientation, where the base moiety interacts with a pair of charged residues (R92/E191) that are specifically conserved in the subfamily of UCPs. E191, among others, interacts with purine but not pyrimidine bases, suggesting a rationale for nucleotide specificity in UCP1 inhibition. We also identify a triplet of uncharged residues involved in hydrophobic contacts with GDP. Site-directed mutagenesis of either I187 or W281 to alanine increases lauric acid-induced uncoupling activity of UCP1 and partially suppresses inhibition of UCP1 activity by GDP in yeast spheroplasts. The triple mutant (F88, I187, W281) to alanine is overactivated by lauric acid even in a high concentration of purine nucleotides. Variants at these positions may help increase energy expenditure in a cellular and therapeutic context.

Misfolding diseases that result in loss of function represent a considerable burden for both individuals and society. Primary hyperoxaluria type 1 (PH1) is a rare genetic disorder caused by mutations in the alanine:glyoxylate aminotransferase 1 (AGT) enzyme. The underlying molecular mechanisms causing PH1 are associated with protein misfolding (enhanced aggregation and mitochondrial mistargeting). The main therapeutic approach to increase patients lifespan and quality of life is a double kidney and liver transplantation. Alternative treatments such as gene and enzyme replacement and pharmacological chaperones are currently being introduced, but other alternatives are necessary. In this work, we developed and characterized a novel biotechnological approach using six single-domain nanobodies (NB-AGT-1 to -6) as potential therapeutics for PH1 misfolding. We show that NB-AGTs are very stable proteins and bind to pathogenic and non-pathogenic variants of AGT with extreme affinities (with Kd values from low nM to low pM). Structural studies showed that NB-AGTs bind to different epitopes of AGT with selectivity for different AGT variants. Experiments in cellular PH1 models showed that internalization of engineered NB-AGT-3 enhanced the specific activity of disease-associated variants. Overall, we show that NBs are a novel and promising approach to treat PH1 and other loss-of-function misfolding diseases.

Schnyder corneal dystrophy (SCD) is an autosomal dominant disease characterized by abnormal deposition of cholesterol and lipid in the cornea. The molecular mechanism underlying this process, which involves the interaction between UBAID1 and HMGCR, remains unclear. Here we investigate these events with in silico approaches. We built the homology models of UBIAD1 and HMGCR based on the existing crystal and cryo-EM structures. The UBIAD1 and HMGCR models are docked and their binding interactions are interrogated by MD simulation. We find that the transmembrane helices of UBIAD1 bind to sterol sensing domain of HMGCR. Upon binding of the GGPP substrate, UBIAD1 shows lower structural flexibility in the TM regions binding to HMGCR. The N102S and G177R mutations disrupts GGPP binding, thereby lowering the binding affinity of HMGCR. Overall, our modeling suggests that SCD mutations in UBIAD1 or lower GGPP concentration increase the structural flexibility of UBIAD1, thereby facilitating its association with HMGCR.Competing Interest StatementThe authors have declared no competing interest.View Full Text

The N34S variant of the trypsin inhibitor SPINK1 is the clinically most significant risk factor for chronic pancreatitis, but the underlying molecular mechanism could not be identified. Molecular dynamics simulations and docking of the generated conformational ensemble of SPINK1 to trypsin show that the mutation reduces the fraction of conformations that can directly participate in productive association, thereby reducing the association rate. The small change is difficult to detect by measuring the kinetics of SPINK1 binding to trypsin. However, kinetic modeling reveals that even a small change in the inhibition rate affects the trypsinogen to trypsin conversion rate at the early stage of the reaction when the trypsin concentration is very low, and the impact is substantially amplified by the autocatalytic mechanism of the conversion. Thus, the slightly reduced inhibition rate shortens the delay in the activation of trypsin release, which is therefore occurs within the pancreas.

Genetic variations expand the conformational landscape of proteins and may underlie cryptic properties that promote environmental adaptability. However, they can also represent modifying factors for disease susceptibility, by changing frustrated regions that in turn affect protein overall intracellular fitness. In this dichotomy between conservation and innovation, understanding at structural level how genetic variations keep the balance to maintain protein fitness represents an unmet need. Herein, we took advantage of known genetic variations of human alanine:glyoxylate aminotransferase (AGT1), which is present as a common major allelic form (AGT-Ma) and a minor polymorphic form (AGT-Mi) expressed in 20% of Caucasian population. By crystallographic studies and molecular dynamics simulations we showed that the polymorphic amino acid substitutions shape the conformational flexibility of AGT1 so that three surface regions that are structured in AGT-Ma become disordered in AGT-Mi, thanks to plasticity effects propagated from the mutation site(s) to the whole structure. In-depth biochemical characterisation of variants from a library encompassing the three regions correlate this plasticity to a fitness window between AGT-Ma and AGT-Mi, and suggest the existence of cryptic functions related to protein-protein interactions. These results establish that naturally-occurring genetic variations tip the balance between stability and frustration to expand the potential innovability of the protein.

BackgroundOnly twice have variants in the ITPR1 gene been described among patients with ataxia and miosis. Functional characterization of these variants is lacking. ObjectiveTo characterize a family affected by congenital ataxia and miosis associated with a novel ITPR1 variant and to provide a functional assessment for it and two previously reported variants. MethodsClinical characterization, genetic investigations, and segregation were performed. A novel variant c.7697T>C in ITPR1 was identified, HEK cells were transfected with vectors carrying our variant and two other previously published variants associated with ataxia and miosis. ResultsAtaxia was non-progressive in the reported family, the c.7697T>c ITPR1 variant segregated with disease. Functional validation showed that all the three ITPR1 variants were associated with reduced intracellular calcium release. ConclusionsHere, we present for the first time evidence of pathogenicity for 3 heterozygous ITPR1 variants in association with ataxia and miosis. Despite being localized in different ITPR1 protein domains, these variants converged on common functional defects.

With the expanding catalogue of novel disease-genes, there is an increasing need to establish the significance of potential disease-causing variants. Based on the idea that pathogenic variants in structured protein domains disturb folding and association with macromolecular assemblies, we employed Fluorescence Correlation and Cross-Correlation Spectroscopy (FCS and FCCS) to assess in vivo protein complex formation. Since the molecular underpinning of BRCA-associated breast and ovarian cancers is well defined and data from a recent genome editing screening allowed us to compare variant binding data with a reliable functional HRD test in addition to ClinVar and AlphaMissense data, we examined the binding of mutated wild-type BRCA1 or isolated RING and BRCT domains to BARD1 and RBBP8, respectively. The results demonstrate that FCCS, whether applied to full-length BRCA1 in live cells and/or to isolated domains in cellular lysates identified pathogenic BRCA1 RING or BRCT domain variants. We moreover demonstrate the feasibility of employing FCCS for analysis of HNPCC-related factor MSH2 and MEN1 factor Menin variants in combination with DNA mismatch repair factor MSH6 and transcription factor JUND, respectively. We propose that FCCS may be an appealing complement to current clinical procedures for classifying variants, for many monogenic diseases given its generic nature and ease of use.

ABCA4 is an ATP-binding cassette (ABC) transporter predominantly expressed in photoreceptors where it transports the substrate N-retinylidene-phosphatidylethanolamine across disc membranes thereby facilitating the clearance of retinal compounds from photoreceptor outer segments. Loss of function mutations in ABCA4 cause the accumulation of bisretinoids leading to Stargardt disease (STGD1) and other retinopathies. In this study, we examined the expression and functional properties of ABCA4 harboring disease-causing missense mutations in the two transmembrane domains (TMDs) of ABCA4. Our results indicate that these mutations lead to protein misfolding, loss in substrate binding, decreased ATPase activity or a combination of these properties. Additionally, we identified an arginine (R653) in transmembrane segment 2 of ABCA4 as a residue essential for substrate binding and substrate-stimulated ATPase activity. The expression and functional activity of the TMD variants correlate well with the severity of STGD1. Our studies provide a basis for developing and evaluating novel treatments for STGD1.

Aldo-keto oxidoreductase (AKR) inhibitors could reverse several cancer cells resistance to Cis-platin, but their role in resistance remains unclear. Our RNA-seq results showed de novo NAD biosynthesis-related genes, and NAD(P)H-dependent oxidoreductases were significantly upregulated in Cis-platin-resistant HepG2 hepatic cancer cells (HepG2-RC cells) compared with HepG2 cells. Knockdown of AKR1Cs could increase Cis-platin sensitivity in HepG2-RC cells about two-fold. Interestingly, the AKR1C inhibitor meclofenamic acid could increase Cis-platin sensitivity of HepG2-RC cells about eight-fold, indicating that knockdown of AKR1Cs only partially reversed the resistance. Meanwhile, the amount of total NAD and the ratio of NADH/NAD+ were increased in HepG2-RC cells compared with HepG2 cells. The increased NADH could be explained as a directly operating antioxidant to scavenge radicals induced by Cis-platin. We report here that NADH, which is produced by NAD(P)H-dependent oxidoreductases, plays a key role in the AKR-associated Cis-platin resistance of HepG2 hepatic cancer cells.