Modeling the Orthosteric Binding Site of the G Protein-Coupled Odorant Receptor OR5K1.

With approximately 400 encoding genes in humans, odorant receptors (ORs) are the largest subfamily of class A G protein-coupled receptors (GPCRs). Despite its high relevance and representation, the odorant-GPCRome is structurally poorly characterized: no experimental structures are available, and the low sequence identity of ORs to experimentally solved GPCRs is a significant challenge for their modeling. Moreover, the receptive range of most ORs is unknown. The odorant receptor OR5K1 was recently and comprehensively characterized in terms of cognate agonists. Here, we report two additional agonists and functional data of the most potent compound on two mutants, L1043.32 and L2556.51. Experimental data was used to guide the investigation of the binding modes of OR5K1 ligands into the orthosteric binding site using structural information from AI-driven modeling, as recently released in the AlphaFold Protein Structure Database, and from homology modeling. Induced-fit docking simulations were used to sample the binding site conformational space for ensemble docking. Mutagenesis data guided side chain residue sampling and model selection. We obtained models that could better rationalize the different activity of active (agonist) versus inactive molecules with respect to starting models and also capture differences in activity related to minor structural differences. Therefore, we provide a model refinement protocol that can be applied to model the orthosteric binding site of ORs as well as that of GPCRs with low sequence identity to available templates.

MERMAID: dedicated web server to prepare and run coarse-grained membrane protein dynamics.

Atomistic molecular dynamics simulations of membrane proteins have been shown to be extremely useful for characterizing the molecular features underlying their function, but require high computational power, limiting the understanding of complex events in membrane proteins, e.g. ion channels gating, GPCRs activation. To overcome this issue, it has been shown that coarse-grained approaches, although requiring less computational power, are still capable of correctly describing molecular events underlying big conformational changes in biological systems. Here, we present the Martini coarse-grained membrane protein dynamics (MERMAID), a publicly available web interface that allows the user to prepare and run coarse-grained molecular dynamics (CGMD) simulations and to analyse the trajectories.

Genetic and Clinical Heterogeneity in Thirteen New Cases with Aceruloplasminemia. Atypical Anemia as a Clue for an Early Diagnosis.

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.

Impact of Cholesterol on the Stability of Monomeric and Dimeric Forms of the Translocator Protein TSPO: A Molecular Simulation Study.

The translocator protein (TSPO) is a transmembrane protein present across the three domains of life. Its functional quaternary structure consists of one or more subunits. In mice, the dimer-to-monomer equilibrium is shifted in vitro towards the monomer by adding cholesterol, a natural component of mammalian membranes. Here, we present a coarse-grained molecular dynamics study on the mouse protein in the presence of a physiological content and of an excess of cholesterol. The latter turns out to weaken the interfaces of the dimer by clusterizing mostly at the inter-monomeric space and pushing the contact residues apart. It also increases the compactness and the rigidity of the monomer. These two factors might play a role for the experimentally observed incremented stability of the monomeric form with increased content of cholesterol. Comparison with simulations on bacterial proteins suggests that the effect of cholesterol is much less pronounced for the latter than for the mouse protein.

CGMD Platform: Integrated Web Servers for the Preparation, Running, and Analysis of Coarse-Grained Molecular Dynamics Simulations.

Advances in coarse-grained molecular dynamics (CGMD) simulations have extended the use of computational studies on biological macromolecules and their complexes, as well as the interactions of membrane protein and lipid complexes at a reduced level of representation, allowing longer and larger molecular dynamics simulations. Here, we present a computational platform dedicated to the preparation, running, and analysis of CGMD simulations. The platform is built on a completely revisited version of our Martini coarsE gRained MembrAne proteIn Dynamics (MERMAID) web server, and it integrates this with other three dedicated services. In its current version, the platform expands the existing implementation of the Martini force field for membrane proteins to also allow the simulation of soluble proteins using the Martini and the SIRAH force fields. Moreover, it offers an automated protocol for carrying out the backmapping of the coarse-grained description of the system into an atomistic one.

Towards design of drugs and delivery systems with the Martini coarse-grained model.

Coarse-grained (CG) modelling with the Martini force field has come of age. By combining a variety of bead types and sizes with a new mapping approach, the newest version of the model is able to accurately simulate large biomolecular complexes at millisecond timescales. In this perspective, we discuss possible applications of the Martini 3 model in drug discovery and development pipelines and highlight areas for future development. Owing to its high simulation efficiency and extended chemical space, Martini 3 has great potential in the area of drug design and delivery. However, several aspects of the model should be improved before Martini 3 CG simulations can be routinely employed in academic and industrial settings. These include the development of automatic parameterisation protocols for a variety of molecule types, the improvement of backmapping procedures, the description of protein flexibility and the development of methodologies enabling efficient sampling. We illustrate our view with examples on key areas where Martini could give important contributions such as drugs targeting membrane proteins, cryptic pockets and protein-protein interactions and the development of soft drug delivery systems.

A Novel ALAS2 Missense Mutation in Two Brothers With Iron Overload and Associated Alterations in Serum Hepcidin/Erythroferrone Levels.

Iron loading anemias are characterized by ineffective erythropoiesis and iron overload. The prototype is non-transfusion dependent ß-thalassemia (NTDT), with other entities including congenital sideroblastic anemias, congenital dyserythropoietic anemias, some hemolytic anemias, and myelodysplastic syndromes. Differential diagnosis of iron loading anemias may be challenging due to heterogeneous genotype and phenotype. Notwithstanding the recent advances in linking ineffective erythropoiesis to iron overload, many pathophysiologic aspects are still unclear. Moreover, measurement of hepcidin and erythroferrone (ERFE), two key molecules in iron homeostasis and erythropoiesis, is scarcely used in clinical practice and of uncertain utility. Here, we describe a comprehensive diagnostic approach, including next-generation sequencing (NGS), in silico modeling, and measurement of hepcidin and erythroferrone (ERFE), in two brothers eventually diagnosed as X-linked sideroblastic anemia (XLSA). A novel pathogenic ALAS2 missense mutation (c.1382T>A, p.Leu461His) is described. Hyperferritinemia with high hepcidin-25 levels (but decreased hepcidin:ferritin ratio) and mild-to-moderate iron overload were detected in both patients. ERFE levels were markedly elevated in both patients, especially in the proband, who had a more expressed phenotype. Our study illustrates how new technologies, such as NGS, in silico modeling, and measurement of serum hepcidin-25 and ERFE, may help in diagnosing and studying iron loading anemias. Further studies on the hepcidin-25/ERFE axis in additional patients with XLSA and other iron loading anemias may help in establishing its usefulness in differential diagnosis, and it may also aid our understanding of the pathophysiology of these genetically and phenotypically heterogeneous entities.