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.

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 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.

Interactions of hydrophobic peptides with lipid bilayers: Monte Carlo simulations with M2delta.

We introduce here a novel Monte Carlo simulation method for studying the interactions of hydrophobic peptides with lipid membranes. Each of the peptide's amino acids is represented as two interaction sites: one corresponding to the backbone alpha-carbon and the other to the side chain, with the membrane represented as a hydrophobic profile. Peptide conformations and locations in the membrane and changes in the membrane width are sampled using the Metropolis criterion, taking into account the underlying energetics. Using this method we investigate the interactions between the hydrophobic peptide M2delta and a model membrane. The simulations show that starting from an extended conformation in the aqueous phase, the peptide first adsorbs onto the membrane surface, while acquiring an ordered helical structure. This is followed by formation of a helical-hairpin and insertion into the membrane. The observed path is in agreement with contemporary understanding of peptide insertion into biological membranes. Two stable orientations of membrane-associated M2delta were obtained: transmembrane (TM) and surface, and the value of the water-to-membrane transfer free energy of each of them is in agreement with calculations and measurements on similar cases. M2delta is most stable in the TM orientation, where it assumes a helical conformation with a tilt of 14 degrees between the helix principal axis and the membrane normal. The peptide conformation agrees well with the experimental data; average root-mean-square deviations of 2.1 A compared to nuclear magnetic resonance structures obtained in detergent micelles and supported lipid bilayers. The average orientation of the peptide in the membrane in the most stable configurations reported here, and in particular the value of the tilt angle, are in excellent agreement with the ones calculated using the continuum-solvent model and the ones observed in the nuclear magnetic resonance studies. This suggests that the method may be used to predict the three-dimensional structure of TM peptides.

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.

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.