Vibrational spectroscopy coupled with machine learning sheds light on the cellular effects induced by rationally designed TLR4 agonists.

In this work, we present the potential of Fourier transform infrared (FTIR) microspectroscopy to compare on whole cells, in an unbiased and untargeted way, the capacity of bacterial lipopolysaccharide (LPS) and two rationally designed molecules (FP20 and FP20Rha) to activate molecular circuits of innate immunity. These compounds are important drug hits in the development of vaccine adjuvants and tumor immunotherapeutics. The biological assays indicated that FP20Rha was more potent than FP20 in inducing cytokine production in cells and in stimulating IgG antibody production post-vaccination in mice. Accordingly, the overall significant IR spectral changes induced by the treatment with LPS and FP20Rha were similar, lipids and glycans signals being the most diagnostic, while the effect of the less potent molecule FP20 on cells resulted to be closer to control untreated cells. We propose here the use of FTIR spectroscopy supported by artificial intelligence (AI) to achieve a more holistic understanding of the cell response to new drug candidates while screening them in cells.

A Global Picture of Molecular Changes Associated to LPS Treatment in THP-1 Derived Human Macrophages by Fourier Transform Infrared Microspectroscopy.

Macrophages are among the first immune cells involved in the initiation of the inflammatory response to protect the host from pathogens. THP-1 derived macrophages (TDM) are used as a model to study the pro-inflammatory effects of lipopolysaccharide (LPS) exposure. Intact TDM cells were analysed by Fourier transform infrared (FTIR) microspectroscopy, supported by multivariate analysis, to obtain a snapshot of the molecular events sparked by LPS stimulation in macrophage-like cells. This spectroscopic analysis enabled the untargeted identification of the most significant spectral components affected by the treatment, ascribable mainly to lipid, protein, and sulfated sugar bands, thus stressing the fundamental role of these classes of molecules in inflammation and in immune response. Our study, therefore, shows that FTIR microspectroscopy enabled the identification of spectroscopic markers of LPS stimulation and has the potential to become a tool to assess those global biochemical changes related to inflammatory and anti-inflammatory stimuli of synthetic and natural immunomodulators different from LPS.

Contribution of Infrared Spectroscopy to the Understanding of Amyloid Protein Aggregation in Complex Systems.

Infrared (IR) spectroscopy is a label-free and non-invasive technique that probes the vibrational modes of molecules, thus providing a structure-specific spectrum. The development of infrared spectroscopic approaches that enable the collection of the IR spectrum from a selected sample area, from micro- to nano-scale lateral resolutions, allowed to extend their application to more complex biological systems, such as intact cells and tissues, thus exerting an enormous attraction in biology and medicine. Here, we will present recent works that illustrate in particular the applications of IR spectroscopy to the in situ characterization of the conformational properties of protein aggregates and to the investigation of the other biomolecules surrounding the amyloids. Moreover, we will discuss the potential of IR spectroscopy to the monitoring of cell perturbations induced by protein aggregates. The essential support of multivariate analyses to objectively pull out the significant and non-redundant information from the spectra of highly complex systems will be also outlined.

Tear-Based Vibrational Spectroscopy Applied to Amyotrophic Lateral Sclerosis.

Biofluid analysis by optical spectroscopy techniques is attracting considerable interest due to its potential to revolutionize diagnostics and precision medicine, particularly for neurodegenerative diseases. However, the lack of effective biomarkers combined with the unaccomplished identification of convenient biofluids has drastically hampered optical advancements in clinical diagnosis and monitoring of neurodegenerative disorders. Here, we show that vibrational spectroscopy applied to human tears opens a new route, offering a non-invasive, label-free identification of a devastating disease such as amyotrophic lateral sclerosis (ALS). Our proposed approach has been validated using two widespread techniques, namely, Fourier transform infrared (FTIR) and Raman microspectroscopies. In conjunction with multivariate analysis, this vibrational approach made it possible to discriminate between tears from ALS patients and healthy controls (HCs) with high specificity (∼97% and ∼100% for FTIR and Raman spectroscopy, respectively) and sensitivity (∼88% and ∼100% for FTIR and Raman spectroscopy, respectively). Additionally, the investigation of tears allowed us to disclose ALS spectroscopic markers related to protein and lipid alterations, as well as to a reduction of the phenylalanine level, in comparison with HCs. Our findings show that vibrational spectroscopy is a new potential ALS diagnostic approach and indicate that tears are a reliable and non-invasive source of ALS biomarkers.

Pathological ATX3 Expression Induces Cell Perturbations in E. coli as Revealed by Biochemical and Biophysical Investigations.

Amyloid aggregation of human ataxin-3 (ATX3) is responsible for spinocerebellar ataxia type 3, which belongs to the class of polyglutamine neurodegenerative disorders. It is widely accepted that the formation of toxic oligomeric species is primarily involved in the onset of the disease. For this reason, to understand the mechanisms underlying toxicity, we expressed both a physiological (ATX3-Q24) and a pathological ATX3 variant (ATX3-Q55) in a simplified cellular model, Escherichia coli. It has been observed that ATX3-Q55 expression induces a higher reduction of the cell growth compared to ATX3-Q24, due to the bacteriostatic effect of the toxic oligomeric species. Furthermore, the Fourier transform infrared microspectroscopy investigation, supported by multivariate analysis, made it possible to monitor protein aggregation and the induced cell perturbations in intact cells. In particular, it has been found that the toxic oligomeric species associated with the expression of ATX3-Q55 are responsible for the main spectral changes, ascribable mainly to the cell envelope modifications. A structural alteration of the membrane detected through electron microscopy analysis in the strain expressing the pathological form supports the spectroscopic results.

ATR-FTIR Spectroscopy Supported by Multivariate Analysis for the Characterization of Adipose Tissue Aspirates from Patients Affected by Systemic Amyloidosis.

Deposition of misfolded proteins as extracellular amyloid aggregates is the pathological hallmark of systemic amyloidoses. Subcutaneous fat acquired by fine needle aspiration is the preferred screening tissue in suspected patients. In this study we employed Fourier transform infrared (FTIR) spectroscopy in attenuated total reflection (ATR) to investigate human abdominal fat aspirates with the aim of detecting disease-related changes in the molecular structure and composition of the tissue and exploiting the potentiality of the method to discriminate between amyloid-positive and -negative samples. The absorption and second-derivative spectra of Congo Red (CR) positive and CR-negative specimens were analyzed by three multivariate methods in four spectral regions. The proposed ATR-FTIR method is label-free, rapid, and relatively inexpensive and requires minimal sample preparation. We found that the ATR-FTIR approach can differentiate fat aspirates containing amyloid deposits from control specimens with high sensitivity and specificity, both at 100 [89-100]%. It is worth noting that the wavenumbers most important for discrimination indicate that changes both in the protein conformation and in resident lipids are intrinsic features of affected subcutaneous fat in comparison with the CR-negative controls. In this proof of concept study, we show that this approach could be useful for assessing tissue amyloid aggregates and for acquiring novel knowledge of the molecular bases of the disease.

Correction to: Long range Debye-Hückel correction for computation of grid-based electrostatic forces between biomacromolecules.

[This corrects the article DOI: 10.1186/2046-1682-7-4.].

A FTIR microspectroscopy study of the structural and biochemical perturbations induced by natively folded and aggregated transthyretin in HL-1 cardiomyocytes.

Protein misfolding and aggregation are associated with a number of human degenerative diseases. In spite of the enormous research efforts to develop effective strategies aimed at interfering with the pathogenic cascades induced by misfolded/aggregated peptides/proteins, the necessary detailed understanding of the molecular bases of amyloid formation and toxicity is still lacking. To this aim, approaches able to provide a global insight in amyloid-mediated physiological alterations are of importance. In this study, we exploited Fourier transform infrared microspectroscopy, supported by multivariate analysis, to investigate in situ the spectral changes occurring in cultured intact HL-1 cardiomyocytes exposed to wild type (WT) or mutant (L55P) transthyretin (TTR) in native, or amyloid conformation. The presence of extracellular deposits of amyloid aggregates of WT or L55P TTR, respectively, is a key hallmark of two pathological conditions, known as senile systemic amyloidosis and familial amyloid polyneuropathy. We found that the major effects, associated with modifications in lipid properties and in the cell metabolic/phosphorylation status, were observed when natively folded WT or L55P TTR was administered to the cells. The effects induced by aggregates of TTR were milder and in some cases displayed a different timing compared to those elicited by the natively folded protein.

Rational Design of a Transferrin-Binding Peptide Sequence Tailored to Targeted Nanoparticle Internalization.

The transferrin receptor (TfR) is a promising target in cancer therapy owing to its overexpression in most solid tumors and on the blood-brain barrier. Nanostructures chemically derivatized with transferrin are employed in TfR targeting but often lose their functionality upon injection in the bloodstream. As an alternative strategy, we rationally designed a peptide coating able to bind transferrin on suitable pockets not involved in binding to TfR or iron by using an iterative multiscale-modeling approach coupled with quantitative structure-activity and relationship (QSAR) analysis and evolutionary algorithms. We tested that selected sequences have low aspecific protein adsorption and high binding energy toward transferrin, and one of them is efficiently internalized in cells with a transferrin-dependent pathway. Furthermore, it promotes transferrin-mediated endocytosis of gold nanoparticles by modifying their protein corona and promoting oriented adsorption of transferrin. This strategy leads to highly effective nanostructures, potentially useful in diagnostic and therapeutic applications, which exploit (and do not suffer) the protein solvation for achieving a better targeting.

Induced Fit in Protein Multimerization: The HFBI Case.

Hydrophobins, produced by filamentous fungi, are small amphipathic proteins whose biological functions rely on their unique surface-activity properties. Understanding the mechanistic details of the multimerization process is of primary importance to clarify the interfacial activity of hydrophobins. We used free energy calculations to study the role of a flexible ß-hairpin in the multimerization process in hydrophobin II from Trichoderma reesei (HFBI). We characterized how the displacement of this ß-hairpin controls the stability of the monomers/dimers/tetramers in solution. The regulation of the oligomerization equilibrium of HFBI will necessarily affect its interfacial properties, fundamental for its biological function and for technological applications. Moreover, we propose possible routes for the multimerization process of HFBI in solution. This is the first case where a mechanism by which a flexible loop flanking a rigid patch controls the protein-protein binding equilibrium, already known for proteins with charged binding hot-spots, is described within a hydrophobic patch.

Assessing an effective feeding strategy to optimize crude glycerol utilization as sustainable carbon source for lipid accumulation in oleaginous yeasts.

BACKGROUND: Microbial lipids can represent a valuable alternative feedstock for biodiesel production in the context of a viable bio-based economy. This production can be driven by cultivating some oleaginous microorganisms on crude-glycerol, a 10% (w/w) by-product produced during the transesterification process from oils into biodiesel. Despite attractive, the perspective is still economically unsustainable, mainly because impurities in crude glycerol can negatively affect microbial performances. In this view, the selection of the best cell factory, together with the development of a robust and effective production process are primary requirements. RESULTS: The present work compared crude versus pure glycerol as carbon sources for lipid production by three different oleaginous yeasts: Rhodosporidium toruloides (DSM 4444), Lipomyces starkeyi (DSM 70295) and Cryptococcus curvatus (DSM 70022). An efficient yet simple feeding strategy for avoiding the lag phase caused by growth on crude glycerol was developed, leading to high biomass and lipid production for all the tested yeasts. Flow-cytometry and fourier transform infrared (FTIR) microspectroscopy, supported by principal component analysis (PCA), were used as non-invasive and quick techniques to monitor, compare and analyze the lipid production over time. Gas chromatography (GC) analysis completed the quali-quantitative description. Under these operative conditions, the highest lipid content (up to 60.9% wt/wt) was measured in R. toruloides, while L. starkeyi showed the fastest glycerol consumption rate (1.05 g L(-1) h(-1)). Being productivity the most industrially relevant feature to be pursued, under the presented optimized conditions R. toruloides showed the best lipid productivity (0.13 and 0.15 g L(-1) h(-1) on pure and crude glycerol, respectively). CONCLUSIONS: Here we demonstrated that the development of an efficient feeding strategy is sufficient in preventing the inhibitory effect of crude glycerol, and robust enough to ensure high lipid accumulation by three different oleaginous yeasts. Single cell and in situ analyses allowed depicting and comparing the transition between growth and lipid accumulation occurring differently for the three different yeasts. These data provide novel information that can be exploited for screening the best cell factory, moving towards a sustainable microbial biodiesel production.

Optimization of Analytical Potentials for Coarse-Grained Biopolymer Models.

The increasing trend in the recent literature on coarse grained (CG) models testifies their impact in the study of complex systems. However, the CG model landscape is variegated: even considering a given resolution level, the force fields are very heterogeneous and optimized with very different parametrization procedures. Along the road for standardization of CG models for biopolymers, here we describe a strategy to aid building and optimization of statistics based analytical force fields and its implementation in the software package AsParaGS (Assisted Parameterization platform for coarse Grained modelS). Our method is based on the use and optimization of analytical potentials, optimized by targeting internal variables statistical distributions by means of the combination of different algorithms (i.e., relative entropy driven stochastic exploration of the parameter space and iterative Boltzmann inversion). This allows designing a custom model that endows the force field terms with a physically sound meaning. Furthermore, the level of transferability and accuracy can be tuned through the choice of statistical data set composition. The method-illustrated by means of applications to helical polypeptides-also involves the analysis of two and three variable distributions, and allows handling issues related to the FF term correlations. AsParaGS is interfaced with general-purpose molecular dynamics codes and currently implements the "minimalist" subclass of CG models (i.e., one bead per amino acid, Cα based). Extensions to nucleic acids and different levels of coarse graining are in the course.

Long range Debye-Hückel correction for computation of grid-based electrostatic forces between biomacromolecules.

BACKGROUND: Brownian dynamics (BD) simulations can be used to study very large molecular systems, such as models of the intracellular environment, using atomic-detail structures. Such simulations require strategies to contain the computational costs, especially for the computation of interaction forces and energies. A common approach is to compute interaction forces between macromolecules by precomputing their interaction potentials on three-dimensional discretized grids. For long-range interactions, such as electrostatics, grid-based methods are subject to finite size errors. We describe here the implementation of a Debye-Hückel correction to the grid-based electrostatic potential used in the SDA BD simulation software that was applied to simulate solutions of bovine serum albumin and of hen egg white lysozyme. RESULTS: We found that the inclusion of the long-range electrostatic correction increased the accuracy of both the protein-protein interaction profiles and the protein diffusion coefficients at low ionic strength. CONCLUSIONS: An advantage of this method is the low additional computational cost required to treat long-range electrostatic interactions in large biomacromolecular systems. Moreover, the implementation described here for BD simulations of protein solutions can also be applied in implicit solvent molecular dynamics simulations that make use of gridded interaction potentials.

A Fourier transform infrared spectroscopy study of cell membrane domain modifications induced by docosahexaenoic acid.

BACKGROUND: Detergent resistant membranes (DRMs) are a useful model system for the in vitro characterization of cell membrane domains. Indeed, DRMs provide a simple model to study the mechanisms underlying several key cell processes based on the interplay between specific cell membrane domains on one hand, and specific proteins and/or lipids on the other. Considering therefore their biological relevance, the development of methods capable to provide information on the composition and structure of membrane domains and to detect their modifications is highly desirable. In particular, Fourier transform infrared (FTIR) spectroscopy is a vibrational tool widely used for the study not only of isolated and purified biomolecules but also of complex biological systems, including intact cells and tissues. One of the main advantages of this non-invasive approach is that it allows obtaining a molecular fingerprint of the sample under investigation in a rapid and label-free way. METHODS: Here we present an FTIR characterization of DRM fractions purified from the human breast cancer cells MCF-7, before and after treatment with the omega 3 fatty acid docosahexaenoic acid (DHA), which was found to promote membrane microdomain reorganization. RESULTS AND CONCLUSIONS: We will show that FTIR spectroscopy coupled with multivariate analysis enables to monitor changes in the composition of DRMs, induced in particular by the incorporation of DHA in cell membrane phospholipids. GENERAL SIGNIFICANCE: This study paves the way for a new label-free characterization of specific membrane domains within intact cells, which could provide complementary information to the fluorescence approaches presently used.

Two interconvertible folds modulate the activity of a DNA aptamer against transferrin receptor.

Thanks to their ability to recognize biomolecular targets with high affinity and specificity, nucleic acid aptamers are increasingly investigated as diagnostic and therapeutic tools, particularly when their targets are cell-surface receptors. Here, we investigate the relationship between the folding of an anti-mouse transferrin receptor DNA aptamer and its interaction with the transferrin receptor both in vitro and in living cells. We identified and purified two aptamer conformers by means of chromatographic techniques. Fluorescence-anisotropy measurements showed that only one fold is able to bind mouse transferrin receptor. Besides displaying enhanced endocytosis in living mouse fibroblasts, the purified active fold is internalized also in human pancreatic cancer cells. Starting from these observations, we rationally designed variations of the parent sequence aimed at stabilizing the active fold, and consequently increase aptamer activity. A truncated version and full-length mutants with higher affinity than the parent sequence are shown.Molecular Therapy-Nucleic Acids (2014) 3, e144; doi:10.1038/mtna.2013.71; published online 28 January 2014.

Fourier transform infrared spectroscopy as a method to study lipid accumulation in oleaginous yeasts.

BACKGROUND: Oleaginous microorganisms, such as different yeast and algal species, can represent a sustainable alternative to plant oil for the production of biodiesel. They can accumulate fatty acids (FA) up to 70% of their dry weight with a predominance of (mono)unsaturated species, similarly to what plants do, but differently from animals. In addition, their growth is not in competition either with food, feed crops, or with agricultural land.Despite these advantages, the exploitation of the single cell oil system is still at an early developmental stage. Cultivation mode and conditions, as well as lipid extraction technologies, represent the main limitations. The monitoring of lipid accumulation in oleaginous microorganisms is consequently crucial to develop and validate new approaches, but at present the majority of the available techniques is time consuming, invasive and, when relying on lipid extraction, can be affected by FA degradation. RESULTS: In this work the fatty acid accumulation of the oleaginous yeasts Cryptococcus curvatus and Rhodosporidium toruloides and of the non-oleaginous yeast Saccharomyces cerevisiae (as a negative control) was monitored in situ by Fourier Transform Infrared Spectroscopy (FTIR). Indeed, this spectroscopic tool can provide complementary information to those obtained by classical techniques, such as microscopy, flow cytometry and gas chromatography. As shown in this work, through the analysis of the absorption spectra of intact oleaginous microorganisms it is possible not only to monitor the progression of FA accumulation but also to identify the most represented classes of the produced lipids. CONCLUSIONS: Here we propose FTIR microspectroscopy - supported by multivariate analysis - as a fast, reliable and non invasive method to monitor and analyze FA accumulation in intact oleaginous yeasts. The results obtained by the FTIR approach were in agreement with those obtained by the other classical methods like flow cytometry and gas chromatography. Moreover, the possibility to track lipid production in real time is highly desirable to support the initial screening of strains and media as well as to optimize the scaling up experiments, which are essential for a viable and successful development of an industrial production process.

The shape of protein crowders is a major determinant of protein diffusion.

As a model for understanding how molecular crowding influences diffusion and transport of proteins in cellular environments, we combined experimental and theoretical approaches to study the diffusion of proteins in highly concentrated protein solutions. Bovine serum albumin and γ-Globulin were chosen as molecular crowders and as tracers. These two proteins are representatives of the main types of plasma protein and have different shapes and sizes. Solutions consisting of one or both proteins were studied. The self-diffusion coefficients of the fluorescently labeled tracer proteins were measured by means of fluorescence correlation spectroscopy at a total protein concentration of up to 400 g/L. γ-Globulin is found to have a stronger influence as a crowder on the tracer self-diffusion coefficient than Bovine serum albumin. Brownian dynamics simulations show that the excluded volume and the shape of the crowding protein have a significantly stronger influence on translational and rotational diffusion coefficients, as well as transient oligomerization, than hydrodynamic or direct interactions. Anomalous subdiffusion, which is not observed at the experimental fluorescence correlation spectroscopy timescales (>100 µs), appears only at very short timescales (<1 µs) in the simulations due to steric effects of the proteins. We envision that the combined experimental and computational approach employed here can be developed to unravel the different biophysical contributions to protein motion and interaction in cellular environments by systematically varying protein properties such as molecular weight, size, shape, and electrostatic interactions.

MELAS-like encephalomyopathy caused by a new pathogenic mutation in the mitochondrial DNA encoded cytochrome c oxidase subunit I.

We report a 35-year-old woman presenting a stroke-like episode with transitory aphasia followed by generalized tonic-clonic seizures. She had severe hearing loss and suffered from frequent episodes of migraine. Although a brain MRI disclosed a T2-hyperintense lesion in the left parietal lobe, she had hardly any long-term sequela. Exercise intolerance, myalgias and limb-girdle muscle weakness indicated a slowly progressive myopathy. Extra-neurological features included short stature, and secondary amenorrhea with low gonadotropin levels, indicating secondary hypogonadism. However, she had three mutation-free, healthy children by ovarian stimulation. A muscle biopsy showed ragged-red, cytochrome c oxidase-negative fibers, and an isolated defect of cytochrome c oxidase activity in muscle mitochondria. Sequence analysis of muscle mtDNA revealed a previously unreported heteroplasmic m.6597C>A transversion in the MTCOI gene, encoding subunit I of cytochrome c oxidase, corresponding to p.Q232K aminoacid change. Analysis on transmitochondrial cybrids demonstrated that the mutation is indeed associated with COX deficiency, i.e. pathogenic.

Atomic detail brownian dynamics simulations of concentrated protein solutions with a mean field treatment of hydrodynamic interactions.

High macromolecular concentrations are a distinguishing feature of living organisms. Understanding how the high concentration of solutes affects the dynamic properties of biological macromolecules is fundamental for the comprehension of biological processes in living systems. In this paper, we describe the implementation of mean field models of translational and rotational hydrodynamic interactions into an atomically detailed many-protein brownian dynamics simulation method. Concentrated solutions (30-40% volume fraction) of myoglobin, hemoglobin A, and sickle cell hemoglobin S were simulated, and static structure factors, oligomer formation, and translational and rotational self-diffusion coefficients were computed. Good agreement of computed properties with available experimental data was obtained. The results show the importance of both solvent mediated interactions and weak protein-protein interactions for accurately describing the dynamics and the association properties of concentrated protein solutions. Specifically, they show a qualitative difference in the translational and rotational dynamics of the systems studied. Although the translational diffusion coefficient is controlled by macromolecular shape and hydrodynamic interactions, the rotational diffusion coefficient is affected by macromolecular shape, direct intermolecular interactions, and both translational and rotational hydrodynamic interactions.

Diffusion of hydrophobin proteins in solution and interactions with a graphite surface.

BACKGROUND: Hydrophobins are small proteins produced by filamentous fungi that have a variety of biological functions including coating of spores and surface adhesion. To accomplish these functions, they rely on unique interface-binding properties. Using atomic-detail implicit solvent rigid-body Brownian dynamics simulations, we studied the diffusion of HFBI, a class II hydrophobin from Trichoderma reesei, in aqueous solution in the presence and absence of a graphite surface. RESULTS: In the simulations, HFBI exists in solution as a mixture of monomers in equilibrium with different types of oligomers. The oligomerization state depends on the conformation of HFBI. When a Highly Ordered Pyrolytic Graphite (HOPG) layer is present in the simulated system, HFBI tends to interact with the HOPG layer through a hydrophobic patch on the protein. CONCLUSIONS: From the simulations of HFBI solutions, we identify a tetrameric encounter complex stabilized by non-polar interactions between the aliphatic residues in the hydrophobic patch on HFBI. After the formation of the encounter complex, a local structural rearrangement at the protein interfaces is required to obtain the tetrameric arrangement seen in HFBI crystals. Simulations performed with the graphite surface show that, due to a combination of a geometric hindrance and the interaction of the aliphatic sidechains with the graphite layer, HFBI proteins tend to accumulate close to the hydrophobic surface.