Profile and potential bioactivity of the miRNome and metabolome expressed in Malva sylvestris L. leaf and flower

Malva sylvestris L. (common mallow) is a plant species widely used in phytotherapy and ethnobotanical practices since time immemorial. Characterizing the components of this herb might promote a better comprehension of its biological effects on the human body but also favour the identification of the molecular processes that occur in the plant tissues. Thus, in the present contribution, the scientific knowledge about the metabolomic profile of the common mallow was expanded. In particular, the phytocomplex of leaves and flowers from this botanical species and the extraction capacity of different concentrations of ethanol (i.e., 95%, 70%, 50%, and 0%; v/v in ddH2O) for it were investigated by spectrophotometric and chromatographic approaches. In detail, 95% ethanol extracts showed the worst capacity in isolating total phenols and flavonoids, while all the hydroalcoholic samples revealed a specific ability in purifying the anthocyanins. HPLC–DAD system detected and quantified 20 phenolic secondary metabolites, whose concentration in the several extracts depended on their own chemical nature and the percentage of ethanol used in the preparation. In addition, the stability of the purified phytochemicals after resuspension in pure ddH2O was also proved, considering a potential employment of them in biological/medical studies which include in vitro and in vivo experiments on mammalian models. Here, for the first time, the expressed miRNome in M. sylvestris was also defined by Next Generation Sequencing, revealing the presence of 33 microRNAs (miRNAs), 10 typical for leaves and 2 for flowers. Then, both plant and human putative mRNA targets for the detected miRNAs were predicted by bioinformatics analyses, with the aim to clarify the possible role of these small nucleic acids in the common mallow plant tissues and to try to understand if they could exert a potential cross-kingdom regulatory activity on the human health. Surprisingly, our investigations revealed that 19 miRNAs out of 33 were putatively able to modulate, in the plant cells, the expression of various chromosome scaffold proteins. In parallel, we found, in the human transcriptome, a total of 383 mRNAs involved in 5 fundamental mammalian cellular processes (i.e., apoptosis, senescence, cell-cycle, oxidative stress, and invasiveness) that theoretically could be bound and regulated by M. sylvestris miRNAs. The evidence collected in this work would suggest that the beneficial properties of the use of M. sylvestris, documented by the folk medicine, are probably linked to their content of miRNAs and not only to the action of phytochemicals (e.g., anthocyanins). This would open new perspectives about the possibility to develop gene therapies based on miRNAs isolated from medicinal plants, including M. sylvestris

Holo-like and Druggable Protein Conformations from Enhanced Sampling of Binding Pocket Volume and Shape

Understanding molecular recognition of small molecules by proteins in atomistic detail is key for drug design. Molecular docking is a widely used computational method to mimic ligand-protein association in silico. However, predicting conformational changes occurring in proteins upon ligand binding is still a major challenge. Ensemble docking approaches address this issue by considering a set of different conformations of the protein obtained either experimentally or from computer simulations, e.g., molecular dynamics. However, holo structures prone to host (the correct) ligands are generally poorly sampled by standard molecular dynamics simulations of the apo protein. In order to address this limitation, we introduce a computational approach based on metadynamics simulations called ensemble docking with enhanced sampling of pocket shape (EDES) that allows holo-like conformations of proteins to be generated by exploiting only their apo structures. This is achieved by defining a set of collective variables that effectively sample different shapes of the binding site, ultimately mimicking the steric effect due to the ligand. We assessed the method on three challenging proteins undergoing different extents of conformational changes upon ligand binding. In all cases our protocol generates a significant fraction of structures featuring a low RMSD from the experimental holo geometry. Moreover, ensemble docking calculations using those conformations yielded in all cases native-like poses among the top-ranked ones

Intestinal Taxa Abundance and Diversity in Inflammatory Bowel Disease Patients: An Analysis including Covariates and Confounders

Intestinal dysbiosis has been widely documented in inflammatory bowel diseases (IBDs) and is thought to influence the onset and perpetuation of gut inflammation. However, it remains unclear whether such bacterial changes rely in part on the modification of an IBD-associated lifestyle (e.g., smoking and physical activity) and diet (e.g., rich in dairy products, cereals, meat and vegetables). In this study, we investigated the impact of these habits, which we defined as confounders and covariates, on the modulation of intestinal taxa abundance and diversity in IBD patients. 16S rRNA gene sequence analysis was performed using genomic DNA extracted from the faecal samples of 52 patients with Crohn's disease (CD) and 58 with ulcerative colitis (UC), which are the two main types of IBD, as well as 42 healthy controls (HC). A reduced microbial diversity was documented in the IBD patients compared with the HC. Moreover, we identified specific confounders and covariates that influenced the association between some bacterial taxa and disease extent (in UC patients) or behaviour (in CD patients) compared with the HC. In particular, a PERMANOVA stepwise regression identified the variables "age", "eat yogurt at least four days per week" and "eat dairy products at least 4 days per week" as covariates when comparing the HC and patients affected by ulcerative proctitis (E1), left-sided UC (distal UC) (E2) and extensive UC (pancolitis) (E3). Instead, the variables "age", "gender", "eat meat at least four days per week" and "eat bread at least 4 days per week" were considered as covariates when comparing the HC with the CD patients affected by non-stricturing, non-penetrating (B1), stricturing (B2) and penetrating (B3) diseases. Considering such variables, our analysis indicated that the UC extent differentially modulated the abundance of the Bifidobacteriaceae, Rikenellaceae, Christensenellaceae, Marinifilaceae, Desulfovibrionaceae, Lactobacillaceae, Streptococcaceae and Peptostreptococcaceae families, while the CD behaviour influenced the abundance of Christensenellaceae, Marinifilaceae, Rikenellaceae, Ruminococcaceae, Barnesiellaceae and Coriobacteriaceae families. In conclusion, our study indicated that some covariates and confounders related to an IBD-associated lifestyle and dietary habits influenced the intestinal taxa diversity and relative abundance in the CD and UC patients compared with the HC. Indeed, such variables should be identified and excluded from the analysis to characterize the bacterial families whose abundance is directly modulated by IBD status, as well as disease extent or behaviour

Ab initio study of the vibrational properties of crystalline TeO2: The alpha, beta, and gamma phases

Based on density functional perturbation theory, we have studied the vibrational properties of three crystalline phases of tellurium dioxide: paratellurite alpha-TeO2, tellurite beta-TeO2, and the new phase gamma-TeO2, recently identified experimentally. Calculated Raman and IR spectra are in good agreement with available experimental data. The vibrational spectra of alpha- and beta-TeO2 can be interpreted in terms of vibrations of TeO2 molecular units

Exploring the Universe of Protein Structures beyond the Protein Data Bank

It is currently believed that the atlas of existing protein structures is faithfully represented in the Protein Data Bank. However, whether this atlas covers the full universe of all possible protein structures is still a highly debated issue. By using a sophisticated numerical approach, we performed an exhaustive exploration of the conformational space of a 60 amino acid polypeptide chain described with an accurate all-atom interaction potential. We generated a database of around 30,000 compact folds with at least of secondary structure corresponding to local minima of the potential energy. This ensemble plausibly represents the universe of protein folds of similar length; indeed, all the known folds are represented in the set with good accuracy. However, we discover that the known folds form a rather small subset, which cannot be reproduced by choosing random structures in the database. Rather, natural and possible folds differ by the contact order, on average significantly smaller in the former. This suggests the presence of an evolutionary bias, possibly related to kinetic accessibility, towards structures with shorter loops between contacting residues. Beside their conceptual relevance, the new structures open a range of practical applications such as the development of accurate structure prediction strategies, the optimization of force fields, and the identification and design of novel folds

Accurate multiple time step in biased molecular simulations

Many recently introduced enhanced sampling techniques are based on biasing coarse descriptors (collective variables) of a molecular system on the fly. Sometimes the calculation of such collective variables is expensive and becomes a bottleneck in molecular dynamics simulations. An algorithm to treat smooth biasing forces within a multiple time step framework is here discussed. The implementation is simple and allows a speed up when expensive collective variables are employed. The gain can be substantial when using massively parallel or GPU-based molecular dynamics software. Moreover, a theoretical framework to assess the sampling accuracy is introduced, which can be used to assess the choice of the integration time step in both single and multiple time step biased simulations

The inverted free energy landscape of an intrinsically disordered peptide by simulations and experiments

The free energy landscape theory has been very successful in rationalizing the folding behaviour of globular proteins, as this representation provides intuitive information on the number of states involved in the folding process, their populations and pathways of interconversion. We extend here this formalism to the case of the A\u3b240 peptide, a 40-residue intrinsically disordered protein fragment associated with Alzheimer's disease. By using an advanced sampling technique that enables free energy calculations to reach convergence also in the case of highly disordered states of proteins, we provide a precise structural characterization of the free energy landscape of this peptide. We find that such landscape has inverted features with respect to those typical of folded proteins. While the global free energy minimum consists of highly disordered structures, higher free energy regions correspond to a large variety of transiently structured conformations with secondary structure elements arranged in several different manners, and are not separated from each other by sizeable free energy barriers. From this peculiar structure of the free energy landscape we predict that this peptide should become more structured and not only more compact, with increasing temperatures, and we show that this is the case through a series of biophysical measurements

Multiple Routes and Milestones in the Folding of HIV–1 Protease Monomer

Proteins fold on a time scale incompatible with a mechanism of random search in conformational space thus indicating that somehow they are guided to the native state through a funneled energetic landscape. At the same time the heterogeneous kinetics suggests the existence of several different folding routes. Here we propose a scenario for the folding mechanism of the monomer of HIV–1 protease in which multiple pathways and milestone events coexist. A variety of computational approaches supports this picture. These include very long all-atom molecular dynamics simulations in explicit solvent, an analysis of the network of clusters found in multiple high-temperature unfolding simulations and a complete characterization of free-energy surfaces carried out using a structure-based potential at atomistic resolution and a combination of metadynamics and parallel tempering. Our results confirm that the monomer in solution is stable toward unfolding and show that at least two unfolding pathways exist. In our scenario, the formation of a hydrophobic core is a milestone in the folding process which must occur along all the routes that lead this protein towards its native state. Furthermore, the ensemble of folding pathways proposed here substantiates a rational drug design strategy based on inhibiting the folding of HIV–1 protease

Strength of Hydrogen Bond Network Takes Crucial Roles in the Dissociation Process of Inhibitors from the HIV-1 Protease Binding Pocket

To understand the underlying mechanisms of significant differences in dissociation rate constant among different inhibitors for HIV-1 protease, we performed steered molecular dynamics (SMD) simulations to analyze the entire dissociation processes of inhibitors from the binding pocket of protease at atomistic details. We found that the strength of hydrogen bond network between inhibitor and the protease takes crucial roles in the dissociation process. We showed that the hydrogen bond network in the cyclic urea inhibitors AHA001/XK263 is less stable than that of the approved inhibitor ABT538 because of their large differences in the structures of the networks. In the cyclic urea inhibitor bound complex, the hydrogen bonds often distribute at the flap tips and the active site. In contrast, there are additional accessorial hydrogen bonds formed at the lateral sides of the flaps and the active site in the ABT538 bound complex, which take crucial roles in stabilizing the hydrogen bond network. In addition, the water molecule W301 also plays important roles in stabilizing the hydrogen bond network through its flexible movement by acting as a collision buffer and helping the rebinding of hydrogen bonds at the flap tips. Because of its high stability, the hydrogen bond network of ABT538 complex can work together with the hydrophobic clusters to resist the dissociation, resulting in much lower dissociation rate constant than those of cyclic urea inhibitor complexes. This study may provide useful guidelines for design of novel potent inhibitors with optimized interactions