Constant charge method or constant potential method:Which is better for molecular modeling of electrical double layers?

In molecular modeling of electrical double layers(EDLs),the constant charge method(CCM)is prized for its computational efficiency but cannot maintain electrode equipotentiality like the more resourceintensive constant potential method(CPM),potentially leading to inaccuracies.In certain scenarios,CCM can yield results identical to CPM.However,there are no clear guidelines to determine when CCM is sufficient and when CPM is required.Here,we conduct a series of molecular simulations across various electrodes and electrolytes to present a comprehensive comparison between CCM and CPM under different charging modes.Results reveal that CCM approximates CPM effectively in capturing equilibrium EDL and current-driven dynamics in open electrode systems featuring ionic liquids or regular concentration aqueous electrolytes,while CPM is indispensable in scenarios involving organic and highly concentrated aqueous electrolytes,nanoconfinement effects,and voltage-driven dynamics.This work helps to select appropriate methods for modeling EDL systems,prioritizing accuracy while considering computationalefficiency.

Predicting liposome formulations by the integrated machine learning and molecular modeling approaches

Liposome is one of the most widely used carriers for drug delivery because of the great biocompatibility and biodegradability.Due to the complex formulation components and preparation process,formulation screening mostly relies on trial-and-error process with low efficiency.Here liposome formulation prediction models have been built by machine learning(ML)approaches.The important parameters of liposomes,including size,polydispersity index(PDI),zeta potential and encapsulation,are predicted individually by optimal ML algorithm,while the formulation features are also ranked to provide important guidance for formulation design.The analysis of key parameter reveals that drug molecules with logS[-3,-6],molecular complexity[500,1000]and XLogP3(≥2)are priority for preparing liposome with higher encapsulation.In addition,naproxen(NAP)and palmatine HCl(PAL)represented the insoluble and water-soluble molecules are prepared as liposome formulations to validate prediction ability.The consistency between predicted and experimental value verifies the satisfied accuracy of ML models.As the drug properties are critical for liposome particles,the molecular interactions and dynamics of NAP and PAL liposome are further investigated by coarse-grained molecular dynamics simulations.The modeling structure reveals that NAP molecules could distribute into lipid layer,while most PAL molecules aggregate in the inner aqueous phase of liposome.The completely different physical state of NAP and PAL confirms the importance of drug properties for liposome formulations.In summary,the general prediction models are built to predict liposome formulations,and the impacts of key factors are analyzed by combing ML with molecular modeling.The availability and rationality of these intelligent prediction systems have been proved in this study,which could be applied for liposome formulation development in the future.

A SARS-CoV-2 neutralizing antibody discovery by single cell sequencing and molecular modeling

Although vaccines have been developed,mutations of SARS-CoV-2,especially the dominant B.1.617.2(delta)and B.1.529(omicron)strains with more than 30 mutations on their spike protein,have caused a significant decline in prophylaxis,calling for the need for drug improvement.Antibodies are drugs preferentially used in infectious diseases and are easy to get from immunized organisms.The current study combined molecular modeling and single memory B cell sequencing to assess candidate sequences before experiments,providing a strategy for the fabrication of SARS-CoV-2 neutralizing antibodies.A total of 128 sequences were obtained after sequencing 196 memory B cells,and 42 sequences were left after merging extremely similar ones and discarding incomplete ones,followed by homology modeling of the antibody variable region.Thirteen candidate sequences were expressed,of which three were tested positive for receptor binding domain recognition but only one was confirmed as having broad neutralization against several SARS-CoV-2 variants.The current study successfully obtained a SARS-CoV-2 antibody with broad neutralizing abilities and provided a strategy for antibody development in emerging infectious diseases using single memory B cell BCR sequencing and computer assistance in antibody fabrication.

Molecular Mechanism and Molecular Design of Lubricating Oil Antioxidants

To overcome the limitations of traditional experimental“trial and error”methods in lubricant additive design,a new molecular design method based on molecular structure parameters is established here.The molecular mechanism of the antioxidant reaction of hindered phenol,diphenylamine,and alkyl sulfide are studied via molecular simulations.Calculation results show that the strong electron-donating ability and high hydrogen-donating activity of the antioxidant molecule and the low hydrogen-abstracting activity of free radicals formed after dehydrogenation are the internal molecular causes of the shielding of phenol and diphenylamine from scavenging peroxy free radicals,and the strong electron-donating ability is the internal molecular cause of the high activity of thioether in decomposing alkyl hydrogen peroxide.Based on this antioxidant molecular mechanism,a molecular design rule of antioxidant is proposed,namely“high EHOMO,large Q(S),low bond dissociation energy BDE(O—H)and BDE(N—H)”.Two new antioxidants,PAS-I and PAS-II,are designed and prepared by chemical bonding of hindered phenol,diphenylamine,and sulfur atoms.Experimental results show that these antioxidants both have excellent antioxidant effects in lubricating oil,and that PAS-II is the superior antioxidant,consistent with theoretical predictions.

Novel umami peptides from two Termitomyces mushrooms and molecular docking to the taste receptor T1R1/T1R3

Wild edible Termitomyces mushrooms are popular in Southwest China and umami is important flavor qualities of edible mushrooms.This study aimed to understand the umami taste of Termitomyces intermedius and Termitomyces aff.bulborhizus.Ten umami peptides from aqueous extracts were separated using a Sephadex G-15 gel filtration chromatography.The intense umami fraction was evaluated by both sensory evaluation and electronic tongue.They were identified as KLNDAQAPK,DSTDEKFLR,VGKGAHLSGEH,MLKKKKLA,SLGFGGPPGY,TVATFSSSTKPDD,AMDDDEADLLLLAM,VEDEDEKPKEK,SPEEKKEEET and PEGADKPNK.Seven peptides,except VEDEDEKPKEK,SPEEKKEEET and PEGADKPNK were selectively synthesized to verify their taste characteristics.All these 10 peptides had umami or salt taste.The 10 peptides were conducted by molecular docking to study their interaction with identified peptides and the umami taste receptor T1R1/T1R3.All these 10 peptides perfectly docked the active residues in the T1R3 subunit.Our results provide theoretical basis for the umami taste and address the umami mechanism of two wild edible Termitomyces mushrooms.

Analysis of CH_(4) and H_(2) Adsorption on Heterogeneous Shale Surfaces Using aMolecular Dynamics Approach

Determining the adsorption of shale gas on complex surfaces remains a challenge in molecular simulation studies.Difficulties essentially stem from the need to create a realistic shale structure model in terms of mineral heterogeneityand multiplicity.Moreover,precise characterization of the competitive adsorption of hydrogen andmethane in shale generally requires the experimental determination of the related adsorptive capacity.In thisstudy,the adsorption of adsorbates,methane(CH_(4)),and hydrogen(H_(2))on heterogeneous shale surface modelsof Kaolinite,Orthoclase,Muscovite,Mica,C_(60),and Butane has been simulated in the frame of a moleculardynamic’s numerical technique.The results show that these behaviors are influenced by pressure and potentialenergy.On increasing the pressure from 500 to 2000 psi,the sorption effect for CH_(4)significantly increasesbut shows a decline at a certain stage(if compared to H_(2)).The research findings also indicate that raw shalehas a higher capacity to adsorb CH_(4)compared to hydrogen.However,in shale,this difference is negligible.

Molecular mechanisms underlying microglial sensing and phagocytosis in synaptic pruning

Microglia are the main non-neuronal cells in the central nervous system that have important roles in brain development and functional connectivity of neural circuits.In brain physiology,highly dynamic microglial processes are facilitated to sense the surrounding environment and stimuli.Once the brain switches its functional states,microglia are recruited to specific sites to exert their immune functions,including the release of cytokines and phagocytosis of cellular debris.The crosstalk of microglia between neurons,neural stem cells,endothelial cells,oligodendrocytes,and astrocytes contributes to their functions in synapse pruning,neurogenesis,vascularization,myelination,and blood-brain barrier permeability.In this review,we highlight the neuron-derived“find-me,”“eat-me,”and“don't eat-me”molecular signals that drive microglia in response to changes in neuronal activity for synapse refinement during brain development.This review reveals the molecular mechanism of neuron-microglia interaction in synaptic pruning and presents novel ideas for the synaptic pruning of microglia in disease,thereby providing important clues for discovery of target drugs and development of nervous system disease treatment methods targeting synaptic dysfunction.

Risk stratification for radioactive iodine refractoriness using molecular alterations in distant metastatic differentiated thyroid cancer

Objective: Patients with radioactive iodine-refractory differentiated thyroid cancer(RAIR-DTC) are often diagnosed with delay and constrained to limited treatment options. The correlation between RAI refractoriness and the underlying genetic characteristics has not been extensively studied.Methods: Adult patients with distant metastatic DTC were enrolled and assigned to undergo next-generation sequencing of a customized 26-gene panel(Thyro Lead). Patients were classified into RAIR-DTC or non-RAIR groups to determine the differences in clinicopathological and molecular characteristics. Molecular risk stratification(MRS) was constructed based on the association between molecular alterations identified and RAI refractoriness, and the results were classified as high, intermediate or low MRS.Results: A total of 220 patients with distant metastases were included, 63.2% of whom were identified as RAIRDTC. Genetic alterations were identified in 90% of all the patients, with BRAF(59.7% vs. 17.3%), TERT promoter(43.9% vs. 7.4%), and TP53 mutations(11.5% vs. 3.7%) being more prevalent in the RAIR-DTC group than in the non-RAIR group, except for RET fusions(15.8% vs. 39.5%), which had the opposite pattern. BRAF and TERT promoter are independent predictors of RAIR-DTC, accounting for 67.6% of patients with RAIR-DTC. MRS was strongly associated with RAI refractoriness(P<0.001), with an odds ratio(OR) of high to low MRS of 7.52 [95%confidence interval(95% CI), 3.96-14.28;P<0.001] and an OR of intermediate to low MRS of 3.20(95% CI,1.01-10.14;P=0.041).Conclusions: Molecular alterations were associated with RAI refractoriness, with BRAF and TERT promoter mutations being the predominant contributors, followed by TP53 and DICER1 mutations. MRS might serve as a valuable tool for both prognosticating clinical outcomes and directing precision-based therapeutic interventions.

Development of in situ characterization techniques in molecular beam epitaxy

Ex situ characterization techniques in molecular beam epitaxy(MBE)have inherent limitations,such as being prone to sample contamination and unstable surfaces during sample transfer from the MBE chamber.In recent years,the need for improved accuracy and reliability in measurement has driven the increasing adoption of in situ characterization techniques.These techniques,such as reflection high-energy electron diffraction,scanning tunneling microscopy,and X-ray photoelectron spectroscopy,allow direct observation of film growth processes in real time without exposing the sample to air,hence offering insights into the growth mechanisms of epitaxial films with controlled properties.By combining multiple in situ characterization techniques with MBE,researchers can better understand film growth processes,realizing novel materials with customized properties and extensive applications.This review aims to overview the benefits and achievements of in situ characterization techniques in MBE and their applications for material science research.In addition,through further analysis of these techniques regarding their challenges and potential solutions,particularly highlighting the assistance of machine learning to correlate in situ characterization with other material information,we hope to provide a guideline for future efforts in the development of novel monitoring and control schemes for MBE growth processes with improved material properties.

Absorption characteristics,model,and molecular mechanism of hydrogen sulfide in morpholine acetate aqueous solution

The solubility of H_(2)S was measured in solutions of N-butyl-N-methylmorpholine acetate([Bmmorp][Ac])containing 20%-40%(mass)water at experimental temperatures ranged from 298.15 to 328.15 K and pressures up to 320 k Pa.The total solubility of H_(2)S increased with higher temperatures,lower pressures,and reduced water content.The reaction equilibrium thermodynamic model was used to correlate the solubility data.The results indicate that the chemical reaction equilibrium constant decrease with increasing water content and temperature,whereas Henry constant increase with increasing water content and temperature.Compared with other ionic liquids,H_(2)S exhibits a higher physical absorption enthalpy and a lower chemical absorption enthalpy in[Bmmorp][Ac]aqueous solution.This suggests that[Bmmorp][Ac]has a strong physical affinity for H_(2)S and low energy requirement for desorption.Quantum chemical methods were used to investigate the molecular mechanism of H_(2)S absorption in ionic liquids.The interaction energy analysis revealed that the binding of H_(2)S with the ionic liquid in a1:2 ratio is more stable.Detailed analyses by the methods of the interaction region indicator and the atoms in molecules were conducted to the interactions between H_(2)S and the ionic liquid.

Lithium Ion Transport Environment by Molecular Vibrations in Ion-Conducting Glasses

Controlling Li ion transport in glasses at atomic and molecular levels is key to realizing all-solid-state batteries,a promising technology for electric vehicles.In this context,Li_(3)PS_(4)glass,a promising solid electrolyte candidate,exhibits dynamic coupling between the Li^(+)cation mobility and the PS_(4)^(3-)anion libration,which is commonly referred to as the paddlewheel effect.In addition,it exhibits a concerted cation diffusion effect(i.e.,a cation-cation interaction),which is regarded as the essence of high Li ion transport.However,the correlation between the Li^(+)ions within the glass structure can only be vaguely determined,due to the limited experimental information that can be obtained.Here,this study reports that the Li ions present in glasses can be classified by evaluating their valence oscillations via Bader analysis to topologically analyze the chemical bonds.It is found that three types of Li ions are present in Li_(3)PS_(4)glass,and that the more mobile Li ions(i.e.,the Li3-type ions)exhibit a characteristic correlation at relatively long distances of 4.0-5.0A.Furthermore,reverse Monte Carlo simulations combined with deep learning potentials that reproduce X-ray,neutron,and electron diffraction pair distribution functions showed an increase in the number of Li3-type ions for partially crystallized glass structures with improved Li ion transport properties.Our results show order within the disorder of the Li ion distribution in the glass by a topological analysis of their valences.Thus,considering the molecular vibrations in the glass during the evaluation of the Li ion valences is expected to lead to the development of new solid electrolytes.

Comparative study on phase transition behaviors of fractional molecular field theory and random-site Ising model

Fractional molecular field theory(FMFT)is a phenomenological theory that describes phase transitions in crystals with randomly distributed components,such as the relaxor-ferroelectrics and spin glasses.In order to verify the feasibility of this theory,this paper fits it to the Monte Carlo simulations of specific heat and susceptibility versus temperature of two-dimensional(2D)random-site Ising model(2D-RSIM).The results indicate that the FMFT deviates from the 2D-RSIM significantly.The main reason for the deviation is that the 2D-RSIM is a typical system of component random distribution,where the real order parameter is spatially heterogeneous and has no symmetry of space translation,but the basic assumption of FMFT means that the parameter is spatially uniform and has symmetry of space translation.

Predictive modeling for postoperative delirium in elderly patients with abdominal malignancies using synthetic minority oversampling technique

BACKGROUND Postoperative delirium,particularly prevalent in elderly patients after abdominal cancer surgery,presents significant challenges in clinical management.AIM To develop a synthetic minority oversampling technique(SMOTE)-based model for predicting postoperative delirium in elderly abdominal cancer patients.METHODS In this retrospective cohort study,we analyzed data from 611 elderly patients who underwent abdominal malignant tumor surgery at our hospital between September 2020 and October 2022.The incidence of postoperative delirium was recorded for 7 d post-surgery.Patients were divided into delirium and non-delirium groups based on the occurrence of postoperative delirium or not.A multivariate logistic regression model was used to identify risk factors and develop a predictive model for postoperative delirium.The SMOTE technique was applied to enhance the model by oversampling the delirium cases.The model’s predictive accuracy was then validated.RESULTS In our study involving 611 elderly patients with abdominal malignant tumors,multivariate logistic regression analysis identified significant risk factors for postoperative delirium.These included the Charlson comorbidity index,American Society of Anesthesiologists classification,history of cerebrovascular disease,surgical duration,perioperative blood transfusion,and postoperative pain score.The incidence rate of postoperative delirium in our study was 22.91%.The original predictive model(P1)exhibited an area under the receiver operating characteristic curve of 0.862.In comparison,the SMOTE-based logistic early warning model(P2),which utilized the SMOTE oversampling algorithm,showed a slightly lower but comparable area under the curve of 0.856,suggesting no significant difference in performance between the two predictive approaches.CONCLUSION This study confirms that the SMOTE-enhanced predictive model for postoperative delirium in elderly abdominal tumor patients shows performance equivalent to that of traditional methods,effectively addressing data imbalance.

Volume-averaged modeling of multiphase solidification with equiaxed crystal sedimentation in a steel ingot

Macrosegregation is a critical factor that limits the mechanical properties of materials.The impact of equiaxed crystal sedimentation on macrosegregation has been extensively studied,as it plays a significant role in determining the distribution of alloying elements and impurities within a material.To improve macrosegregation in steel connecting shafts,a multiphase solidification model that couples melt flow,heat transfer,microstructure evolution,and solute transport was established based on the volume-averaged Eulerian-Eulerian approach.In this model,the effects of liquid phase,equiaxed crystals,columnar dendrites,and columnar-to-equiaxed transition(CET)during solidification and evolution of microstructure can be considered simultaneously.The sedimentation of equiaxed crystals contributes to negative macrosegregation,where regions between columnar dendrites and equiaxed crystals undergo significant A-type positive macrosegregation due to the CET.Additionally,noticeable positive macrosegregation occurs in the area of final solidification in the ingot.The improvement in macrosegregation is beneficial for enhancing the mechanical properties of connecting shafts.To mitigate the thermal convection of molten steel resulting from excessive superheating,reducing the superheating during casting without employing external fields or altering the design of the ingot mold is indeed an effective approach to control macrosegregation.

Structure-based multi-ligand molecular modeling to predict the synergistic effects of limonin and obacunone from Simiao pill against nitric oxide synthase 3 associated with hyperuricemia

Hyperuricemia(HUA)mainly occurs because of purine metabolism disorders.We recently proposed that limonin from Simiao pill may have therapeutic effects on nitric oxide synthase 3(NOS3)that is related to HUA.Concurrently,our previous work employed a new method,structure-based multi-ligand molecular modeling,to identify potential agents from a herbal formula that may produce synergistic effects and may have the potential to develop combination drugs.Therefore,we employed multi-ligand modeling to seek compounds with potential synergistic effects with limonin against NOS3.We obtained 403 multi-ligand docking results between 403 compounds and the limonin-NOS3 complex(average affinity–8.297 kcal/mol).Then we selected the top 10 highest binding affinity compounds for virtual pharmacokinetic and toxicity screening and we found that only obacunone passed all filters.We further subjected obacunone,bound to limonin and NOS3,to molecular dynamics simulations.We found that the NOS3-limonin-obacunone complex was more stable than the NOS3-limonin complex,based on the root mean square deviation of backbone Cαatoms and root mean square fluctuation,which suggests that synergistic effects may exist between limonin and obacunone.Further cell and animal experimental research is required to verify our results.

Background removal from global auroral images:Data-driven dayglow modeling

Global images of auroras obtained by cameras on spacecraft are a key tool for studying the near-Earth environment.However,the cameras are sensitive not only to auroral emissions produced by precipitating particles,but also to dayglow emissions produced by photoelectrons induced by sunlight.Nightglow emissions and scattered sunlight can contribute to the background signal.To fully utilize such images in space science,background contamination must be removed to isolate the auroral signal.Here we outline a data-driven approach to modeling the background intensity in multiple images by formulating linear inverse problems based on B-splines and spherical harmonics.The approach is robust,flexible,and iteratively deselects outliers,such as auroral emissions.The final model is smooth across the terminator and accounts for slow temporal variations and large-scale asymmetries in the dayglow.We demonstrate the model by using the three far ultraviolet cameras on the Imager for Magnetopause-to-Aurora Global Exploration(IMAGE)mission.The method can be applied to historical missions and is relevant for upcoming missions,such as the Solar wind Magnetosphere Ionosphere Link Explorer(SMILE)mission.

Surrogate modeling for unsaturated infiltration via the physics and equality-constrained artificial neural networks

Machine learning(ML)provides a new surrogate method for investigating groundwater flow dynamics in unsaturated soils.Traditional pure data-driven methods(e.g.deep neural network,DNN)can provide rapid predictions,but they do require sufficient on-site data for accurate training,and lack interpretability to the physical processes within the data.In this paper,we provide a physics and equalityconstrained artificial neural network(PECANN),to derive unsaturated infiltration solutions with a small amount of initial and boundary data.PECANN takes the physics-informed neural network(PINN)as a foundation,encodes the unsaturated infiltration physical laws(i.e.Richards equation,RE)into the loss function,and uses the augmented Lagrangian method to constrain the learning process of the solutions of RE by adding stronger penalty for the initial and boundary conditions.Four unsaturated infiltration cases are designed to test the training performance of PECANN,i.e.one-dimensional(1D)steady-state unsaturated infiltration,1D transient-state infiltration,two-dimensional(2D)transient-state infiltration,and 1D coupled unsaturated infiltration and deformation.The predicted results of PECANN are compared with the finite difference solutions or analytical solutions.The results indicate that PECANN can accurately capture the variations of pressure head during the unsaturated infiltration,and present higher precision and robustness than DNN and PINN.It is also revealed that PECANN can achieve the same accuracy as the finite difference method with fewer initial and boundary training data.Additionally,we investigate the effect of the hyperparameters of PECANN on solving RE problem.PECANN provides an effective tool for simulating unsaturated infiltration.

Identification of prognostic molecular subtypes and model based on CD8+ T cells for lung adenocarcinoma

Background:Cytotoxic T lymphocytes(CD8+T)cells function critically in mediating anti-tumor immune response in cancer patients.Characterizing the specific functions of CD8+T cells in lung adenocarcinoma(LUAD)could help better understand local anti-tumor immune responses and estimate the effect of immunotherapy.Methods:Gens related to CD8+T cells were identified by cluster analysis based on the single-cell sequencing data of three LUAD tissues and their paired normal tissues.Weighted gene co-expression network analysis(WGCNA),consensus clustering,differential expression analysis,least absolute shrinkage and selection operator(LASSO)and Cox regression analysis were conducted to classify molecular subtypes for LUAD and to develop a risk model using prognostic genes related to CD8+T cells.Expression of the genes in the prognostic model,their effects on tumor cell invasion,and interactions with CD8+T cells were verified by cell experiments.Results:This study defined two LUAD clusters(CD8+0 and CD8+1)based on CD8+T cells,with cluster CD8+0 being significantly associated with the prognosis of LUAD.Three heterogeneous subtypes(clusters 1,2,and 3)differing in prognosis,genome mutation events,and immune status were categorized using 42 prognostic genes.A prognostic model created based on 11 significant genes(including CD200R1,CLEC17A,ZC3H12D,GNG7,SNX30,CDCP1,NEIL3,IGF2BP1,RHOV,ABCC2,and KRT81)was able to independently estimate the death risk for patients in different LUAD cohorts.Moreover,the model also showed general applicability in external validation cohorts.Low-risk patients could benefit more from taking immunotherapy and were significantly related to the resistance to anticancer drugs.The results from cell experiments demonstrated that the expression of CD200R1,CLEC17A,ZC3H12D,GNG7,and SNX30 was significantly downregulated,while that of CDCP1,NEIL3,IGF2BP1,RHOV,ABCC2 and KRT81 was upregulated in LUAD cells.Inhibition of CD200R1 greatly increased the invasiveness of the LUAD cells,but inhibiting CDCP1 expression weakened the invasion ability of LUAD cells.Conclusion:This study defined two prognostic CD8+T cell clusters and classified three heterogeneous molecular subtypes for LUAD.A prognostic model predictive of the potential effects of immunotherapy on LUAD patients was developed.

A methodology for damage evaluation of underground tunnels subjected to static loading using numerical modeling

We have proposed a methodology to assess the robustness of underground tunnels against potential failure.This involves developing vulnerability functions for various qualities of rock mass and static loading intensities.To account for these variations,we utilized a Monte Carlo Simulation(MCS)technique coupled with the finite difference code FLAC^(3D),to conduct two thousand seven hundred numerical simulations of a horseshoe tunnel located within a rock mass with different geological strength index system(GSIs)and subjected to different states of static loading.To quantify the severity of damage within the rock mass,we selected one stress-based(brittle shear ratio(BSR))and one strain-based failure criterion(plastic damage index(PDI)).Based on these criteria,we then developed fragility curves.Additionally,we used mathematical approximation techniques to produce vulnerability functions that relate the probabilities of various damage states to loading intensities for different quality classes of blocky rock mass.The results indicated that the fragility curves we obtained could accurately depict the evolution of the inner and outer shell damage around the tunnel.Therefore,we have provided engineers with a tool that can predict levels of damages associated with different failure mechanisms based on variations in rock mass quality and in situ stress state.Our method is a numerically developed,multi-variate approach that can aid engineers in making informed decisions about the robustness of underground tunnels.

Artificial intelligence-driven radiomics study in cancer:the role of feature engineering and modeling

Modern medicine is reliant on various medical imaging technologies for non-invasively observing patients’anatomy.However,the interpretation of medical images can be highly subjective and dependent on the expertise of clinicians.Moreover,some potentially useful quantitative information in medical images,especially that which is not visible to the naked eye,is often ignored during clinical practice.In contrast,radiomics performs high-throughput feature extraction from medical images,which enables quantitative analysis of medical images and prediction of various clinical endpoints.Studies have reported that radiomics exhibits promising performance in diagnosis and predicting treatment responses and prognosis,demonstrating its potential to be a non-invasive auxiliary tool for personalized medicine.However,radiomics remains in a developmental phase as numerous technical challenges have yet to be solved,especially in feature engineering and statistical modeling.In this review,we introduce the current utility of radiomics by summarizing research on its application in the diagnosis,prognosis,and prediction of treatment responses in patients with cancer.We focus on machine learning approaches,for feature extraction and selection during feature engineering and for imbalanced datasets and multi-modality fusion during statistical modeling.Furthermore,we introduce the stability,reproducibility,and interpretability of features,and the generalizability and interpretability of models.Finally,we offer possible solutions to current challenges in radiomics research.