Wild Patagonian yeast improve the evolutionary potential of novel interspecific hybrid strains for lager brewing.

Lager yeasts are limited to a few strains worldwide, imposing restrictions on flavour and aroma diversity and hindering our understanding of the complex evolutionary mechanisms during yeast domestication. The recent finding of diverse S. eubayanus lineages from Patagonia offers potential for generating new lager yeasts with different flavour profiles. Here, we leverage the natural genetic diversity of S. eubayanus and expand the lager yeast repertoire by including three distinct Patagonian S. eubayanus lineages. We used experimental evolution and selection on desirable traits to enhance the fermentation profiles of novel S. cerevisiae x S. eubayanus hybrids. Our analyses reveal an intricate interplay of pre-existing diversity, selection on species-specific mitochondria, de-novo mutations, and gene copy variations in sugar metabolism genes, resulting in high ethanol production and unique aroma profiles. Hybrids with S. eubayanus mitochondria exhibited greater evolutionary potential and superior fitness post-evolution, analogous to commercial lager hybrids. Using genome-wide screens of the parental subgenomes, we identified genetic changes in IRA2, IMA1, and MALX genes that influence maltose metabolism, and increase glycolytic flux and sugar consumption in the evolved hybrids. Functional validation and transcriptome analyses confirmed increased maltose-related gene expression, influencing greater maltotriose consumption in evolved hybrids. This study demonstrates the potential for generating industrially viable lager yeast hybrids from wild Patagonian strains. Our hybridization, evolution, and mitochondrial selection approach produced hybrids with high fermentation capacity and expands lager beer brewing options.

Nakazawaea atacamensis f.a., sp. nov. a novel nonconventional fermentative ascomycetous yeast species from the Atacama Desert.

In this study, we describe Nakazawaea atacamensis f. a., sp. nov., a novel species obtained from Neltuma chilensis plant samples in Chile's hyperarid Atacama Desert. In total, three strains of N. atacamensis were obtained from independent N. chilensis samples (synonym Prosopis chilensis, Algarrobo). Two strains were obtained from bark samples, while the third strain was obtained from bark-exuded gum from another tree. The novel species was defined using molecular characteristics and subsequently characterized with respect to morphological, physiological, and biochemical properties. A neighbor-joining analysis using the sequences of the D1/D2 domains of the large subunit ribosomal RNA gene revealed that N. atacamensis clustered with Nakazawaea pomicola. The sequence of N. atacamensis differed from closely related species by 1.3%-5.2% in the D1/D2 domains. A phylogenomic analysis based on single-nucleotide polymorphism's data confirms that the novel species belongs to the genus Nakazawaea, where N. atacamensis clustered with N. peltata. Phenotypic comparisons demonstrated that N. atacamensis exhibited distinct carbon assimilation patterns compared to its related species. Genome sequencing of the strain ATA-11A-BT revealed a genome size of approximately 12.4 Mbp, similar to other Nakazawaea species, with 5116 protein-coding genes annotated using InterProScan. In addition, N. atacamensis exhibited the capacity to ferment synthetic wine must, representing a potential new yeast for mono or co-culture wine fermentations. This comprehensive study expands our understanding of the genus Nakazawaea and highlights the ecological and industrial potential of N. atacamensis in fermentation processes. The holotype of N. atacamensis sp. nov. is CBS 18375T . The Mycobank number is MB 849680.

Late Pleistocene-dated divergence between South Hemisphere populations of the non-conventional yeast L. cidri.

Most organisms belonging to the Saccharomycotina subphylum have high genetic diversity and a vast repertoire of metabolisms and lifestyles. Lachancea cidri is an ideal yeast model for exploring the interplay between genetics, ecological function and evolution. Lachancea cidri diverged from the Saccharomyces lineage before the whole-genome duplication and is distributed across the South Hemisphere, displaying an important ecological success. We applied phylogenomics to investigate the genetic variation of L. cidri isolates obtained from Australia and South America. Our approach revealed the presence of two main lineages according to their geographic distribution (Aus and SoAm). Estimation of the divergence time suggests that SoAm and Aus lineages diverged near the last glacial maximum event during the Pleistocene (64-8 KYA). Interestingly, we found that the French reference strain is closely related to the Australian strains, with a recent divergence (405-51 YA), likely associated to human movements. Additionally, we identified different lineages within the South American population, revealing that Patagonia contains a similar genetic diversity comparable to that of other lineages in S. cerevisiae. These findings support the idea of a Pleistocene-dated divergence between South Hemisphere lineages, where the Nothofagus and Araucaria ecological niches likely favoured the extensive distribution of L. cidri in Patagonia.

Identification of new ethanol-tolerant yeast strains with fermentation potential from central Patagonia.

The quest for new wild yeasts has increasingly gained attention because of their potential ability to provide unique organoleptic characters to fermented beverages. In this sense, Patagonia offers a wide diversity of ethanol-tolerant yeasts and stands out as a bioprospecting alternative. This study characterized the genetic and phenotypic diversity of yeast isolates obtained from Central Chilean Patagonia and analyzed their fermentation potential under different fermentative conditions. We recovered 125 colonies from Nothofagus spp. bark samples belonging to five yeast species: Saccharomyces eubayanus, Saccharomyces uvarum, Lachancea cidri, Kregervanrija delftensis, and Hanseniaspora valbyensis. High-throughput microcultivation assays demonstrated the extensive phenotypic diversity among Patagonian isolates, where Saccharomyces spp and L. cidri isolates exhibited the most outstanding fitness scores across the conditions tested. Fermentation performance assays under wine, mead, and beer conditions demonstrated the specific potential of the different species for each particular beverage. Saccharomyces spp. were the only isolates able to ferment beer wort. Interestingly, we found that L. cidri is a novel candidate species to ferment wine and mead, exceeding the fermentation capacity of a commercial strain. Unlike commercial strains, we found that L. cidri does not require nutritional supplements for efficient mead fermentation. In addition, L. cidri produces succinic and acetic acids, providing a distinct profile to the final fermented product. This work demonstrates the importance of bioprospecting efforts in Patagonia to isolate novel wild yeast strains with extraordinary biotechnological potential for the fermentation industry.

Yeasts from temperate forests.

Yeasts are ubiquitous in temperate forests. While this broad habitat is well-defined, the yeasts inhabiting it and their life cycles, niches, and contributions to ecosystem functioning are less understood. Yeasts are present on nearly all sampled substrates in temperate forests worldwide. They associate with soils, macroorganisms, and other habitats and no doubt contribute to broader ecosystem-wide processes. Researchers have gathered information leading to hypotheses about yeasts' niches and their life cycles based on physiological observations in the laboratory as well as genomic analyses, but the challenge remains to test these hypotheses in the forests themselves. Here, we summarize the habitat and global patterns of yeast diversity, give some information on a handful of well-studied temperate forest yeast genera, discuss the various strategies to isolate forest yeasts, and explain temperate forest yeasts' contributions to biotechnology. We close with a summary of the many future directions and outstanding questions facing researchers in temperate forest yeast ecology. Yeasts present an exciting opportunity to better understand the hidden world of microbial ecology in this threatened and global habitat.

Helium structures around SF5+ and SF6+: novel intermolecular potential and mass spectrometry experiments.

Helium clusters around the recently experimentally observed sulphur hexafluoride SF6+ and sulphur pentafluoride SF5+ ions are investigated using a combined experimental and theoretical effort. Mass spectrometry ion yields are obtained and the energetics and structure of the corresponding HeN-SF6+ and HeN-SF5+ clusters are analyzed using path integral molecular dynamics calculations as a function of N, the number of He atoms, employing a new intermolecular potential describing the interaction between the dopant and the surrounding helium. The new force field is optimized on benchmark potential energy ab initio calculations and represented by improved Lennard-Jonnes analytical expressions. This procedure improves the previous potentials employed in similar simulations for neutral SF6 attached to helium nanodroplets. The theoretical analysis explains the characteristic features observed in the experimental ion yields which suggest the existence of stable configurations at specific sizes.

A Path Integral Molecular Dynamics Simulation of a Harpoon-Type Redox Reaction in a Helium Nanodroplet.

We present path integral molecular dynamics (PIMD) calculations of an electron transfer from a heliophobic Cs2 dimer in its (3Σu) state, located on the surface of a He droplet, to a heliophilic, fully immersed C60 molecule. Supported by electron ionization mass spectroscopy measurements (Renzler et al., J. Chem. Phys.2016, 145, 181101), this spatially quenched reaction was characterized as a harpoon-type or long-range electron transfer in a previous high-level ab initio study (de Lara-Castells et al., J. Phys. Chem. Lett.2017, 8, 4284). To go beyond the static approach, classical and quantum PIMD simulations are performed at 2 K, slightly below the critical temperature for helium superfluidity (2.172 K). Calculations are executed in the NVT ensemble as well as the NVE ensemble to provide insights into real-time dynamics. A droplet size of 2090 atoms is assumed to study the impact of spatial hindrance on reactivity. By changing the number of beads in the PIMD simulations, the impact of quantization can be studied in greater detail and without an implicit assumption of superfluidity. We find that the reaction probability increases with higher levels of quantization. Our findings confirm earlier, static predictions of a rotational motion of the Cs2 dimer upon reacting with the fullerene, involving a substantial displacement of helium. However, it also raises the new question of whether the interacting species are driven out-of-equilibrium after impurity uptake, since reactivity is strongly quenched if a full thermal equilibration is assumed. More generally, our work points towards a novel mechanism for long-range electron transfer through an interplay between nuclear quantum delocalization within the confining medium and delocalized electronic dispersion forces acting on the two reactants.

Theoretical Study of Cationic Alkali Dimers Interacting with He: Li2+-He and Na2+-He van der Waals Complexes.

We present a theoretical study on the potential energy surface and bound states of He-A2+ complexes, where A is one of the alkali Li or Na atoms. The intermolecular interactions were systematically investigated by high-level ab initio electronic structure computations, and the corresponding raw data were then employed to reproduce accurate analytical expressions of the potential surfaces. In turn, we used these potentials to evaluate bound configurations of the trimers from nuclear quantum calculations and to extract information on the effect of orientational anisotropy of the forces and the interplay between repulsive and attractive interaction within the potential surfaces. The spatial features of the bound states are analyzed and discussed in detail. We found that both systems are going under large amplitude stretching and bending motions with high zero-point energies. Despite the large differences in the potential well-depths, the correct treatment of nuclear quantum effects provides insights on the effect of different strength of the ionic interaction on the spectral structure of such cationic alkali van der Waals complexes, related to the mobility of ions and the formation of cold-molecules in He-controlled environments.

A combined experimental and theoretical investigation of Cs+ ions solvated in HeN clusters.

Solvation of Cs+ ions inside helium droplets has been investigated both experimentally and theoretically. On the one hand, mass spectra of doped helium clusters ionized with a crossed electron beam, HeNCs+, have been recorded for sizes up to N = 60. The analysis of the ratio between the observed peaks for each size N reveals evidences of the closure of the first solvation shell when 17 He atoms surround the alkali ion. On the other hand, we have obtained energies and geometrical structures of the title clusters by means of basin-hopping, diffusion Monte Carlo (DMC), and path integral Monte Carlo (PIMC) methods. The analytical He-Cs+ interaction potential employed in our calculations is represented by the improved Lennard-Jones expression optimized on high level ab initio energies. The weakness of the existing interaction between helium and Cs+ in comparison with some other alkali ions such as Li+ is found to play a crucial role. Our theoretical findings confirm that the first solvation layer is completed at N = 17 and both evaporation and second difference energies obtained with the PIMC calculation seem to reproduce a feature observed at N = 12 for the experimental ion abundance. The analysis of the DMC probability distributions reveals the important contribution from the icosahedral structure to the overall configuration for He12Cs+.

Antarctic yeasts: analysis of their freeze-thaw tolerance and production of antifreeze proteins, fatty acids and ergosterol.

BACKGROUND: Microorganisms have evolved a number of mechanisms to thrive in cold environments, including the production of antifreeze proteins, high levels of polyunsaturated fatty acids, and ergosterol. In this work, several yeast species isolated from Antarctica were analyzed with respect to their freeze-thaw tolerance and production of the three abovementioned compounds, which may also have economic importance. RESULTS: The freeze-thaw tolerance of yeasts was widely variable among species, and a clear correlation with the production of any of the abovementioned compounds was not observed. Antifreeze proteins that were partially purified from Goffeauzyma gastrica maintained their antifreeze activities after several freeze-thaw cycles. A relatively high volumetric production of ergosterol was observed in the yeasts Vishniacozyma victoriae, G. gastrica and Leucosporidium creatinivorum, i.e., 19, 19 and 16 mg l- 1, respectively. In addition, a high percentage of linoleic acid with respect to total fatty acids was observed in V. victoriae (10%), Wickerhamomyces anomalus (12%) and G. gastrica (13%), and a high percentage of alpha linoleic acid was observed in L. creatinivorum (3.3%). CONCLUSIONS: Given these results, the abovementioned yeasts are good candidates to be evaluated for use in the production of antifreeze proteins, fatty acids, and ergosterol at the industrial scale.

Cold tolerance mechanisms of two arthropods from the Andean Range of Central Chile: Agathemera crassa (Insecta: Agathemeridae) and Euathlus condorito (Arachnida: Theraphosidae).

Two strategies have been described for cold tolerance in arthropods: (1) freeze-tolerant organisms, which can survive the formation of ice crystals and (2) freeze-avoidant organisms, which prevent the ice crystal formation by super cooling their internal fluids. We studied two arthropods from the Andean Range in central Chile (2400 m a.s.l.), the stick insect Agathemera crassa commonly named as "Chinchemolle", and the tarantula spider Euathlus condorito commonly named as "Araña pollito", in order to evaluate how they respond to low temperatures at the physiological and molecular levels. We sampled the soil temperature during one year to track the temperature changes that these organisms must overcome. We found minimum temperatures around -6 °C in autumn, while the temperature were stable at 0 °C in winter due to the snow. The average field-cooling rate was 0.01 ±â€¯0.006 °C min-1. For both arthropods we determined the super cooling point (SCP) at a cooling rate of 1 °C min-1 and its subsequent survival, finding that A. crassa is a freezing tolerant organism with a SCP of -3.8 ±â€¯1.8 °C and 100% survival, while E. condorito is a freezing avoidant organism with a SCP of -3.0 ±â€¯1.3 °C and 0% survival. The SCP and survival were not affected by the season in which individuals were collected, the SCP was significantly affected by the cooling rate of the experiment. Both species had low molecular weight cryoprotective in their hemolymph that could explain their cold-tolerance behavior. Glucose, glycerol, and trehalose were found in A. crassa's hemolymph, only glucose and glycerol were found in E. condorito's. We analyzed the hemolymph proteins and found no seasonal differences in composition for either species and also we detected protein antifreeze activity in the hemolymph from both arthropods.

Adsorption of molecular hydrogen on coronene with a new potential energy surface.

Benchmark interaction energies between coronene, C24H12, and molecular hydrogen, H2, have been computed by means of high level electronic structure calculations. Binding energies, equilibrium distances and strengths of the long range attraction, evaluated for the basic configurations of the H2-C24H12 complex, indicate that the system is not too affected by the relative orientations of the diatom, suggesting that its behavior can be approximated to that of a pseudoatom. The obtained energy profiles have confirmed the noncovalent nature of the bonding and serve to tune-up the parameters of a new force field based on the atom-bond approach which correctly describes the main features of the H2-coronene interaction. The structure and binding energies of (para-H2)N-coronene clusters have been investigated with an additive model for the above mentioned interactions and exploiting basin-hopping and path integral Monte Carlo calculations for N = 1-16 at T = 2 K. Differences with respect to the prototypical (rare gas)N-coronene aggregates have been discussed.

Comparative investigation of pure and mixed rare gas atoms on coronene molecules.

Clusters formed by the combination of rare gas (RG) atoms of He, Ne, Ar, and Kr on coronene have been investigated by means of a basin-hopping algorithm and path integral Monte Carlo calculations at T = 2 K. Energies and geometries have been obtained and the role played by the specific RG-RG and RG-coronene interactions on the final results is analysed in detail. Signatures of diffuse behavior of the He atoms on the surface of the coronene are in contrast with the localization of the heavier species, Ar and Kr. The observed coexistence of various geometries for Ne suggests the motion of the RG atoms on the multi-well potential energy surface landscape offered by the coronene. Therefore, the investigation of different clusters enables a comparative analysis of localized versus non-localized features. Mixed Ar-He-coronene clusters have also been considered and the competition of the RG atoms to occupy the docking sites on the molecule is discussed. All the obtained information is crucial to assess the behavior of coronene, a prototypical polycyclic aromatic hydrocarbon clustering with RG atoms at a temperature close to that of interstellar medium, which arises from the critical balance of the interactions involved.

Screening and characterization of amylase and cellulase activities in psychrotolerant yeasts.

BACKGROUND: Amylases and cellulases have great potential for application in industries such as food, detergent, laundry, textile, baking and biofuels. A common requirement in these fields is to reduce the temperatures of the processes, leading to a continuous search for microorganisms that secrete cold-active amylases and cellulases. Psychrotolerant yeasts are good candidates because they inhabit cold-environments. In this work, we analyzed the ability of yeasts isolated from the Antarctic region to grow on starch or carboxymethylcellulose, and their potential extracellular amylases and cellulases. RESULT: All tested yeasts were able to grow with soluble starch or carboxymethylcellulose as the sole carbon source; however, not all of them produced ethanol by fermentation of these carbon sources. For the majority of the yeast species, the extracellular amylase or cellulase activity was higher when cultured in medium supplemented with glucose rather than with soluble starch or carboxymethylcellulose. Additionally, higher amylase activities were observed when tested at pH 5.4 and 6.2, and at 30-37 °C, except for Rhodotorula glacialis that showed elevated activity at 10-22 °C. In general, cellulase activity was high until pH 6.2 and between 22-37 °C, while the sample from Mrakia blollopis showed high activity at 4-22 °C. Peptide mass fingerprinting analysis of a potential amylase from Tetracladium sp. of about 70 kDa, showed several peptides with positive matches with glucoamylases from other fungi. CONCLUSIONS: Almost all yeast species showed extracellular amylase or cellulase activity, and an inducing effect by the respective substrate was observed in a minor number of yeasts. These enzymatic activities were higher at 30 °C in most yeast, with highest amylase and cellulase activity in Tetracladium sp. and M. gelida, respectively. However, Rh. glacialis and M. blollopis displayed high amylase or cellulase activity, respectively, under 22 °C. In this sense, these yeasts are interesting candidates for industrial processes that require lower temperatures.

Examination of the Feynman-Hibbs Approach in the Study of NeN-Coronene Clusters at Low Temperatures.

Feynman-Hibbs (FH) effective potentials constitute an appealing approach for investigations of many-body systems at thermal equilibrium since they allow us to easily include quantum corrections within standard classical simulations. In this work we apply the FH formulation to the study of NeN-coronene clusters (N = 1-4, 14) in the 2-14 K temperature range. Quadratic (FH2) and quartic (FH4) contributions to the effective potentials are built upon Ne-Ne and Ne-coronene analytical potentials. In particular, a new corrected expression for the FH4 effective potential is reported. FH2 and FH4 cluster energies and structures-obtained from energy optimization through a basin-hopping algorithm as well as classical Monte Carlo simulations-are reported and compared with reference path integral Monte Carlo calculations. For temperatures T > 4 K, both FH2 and FH4 potentials are able to correct the purely classical calculations in a consistent way. However, the FH approach fails at lower temperatures, especially the quartic correction. It is thus crucial to assess the range of applicability of this formulation and, in particular, to apply the FH4 potentials with great caution. A simple model of N isotropic harmonic oscillators allows us to propose a means of estimating the cutoff temperature for the validity of the method, which is found to increase with the number of atoms adsorbed on the coronene molecule.

Tolerance to Ultraviolet Radiation of Psychrotolerant Yeasts and Analysis of Their Carotenoid, Mycosporine, and Ergosterol Content.

Yeasts colonizing the Antarctic region are exposed to a high ultraviolet radiation evolving mechanisms to minimize the UV radiation damages, such as the production of UV-absorbing or antioxidant compounds like carotenoid pigments and mycosporines. Ergosterol has also been suggested to play a role in this response. These compounds are also economically attractive for several industries such as pharmaceutical and food, leading to a continuous search for biological sources of them. In this work, the UV-C radiation tolerance of yeast species isolated from the sub-Antarctic region and their production of carotenoids, mycosporines, and ergosterol were evaluated. Dioszegia sp., Leuconeurospora sp. (T27Cd2), Rhodotorula laryngis, Rhodotorula mucilaginosa, and Cryptococcus gastricus showed the highest UV-C radiation tolerance. The yeasts with the highest content of carotenoids were Dioszegia sp. (OHK torulene), Rh. laryngis (torulene and lycopene), Rh. mucilaginosa, (torulene, gamma carotene, and lycopene), and Cr. gastricus (2-gamma carotene). Probable mycosporine molecules and biosynthesis intermediates were found in Rh. laryngis, Dioszegia sp., Mrakia sp., Le. creatinivora, and Leuconeurospora sp. (T27Cd2). Ergosterol was the only sterol detected in all yeasts, and M. robertii and Le. creatinivora showed amounts higher than 4 mg g−1. Although there was not a well-defined relation between UV-C tolerance and the production of these three kinds of compounds, the majority of the yeasts with lower amounts of carotenoids showed lower UV-C tolerance. Dioszegia sp., M. robertii, and Le. creatinivora were the greatest producers of carotenoids, ergosterol, and mycosporines, respectively, representing good candidates for future studies intended to increase their production for large-scale applications.

Quantum Features of Anionic Species He*⁻ and He2*⁻ in Small He(N) Clusters.

We present variational calculations on systems containing a few boson helium atoms attached to electronically excited atomic and molecular helium anions He*⁻ and He2*⁻ and characterize their structures and energetics. Previously reported high-level ab initio results [Huber, S. E.; Mauracher, A. Mol. Phys. 2014, 112, 794] to describe the interactions between excited (metastable) anions and a neutral He atom have been employed. For the case of the atomic species He*⁻, the corresponding interaction with He suggests large anharmonicity effects due to the presence of a deep well of ∼17,500 cm⁻¹ at short distances, together with a more external shallow secondary well of ∼4 cm⁻¹, both supporting bound levels. Moreover, when a sum of pairwise interactions is assumed to describe the full PES corresponding to the presence of several neutral He atoms, geometrical constraints already predict the complete solvation of the anionic impurity by six helium atoms, giving rise to a bipyramidal structure. In turn, for the anisotropic weak interaction He-He2*⁻, where the anionic dimer is considered as a rigid rotor, the obtained structures show the tendency of the helium atoms to pack themselves together and largely far away from the dopant, thereby confirming the heliophobic character of He2*⁻.

A configurational study of helium clusters doped with He(∗-) and He2(∗-).

Helium clusters doped with electronically excited atomic and molecular helium anions He(∗-) and He2(∗-) at T = 0.4 K are studied by means of path integral Monte Carlo calculations. Geometry and energetics of the systems with up to 32 solvating He atoms are characterised. The interactions between the anions and the neutral He atoms have been described by fitting previously reported ab initio points to analytical expressions. The HeN-He(∗-) clusters with N > 6 display a structure defined by a bipyramid which completely solvates the atomic anion, whereas the rest of surrounding He atoms form a dimple around that initial cage. On the contrary, the structures observed for the HeN-He2(∗-) clusters clearly show the dopant located outside the helium droplet, thereby confirming the heliophobic character of He2(∗-).

Reactive scattering calculations for (87)Rb+(87)RbHe→Rb2((3)Σ(u)(+),v)+He from ultralow to intermediate energies.

We investigate atom-diatom reactive collisions, as a preliminary step,in order to assess the possibility of forming Rb(2) molecules in their lowest triplet electronic state by cold collisions of rubidium atoms on the surface of helium nanodroplets [corrected]. A simple model related to the well-known Rosen treatment of linear triatomic molecules [N. Rosen, J. Chem. Phys. 1, 319 (1933)] in relative coordinates is used, allowing to estimate reactive probabilities for different values of the total angular momentum. The best available full dimensional potential energy surface [Guillon et al., J. Chem. Phys. 136, 174307 (2012)] is employed through the calculations. Noticeable values of the probabilities in the ultracold regime, which numerically fulfill the Wigner threshold law, support the feasibility of the process. The rubidium dimer is mainly produced at high vibrational states, and the reactivity is more efficient for a bosonic helium partner than when the fermion species is considered.

Coronene molecules in helium clusters: Quantum and classical studies of energies and configurations.

Coronene-doped helium clusters have been studied by means of classical and quantum mechanical (QM) methods using a recently developed He-C24H12 global potential based on the use of optimized atom-bond improved Lennard-Jones functions. Equilibrium energies and geometries at global and local minima for systems with up to 69 He atoms were calculated by means of an evolutive algorithm and a basin-hopping approach and compared with results from path integral Monte Carlo (PIMC) calculations at 2 K. A detailed analysis performed for the smallest sizes shows that the precise localization of the He atoms forming the first solvation layer over the molecular substrate is affected by differences between relative potential minima. The comparison of the PIMC results with the predictions from the classical approaches and with diffusion Monte Carlo results allows to examine the importance of both the QM and thermal effects.