Metabolomics, similarly to other high-throughput “-omics” techniques, generates large arrays of data, whose analysis and interpretation can be difficult and not always straightforward. Several software for the detailed metabolomics statistical analysis are available, however there is a lack of simple protocols guiding the user through a standard statistical analysis of the data. Herein we present “muma”, an R package providing a simple step-wise pipeline for metabolomics univariate and multivariate statistical analyses. Based on published statistical algorithms and techniques, muma provides user-friendly tools for the whole process of data analysis, ranging from data imputation and preprocessing, to dataset exploration, to data interpretation through unsupervised/supervised multivariate and/or univariate techniques. Of note, specific tools and graphics aiding the explanation of statistical outcomes have been developed. Finally, a section dedicated to metabolomics data interpretation has been implemented, providing specific techniques for molecular assignments and biochemical interpretation of metabolic patterns. muma is a free, user-friendly and versatile tool suite tailored to assist the user in the interpretation of metabolomics data in the identification of biomarkers and in the analysis of metabolic patterns.

PHD fingers represent one of the largest families of epigenetic readers capable of decoding post-translationally modified or unmodified histone H3 tails. Because of their direct involvement in human pathologies they are increasingly considered as a potential therapeutic target. Several PHD/histone-peptide structures have been determined, however relatively little information is available on their dynamics. Studies aiming to characterize the dynamic and energetic determinants driving histone peptide recognition by epigenetic readers would strongly benefit from computational studies. Herein we focus on the dynamic and energetic characterization of the PHD finger subclass specialized in the recognition of histone H3 peptides unmodified in position K4 (H3K4me0). As a case study we focused on the first PHD finger of autoimmune regulator protein (AIRE-PHD1) in complex with H3K4me0. PCA analysis of the covariance matrix of free AIRE-PHD1 highlights the presence of a “flapping” movement, which is blocked in an open conformation upon binding to H3K4me0. Moreover, binding free energy calculations obtained through Molecular Mechanics/Poisson-Boltzmann Surface Area (MM/PBSA) methodology are in good qualitative agreement with experiments and allow dissection of the energetic terms associated with native and alanine mutants of AIRE-PHD1/H3K4me0 complexes. MM/PBSA calculations have also been applied to the energetic analysis of other PHD fingers recognizing H3K4me0. In this case we observe excellent correlation between computed and experimental binding free energies. Overall calculations show that H3K4me0 recognition by PHD fingers relies on compensation of the electrostatic and polar solvation energy terms and is stabilized by non-polar interactions.

Asparagine deamidation at the NGR sequence in the 5th type I repeat of fibronectin (FN-I5) generates isoDGR, an alphavbeta3 integrin-binding motif regulating endothelial cell adhesion and proliferation. By NMR and molecular dynamics studies, we analyzed the structure of CisoDGRC (isoDGR-2C), a cyclic beta-peptide mimicking the FN-I5 site, and compared it with NGR, RGD, or DGR-containing cyclopeptides. Docking experiments show that isoDGR, exploiting an inverted orientation as compared with RGD, favorably interacts with the RGD-binding site of alphavbeta3, both recapitulating canonical RGD-alphavbeta3 contacts and establishing additional polar interactions. Conversely, NGR and DGR motifs lack the fundamental pharmacophoric requirements for high receptor affinity. Therefore, unlike NGR and DGR, isoDGR is a new natural recognition motif of the RGD-binding pocket of alphavbeta3. These findings contribute to explain the different functional properties of FN-I5 before and after deamidation, and provide support for the hypothesis that NGR --> isoDGR transition can work as a molecular timer for activating latent integrin-binding sites in proteins, thus regulating protein function.