Rare Diseases: Implementation of Molecular Diagnosis, Pathogenesis Insights and Precision Medicine Treatment

: Rare Diseases (RD) do not have an exact definition since local authorities define the criteria in different ways, from fewer than 5 people in 10,000, according to the European Union, to the standard world average of 40 cases per 100,000 people [...]

Secondary structure assignment that accurately reflects physical and evolutionary characteristics.

RIGHTS : This article is licensed under the BioMed Central licence at http://www.biomedcentral.com/about/license which is similar to the 'Creative Commons Attribution Licence'. In brief you may : copy, distribute, and display the work; make derivative works; or make commercial use of the work - under the following conditions: the original author must be given credit; for any reuse or distribution, it must be made clear to others what the license terms of this work are.BACKGROUND: Secondary structure is used in hierarchical classification of protein structures, identification of protein features, such as helix caps and loops, for fold recognition, and as a precursor to ab initio structure prediction. There are several methods available for assigning secondary structure if the three-dimensional structure of the protein is known. Unfortunately they differ in their definitions, particularly in the exact positions of the termini. Additionally, most existing methods rely on hydrogen bonding, which means that important secondary structural classes, such as isolated beta-strands and poly-proline helices cannot be identified as they do not have characteristic hydrogen-bonding patterns. For this reason we have developed a more accurate method for assigning secondary structure based on main chain geometry, which also allows a more comprehensive assignment of secondary structure. RESULTS: We define secondary structure based on a number of geometric parameters. Helices are defined based on whether they fit inside an imaginary cylinder: residues must be within the correct radius of a central axis. Different types of helices (alpha, 3(10) or pi) are assigned on the basis of the angle between successive peptide bonds. beta-strands are assigned based on backbone dihedrals and with alternating peptide bonds. Thus hydrogen bonding is not required and beta-strands can be within a parallel sheet, antiparallel sheet, or can be isolated. Poly-proline helices are defined similarly, although with three-fold symmetry. CONCLUSION: We find that our method better assigns secondary structure than existing methods. Specifically, we find that comparing our methods with those of others, amino-acid trends at helix caps are stronger, secondary structural elements less likely to be concatenated together and secondary structure guided sequence alignment is improved. We conclude, therefore, that secondary structure assignments using our method better reflects physical and evolutionary characteristics of proteins. The program is available from http://www.bioinf.man.ac.uk/~lovell/segno.shtml

ReBaTSA: A simplified CeTSA protocol for studying recombinant mutant proteins in bacterial extracts

Introduction The study of protein stability is crucial to biochemistry and relies on different methodologies. Recently, the Cellular Thermal Shift Assay has been introduced to study protein stability in whole cells. Methods We report a novel application of CeTSA named ReBaTSA. This Recombinant Bacterial TSA was performed using clear extracts from bacteria expressing a recombinant protein, incubated at different temperatures, centrifuged and analyzed via SDS-PAGE. Results and conclusions We demonstrated the feasibility and reliability of this simplified approach. We validated the method using the protein phosphomannomutase-2 and its common mutants, which were compared in the presence or the absence of a known ligand

Identification and analysis of conserved pockets on protein surfaces

BACKGROUND: The interaction between proteins and ligands occurs at pockets that are often lined by conserved amino acids. These pockets can represent the targets for low molecular weight drugs. In order to make the research for new medicines as productive as possible, it is necessary to exploit "in silico" techniques, high throughput and fragment-based screenings that require the identification of druggable pockets on the surface of proteins, which may or may not correspond to active sites. RESULTS: We developed a tool to evaluate the conservation of each pocket detected on the protein surface by CastP. This tool was named DrosteP because it recursively searches for optimal input sequences to be used to calculate conservation. DrosteP uses a descriptor of statistical significance, Poisson p-value, as a target to optimize the choice of input sequences. To benchmark DrosteP we used monomeric or homodimer human proteins with known 3D-structure whose active site had been annotated in UniProt. DrosteP is able to detect the active site with high accuracy because in 81% of the cases it coincides with the most conserved pocket. Comparing DrosteP with analogous programs is difficult because the outputs are different. Nonetheless we could assess the efficacy of the recursive algorithm in the identification of active site pockets by calculating conservation with the same input sequences used by other programs. We analyzed the amino-acid composition of conserved pockets identified by DrosteP and we found that it differs significantly from the amino-acid composition of non conserved pockets. CONCLUSIONS: Several methods for predicting ligand binding sites on protein surfaces, that combine 3D-structure and evolutionary sequence conservation, have been proposed. Any method relying on conservation mainly depends on the choice of the input sequences. DrosteP chooses how deeply distant homologs must be collected to evaluate conservation and thus optimizes the identification of active site pockets. Moreover it recognizes conserved pockets other than those coinciding with the sites annotated in UniProt that might represent useful druggable sites. The distinctive amino-acid composition of conserved pockets provides useful hints on the fundamental principles underlying protein-ligand interaction. AVAILABILITY: http://www.icb.cnr.it/project/drosteppy

A thermodynamic assay to test pharmacological chaperones for Fabry disease

AbstractBackgroundThe majority of the disease-causing mutations affect protein stability, but not functional sites and are amenable, in principle, to be treated with pharmacological chaperones. These drugs enhance the thermodynamic stability of their targets. Fabry disease, a disorder caused by mutations in the gene encoding lysosomal alpha-galactosidase, represents an excellent model system to develop experimental protocols to test the efficiency of such drugs.MethodsThe stability of lysosomal alpha-galactosidase under different conditions was studied by urea-induced unfolding followed by limited proteolysis and Western blotting.ResultsWe measured the concentration of urea needed to obtain half-maximal unfolding because this parameter represents an objective indicator of protein stability.ConclusionsUrea-induced unfolding is a versatile technique that can be adapted to cell extracts containing tiny amounts of wild-type or mutant proteins. It allows testing of protein stability as a function of pH, in the presence or in the absence of drugs. Results are not influenced by the method used to express the protein in transfected cells.General significanceScarce and dispersed populations pose a problem for the clinical trial of drugs for rare diseases. This is particularly true for pharmacological chaperones that must be tested on each mutation associated with a given disease. Diverse in vitro tests are needed. We used a method based on chemically induced unfolding as a tool to assess whether a particular Fabry mutation is responsive to pharmacological chaperones, but, by no means is our protocol limited to this disease

In silico docking of urokinase plasminogen activator and integrins

Background: Urokinase, its receptor and the integrins are functionally associated and involved in regulation of cell signaling, migration, adhesion and proliferation. No structural information is available on this potential multimolecular complex. However, the tri-dimensional structure of urokinase, urokinase receptor and integrins is known. Results: We have modeled the interaction of urokinase on two integrins, alpha IIb beta 3 in the open configuration and alpha v beta 3 in the closed configuration. We have found that multiple lowest energy solutions point to an interaction of the kringle domain of uPA at the boundary between alpha and beta chains on the surface of the integrins. This region is not far away from peptides that have been previously shown to have a biological role in urokinase receptor/integrins dependent signaling. Conclusions: We demonstrated that in silico docking experiments can be successfully carried out to identify the binding mode of the kringle domain of urokinase on the scaffold of integrins in the open and closed conformation. Importantly we found that the binding mode was the same on different integrins and in both configurations. To get a molecular view of the system is a prerequisite to unravel the complex protein-protein interactions underlying urokinase/urokinase receptor/integrin mediated cell motility, adhesion and proliferation and to design rational in vitro experiments

Prediction of the responsiveness to pharmacological chaperones: lysosomal human alpha-galactosidase, a case of study

<p>Abstract</p> <p>Background</p> <p>The pharmacological chaperones therapy is a promising approach to cure genetic diseases. It relies on substrate competitors used at sub-inhibitory concentration which can be administered orally, reach difficult tissues and have low cost. Clinical trials are currently carried out for Fabry disease, a lysosomal storage disorder caused by inherited genetic mutations of alpha-galactosidase. Regrettably, not all genotypes respond to these drugs.</p> <p>Results</p> <p>We collected the experimental data available in literature on the enzymatic activity of ninety-six missense mutants of lysosomal alpha-galactosidase measured in the presence of pharmacological chaperones. We associated with each mutation seven features derived from the analysis of 3D-structure of the enzyme, two features associated with their thermo-dynamic stability and four features derived from sequence alone. Structural and thermodynamic analysis explains why some mutants of human lysosomal alpha-galactosidase cannot be rescued by pharmacological chaperones: approximately forty per cent of the non responsive cases examined can be correctly associated with a negative prognostic feature. They include mutations occurring in the active site pocket, mutations preventing disulphide bridge formation and severely destabilising mutations. Despite this finding, prediction of mutations responsive to pharmacological chaperones cannot be achieved with high accuracy relying on combinations of structure- and thermodynamic-derived features even with the aid of classical and state of the art statistical learning methods.</p> <p>We developed a procedure to predict responsive mutations with an accuracy as high as 87%: the method scores the mutations by using a suitable position-specific substitution matrix. Our approach is of general applicability since it does not require the knowledge of 3D-structure but relies only on the sequence.</p> <p>Conclusions</p> <p>Responsiveness to pharmacological chaperones depends on the structural/functional features of the disease-associated protein, whose complex interplay is best reflected on sequence conservation by evolutionary pressure. We propose a predictive method which can be applied to screen novel mutations of alpha galactosidase. The same approach can be extended on a genomic scale to find candidates for therapy with pharmacological chaperones among proteins with unknown tertiary structures.</p

The polymorphism L412F in TLR3 inhibits autophagy and is a marker of severe COVID-19 in males

The polymorphism L412F in TLR3 has been associated with several infectious diseases. However, the mechanism underlying this association is still unexplored. Here, we show that the L412F polymorphism in TLR3 is a marker of severity in COVID-19. This association increases in the sub-cohort of males. Impaired macroautophagy/autophagy and reduced TNF/TNFα production was demonstrated in HEK293 cells transfected with TLR3L412F-encoding plasmid and stimulated with specific agonist poly(I:C). A statistically significant reduced survival at 28 days was shown in L412F COVID-19 patients treated with the autophagy-inhibitor hydroxychloroquine (p = 0.038). An increased frequency of autoimmune disorders such as co-morbidity was found in L412F COVID-19 males with specific class II HLA haplotypes prone to autoantigen presentation. Our analyses indicate that L412F polymorphism makes males at risk of severe COVID-19 and provides a rationale for reinterpreting clinical trials considering autophagy pathways.publishedVersio

The KCNQ1OT1 imprinting control region and non-coding RNA: new properties derived from the study of Beckwith–Wiedemann syndrome and Silver–Russell syndrome cases

A cluster of imprinted genes at chromosome 11p15.5 is associated with the growth disorders, Silver–Russell syndrome (SRS) and Beckwith–Wiedemann syndrome (BWS). The cluster is divided into two domains with independent imprinting control regions (ICRs). We describe two maternal 11p15.5 microduplications with contrasting phenotypes. The first is an inverted and in cis duplication of the entire 11p15.5 cluster associated with the maintenance of genomic imprinting and with the SRS phenotype. The second is a 160 kb duplication also inverted and in cis, but resulting in the imprinting alteration of the centromeric domain. It includes the centromeric ICR (ICR2) and the most 5′ 20 kb of the non-coding KCNQ1OT1 gene. Its maternal transmission is associated with ICR2 hypomethylation and the BWS phenotype. By excluding epigenetic mosaicism, cell clones analysis indicated that the two closely located ICR2 sequences resulting from the 160 kb duplication carried discordant DNA methylation on the maternal chromosome and supported the hypothesis that the ICR2 sequence is not sufficient for establishing imprinted methylation and some other property, possibly orientation-dependent, is needed. Furthermore, the 1.2 Mb duplication demonstrated that all features are present for correct imprinting at ICR2 when this is duplicated and inverted within the entire cluster. In the individuals maternally inheriting the 160 kb duplication, ICR2 hypomethylation led to the expression of a truncated KCNQ1OT1 transcript and to down-regulation of CDKN1C. We demonstrated by chromatin RNA immunopurification that the KCNQ1OT1 RNA interacts with chromatin through its most 5′ 20 kb sequence, providing a mechanism likely mediating the silencing activity of this long non-coding RNA