Multiphysical graph neural network (MP-GNN) for COVID-19 drug design.

Graph neural networks (GNNs) are the most promising deep learning models that can revolutionize non-Euclidean data analysis. However, their full potential is severely curtailed by poorly represented molecular graphs and features. Here, we propose a multiphysical graph neural network (MP-GNN) model based on the developed multiphysical molecular graph representation and featurization. All kinds of molecular interactions, between different atom types and at different scales, are systematically represented by a series of scale-specific and element-specific graphs with distance-related node features. From these graphs, graph convolution network (GCN) models are constructed with specially designed weight-sharing architectures. Base learners are constructed from GCN models from different elements at different scales, and further consolidated together using both one-scale and multi-scale ensemble learning schemes. Our MP-GNN has two distinct properties. First, our MP-GNN incorporates multiscale interactions using more than one molecular graph. Atomic interactions from various different scales are not modeled by one specific graph (as in traditional GNNs), instead they are represented by a series of graphs at different scales. Second, it is free from the complicated feature generation process as in conventional GNN methods. In our MP-GNN, various atom interactions are embedded into element-specific graph representations with only distance-related node features. A unique GNN architecture is designed to incorporate all the information into a consolidated model. Our MP-GNN has been extensively validated on the widely used benchmark test datasets from PDBbind, including PDBbind-v2007, PDBbind-v2013 and PDBbind-v2016. Our model can outperform all existing models as far as we know. Further, our MP-GNN is used in coronavirus disease 2019 drug design. Based on a dataset with 185 complexes of inhibitors for severe acute respiratory syndrome coronavirus (SARS-CoV/SARS-CoV-2), we evaluate their binding affinities using our MP-GNN. It has been found that our MP-GNN is of high accuracy. This demonstrates the great potential of our MP-GNN for the screening of potential drugs for SARS-CoV-2. Availability: The Multiphysical graph neural network (MP-GNN) model can be found in https://github.com/Alibaba-DAMO-DrugAI/MGNN. Additional data or code will be available upon reasonable request.

Molecular persistent spectral image (Mol-PSI) representation for machine learning models in drug design.

Artificial intelligence (AI)-based drug design has great promise to fundamentally change the landscape of the pharmaceutical industry. Even though there are great progress from handcrafted feature-based machine learning models, 3D convolutional neural networks (CNNs) and graph neural networks, effective and efficient representations that characterize the structural, physical, chemical and biological properties of molecular structures and interactions remain to be a great challenge. Here, we propose an equal-sized molecular 2D image representation, known as the molecular persistent spectral image (Mol-PSI), and combine it with CNN model for AI-based drug design. Mol-PSI provides a unique one-to-one image representation for molecular structures and interactions. In general, deep models are empowered to achieve better performance with systematically organized representations in image format. A well-designed parallel CNN architecture for adapting Mol-PSIs is developed for protein-ligand binding affinity prediction. Our results, for the three most commonly used databases, including PDBbind-v2007, PDBbind-v2013 and PDBbind-v2016, are better than all traditional machine learning models, as far as we know. Our Mol-PSI model provides a powerful molecular representation that can be widely used in AI-based drug design and molecular data analysis.

Ultraviolet photodetector based on RbCu2I3microwire.

Copper-based halide perovskites have shown great potential in lighting and photodetection due to their excellent photoelectric properties, good stability and lead-free nature. However, as an important piece of copper-based perovskites, the synthesis and application of RbCu2I3have never been reported. Here, we demonstrate the synthesis of high-quality RbCu2I3microwires (MWs) by a fast-cooling hot saturated solution method. The prepared MWs exhibit an orthorhombic structure with a smooth surface. Optical measurements show the RbCu2I3MWs have a sharp ultraviolet absorption edge with 3.63 eV optical band gap and ultra-large stokes shift (300 nm) in photoluminescence. The subsequent photodetector based on a single RbCu2I3MW shows excellent ultraviolet detection performance. Under the 340 nm illumination, the device shows a specific detectivity of 5.0 × 109Jones and a responsivity of 380 mA·W-1. The synthesis method and physical properties of RbCu2I3could be a guide to the future optoelectronic application of the new material.

Association of Polymorphisms of Mismatch Repair Genes hMLHI and hMSH2 with Breast Cancer Susceptibility: A Meta-Analysis.

This article serves to evaluate the association of polymorphisms of mismatch repair genes (hMLH1 and hMSH2) with breast cancer (BC) susceptibility through a meta-analysis. Our methods involved extensive research in Chinese and English databases that examined the association of hMLH1 and hMSH2 polymorphisms with susceptibility to BC, strictly abiding by established inclusion and exclusion criteria. Software Stata 12.0 was used for statistical data analysis. A total of 12 studies were available for meta-analysis, published between 2014 and 2017, of which respectively 9 studies explored the association of hMLH1 (rs1799977 A > G and rs63750447 T > A) and 3 studies explored the association of hMSH2 (rs4987188 [Gly322Asp] and rs17217772 [Asn127Ser]) with patients' susceptibility to BC. The results showed that both the rs1799977 A > G polymorphism GA + GG genotype (especially in the Caucasian population) and the rs63750447 T > A polymorphism TA + AA genotype in the hMLH1 gene increased patients' susceptibility to BC. The genotype detection method was selected as a target for subgroup analysis. According to studies where MassARRAY assay was conducted, the rs1799977 A > G polymorphism was correlated with BC susceptibility in the dominant model, while rs4987188 (Gly322Asp) and rs17217772 (Asn127Ser) of the hMSH2 gene presented no observable correlation with the risk for BC. Both the rs1799977 A > G and rs63750447 T > A polymorphisms in the hMLH1 gene showed a significant association with a markedly increased risk for BC, while rs4987188 (Gly322Asp) and rs17217772 (Asn127Ser) of the hMSH2 gene were not clearly correlated with BC susceptibility.

Thermal tolerance of dried shoots of the moss Bryum argenteum.

We conducted laboratory experiments to determine the lethal temperatures of the shoots of dried Bryum argenteum and to determine how this restoration species responds to extreme environments. We specifically assessed changes in gene expression levels in the shoots of dried B. argenteum plants that were subjected to sudden heat shock (control (20 ± 2°C), 80°C, 100°C, 110°C or 120°C) followed by exposure to heat for an additional 10, 20, 30 or 60 min. After they were exposed to heat, the samples were placed in wet sand medium, and their survival and regeneration abilities were evaluated daily for 56 days. The results showed that lethal temperatures significantly reduced the shoot regeneration potential, delayed both shoot and protonemal emergence times and reduced the protonemal emergence area. In addition, the expression of nine genes (HSF3, HSP70, ERF, LEA, ELIP, LHCA, LHCB, Tr288 and DHN) was induced by temperature stress, as assessed after 30 min of exposure. Additionally, a new thermal tolerance level for dried B. argenteum - 120°C for 20 min - was determined, which was the highest temperature recorded for this moss; this tolerance exceeded the previous record of 110°C for 10 min. These findings help elucidate the survival mechanism of this species under heat shock stress and facilitate the recovery and restoration of destroyed ecosystems.

Metal-organic framework MIL-101 as sorbent based on double-pumps controlled on-line solid-phase extraction coupled with high-performance liquid chromatography for the determination of flavonoids in environmental water samples.

A novel method with high sensitivity for the rapid determination of chrysin, apigenin and luteolin in environment water samples was developed by double-pumps controlled on-line solid-phase extraction (SPE) coupled with high-performance liquid chromatography (HPLC). In the developed technique, metal organic framework MIL-101 was synthesized and applied as a sorbent for SPE. The as-synthesized MIL-101 was characterized by scanning electron microscope, X-ray diffraction spectrometry, thermal gravimetric analysis and micropore physisorption analysis. The MIL-101 behaved as a fast kinetics in the adsorption of chrysin, apigenin and luteolin. On-line SPE of chrysin, apigenin and luteolin was processed by loading a sample solution at a flow rate of 1.0 mL/min for 10 min. The extracted analytes were subsequently eluted into a ZORBAX Bonus-RP analytical column (25 cm long × 4.6 mm i.d.) for HPLC separation under isocratic condition with a mobile phase (MeOH: ACN: 0.02 M H3 PO4 = 35:35:30) at a flow rate of 1.0 mL/min. Experimental conditions, including ionic strength, sample pH, sample loading rates, sample loading time and desorption analytes time, were further optimized to obtain efficient preconcentration and high-precision determination of the analytes mentioned above. The method achieved the merits of simplicity, rapidity, sensitivity, wide linear range and high sample throughput. The possible mechanism for the adsorption of flavonoids on MIL-101 was proposed. The developed method has been applied to determine trace chrysin, apigenin and luteolin in a variety of environmental water samples.

Accompanying mild hypothermia significantly improved the prognosis of septic mice than artificial mild hypothermia.

BACKGROUND: Some patients with sepsis are found with accompanying mild hypothermia (ACMH); however, the effects of ACMH on the patients with sepsis are poorly understood. OBJECTIVE: To compare the impacts of ACMH and artificial mild hypothermia (ATMH) on mortality, systemic inflammatory reactions, and organ functions in mice with sepsis. METHODS: Septic mouse models were induced and divided into ACMH, un-hypothermia, keep normothermia, and ATMH groups, according to the anal temperature and the thermic intervention strategy. The mortality rate, serum levels of tumor necrosis factor α (TNF-α), interferon γ (IFN-γ), and interleukin (IL)-4 and liver and renal functions of the mice in each group were recorded. Liver, lung, and renal tissues of the mice were stained and examined under optic microscope. RESULTS: The mortality rate in the ACMH group was the lowest among all the sepsis groups. Increased serum levels of TNF-α, IFN-γ, and IL-4 and impairments of the liver and renal functions were found in the septic mice. The serum levels of TNF-α, IFN-γ, and IL-4 were significantly lower and the liver and renal functions of ACMH group were not impaired significantly as compared with other sepsis groups. Pathological examinations of the lung, liver, and renal tissues showed that the ACMH group were with the lowest pathological score among all the sepsis groups. CONCLUSION: Accompanying mild hypothermia and ATMH could both reduce mortalities in mice with sepsis, and ACMH could reduce mortality even lower, and more alleviate systemic inflammatory responses and the damages in lung, kidney, and other organs were lighter.

Activating PV-positive neurons in ventral thalamic reticular nucleus reduces pain sensitivity in mice.

Previous studies have demonstrated that thalamic reticular nucleus (TRN) and the sub-nuclei play important roles in pain sensation. Our previous findings showed that activating parvalbumin-positive (PV+) neurons in dorsal sector of TRN (dTRN) could reduce the pain threshold and consequently increase the pain sensitivity of mice. Recent studies have shown that activation of GABAergic projection of TRN to ventrobasal thalamus (VB) alleviated pathological pain. GABAergic neurons in TRN are mainly PV+ neurons. However, the exact roles of ventral TRN (vTRN) PV+ neurons in pain sensation remain unclear. In this study, the designer receptors exclusively activated by designer drugs (DREADD) method was used to activate the PV+ neurons in vTRN of PV-Cre transgenic mice, and the mechanical threshold and thermal latency were measured to investigate the regulatory effects of vTRN on pain sensitivity in mice. Thereafter, PV-Cre transgenic mice, conditional anterograde axonal tract tracing, and immunohistochemistry were used to investigate the distribution of PV+ neurons fibers in vTRN. The results showed that the activation of PV+ neurons in vTRN increased the mechanical threshold and thermal latency, which indicated reduction of pain sensitivity. The fibers of these neurons mainly projected to ventral posterolateral thalamic nucleus (VPL), ventral posteromedial thalamic nucleus (VPM), ventrolateral thalamic nucleus (VL), centrolateral thalamic nucleus (CL) and various other brain regions. These findings indicated that activation of PV+ neurons in the vTRN decreased pain sensitivity in mice, which provided additional evidence on the mechanisms of PV+ neurons of TRN in regulating neuralgia.

The distribution of HCN2-positive cells in the gastrointestinal tract of mice.

HCN2 channels are involved in the spontaneous rhythmic activities of some CNS neurons and act by generating I(f) current. The gastrointestinal (GI) tract is known to be capable of spontaneous rhythmic activity; however, the possible role of HCN2 channels in this organ has not yet been elucidated. This study investigated the distribution of HCN2-positive cells in the mouse GI tract using immunohistochemistry. To identify the nature of these HCN2 cells, anti-ChAT and anti-Kit antibodies were used to co-label neurons and the interstitial cells of Cajal (ICCs), respectively. Additionally, differences in the distribution of HCN2-positive cells within the GI tract were also analyzed. Our results showed that HCN2 channels were mainly located within the myenteric neurons of the enteric nervous system in the GI tract. Double-staining revealed that HCN2-positive neurons were labeled by ChAT, indicating that these HCN2-positive cells are also cholinergic neurons. Although the HCN2-positive cells were not stained by the anti-Kit antibody, their processes were in close proximity to ICCs around the myenteric plexus region. Moreover, several differences in the distribution of HCN2 in the stomach, small intestine and colon were partly consistent with the regional differences in the spontaneous rhythmic activities of these organs. Basing on the role HCN2, we suggested that HCN2 channels facilitate the release of Ach from cholinergic neurons to affect the GI peristalsis by acting on M receptors on the ICCs. However, the HCN2 channels are not directly involved in spontaneous slow-wave initiation by ICCs.

Reference gene selection in the desert plant Eremosparton songoricum.

Eremosparton songoricum (Litv.) Vass. (E. songoricum) is a rare and extremely drought-tolerant desert plant that holds promise as a model organism for the identification of genes associated with water deficit stress. Here, we cloned and evaluated the expression of eight candidate reference genes using quantitative real-time reverse transcriptase polymerase chain reactions. The expression of these candidate reference genes was analyzed in a diverse set of 20 samples including various E. songoricum plant tissues exposed to multiple environmental stresses. GeNorm analysis indicated that expression stability varied between the reference genes in the different experimental conditions, but the two most stable reference genes were sufficient for normalization in most conditions. EsEF and Esα-TUB were sufficient for various stress conditions, EsEF and EsACT were suitable for samples of differing germination stages, and EsGAPDHand EsUBQ were most stable across multiple adult tissue samples. The Es18S gene was unsuitable as a reference gene in our analysis. In addition, the expression level of the drought-stress related transcription factor EsDREB2 verified the utility of E. songoricum reference genes and indicated that no single gene was adequate for normalization on its own. This is the first systematic report on the selection of reference genes in E. songoricum, and these data will facilitate future work on gene expression in this species.

Steady and dynamic photosynthetic responses of seedlings from contrasting successional groups under low-light growth conditions.

To test the hypothesis that leaf-level photosynthetic-related traits might confer late successionals a competitive advantage over early successionals in low-light growth conditions, steady photosynthetic assimilation- and dynamic photosynthetic induction-related traits were examined in low-light-grown seedlings with contrasting successional status. Compared with the early successionals, late successionals as a group significantly exhibited lower leaf gas exchange rates. While late successionals required a longer time to respond to simulated sunflecks, they had lower rates of induction losses after sunflecks. Such photosynthetic induction traits allowed late successionals to more effectively utilize subsequent sunflecks. It was observed that plants with lower gas exchange rates responded more slowly to simulated sunfelcks, but they had lower rates of induction losses after sunflecks. In addition, the rate of response to sunflecks was positively correlated with the rate of induction loss after sunflecks across the successional status of species. A principal components analysis (PCA) demonstrated that early and late successionals were separated along the first axis of the PCA, and that early successionals were grouped on the right and were associated with higher gas exchange rates, fast responses to sunflecks, and rapid rates of induction loss after sunflecks; late successionals held an opposite pattern. Overall, our results suggest that smaller respiratory carbon losses and lower metabolic costs give late successionals a competitive advantage in low-light growth conditions, that late successionals have an advantage over early successionals in utilizing sunflecks, and thus that the successional status of species are mainly associated with the leaf-level photosynthetic-related traits.

Gas exchange and hydraulics in seedlings of Hevea brasiliensis during water stress and recovery.

The response of plants to drought has received significant attention, but far less attention has been given to the dynamic response of plants during recovery from drought. Photosynthetic performance and hydraulic capacity were monitored in seedlings of Hevea brasiliensis under water stress and during recovery following rewatering. Leaf water relation, gas exchange rate and hydraulic conductivity decreased gradually after water stress fell below a threshold, whereas instantaneous water use efficiency and osmolytes increased significantly. After 5 days of rewatering, leaf water relation, maximum stomatal conductance (g(s-max)) and plant hydraulic conductivity had recovered to the control levels except for sapwood area-specific hydraulic conductivity, photosynthetic assimilation rate and osmolytes. During the phase of water stress, stomata were almost completely closed before water transport efficiency decreased substantially, and moreover, the leaf hydraulic pathway was more vulnerable to water stress-induced embolism than the stem hydraulic pathway. Meanwhile, g(s-max) was linearly correlated with hydraulic capacity when water stress exceeded a threshold. In addition, a positive relationship was shown to occur between the recovery of g(s-max) and of hydraulic capacity during the phase of rewatering. Our results suggest (i) that stomatal closure effectively reduces the risk of xylem dysfunction in water-stressed plants at the cost of gas exchange, (ii) that the leaf functions as a safety valve to protect the hydraulic pathway from water stress-induced dysfunction to a larger extent than does the stem and (iii) that the full drought recovery of gas exchange is restricted by not only hydraulic factors but also non-hydraulic factors.

Regeneration pattern of primary forest species across forest-field gradients in the subtropical mountains of Southwestern China.

Evergreen broad-leaved forest is now gradually degraded and fragmented, and there is an increase in the amount of habitat edges as a result of long-term human activity. However, the role of edges in the regeneration of primary forest species is poorly understood. After 20 years of the edge creation, we analyzed primary forest species distribution and abundance, and changes in floristic composition, vegetation structure across forest-field gradients in Ailao Mountain, SW China. Our results revealed that there was a higher abundance and richness of primary species, late secondary species and thorny lianas at the distances 0-50 m than at the distances more than 50 m from the edge into the forest exterior. At the distances >50 m, no individuals of dominant canopy trees Lithocarpus xylocarpus, Castanopsis wattii, and L. jingdongensis were found, whereas the abundance of early pioneer shrub species and herbaceous cover was significantly greater. The richness of primary species showed a decrease with increasing distances from the forest edge to the exterior, particularly of medium-seeded primary species showing a drastic decrease. Moreover, no large-seeded primary species occurred at the distances >60 m. This study indicates that the forest edge as a buffer zone may be in favor of primary species regeneration. A dense shrub and herb layer, and seed dispersal may be the major factors limiting the forest regeneration farther from the forest edge. Therefore, to facilitate forest recovery processes, management should give priority to the protection of buffer zones of this forest edge.

Independence of stem and leaf hydraulic traits in six Euphorbiaceae tree species with contrasting leaf phenology.

Hydraulic traits and hydraulic-related structural properties were examined in three deciduous (Hevea brasiliensis, Macaranga denticulate, and Bischofia javanica) and three evergreen (Drypetes indica, Aleurites moluccana, and Codiaeum variegatum) Euphorbiaceae tree species from a seasonally tropical forest in south-western China. Xylem water potential at 50% loss of stem hydraulic conductivity (P50(stem)) was more negative in the evergreen tree, but leaf water potential at 50% loss of leaf hydraulic conductivity (P50(leaf)) did not function as P50(stem) did. Furthermore, P50(stem) was more negative than P50(leaf) in the evergreen tree; contrarily, this pattern was not observed in the deciduous tree. Leaf hydraulic conductivity overlapped considerably, but stem hydraulic conductivity diverged between the evergreen and deciduous tree. Correspondingly, structural properties of leaves overlapped substantially; however, structural properties of stem diverged markedly. Consequently, leaf and stem hydraulic traits were closely correlated with leaf and stem structural properties, respectively. Additionally, stem hydraulic efficiency was significantly correlated with stem hydraulic resistance to embolism; nevertheless, such a hydraulic pattern was not found in leaf hydraulics. Thus, these results suggest: (1) that the evergreen and deciduous tree mainly diverge in stem hydraulics, but not in leaf hydraulics, (2) that regardless of leaf or stem, their hydraulic traits result primarily from structural properties, and not from leaf phenology, (3) that leaves are more vulnerable to drought-induced embolism than stem in the evergreen tree, but not always in the deciduous tree and (4) that there exists a trade-off between hydraulic efficiency and safety for stem hydraulics, but not for leaf hydraulics.

Modulation of out-of-plane reflected waves by using acoustic metasurfaces with tapered corrugated holes.

In this paper, modulation of reflected wavefront out of the incident plane by a tunable acoustic metasurface is investigated based on the fully generalized Snell's law in the three-dimensional space. The metasurface is constructed by a square lattice of circular holes with gradient annular bumps. The phase shift is tuned by changing the volume of water filled in the holes. The acoustic wave steering out of the incident plane and the out-of-plane acoustic focusing with the oblique incidence at the subwavelength scale are demonstrated numerically by selecting suitable distributions of water depth. The numerical results show that the wavefront of the reflected wave can be manipulated over a wide frequency range; and the gradient design of the unit cells can suppress the parasitic reflection. The present work is relevant to the practical design of novel acoustic devices.

Discovery of tetrahydrocarbazoles with potent hypoglycemic and hypolipemic activities.

A series of tetrahydrocarbazole derivatives was designed and synthesized on the basis of the AMP-activated protein kinase activator GY3. All the synthesized compounds were screened in HepG2 cell lines for glucose consumption activity and several of them showed potent glucose decreasing activity. In vivo evaluation of the hypoglycemic and hypolipemic effects indicated that 7a exhibited comparable activity with pioglitazone, but with a weaker body-weight increasing effect. The pharmacokinetic profiles of 7a were also investigated.

Author Correction: Identifying long-term stable refugia for relict plant species in East Asia.

In the original version of this Article, an incorrect sample size was provided for the number of relict species (443 instead of 442) and the number of relict forests (423 instead of 422). These errors have been corrected in both the PDF and HTML versions of the Article.

Identifying long-term stable refugia for relict plant species in East Asia.

Today East Asia harbors many "relict" plant species whose ranges were much larger during the Paleogene-Neogene and earlier. The ecological and climatic conditions suitable for these relict species have not been identified. Here, we map the abundance and distribution patterns of relict species, showing high abundance in the humid subtropical/warm-temperate forest regions. We further use Ecological Niche Modeling to show that these patterns align with maps of climate refugia, and we predict species' chances of persistence given the future climatic changes expected for East Asia. By 2070, potentially suitable areas with high richness of relict species will decrease, although the areas as a whole will probably expand. We identify areas in southwestern China and northern Vietnam as long-term climatically stable refugia likely to preserve ancient lineages, highlighting areas that could be prioritized for conservation of such species.

Potential effects of climate change on geographic distribution of the Tertiary relict tree species Davidia involucrata in China.

This study, using species distribution modeling (involving a new approach that allows for uncertainty), predicts the distribution of climatically suitable areas prevailing during the mid-Holocene, the Last Glacial Maximum (LGM), and at present, and estimates the potential formation of new habitats in 2070 of the endangered and rare Tertiary relict tree Davidia involucrata Baill. The results regarding the mid-Holocene and the LGM demonstrate that south-central and southwestern China have been long-term stable refugia, and that the current distribution is limited to the prehistoric refugia. Given future distribution under six possible climate scenarios, only some parts of the current range of D. involucrata in the mid-high mountains of south-central and southwestern China would be maintained, while some shift west into higher mountains would occur. Our results show that the predicted suitable area offering high probability (0.5‒1) accounts for an average of only 29.2% among the models predicted for the future (2070), making D. involucrata highly vulnerable. We assess and propose priority protected areas in light of climate change. The information provided will also be relevant in planning conservation of other paleoendemic species having ecological traits and distribution ranges comparable to those of D. involucrata.

Fabrication of MPEG-b-PMAA capped YVO4:Eu nanoparticles with biocompatibility for cell imaging.

A novel nanoparticle with multilayer core-shell architecture for cell imaging is designed and synthesized by coating a fluorescent YVO4:Eu core with a diblock copolymer, MPEG-b-PMAA. The synthesis of YVO4:Eu core, which further makes MPEG-b-PMAA-YVO4:Eu NPs adapt for cell imaging, is guided by the model determined upon the evaluation of pH and CEu%. The PMAA block attached tightly on the YVO4:Eu core forms the inner shell and the MPEG block forms the biocompatible outermost shell. Factors including reaction time, reaction temperature, CEu% and pH are optimized for the preparation of the YVO4:Eu NPs. A precise defined model is established according to analyzing the coefficients of pH and CEu% during the synthesis. The MPEG-b-PMAA-YVO4:Eu NPs, with an average diameter of 24 nm, have a tetragonal structure and demonstrate luminescence in the red region, which lies in a biological window (optical imaging). Significant enhancement in luminescence intensity by MPEG-b-PMAA-YVO4:Eu NPs formation is observed. The capping copolymer MPEG-b-PMAA improves the dispersibility of hydrophobic YVO4:Eu NPs in water, making the NPs stable under different conditions. In addition, the biocompatibility MPEG layer reduces the cytotoxicity of the nanoparticles effectively. 95% cell viability can be achieved at the NPs concentration of 800 mgL(-1) after 24h of culture. Cellular uptake of the MPEG-b-PMAA-YVO4:Eu NPs is evaluated by cell imaging assay, indicating that the NPs can be taken up rapidly and largely by cancerous or non-cancerous cells through an endocytosis mechanism.