Human G-protein coupled receptors (hGPCRs) constitute a large and highly pharmaceutically relevant membrane receptor superfamily. About half of the hGPCRs' family members are chemosensory receptors, involved in bitter taste and olfaction, along with a variety of other physiological processes. Hence these receptors constitute promising targets for pharmaceutical intervention. Molecular modeling has been so far the most important tool to get insights on agonist binding and receptor activation. Here we investigate both aspects by bioinformatics-based predictions across all bitter taste and odorant receptors for which site-directed mutagenesis data are available. First, we observe that state-of-the-art homology modeling combined with previously used docking procedures turned out to reproduce only a limited fraction of ligand/receptor interactions inferred by experiments. This is most probably caused by the low sequence identity with available structural templates, which limits the accuracy of the protein model and in particular of the side-chains' orientations. Methods which transcend the limited sampling of the conformational space of docking may improve the predictions. As an example corroborating this, we review here multi-scale simulations from our lab and show that, for the three complexes studied so far, they significantly enhance the predictive power of the computational approach. Second, our bioinformatics analysis provides support to previous claims that several residues, including those at positions 1.50, 2.50, and 7.52, are involved in receptor activation.

Aceruloplasminemia is a rare autosomal recessive genetic disease characterized by mild microcytic anemia, diabetes, retinopathy, liver disease, and progressive neurological symptoms due to iron accumulation in pancreas, retina, liver, and brain. The disease is caused by mutations in the Ceruloplasmin (CP) gene that produce a strong reduction or absence of ceruloplasmin ferroxidase activity, leading to an impairment of iron metabolism. Most patients described so far are from Japan. Prompt diagnosis and therapy are crucial to prevent neurological complications since, once established, they are usually irreversible. Here, we describe the largest series of non-Japanese patients with aceruloplasminemia published so far, including 13 individuals from 11 families carrying 13 mutations in the CP gene (7 missense, 3 frameshifts, and 3 splicing mutations), 10 of which are novel. All missense mutations were studied by computational modeling. Clinical manifestations were heterogeneous, but anemia, often but not necessarily microcytic, was frequently the earliest one. This study confirms the clinical and genetic heterogeneity of aceruloplasminemia, a disease expected to be increasingly diagnosed in the Next-Generation Sequencing (NGS) era. Unexplained anemia with low transferrin saturation and high ferritin levels without inflammation should prompt the suspicion of aceruloplasminemia, which can be easily confirmed by low serum ceruloplasmin levels. Collaborative joint efforts are needed to better understand the pathophysiology of this potentially disabling disease.

Hereditary Hemochromatosis (HH) is a genetically heterogeneous disorder caused by mutations in at least five different genes (HFE, HJV, TFR2, SLC40A1, HAMP) involved in the production or activity of the liver hormone hepcidin, a key regulator of systemic iron homeostasis. Nevertheless, patients with an HH-like phenotype that remains completely/partially unexplained despite extensive sequencing of known genes are not infrequently seen at referral centers, suggesting a role of still unknown genetic factors. A compelling candidate is Bone Morphogenetic Protein 6 (BMP6), which acts as a major activator of the BMP-SMAD signaling pathway, ultimately leading to the upregulation of hepcidin gene transcription. A recent seminal study by French authors has described three heterozygous missense mutations in BMP6 associated with mild to moderate lateonset iron overload (IO). Using an updated next-generation sequencing (NGS)-based genetic test in IO patients negative for the classical HFE p.Cys282Tyr mutation, we found three BMP6 heterozygous missense mutations in four patients from three different families. One mutation (p. Leu96Pro) has already been described and proven to be functional. The other two (p.Glu112Gln, p.Arg257His) were novel, and both were located in the pro-peptide domain known to be crucial for appropriate BMP6 processing and secretion. In silico modeling also showed results consistent with their pathogenetic role. The patients' clinical phenotypes were similar to that of other patients with BMP6-related IO recently described. Our results independently add further evidence to the role of BMP6 mutations as likely contributing factors to late-onset moderate IO unrelated to mutations in the established five HH genes.

Clinical and experimental studies indicate that muscarinic acetylcholine receptors are potential pharmacological targets for the treatment of neurological diseases. Although these receptors have been described in human, bovine and rat cerebral microvascular tissue, a subtype functional characterization in mouse brain endothelium is lacking. Here, we show that all muscarinic acetylcholine receptors (M1-M5) are expressed in mouse brain microvascular endothelial cells. The mRNA expression of M2, M3, and M5 correlates with their respective protein abundance, but a mismatch exists for M1 and M4 mRNA versus protein levels. Acetylcholine activates calcium transients in brain endothelium via muscarinic, but not nicotinic, receptors. Moreover, although M1 and M3 are the most abundant receptors, only a small fraction of M1 is present in the plasma membrane and functions in ACh-induced Ca2+ signaling. Bioinformatic analyses performed on eukaryotic muscarinic receptors demonstrate a high degree of conservation of the orthosteric binding site and a great variability of the allosteric site. In line with previous studies, this result indicates muscarinic acetylcholine receptors as potential pharmacological targets in future translational studies. We argue that research on drug development should especially focus on the allosteric binding sites of the M1 and M3 receptors.

Transmembrane proteins allow cells to extensively communicate with the external world in a very accurate and specific way. They form principal nodes in several signaling pathways and attract large interest in therapeutic intervention, as the majority pharmaceutical compounds target membrane proteins. Thus, according to the current genome annotation methods, a detailed structural/functional characterization at the protein level of each of the elements codified in the genome is also required. The extreme difficulty in obtaining high-resolution three-dimensional structures, calls for computational approaches. Here we review to which extent the efforts made in the last few years, combining the structural characterization of membrane proteins with protein bioinformatics techniques, could help describing membrane proteins at a genome-wide scale. In particular we analyze the use of comparative modeling techniques as a way of overcoming the lack of high-resolution three-dimensional structures in the human membrane proteome.