On the evolution of protein-adenine binding.

Proteins' interactions with ancient ligands may reveal how molecular recognition emerged and evolved. We explore how proteins recognize adenine: a planar rigid fragment found in the most common and ancient ligands. We have developed a computational pipeline that extracts protein-adenine complexes from the Protein Data Bank, structurally superimposes their adenine fragments, and detects the hydrogen bonds mediating the interaction. Our analysis extends the known motifs of protein-adenine interactions in the Watson-Crick edge of adenine and shows that all of adenine's edges may contribute to molecular recognition. We further show that, on the proteins' side, binding is often mediated by specific amino acid segments ("themes") that recur across different proteins, such that different proteins use the same themes when binding the same adenine-containing ligands. We identify numerous proteins that feature these themes and are thus likely to bind adenine-containing ligands. Our analysis suggests that adenine binding has emerged multiple times in evolution.

An angular motion of a conserved four-helix bundle facilitates alternating access transport in the TtNapA and EcNhaA transporters.

There is ongoing debate regarding the mechanism through which cation/proton antiporters (CPAs), like Thermus thermophilus NapA (TtNapA) and Escherichia coli NapA (EcNhaA), alternate between their outward- and inward-facing conformations in the membrane. CPAs comprise two domains, and it is unclear whether the transition is driven by their rocking-bundle or elevator motion with respect to each other. Here we address this question using metadynamics simulations of TtNapA, where we bias conformational sampling along two axes characterizing the two proposed mechanisms: angular and translational motions, respectively. By applying the bias potential for the two axes simultaneously, as well as to the angular, but not the translational, axis alone, we manage to reproduce each of the two known states of TtNapA when starting from the opposite state, in support of the rocking-bundle mechanism as the driver of conformational change. Next, starting from the inward-facing conformation of EcNhaA, we sample what could be its long-sought-after outward-facing conformation and verify it using cross-linking experiments.

The Sigma-1 receptor is an ER-localized type II membrane protein.

The Sigma-1 receptor (S1R) is a transmembrane protein with important roles in cellular homeostasis in normal physiology and in disease. Especially in neurodegenerative diseases, S1R activation has been shown to provide neuroprotection by modulating calcium signaling, mitochondrial function and reducing endoplasmic reticulum (ER) stress. S1R missense mutations are one of the causes of the neurodegenerative Amyotrophic Lateral Sclerosis and distal hereditary motor neuronopathies. Although the S1R has been studied intensively, basic aspects remain controversial, such as S1R topology and whether it reaches the plasma membrane. To address these questions, we have undertaken several approaches. C-terminal tagging with a small biotin-acceptor peptide and BirA biotinylation in cells suggested a type II membrane orientation (cytosolic N-terminus). However, N-terminal tagging gave an equal probability for both possible orientations. This might explain conflicting reports in the literature, as tags may affect the protein topology. Therefore, we studied untagged S1R using a protease protection assay and a glycosylation mapping approach, introducing N-glycosylation sites. Both methods provided unambiguous results showing that the S1R is a type II membrane protein with a short cytosolic N-terminal tail. Assessments of glycan processing, surface fluorescence-activated cell sorting, and cell surface biotinylation indicated ER retention, with insignificant exit to the plasma membrane, in the absence or presence of S1R agonists or of ER stress. These findings may have important implications for S1R-based therapeutic approaches.

A novel biallelic loss-of-function mutation in TMCO1 gene confirming and expanding the phenotype spectrum of cerebro-facio-thoracic dysplasia.

The main clinical features of cerebro-facio-thoracic dysplasia (CFTD) syndrome, which were described over four decades ago, include facial dysmorphism, multiple malformations of the vertebrae and ribs, and intellectual disability. Recently, a TMCO1 gene mutation was shown to be responsible for an autosomal recessive CFTD syndrome characterized by craniofacial dysmorphism, skeletal anomalies, and intellectual disability. In the current report, we describe two members of a consanguineous family from an Arab community in Israel who were clinically diagnosed as suffering from craniofacial dysmorphism, skeletal anomalies, intellectual disability, and epilepsy. Both affected siblings had behavioral difficulties such as anxiety and emotional instability with impulsive behaviors. Whole-exome sequencing revealed a homozygous stop-gain mutation NM_019026.4: c.616C > T; p.(Arg206*) in exon 6 of the TMCO1 gene. Bioinformatics analysis suggested a structural model for the TMCO1 protein and its homologues. The clinical features of our patients were compared with those of the only other five studies available in the literature. We conclude that this mutation in the TMCO1 gene is responsible for the various clinical manifestations of CFTD syndrome exhibited by the patients studied that expand the phenotypic spectrum of the disease to include epilepsy as a characteristic feature of this syndrome.

Clinical, radiological, and genetic characteristics of 16 patients with ACO2 gene defects: Delineation of an emerging neurometabolic syndrome.

Mitochondrial aconitase is the second enzyme in the tricarboxylic acid (TCA) cycle catalyzing the interconversion of citrate into isocitrate and encoded by the nuclear gene ACO2. A homozygous pathogenic variant in the ACO2 gene was initially described in 2012 resulting in a novel disorder termed "infantile cerebellar retinal degeneration" (ICRD, OMIM#614559). Subsequently, additional studies reported patients with pathogenic ACO2 variants, further expanding the genetic and clinical spectrum of this disorder to include milder and later onset manifestations. Here, we report an international multicenter cohort of 16 patients (of whom 7 are newly diagnosed) with biallelic pathogenic variants in ACO2 gene. Most patients present in early infancy with severe truncal hypotonia, truncal ataxia, variable seizures, evolving microcephaly, and ophthalmological abnormalities of which the most dominant are esotropia and optic atrophy with later development of retinal dystrophy. Most patients remain nonambulatory and do no acquire any language, but a subgroup of patients share a more favorable course. Brain magnetic resonance imaging (MRI) is typically normal within the first months but global atrophy gradually develops affecting predominantly the cerebellum. Ten of our patients were homozygous to the previously reported c.336C>G founder mutation while the other six patients were all compound heterozygotes displaying 10 novel mutations of whom 2 were nonsense predicting a deleterious effect on enzyme function. Structural protein modeling predicted significant impairment in aconitase substrate binding in the additional missense mutations. This study provides the most extensive cohort of patients and further delineates the clinical, radiological, biochemical, and molecular features of ACO2 deficiency.

Arbitrium communication controls phage lysogeny through non-lethal modulation of a host toxin-antitoxin defence system.

Temperate Bacillus phages often utilize arbitrium communication to control lysis/lysogeny decisions, but the mechanisms by which this control is exerted remains largely unknown. Here we find that the arbitrium system of Bacillus subtilis phage ϕ3T modulates the host-encoded MazEF toxin-antitoxin system to this aim. Upon infection, the MazF ribonuclease is activated by three phage genes. At low arbitrium signal concentrations, MazF is inactivated by two phage-encoded MazE homologues: the arbitrium-controlled AimX and the later-expressed YosL proteins. At high signal, MazF remains active, promoting lysogeny without harming the bacterial host. MazF cleavage sites are enriched on transcripts of phage lytic genes but absent from the phage repressor in ϕ3T and other Spß-like phages. Combined with low activation levels of MazF during infections, this pattern explains the phage-specific effect. Our results show how a bacterial toxin-antitoxin system has been co-opted by a phage to control lysis/lysogeny decisions without compromising host viability.

SCAPER-Related Autosomal Recessive Retinitis Pigmentosa with Intellectual Disability: Confirming and Extending the Phenotypic Spectrum and Bioinformatics Analyses.

Mutations in the gene SCAPER (S phase Cyclin A-Associated Protein residing in the Endoplasmic Reticulum) have recently been associated with retinitis pigmentosa (RP) and intellectual disability (ID). In 2011, a possible involvement of SCAPER in human diseases was discovered for the first time due to the identification of a homozygous mutation causing ID in an Iranian family. Later, five studies were published in 2019 that described patients with autosomal recessive syndromic retinitis pigmentosa (arRP) accompanied by ID and attention-deficit/hyperactivity disorder (ADHD). This present study describes three patients from an Arab consanguineous family in Israel with similar clinical features of the SCAPER syndrome. In addition, new manifestations of ocular symptoms, nystagmus, glaucoma, and elevator palsy, were observed. Genetic testing of the patients and both parents via whole-exome sequencing revealed the homozygous mutation c.2023-2A>G in SCAPER. Phenotypic and genotypic descriptions for all available cases described in the literature including our current three cases (37 cases) were carried out, in addition to a bioinformatics analysis for all the genetic variants that was undertaken. Our study confirms and extends the clinical manifestations of SCAPER-related disorders.

Slow release of a synthetic auxin induces formation of adventitious roots in recalcitrant woody plants.

Clonal propagation of plants by induction of adventitious roots (ARs) from stem cuttings is a requisite step in breeding programs. A major barrier exists for propagating valuable plants that naturally have low capacity to form ARs. Due to the central role of auxin in organogenesis, indole-3-butyric acid is often used as part of commercial rooting mixtures, yet many recalcitrant plants do not form ARs in response to this treatment. Here we describe the synthesis and screening of a focused library of synthetic auxin conjugates in Eucalyptus grandis cuttings and identify 4-chlorophenoxyacetic acid-L-tryptophan-OMe as a competent enhancer of adventitious rooting in a number of recalcitrant woody plants, including apple and argan. Comprehensive metabolic and functional analyses reveal that this activity is engendered by prolonged auxin signaling due to initial fast uptake and slow release and clearance of the free auxin 4-chlorophenoxyacetic acid. This work highlights the utility of a slow-release strategy for bioactive compounds for more effective plant growth regulation.

Cation/proton antiporters: novel structure-driven pharmaceutical opportunities.

Cation/proton antiporters (CPAs) regulate cells' salt concentration and pH. Their malfunction is associated with a range of human pathologies, yet only a handful of CPA-targeting therapeutics are presently in clinical development. Here, we discuss how recently published mammalian protein structures and emerging computational technologies may help to bridge this gap.

Using evolutionary data to make sense of macromolecules with a "face-lifted" ConSurf.

The ConSurf web-sever for the analysis of proteins, RNA, and DNA provides a quick and accurate estimate of the per-site evolutionary rate among homologues. The analysis reveals functionally important regions, such as catalytic and ligand-binding sites, which often evolve slowly. Since the last report in 2016, ConSurf has been improved in multiple ways. It now has a user-friendly interface that makes it easier to perform the analysis and to visualize the results. Evolutionary rates are calculated based on a set of homologous sequences, collected using hidden Markov model-based search tools, recently embedded in the pipeline. Using these, and following the removal of redundancy, ConSurf assembles a representative set of effective homologues for protein and nucleic acid queries to enable informative analysis of the evolutionary patterns. The analysis is particularly insightful when the evolutionary rates are mapped on the macromolecule structure. In this respect, the availability of AlphaFold model structures of essentially all UniProt proteins makes ConSurf particularly relevant to the research community. The UniProt ID of a query protein with an available AlphaFold model can now be used to start a calculation. Another important improvement is the Python re-implementation of the entire computational pipeline, making it easier to maintain. This Python pipeline is now available for download as a standalone version. We demonstrate some of ConSurf's key capabilities by the analysis of caveolin-1, the main protein of membrane invaginations called caveolae.

Computational identification of a multi-peptide vaccine candidate in E2 glycoprotein against diverse Hepatitis C virus genotypes.

Hepatitis C Virus (HCV) is estimated to affect nearly 180 million people worldwide, culminating in ∼0.7 million yearly casualties. However, a safe vaccine against HCV is not yet available. This study endeavored to identify a multi-genotypic, multi-epitopic, safe, and globally competent HCV vaccine candidate. We employed a consensus epitope prediction strategy to identify multi-epitopic peptides in all known envelope glycoprotein (E2) sequences, belonging to diverse HCV genotypes. The obtained peptides were screened for toxicity, allergenicity, autoimmunity and antigenicity, resulting in two favorable peptides viz., P2 (VYCFTPSPVVVG) and P3 (YRLWHYPCTV). Evolutionary conservation analysis indicated that P2 and P3 are highly conserved, supporting their use as part of a designed multi-genotypic vaccine. Population coverage analysis revealed that P2 and P3 are likely to be presented by >89% Human Leukocyte Antigen (HLA) molecules from six geographical regions. Indeed, molecular docking predicted the physical binding of P2 and P3 to various representative HLAs. We designed a vaccine construct using these peptides and assessed its binding to toll-like receptor 4 (TLR-4) by molecular docking and simulation. Subsequent analysis by energy-based and machine learning tools predicted high binding affinity and pinpointed the key binding residues (i.e. hotspots) in P2 and P3. Also, a favorable immunogenic profile of the construct was predicted by immune simulations. We encourage the scientific community to validate our vaccine construct in vitro and in vivo.Communicated by Ramaswamy H. Sarma.

CLN8 Gene Compound Heterozygous Variants: A New Case and Protein Bioinformatics Analyses.

The CLN8 disease type refers to one of the neuronal ceroid lipofuscinoses (NCLs) which are the most common group of neurodegenerative diseases in childhood. The clinical phenotypes of this disease are progressive neurological deterioration that could lead to seizures, dementia, ataxia, visual failure, and various forms of abnormal movement. In the current study, we describe two patients who presented with atypical phenotypic manifestation and protracted clinical course of CLN8 carrying a novel compound heterozygous variant at the CLN8 gene. Our patients developed a mild phenotype of CLN8 disease: as they presented mild epilepsy, cognitive decline, mild learning disability, attention-deficit/hyperactivity disorder (ADHD), they developed a markedly protracted course of motor decline. Bioinformatic analyses of the compound heterozygous CLN8 gene variants were carried out. Most of the variants seem likely to act by compromising the structural integrity of regions within the protein. This in turn is expected to reduce the overall stability of the protein and render the protein less active to various degrees. The cases in our study confirmed and expanded the effect of compound heterozygous variants in CLN8 disease.

Dimeric and high-resolution structures of Chlamydomonas Photosystem I from a temperature-sensitive Photosystem II mutant.

Water molecules play a pivotal functional role in photosynthesis, primarily as the substrate for Photosystem II (PSII). However, their importance and contribution to Photosystem I (PSI) activity remains obscure. Using a high-resolution cryogenic electron microscopy (cryo-EM) PSI structure from a Chlamydomonas reinhardtii temperature-sensitive photoautotrophic PSII mutant (TSP4), a conserved network of water molecules - dating back to cyanobacteria - was uncovered, mainly in the vicinity of the electron transport chain (ETC). The high-resolution structure illustrated that the water molecules served as a ligand in every chlorophyll that was missing a fifth magnesium coordination in the PSI core and in the light-harvesting complexes (LHC). The asymmetric distribution of the water molecules near the ETC branches modulated their electrostatic landscape, distinctly in the space between the quinones and FX. The data also disclosed the first observation of eukaryotic PSI oligomerisation through a low-resolution PSI dimer that was comprised of PSI-10LHC and PSI-8LHC.

ConSurf-DB: An accessible repository for the evolutionary conservation patterns of the majority of PDB proteins.

Patterns observed by examining the evolutionary relationships among proteins of common origin can reveal the structural and functional importance of specific residue positions. In particular, amino acids that are highly conserved (i.e., their positions evolve at a slower rate than other positions) are particularly likely to be of biological importance, for example, for ligand binding. ConSurf is a bioinformatics tool for accurately estimating the evolutionary rate of each position in a protein family. Here we introduce a new release of ConSurf-DB, a database of precalculated ConSurf evolutionary conservation profiles for proteins of known structure. ConSurf-DB provides high-accuracy estimates of the evolutionary rates of the amino acids in each protein. A reliable estimate of a query protein's evolutionary rates depends on having a sufficiently large number of effective homologues (i.e., nonredundant yet sufficiently similar). With current sequence data, ConSurf-DB covers 82% of the PDB proteins. It will be updated on a regular basis to ensure that coverage remains high-and that it might even increase. Much effort was dedicated to improving the user experience. The repository is available at https://consurfdb.tau.ac.il/. BROADER AUDIENCE: By comparing a protein to other proteins of similar origin, it is possible to determine the extent to which each amino acid position in the protein evolved slowly or rapidly. A protein's evolutionary profile can provide valuable insights: For example, amino acid positions that are highly conserved (i.e., evolved slowly) are particularly likely to be of structural and/or functional importance, for example, for ligand binding and catalysis. We introduce here a new and improved version of ConSurf-DB, a continually updated database that provides precalculated evolutionary profiles of proteins with known structure.

Haploinsufficiency due to a novel ACO2 deletion causes mitochondrial dysfunction in fibroblasts from a patient with dominant optic nerve atrophy.

ACO2 is a mitochondrial protein, which is critically involved in the function of the tricarboxylic acid cycle (TCA), the maintenance of iron homeostasis, oxidative stress defense and the integrity of mitochondrial DNA (mtDNA). Mutations in the ACO2 gene were identified in patients suffering from a broad range of symptoms, including optic nerve atrophy, cortical atrophy, cerebellar atrophy, hypotonia, seizures and intellectual disabilities. In the present study, we identified a heterozygous 51 bp deletion (c.1699_1749del51) in ACO2 in a family with autosomal dominant inherited isolated optic atrophy. A complementation assay using aco1-deficient yeast revealed a growth defect for the mutant ACO2 variant substantiating a pathogenic effect of the deletion. We used patient-derived fibroblasts to characterize cellular phenotypes and found a decrease of ACO2 protein levels, while ACO2 enzyme activity was not affected compared to two age- and gender-matched control lines. Several parameters of mitochondrial function, including mitochondrial morphology, mitochondrial membrane potential or mitochondrial superoxide production, were not changed under baseline conditions. However, basal respiration, maximal respiration, and spare respiratory capacity were reduced in mutant cells. Furthermore, we observed a reduction of mtDNA copy number and reduced mtDNA transcription levels in ACO2-mutant fibroblasts. Inducing oxidative stress led to an increased susceptibility for cell death in ACO2-mutant fibroblasts compared to controls. Our study reveals that a monoallelic mutation in ACO2 is sufficient to promote mitochondrial dysfunction and increased vulnerability to oxidative stress as main drivers of cell death related to optic nerve atrophy.

A novel variant of the human mitochondrial DnaJ protein, Tid1, associates with a human disease exhibiting developmental delay and polyneuropathy.

Here, we describe a single patient from a consanguineous family, who suffers from developmental delay, intellectual disability, hypermetropia, moderate alternating esotropia, unsteady gait, and peripheral polyneuropathy. Brain MRI revealed basal ganglia disease. Exome analysis disclosed a homozygous variant, c.452G>C (p.(Arg151Thr)), in TID1, encoding a mitochondrial J-protein chaperone that is known for its function in assisting the Hsp70 chaperone, mortalin, in mediating the refolding of denatured protein and dissolving protein aggregates. Results from in vitro import assays showed that both wild type and c.452G>C (p.(Arg151Thr)) are efficiently imported into isolated mitochondria. However, the import rate of the c.452G>C (p.(Arg151Thr)) variant was less than that of the wild-type protein. In the second part of this study, we demonstrated, in vitro, that the disaggregation function of the mortalin/Tid1 team is compromised in the TID1 c.452G>C (p.(Arg151Thr)) variant, as its chaperone activity has a level similar to that of the non-functional H→Q HPD domain variant. The results shed light on the essential function played by Tid1 during neuronal development.

Broad phylogenetic analysis of cation/proton antiporters reveals transport determinants.

Cation/proton antiporters (CPAs) play a major role in maintaining living cells' homeostasis. CPAs are commonly divided into two main groups, CPA1 and CPA2, and are further characterized by two main phenotypes: ion selectivity and electrogenicity. However, tracing the evolutionary relationships of these transporters is challenging because of the high diversity within CPAs. Here, we conduct comprehensive evolutionary analysis of 6537 representative CPAs, describing the full complexity of their phylogeny, and revealing a sequence motif that appears to determine central phenotypic characteristics. In contrast to previous suggestions, we show that the CPA1/CPA2 division only partially correlates with electrogenicity. Our analysis further indicates two acidic residues in the binding site that carry the protons in electrogenic CPAs, and a polar residue in the unwound transmembrane helix 4 that determines ion selectivity. A rationally designed triple mutant successfully converted the electrogenic CPA, EcNhaA, to be electroneutral.

On the regulatory role of dipeptidyl peptidase IV (=CD=adenosine deaminase complexing protein) on adenosine deaminase activity.

The molecular mechanism controlling the variable activity of the malignancy marker adenosine deaminase (ADA) is enigmatic. ADA activity was found to be modulated by the membrane-bound adenosine deaminase complexing protein (CP=DPPIV=CD26). The role of lipid-protein interactions in this modulation was sought. While direct solubilization of ADA in vesicles resulted in loss of ADA activity, the binding of ADA to CP reconstituted in vesicles restored the specific activity. The activity of ADA, free or bound to CP in solution, resulted in continuous linear Arrhenius plots. However, ADA bound to reconstituted CP exhibited two breaks associated with approximately 30% increased activity, at 25 and 13 degrees C, yielding three lines with similar apparent activation energies (E(a)). Continuum solvent model calculations of the free energy of transfer of the transmembrane helix of CP from the aqueous phase into membranes of various widths show that the most favorable orientations of the helix above and below the main phase transition may be different. We suggest that the 20% change in the thickness of the bilayer below and above the main phase transition may modify the orientation of CP in the membrane, thereby affecting substrate accessibility of ADA. This could account for ADA's reduced activity associated with increased membrane fluidity in transformed vs. normal fibroblasts.

The ConSurf-HSSP database: the mapping of evolutionary conservation among homologs onto PDB structures.

The HSSP (Homology-Derived Secondary Structure of Proteins) database provides multiple sequence alignments (MSAs) for proteins of known three-dimensional (3D) structure in the Protein Data Bank (PDB). The database also contains an estimate of the degree of evolutionary conservation at each amino acid position. This estimate, which is based on the relative entropy, correlates with the functional importance of the position; evolutionarily conserved positions (i.e., positions with limited variability and low entropy) are occasionally important to maintain the 3D structure and biological function(s) of the protein. We recently developed the Rate4Site algorithm for scoring amino acid conservation based on their calculated evolutionary rate. This algorithm takes into account the phylogenetic relationships between the homologs and the stochastic nature of the evolutionary process. Here we present the ConSurf-HSSP database of Rate4Site estimates of the evolutionary rates of the amino acid positions, calculated using HSSP's MSAs. The database provides precalculated evolutionary rates for nearly all of the PDB. These rates are projected, using a color code, onto the protein structure, and can be viewed online using the ConSurf server interface. To exemplify the database, we analyzed in detail the conservation pattern obtained for pyruvate kinase and compared the results with those observed using the relative entropy scores of the HSSP database. It is reassuring to know that the main functional region of the enzyme is detectable using both conservation scores. Interestingly, the ConSurf-HSSP calculations mapped additional functionally important regions, which are moderately conserved and were overlooked by the original HSSP estimate. The ConSurf-HSSP database is available online (http://consurf-hssp.tau.ac.il).

Interactions of the M2delta segment of the acetylcholine receptor with lipid bilayers: a continuum-solvent model study.

M2delta, one of the transmembrane segments of the nicotinic acetylcholine receptor, is a 23-amino-acid peptide, frequently used as a model for peptide-membrane interactions. In this and the companion article we describe studies of M2delta-membrane interactions, using two different computational approaches. In the present work, we used continuum-solvent model calculations to investigate key thermodynamic aspects of its interactions with lipid bilayers. M2delta was represented in atomic detail and the bilayer was represented as a hydrophobic slab embedded in a structureless aqueous phase. Our calculations show that the transmembrane orientation is the most favorable orientation of the peptide in the bilayer, in good agreement with both experimental and computational data. Moreover, our calculations produced the free energy of association of M2delta with the lipid bilayer, which, to our knowledge, has not been reported to date. The calculations included 10 structures of M2delta, determined by nuclear magnetic resonance in dodecylphosphocholine micelles. All the structures were found to be stable inside the lipid bilayer, although their water-to-membrane transfer free energies differed by as much as 12 kT. Although most of the structures were roughly linear, a single structure had a kink in its central region. Interestingly, this structure was found to be the most stable inside the lipid bilayer, in agreement with molecular dynamics simulations of the peptide and with the recently determined structure of the intact receptor. Our analysis showed that the kink reduced the polarity of the peptide in its central region by allowing the electrostatic masking of the Gln13 side chain in that area. Our calculations also showed a tendency for the membrane to deform in response to peptide insertion, as has been previously found for the membrane-active peptides alamethicin and gramicidin. The results are compared to Monte Carlo simulations of the peptide-membrane system, as presented in the accompanying article.