Nucleotide substitutions at the p.Gly117 and p.Thr180 mutational hot-spots of SKI alter molecular dynamics and may affect cell cycle.

Heterozygous deleterious variants in SKI cause Shprintzen-Goldberg Syndrome, which is mainly characterized by craniofacial features, neurodevelopmental disorder and thoracic aorta dilatations/aneurysms. The encoded protein is a member of the transforming growth factor beta signaling. Paucity of reported studies exploring the SGS molecular pathogenesis hampers disease recognition and clinical interpretation of private variants. Here, the unpublished c.349G>A, p.[Gly117Ser] and the recurrent c.539C>T, p.[Thr180Met] SKI variants were studied combining in silico and in vitro approach. 3D comparative modeling and calculation of the interaction energy predicted that both variants alter the SKI tertiary protein structure and its interactions. Computational data were functionally corroborated by the demonstration of an increase of MAPK phosphorylation levels and alteration of cell cycle in cells expressing the mutant SKI. Our findings confirmed the effects of SKI variants on MAPK and opened the path to study the role of perturbations of the cell cycle in SGS.

Deconstructing SARS-CoV-2 neutralization: A modular molecular framework for computational design and comparison of antibodies and nanobodies targeting the spike RBD.

Since 2020 the COVID-19 pandemic has led scientists to search for strategies to predict the transmissibility and virulence of new severe acute respiratory syndrome coronavirus 2 variants based on the estimation of the affinity of the spike receptor binding domain (RBD) for the human angiotensin-converting enzyme 2 (ACE2) receptor and/or neutralizing antibodies. In this context, our lab developed a computational pipeline to quickly quantify the free energy of interaction at the spike RBD/ACE2 protein-protein interface, reflecting the incidence trend observed in the transmissibility/virulence of the investigated variants. In this new study, we used our pipeline to estimate the free energy of interaction between the RBD from 10 variants, and 14 antibodies (ab), or 5 nanobodies (nb), highlighting the RBD regions preferentially targeted by the investigated ab/nb. Our structural comparative analysis and interaction energy calculations allowed us to propose the most promising RBD regions to be targeted by future ab/nb to be designed by site-directed mutagenesis of existing high-affinity ab/nb, to increase their affinity for the target RBD region, for preventing spike-RBD/ACE2 interactions and virus entry in host cells. Furthermore, we evaluated the ability of the investigated ab/nb to simultaneously interact with the three RBD located on the surface of the trimeric spike protein, which can alternatively be in up- or down- (all-3-up-, all-3-down-, 1-up-/2-down-, 2-up-/1-down-) conformations.

TNFRSF13B gene mutation in familial acute myeloid leukemia: A new piece in the complex scenario of hereditary predisposition?

TNFRSF13B mutations are widely associated with common variable immunodeficiency. TNFRSF13B was recently counted among relevant genes associated with childhood-onset of hematological malignancies; nonetheless, its role in acute myeloid leukemia (AML) remains unexplored. We report the study of a family with two cases of AML, sharing a germline TNFRSF13B mutation favoring the formation of a more stable complex with its ligand TNFSF13: a positive regulator of AML-initiating cells. Our data turn the spotlight onto the TNFRSF13B role in AML onset, inserting a new fragment into the complex scenario of a hereditary predisposition to myeloid neoplasms.

The Mycotoxin Patulin Inhibits the Mitochondrial Carnitine/Acylcarnitine Carrier (SLC25A20) by Interaction with Cys136 Implications for Human Health.

The effect of mycotoxin patulin (4-hydroxy-4H-furo [3,2c] pyran-2 [6H] -one) on the mitochondrial carnitine/acylcarnitine carrier (CAC, SLC25A20) was investigated. Transport function was measured as [3H]-carnitineex/carnitinein antiport in proteoliposomes reconstituted with the native protein extracted from rat liver mitochondria or with the recombinant CAC over-expressed in E. coli. Patulin (PAT) inhibited both the mitochondrial native and recombinant transporters. The inhibition was not reversed by physiological and sulfhydryl-reducing reagents, such as glutathione (GSH) or dithioerythritol (DTE). The IC50 derived from the dose-response analysis indicated that PAT inhibition was in the range of 50 µM both on the native and on rat and human recombinant protein. The kinetics process revealed a competitive type of inhibition. A substrate protection experiment confirmed that the interaction of PAT with the protein occurred within a protein region, including the substrate-binding area. The mechanism of inhibition was identified using the site-directed mutagenesis of CAC. No inhibition was observed on Cys mutants in which only the C136 residue was mutated. Mass spectrometry studies and in silico molecular modeling analysis corroborated the outcomes derived from the biochemical assays.

Heterozygous RNF13 Gain-of-Function Variants Are Associated with Congenital Microcephaly, Epileptic Encephalopathy, Blindness, and Failure to Thrive.

Accumulation of unfolded proteins in the endoplasmic reticulum (ER) initiates a stress response mechanism to clear out the unfolded proteins by either facilitating their re-folding or inducing their degradation. When this fails, an apoptotic cascade is initiated so that the affected cell is eliminated. IRE1α is a critical sensor of the unfolded-protein response, essential for initiating the apoptotic signaling. Here, we report an infantile neurodegenerative disorder associated with enhanced activation of IRE1α and increased apoptosis. Three unrelated affected individuals with congenital microcephaly, infantile epileptic encephalopathy, and profound developmental delay were found to carry heterozygous variants (c.932T>C [p.Leu311Ser] or c.935T>C [p.Leu312Pro]) in RNF13, which codes for an IRE1α-interacting protein. Structural modeling predicted that the variants, located on the surface of the protein, would not alter overall protein folding. Accordingly, the abundance of RNF13 and IRE1α was not altered in affected individuals' cells. However, both IRE1α-mediated stress signaling and stress-induced apoptosis were increased in affected individuals' cells. These results indicate that the RNF13 variants confer gain of function to the encoded protein and thereby lead to altered signaling of the ER stress response associated with severe neurodegeneration in infancy.

Pathogenic Variants in NUP214 Cause "Plugged" Nuclear Pore Channels and Acute Febrile Encephalopathy.

We report biallelic missense and frameshift pathogenic variants in the gene encoding human nucleoporin NUP214 causing acute febrile encephalopathy. Clinical symptoms include neurodevelopmental regression, seizures, myoclonic jerks, progressive microcephaly, and cerebellar atrophy. NUP214 and NUP88 protein levels were reduced in primary skin fibroblasts derived from affected individuals, while the total number and density of nuclear pore complexes remained normal. Nuclear transport assays exhibited defects in the classical protein import and mRNA export pathways in affected cells. Direct surface imaging of fibroblast nuclei by scanning electron microscopy revealed a large increase in the presence of central particles (known as "plugs") in the nuclear pore channels of affected cells. This observation suggests that large transport cargoes may be delayed in passage through the nuclear pore channel, affecting its selective barrier function. Exposure of fibroblasts from affected individuals to heat shock resulted in a marked delay in their stress response, followed by a surge in apoptotic cell death. This suggests a mechanistic link between decreased cell survival in cell culture and severe fever-induced brain damage in affected individuals. Our study provides evidence by direct imaging at the single nuclear pore level of functional changes linked to a human disease.

Protein structure analysis of the interactions between SARS-CoV-2 spike protein and the human ACE2 receptor: from conformational changes to novel neutralizing antibodies.

The recent severe acute respiratory syndrome, known as Coronavirus Disease 2019 (COVID-19) has spread so much rapidly and severely to induce World Health Organization (WHO) to declare a state of emergency over the new coronavirus SARS-CoV-2 pandemic. While several countries have chosen the almost complete lock-down for slowing down SARS-CoV-2 spread, the scientific community is called to respond to the devastating outbreak by identifying new tools for diagnosis and treatment of the dangerous COVID-19. With this aim, we performed an in silico comparative modeling analysis, which allows gaining new insights into the main conformational changes occurring in the SARS-CoV-2 spike protein, at the level of the receptor-binding domain (RBD), along interactions with human cells angiotensin-converting enzyme 2 (ACE2) receptor, that favor human cell invasion. Furthermore, our analysis provides (1) an ideal pipeline to identify already characterized antibodies that might target SARS-CoV-2 spike RBD, aiming to prevent interactions with the human ACE2, and (2) instructions for building new possible neutralizing antibodies, according to chemical/physical space restraints and complementary determining regions (CDR) mutagenesis of the identified existing antibodies. The proposed antibodies show in silico high affinity for SARS-CoV-2 spike RBD and can be used as reference antibodies also for building new high-affinity antibodies against present and future coronaviruses able to invade human cells through interactions of their spike proteins with the human ACE2. More in general, our analysis provides indications for the set-up of the right biological molecular context for investigating spike RBD-ACE2 interactions for the development of new vaccines, diagnostic kits, and other treatments based on the targeting of SARS-CoV-2 spike protein.

The Non-Gastric H+/K+ ATPase (ATP12A) Is Expressed in Mammalian Spermatozoa.

H+/K+ ATPase Type 2 is an heteromeric membrane protein involved in cation transmembrane transport and consists of two subunits: a specific α subunit (ATP12A) and a non-specific ß subunit. The aim of this study was to demonstrate the presence and establish the localization of ATP12A in spermatozoa from Bubalus bubalis, Bos taurus and Ovis aries. Immunoblotting revealed, in all three species, a major band (100 kDa) corresponding to the expected molecular mass. The ATP12A immunolocalization pattern showed, consistently in the three species, a strong signal at the acrosome. These results, described here for the first time in spermatozoa, are consistent with those observed for the ß1 subunit of Na+/K+ ATPase, suggesting that the latter may assemble with the α subunit to produce a functional ATP12A dimer in sperm cells. The above scenario appeared to be nicely supported by 3D comparative modeling and interaction energy calculations. The expression of ATP12A during different stages of bovine sperm maturation progressively increased, moving from epididymis to deferent ducts. Based on overall results, we hypothesize that ATP12A may play a role in acrosome reactions. Further studies will be required in order to address the functional role of this target protein in sperm physiology.

Personalized Medicine in Mitochondrial Health and Disease: Molecular Basis of Therapeutic Approaches Based on Nutritional Supplements and Their Analogs.

Mitochondrial diseases (MDs) may result from mutations affecting nuclear or mitochondrial genes, encoding mitochondrial proteins, or non-protein-coding mitochondrial RNA. Despite the great variability of affected genes, in the most severe cases, a neuromuscular and neurodegenerative phenotype is observed, and no specific therapy exists for a complete recovery from the disease. The most used treatments are symptomatic and based on the administration of antioxidant cocktails combined with antiepileptic/antipsychotic drugs and supportive therapy for multiorgan involvement. Nevertheless, the real utility of antioxidant cocktail treatments for patients affected by MDs still needs to be scientifically demonstrated. Unfortunately, clinical trials for antioxidant therapies using α-tocopherol, ascorbate, glutathione, riboflavin, niacin, acetyl-carnitine and coenzyme Q have met a limited success. Indeed, it would be expected that the employed antioxidants can only be effective if they are able to target the specific mechanism, i.e., involving the central and peripheral nervous system, responsible for the clinical manifestations of the disease. Noteworthily, very often the phenotypes characterizing MD patients are associated with mutations in proteins whose function does not depend on specific cofactors. Conversely, the administration of the antioxidant cocktails might determine the suppression of endogenous oxidants resulting in deleterious effects on cell viability and/or toxicity for patients. In order to avoid toxicity effects and before administering the antioxidant therapy, it might be useful to ascertain the blood serum levels of antioxidants and cofactors to be administered in MD patients. It would be also worthwhile to check the localization of mutations affecting proteins whose function should depend (less or more directly) on the cofactors to be administered, for estimating the real need and predicting the success of the proposed cofactor/antioxidant-based therapy.

Inverse Virtual Screening for the rapid re-evaluation of the presumed biological safe profile of natural products. The case of steviol from Stevia rebaudiana glycosides on farnesoid X receptor (FXR).

Nonnutritive sweeteners (NNSs) are widely employed as dietary substitutes for classical sugars thanks to their safety profile and low toxicity. In this study, a re-evaluation of the biological effects of steviol (1), the main metabolite from Stevia rebaudiana glycosides, was performed using the Inverse Virtual Screening (IVS) target fishing computational approach. Starting from well-known pharmacological properties of Stevia rebaudiana glycosides, this computational tool was employed for predicting the putative interacting targets of 1 and, afterwards, of its five synthetic ester derivatives 2-6, accounting a large panel of proteins involved in cancer and inflammation events. Applying this methodology, the farnesoid X receptor (FXR) was identified as the putative target partner of 1-6. The predicted ligand-protein interactions were corroborated by transactivation assays, specifically disclosing the agonistic activity of 1 and the antagonistic activities of 2-6 on FXR. The reported results highlight the feasibility of IVS as a fast and potent tool for predicting the interacting targets of query compounds, addressing the re-evaluation of their bioactivity. In light of the obtained results, the presumably safe profile of known compounds, such as the case of steviol (1), is critically discussed.

Mitochondrial Membranes of Human SH-SY5Y Neuroblastoma Cells Express Serotonin 5-HT7 Receptor.

Mitochondria in neurons contribute to energy supply, the regulation of synaptic transmission, Ca2+ homeostasis, neuronal excitability, and stress adaptation. In recent years, several studies have highlighted that the neurotransmitter serotonin (5-HT) plays an important role in mitochondrial biogenesis in cortical neurons, and regulates mitochondrial activity and cellular function in cardiomyocytes. 5-HT exerts its diverse actions by binding to cell surface receptors that are classified into seven distinct families (5-HT1 to 5-HT7). Recently, it was shown that 5-HT3 and 5-HT4 receptors are located on the mitochondrial membrane and participate in the regulation of mitochondrial function. Furthermore, it was observed that activation of brain 5-HT7 receptors rescued mitochondrial dysfunction in female mice from two models of Rett syndrome, a rare neurodevelopmental disorder characterized by severe behavioral and physiological symptoms. Our Western blot analyses performed on cell-lysate and purified mitochondria isolated from neuronal cell line SH-SY5Y showed that 5-HT7 receptors are also expressed into mitochondria. Maximal binding capacity (Bmax) obtained by Scatchard analysis on purified mitochondrial membranes was 0.081 pmol/mg of 5-HT7 receptor protein. Lastly, we evaluated the effect of selective 5-HT7 receptor agonist LP-211 and antagonist (inverse agonist) SB-269970 on mitochondrial respiratory chain (MRC) cytochrome c oxidase activity on mitochondria from SH-SY5Y cells. Our findings provide the first evidence that 5-HT7 receptor is also expressed in mitochondria.

Effect of cariporide on ram sperm pH regulation and motility: possible role of NHE1.

Sperm motility, a feature essential for in vivo fertilization, is influenced by intracellular pH (pHi) homeostasis. Several mechanisms are involved in pHi regulation, among which sodium-hydrogen exchangers (NHEs), a family of integral transmembrane proteins that catalyze the exchange of Na+ for H+ across lipid bilayers. A preliminary characterization of NHE activity and kinetic parameters, followed by analysis of the expression and localization of the protein in ram spermatozoa was performed. NHE activity showed an apparent Km for external Na+ of 17.61 mM. Immunoblotting revealed a molecular mass of 85 kDa. Immunolocalization pattern showed some species-specific aspects, such as positive labeling at the equatorial region of the sperm head. Cariporide, a selective NHE1 inhibitor, significantly reduced pHi recovery (85%). Similarly, exposure to cariporide significantly inhibited different motility parameters, including those related to sperm capacitation. In vitro fertilization (IVF) was not affected by cariporide, possibly due to the non-dramatic, although significant, drop in motility and velocity parameters or due to prolonged exposure during IVF, which may have caused progressive loss of its inhibitory effect. In conclusion, this is the first study documenting, in a large animal model (sheep) of well-known translational relevance, a direct functional role of NHE on sperm pHi and motility. The postulated specificity of cariporide toward isoform 1 of the Na+/H+ exchanger seems to suggest that NHE1 may contribute to the observed effects on sperm cell functionality.

The Bacterial Protein CNF1 as a Potential Therapeutic Strategy against Mitochondrial Diseases: A Pilot Study.

The Escherichia coli protein toxin cytotoxic necrotizing factor 1 (CNF1), which acts on the Rho GTPases that are key regulators of the actin cytoskeleton, is emerging as a potential therapeutic tool against certain neurological diseases characterized by cellular energy homeostasis impairment. In this brief communication, we show explorative results on the toxin's effect on fibroblasts derived from a patient affected by myoclonic epilepsy with ragged-red fibers (MERRF) that carries a mutation in the m.8344A>G gene of mitochondrial DNA. We found that, in the patient's cells, besides rescuing the wild-type-like mitochondrial morphology, CNF1 administration is able to trigger a significant increase in cellular content of ATP and of the mitochondrial outer membrane marker Tom20. These results were accompanied by a profound F-actin reorganization in MERRF fibroblasts, which is a typical CNF1-induced effect on cell cytoskeleton. These results point at a possible role of the actin organization in preventing or limiting the cell damage due to mitochondrial impairment and at CNF1 treatment as a possible novel strategy against mitochondrial diseases still without cure.

Biochemical characterization of a new mitochondrial transporter of dephosphocoenzyme A in Drosophila melanogaster.

CoA is an essential cofactor that holds a central role in cell metabolism. Although its biosynthetic pathway is conserved across the three domains of life, the subcellular localization of the eukaryotic biosynthetic enzymes and the mechanism behind the cytosolic and mitochondrial CoA pools compartmentalization are still under debate. In humans, the transport of CoA across the inner mitochondrial membrane has been ascribed to two related genes, SLC25A16 and SLC25A42 whereas in D. melanogaster genome only one gene is present, CG4241, phylogenetically closer to SLC25A42. CG4241 encodes two alternatively spliced isoforms, dPCoAC-A and dPCoAC-B. Both isoforms were expressed in Escherichia coli, but only dPCoAC-A was successfully reconstituted into liposomes, where transported dPCoA and, to a lesser extent, ADP and dADP but not CoA, which was a powerful competitive inhibitor. The expression of both isoforms in a Saccharomyces cerevisiae strain lacking the endogenous putative mitochondrial CoA carrier restored the growth on respiratory carbon sources and the mitochondrial levels of CoA. The results reported here and the proposed subcellular localization of some of the enzymes of the fruit fly CoA biosynthetic pathway, suggest that dPCoA may be synthesized and phosphorylated to CoA in the matrix, but it can also be transported by dPCoAC to the cytosol, where it may be phosphorylated to CoA by the monofunctional dPCoA kinase. Thus, dPCoAC may connect the cytosolic and mitochondrial reactions of the CoA biosynthetic pathway without allowing the two CoA pools to get in contact.

The switching mechanism of the mitochondrial ADP/ATP carrier explored by free-energy landscapes.

The ADP/ATP carrier (AAC) of mitochondria has been an early example for elucidating the transport mechanism alternating between the external (c-) and internal (m-) states (M. Klingenberg, Biochim. Biophys. Acta 1778 (2008) 1978-2021). An atomic resolution crystal structure of AAC is available only for the c-state featuring a three repeat transmembrane domain structure. Modeling of transport mechanism remained hypothetical for want of an atomic structure of the m-state. Previous molecular dynamics studies simulated the binding of ADP or ATP to the AAC remaining in the c-state. Here, a full description of the AAC switching from the c- to the m-state is reported using well-tempered metadynamics simulations. Free-energy landscapes of the entire translocation from the c- to the m-state, based on the gyration radii of the c- and m-gates and of the center of mass, were generated. The simulations revealed three free-energy basins attributed to the c-, intermediate- and m-states separated by activation barriers. These simulations were performed with the empty and with the ADP- and ATP-loaded AAC as well as with the poorly transported AMP and guanine nucleotides, showing in the free energy landscapes that ADP and ATP lowered the activation free-energy barriers more than the other substrates. Upon binding AMP and guanine nucleotides a deeper free-energy level stabilized the intermediate-state of the AAC2 hampering the transition to the m-state. The structures of the substrate binding sites in the different states are described producing a full picture of the translocation events in the AAC.

Nitro-substituted tetrahydroindolizines and homologs: Design, kinetics, and mechanism of α-glucosidase inhibition.

A series of 1-[(methylsulfonyl)methyl]-2-nitro-5,6,7,8-tetrahydroindolizines and homologs were designed, prepared, and evaluated as non-sugar-type α-glucosidase inhibitors. The inhibitory activity appeared to be related to cyclo homologation with the best congeners being tetrahydroindolizines. The introduction of a methoxycarbonyl group as an additional hydrogen bond acceptor into the exocyclic methylene group was beneficial affording the most potent congener 3e (half maximal inhibitory concentration, IC50=8.0±0.1µM) which displayed 25-fold higher inhibitory activity than 1-deoxynojirimycin (2, IC50=203±9µM)-the reference compound. Kinetic analysis indicated that compound 3e is a mixed inhibitor with preference for the free enzyme over the α-glucosidase-substrate complex (Ki,free=3.6µM; Ki,bound=7.6µM). Molecular docking experiments were in agreement with kinetic results indicating reliable interactions with both the catalytic cleft and other sites. Circular dichroism spectroscopy studies suggested that the inhibition exerted by 3e may involve changes in the secondary structure of the enzyme. Considering the relatively low molecular weight of 3e together with its high fraction of sp3 hybridized carbon atoms, this nitro-substituted tetrahydroindolizine may be considered as a good starting point towards new leads in the area of α-glucosidase inhibitors.

Truncating mutation in the nitric oxide synthase 1 gene is associated with infantile achalasia.

Nitric oxide is thought to have a role in the pathogenesis of achalasia. We performed a genetic analysis of 2 siblings with infant-onset achalasia. Exome analysis revealed that they were homozygous for a premature stop codon in the gene encoding nitric oxide synthase 1. Kinetic analyses and molecular modeling showed that the truncated protein product has defects in folding, nitric oxide production, and binding of cofactors. Heller myotomy had no effect in these patients, but sildenafil therapy increased their ability to drink. The finding recapitulates the previously reported phenotype of nitric oxide synthase 1-deficient mice, which have achalasia. Nitric oxide signaling appears to be involved in the pathogenesis of achalasia in humans.

Single-nucleotide evolution quantifies the importance of each site along the structure of mitochondrial carriers.

Mitochondrial carriers are membrane-embedded proteins consisting of a tripartite structure, a three-fold pseudo-symmetry, related sequences, and similar folding whose main function is to catalyze the transport of various metabolites, nucleotides, and coenzymes across the inner mitochondrial membrane. In this study, the evolutionary rate in vertebrates was screened at each of the approximately 50,000 nucleotides corresponding to the amino acids of the 53 human mitochondrial carriers. Using this information as a starting point, a scoring system was developed to quantify the evolutionary pressure acting on each site of the common mitochondrial carrier structure and estimate its functional or structural relevance. The degree of evolutionary selection varied greatly among all sites, but it was highly similar among the three symmetric positions in the tripartite structure, known as symmetry-related sites or triplets, suggesting that each triplet constitutes an evolutionary unit. Based on evolutionary selection, 111 structural sites (37 triplets) were found to be important. These sites play a key role in structure/function of mitochondrial carriers and are involved in either conformational changes (sites of the gates, proline-glycine levels, and aromatic belts) or in binding and specificity of the transported substrates (sites of the substrate-binding area in between the two gates). Furthermore, the evolutionary pressure analysis revealed that the matrix short helix sites underwent different degrees of selection with high inter-paralog variability. Evidence is presented that these sites form a new sequence motif in a subset of mitochondrial carriers, including the ADP/ATP translocator, and play a regulatory function by interacting with ligands and/or proteins of the mitochondrial matrix.

Prediction of high- and low-affinity quinol-analogue-binding sites in the aa3 and bo3 terminal oxidases from Bacillus subtilis and Escherichia coli1.

Haem-copper oxidases are the terminal enzymes in both prokaryotic and eukaryotic respiratory chains. They catalyse the reduction of dioxygen to water and convert redox energy into a transmembrane electrochemical proton gradient during their catalytic activity. Haem-copper oxidases show substantial structure similarity, but spectroscopic and biochemical analyses indicate that these enzymes contain diverse prosthetic groups and use different substrates (i.e. cytochrome c or quinol). Owing to difficulties in membrane protein crystallization, there are no definitive structural data about the quinol oxidase physiological substrate-binding site(s). In the present paper, we propose an atomic structure model for the menaquinol:O2 oxidoreductase of Bacillus subtilis (QOx.aa3). Furthermore, a multistep computational approach is used to predict residues involved in the menaquinol/menaquinone binding within B. subtilis QOx.aa3 as well as those involved in quinol/quinone binding within Escherichia coli QOx.bo3. Two specific sequence motifs, R70GGXDX4RXQX3PX3FX[D/N/E/Q]X2HYNE97 and G159GSPX2GWX2Y169 (B. subtilis numbering), were highlighted within QOx from Bacillales. Specific residues within the first and the second sequence motif participate in the high- and low-affinity substrate-binding sites respectively. Using comparative analysis, two analogous motifs, R71GFXDX4RXQX8[Y/F]XPPHHYDQ101 and G163EFX3GWX2Y173 (E. coli numbering) were proposed to be involved in Enterobacteriales/Rhodobacterales/Rhodospirillales QOx high- and low-affinity quinol-derivative-binding sites. Results and models are discussed in the context of the literature.

Mitochondrial glutamate carriers from Drosophila melanogaster: biochemical, evolutionary and modeling studies.

The mitochondrial carriers are members of a family of transport proteins that mediate solute transport across the inner mitochondrial membrane. Two isoforms of the glutamate carriers, GC1 and GC2 (encoded by the SLC25A22 and SLC25A18 genes, respectively), have been identified in humans. Two independent mutations in SLC25A22 are associated with severe epileptic encephalopathy. In the present study we show that two genes (CG18347 and CG12201) phylogenetically related to the human GC encoding genes are present in the D. melanogaster genome. We have functionally characterized the proteins encoded by CG18347 and CG12201, designated as DmGC1p and DmGC2p respectively, by overexpression in Escherichia coli and reconstitution into liposomes. Their transport properties demonstrate that DmGC1p and DmGC2p both catalyze the transport of glutamate across the inner mitochondrial membrane. Computational approaches have been used in order to highlight residues of DmGC1p and DmGC2p involved in substrate binding. Furthermore, gene expression analysis during development and in various adult tissues reveals that CG18347 is ubiquitously expressed in all examined D. melanogaster tissues, while the expression of CG12201 is strongly testis-biased. Finally, we identified mitochondrial glutamate carrier orthologs in 49 eukaryotic species in order to attempt the reconstruction of the evolutionary history of the glutamate carrier function. Comparison of the exon/intron structure and other key features of the analyzed orthologs suggests that eukaryotic glutamate carrier genes descend from an intron-rich ancestral gene already present in the common ancestor of lineages that diverged as early as bilateria and radiata.