Molecular insights and optimization strategies for the competitive binding of engineered ACE2 proteins: a multiple replica molecular dynamics study.

The COVID-19 pandemic caused by severe acute respiratory syndrome coronavirus-2 (SARS-CoV-2) continues to spread globally, and rapid viral evolution and the emergence of new variants pose challenges to pandemic control. During infection, the spike protein of SARS-CoV-2 interacts with the human ACE2 protein via its receptor binding domain (RBD), and it is known that engineered forms of ACE2 can compete with wild-type (WT) ACE2 for binding to inhibit infection. Here, we conducted multiple replica molecular dynamics (MRMD) simulations to study the mechanisms of the engineered ACE2 variants 3N39 and 3N94 and provide directions for optimization. Our findings reveal that engineered ACE2 is notably more efficacious in systems that show weaker binding to WT ACE2 (i.e., WT and BA.1 RBD), but also faces immune escape as the virus evolves. Moreover, by modifying residue types near the binding interface, engineered ACE2 alters the electrostatic potential distribution and reconfigures the hydrogen bonding network, which results in modified binding to the RBD. However, this structural rearrangement does not occur in all RBD variants. In addition, we identified potentially engineerable beneficial residues and potentially engineerable detrimental residues in both ACE2 and RBD. Functional conservation can thus enable the optimization of these residues and improve the binding competitiveness of engineered ACE2, which therefore provides additional immune escape prevention. Finally, we conclude that these findings have implications for understanding the mechanisms responsible for engineered ACE2 and can help us to develop engineered ACE2 proteins that show superior performance.

Deciphering the binding mechanism of inhibitors of the SARS-CoV-2 main protease through multiple replica accelerated molecular dynamics simulations and free energy landscapes.

The pneumonia outbreak caused by the SARS-CoV-2 virus poses a serious threat to human health and the world economy. The development of safe and highly effective antiviral drugs is of great significance for the treatment of COVID-19. The main protease (Mpro) of SARS-CoV-2 is a key enzyme for viral replication and transcription and has no homolog in humans. Therefore, the Mpro is an ideal target for the design of drugs against COVID-19. Insights into the inhibitor-Mpro binding mechanism and conformational changes of the Mpro are essential for the design of potent drugs that target the Mpro. In this study, we analyzed the conformational changes of the Mpro that are induced by the binding of three inhibitors, YTV, YSP and YU4, using multiple replica accelerated molecular dynamics (MR-aMD) simulations, dynamic cross-correlation map (DCCM) calculations, principal component analysis (PCA), and free energy landscape (FEL) analysis. The results from DCCM calculations and PCA show that the binding of inhibitors significantly affects the kinetic behavior of the Mpro and induces a conformational rearrangement of the Mpro. The binding ability and binding mechanism of YTV, YSP and YU4 to the Mpro were investigated using the molecular mechanics Poisson-Boltzmann surface area (MM-PBSA) method. The results indicate that substitution of the tert-butanol group by methylbenzene and trifluoromethyl groups enhances the binding ability of YSP and YU4 to the Mpro compared with YTV; moreover, massive hydrophobic interactions are detected between the inhibitors and the Mpro. Meanwhile, T25, L27, H41, M49, N142, G143, C145, M165, E166 and Q189 are identified as the key residues for inhibitor-Mpro interactions using residue-based free energy decomposition calculations, which can be employed as efficient targets in the design of drugs that inhibit the activity of the Mpro.

Binding mechanism of inhibitors to SARS-CoV-2 main protease deciphered by multiple replica molecular dynamics simulations.

The outbreak caused by SARS-CoV-2 has received extensive worldwide attention. As the main protease (Mpro) in SARS-CoV-2 has no human homologues, it is feasible to reduce the possibility of targeting the host protein by accidental drugs. Thus, Mpro has been an attractive target of efficient drug design for anti-SARS-CoV-2 treatment. In this work, multiple replica molecular dynamics (MRMD) simulations, principal component analysis (PCA), free energy landscapes (FELs), and the molecular mechanics-generalized Born surface area (MM-GBSA) method were integrated together to decipher the binding mechanism of four inhibitors masitinib, O6K, FJC and GQU to Mpro. The results indicate that the binding of four inhibitors clearly affects the structural flexibility and internal dynamics of Mpro along with dihedral angle changes of key residues. The analysis of FELs unveils that the stability in the relative orientation and geometric position of inhibitors to Mpro is favorable for inhibitor binding. Residue-based free energy decomposition reveals that the inhibitor-Mpro interaction networks involving hydrogen bonding interactions and hydrophobic interactions provide significant information for the design of potent inhibitors against Mpro. The hot spot residues including H41, M49, F140, N142, G143, C145, H163, H164, M165, E166 and Q189 identified by computational alanine scanning are considered as reliable targets of clinically available inhibitors inhibiting the activities of Mpro.

Mutation-Induced Impacts on the Switch Transformations of the GDP- and GTP-Bound K-Ras: Insights from Multiple Replica Gaussian Accelerated Molecular Dynamics and Free Energy Analysis.

Mutations yield significant effect on the structural flexibility of two switch domains, SW1 and SW2, in K-Ras, which is considered as an important target of anticancer drug design. To unveil a molecular mechanism with regard to mutation-mediated tuning on the activity of K-Ras, multiple replica Gaussian accelerated molecular dynamics (MR-GaMD) simulations followed by analysis of free energy landscapes (FELs) are performed on the GDP- and GTP-bound wild-type (WT), G12V, and D33E K-Ras. The results suggest that G12V and D33E not only evidently change the flexibility of SW1 and SW2 but also greatly affect correlated motions of SW1 and SW2 separately relative to the P-loop and SW1, which exerts a certain tuning on the activity of K-Ras. The information stemming from the analyses of FELs reveals that the conformations of SW1 and SW2 are in high disorders in the GDP- and GTP-associated WT and mutated K-Ras, possibly producing significant effect on binding of guanine nucleotide exchange factors or effectors to K-Ras. The interaction networks of GDP and GTP with K-Ras are identified and the results uncover that the instability in hydrogen-bonding interactions of SW1 with GDP and GTP is mostly responsible for conformational disorder of SW1 and SW2 as well as tunes the activity of oncogenic K-Ras.

Molecular mechanism with regard to the binding selectivity of inhibitors toward FABP5 and FABP7 explored by multiple short molecular dynamics simulations and free energy analyses.

Recently, fatty acid binding proteins 5 and 7 (FABP5 and FABP7) have been regarded as the prospective targets for clinically treating multiple diseases related to FABPs. In this work, multiple short molecular dynamics (MSMD) simulations followed by binding free energy calculations were performed to investigate the binding selectivity of three inhibitors, namely, 65X, 8KS, and 5M8 toward FABP5 and FABP7. The RMSF analysis suggests that the structural flexibility of FABP5 is stronger than that of FABP7; moreover, the calculated molecular surface area of FABP5 is also larger than that of FABP7. Meanwhile, the results from the cross-correlation analysis show that the inhibitor bindings exert different impacts on the internal dynamics of FABP5 and FABP7. Binding free energies predicted by the molecular mechanics/generalized Born surface area (MM-GBSA) method indicate that the increase in the enthalpy changes caused by the bindings of inhibitors toward FABP7 relative to FABP5 mostly drives the binding selectivity of the inhibitors toward FABP5 versus FABP7. Hierarchical clustering analysis based on the energy contributions of separate residues and calculations of residue-based free energy decompositions were carried out by using the equilibrated MSMD trajectories. The obtained results not only recognize the hot interaction spots of inhibitors with FABP5 and FABP7, but also display that several common residues, namely, (T56, T54), (L60, F58), (E75, E73), (A76, A78), (D79, D77), (R81, R79), (R107, R109), (C120, L118), and (R129, R127) belonging to (FABP5, FABP7) induce obvious binding differences in the inhibitors toward FABP5 and FABP7. Therefore, these residues play significant roles in the binding selectivities of inhibitors toward FABP5 and FABP7.

Fine-mapping of QTLs for individual and total isoflavone content in soybean (Glycine max L.) using a high-density genetic map.

KEY MESSAGE: Fifteen stable QTLs were identified using a high-density soybean genetic map across multiple environments. One major QTL, qIF5-1, contributing to total isoflavone content explained phenotypic variance 49.38, 43.27, 46.59, 45.15 and 52.50%, respectively. Soybeans (Glycine max L.) are a major source of dietary isoflavones. To identify novel quantitative trait loci (QTL) underlying isoflavone content, and to improve the accuracy of marker-assisted breeding in soybean, a valuable mapping population comprised of 196 F7:8-10 recombinant inbred lines (RILs, Huachun 2 × Wayao) was utilized to evaluate individual and total isoflavone content in plants grown in four different environments in Guangdong. A high-density genetic linkage map containing 3469 recombination bin markers based on 0.2 × restriction site-associated DNA tag sequencing (RAD-seq) technology was used to finely map QTLs for both individual and total isoflavone contents. Correlation analyses showed that total isoflavone content, and that of five individual isoflavone, was significantly correlated across the four environments. Based on the high-density genetic linkage map, a total of 15 stable quantitative trait loci (QTLs) associated with isoflavone content across multiple environments were mapped onto chromosomes 02, 05, 07, 09, 10, 11, 13, 16, 17, and 19. Further, one of them, qIF5-1, localized to chromosomes 05 (38,434,171-39,045,620 bp) contributed to almost all isoflavone components across all environments, and explained 6.37-59.95% of the phenotypic variance, especially explained 49.38, 43.27, 46.59, 45.15 and 52.50% for total isoflavone. The results obtained in the present study will pave the way for a better understanding of the genetics of isoflavone accumulation and reveals the scope available for improvement of isoflavone content through marker-assisted selection.

Probing the rate-determining region of the potential energy surface for a prototypical ion-molecule reaction.

We report a joint experimental-theoretical study of the F- + HCl → HF + Cl- reaction kinetics. The experimental measurement of the rate coefficient at several temperatures was made using the selected ion flow tube method. Theoretical rate coefficients are calculated using the quasi-classical trajectory method on a newly developed global potential energy surface, obtained by fitting a large number of high-level ab initio points with augmentation of long-range electrostatic terms. In addition to good agreement between experiment and theory, analyses suggest that the ion-molecule reaction rate is significantly affected by shorter-range interactions, in addition to the traditionally recognized ion-dipole and ion-induced dipole terms. Furthermore, the statistical nature of the reaction is assessed by comparing the measured and calculated HF product vibrational state distributions to that predicted by the phase space theory.This article is part of the theme issue 'Modern theoretical chemistry'.

Molecular Dynamics Exploration of Selectivity of Dual Inhibitors 5M7, 65X, and 65Z toward Fatty Acid Binding Proteins 4 and 5.

Designing highly selective inhibitors of fatty acid binding proteins 4 and 5 (FABP4 and FABP5) is of importance for treatment of some diseases related with inflammation, metabolism, and tumor growth. In this study, molecular dynamics (MD) simulations combined with molecular mechanics generalized Born surface area (MM-GBSA) method were performed to probe binding selectivity of three inhibitors (5M7, 65X, and 65Z) to FABP4/FABP5 with Ki values of 0.022/0.50 µM, 0.011/0.086 µM, and 0.016/0.12 µM, respectively. The results not only suggest that all inhibitors associate more tightly with FABP4 than FABP5, but also prove that the main forces driving the selective bindings of inhibitors to FABP4 and FABP5 stem from the difference in the van der Waals interactions and polar interactions of inhibitors with two proteins. Meanwhile, a residue-based free energy decomposition method was applied to reveal molecular basis that inhibitors selectively interact with individual residues of two different proteins. The calculated results show that the binding difference of inhibitors to the residues (Phe16, Phe19), (Ala33, Gly36), (Phe57, Leu60), (Ala75, Ala78), (Arg126, Arg129), and (Tyr128, Tyr131) in (FABP4, FABP5) drive the selectivity of inhibitors toward FABP4 and FABP5. This study will provide great help for further design of effective drugs to protect against a series of metabolic diseases, arteriosclerosis, and inflammation.

A machine learning-based approach to predict energy layer for each field in spot-scanning proton arc therapy for lung cancer: A feasibility study.

BACKGROUND: Determining the optimal energy layer (EL) for each field, under considering both dose constraints and delivery efficiency, is crucial to promoting the development of proton arc therapy (PAT) technology. PURPOSE: This study aimed to explore the feasibility and potential clinical benefits of utilizing machine learning (ML) technique to automatically select EL for each field in PAT plans of lung cancer. METHODS: Proton Bragg peak position (BPP) was employed to characterize EL. The ground truth BPPs for each field were determined using the modified ELO-SPAT framework. Features in geometric, water-equivalent thicknesses (WET) and beamlet were defined and extracted. By analyzing the relationship between the extracted features and ground truth, a polynomial regression model with L2-norm regularization (Ridge regression) was constructed and trained. The performance of the regression model was reported as an error between the predictions and the ground truth. Besides, the predictions were used to make PAT plans (PAT_PRED). These plans were compared with those using the ground truth BPPs (PAT_TRUTH) and the mid-WET of the target volumes (PAT_MID) in terms of relative biological effectiveness-weighted dose (RWD) distributions. One hundred ten patients with lung cancer, a total of 7920 samples, were enrolled retrospectively, with 5940 cases randomly selected as the training set and the remaining 1980 cases as the testing set. Nine patients (648 samples) were collected additionally to evaluate the regression model in terms of plan quality and robustness. RESULTS: With regard to the prediction errors, the root mean squared errors and mean absolute errors between the ML-predicted and ground truth BPPs for the testing set were 9.165 and 6.572 mm, respectively, indicating differences of approximately two to three ELs. As for plan quality, the PAT_TRUTH and PAT_PRED plans performed similarly in terms of plan robustness, target coverage and organs at risk (OARs) protection, with differences smaller than 0.5 Gy(RBE). This trend was also observed for dose conformity and uniformity. The PAT_MID plans produced the lowest robustness index and lowest doses to OARs, along with the highest heterogeneity index, indicating better protection for OARs, improved plan robustness, but compromised dose homogeneity. Additionally, for relatively small tumor sizes, the PAT_MID plan demonstrated a notably poor dose conformity index. CONCLUSIONS: Within this cohort under investigation, our study demonstrated the feasibility of using ML technique to predict ELs for each field, offering a fast (within 2 s) and memory-efficient reduced way to select ELs for PAT plan.

Nanodosimetric quantity-weighted dose optimization for carbon-ion treatment planning.

Dose verification of treatment plans is an essential step in radiotherapy workflows. In this work, we propose a novel method of treatment planning based on nanodosimetric quantity-weighted dose (NQWD), which could realize biological representation using pure physical quantities for biological-oriented carbon ion-beam treatment plans and their direct verification. The relationship between nanodosimetric quantities and relative biological effectiveness (RBE) was studied with the linear least-squares method for carbon-ion radiation fields. Next, under the framework of the matRad treatment planning platform, NQWD was optimized using the existing RBE-weighted dose (RWD) optimization algorithm. The schemes of NQWD-based treatment planning were compared with the RWD treatment plans in term of the microdosimetric kinetic model (MKM). The results showed that the nanodosimetric quantity F3 - 10 had a good correlation with the radiobiological effect reflected by the relationship between RBE and F3 - 10. Moreover, the NQWD-based treatment plans reproduced the RWD plans generally. Therefore, F3 - 10 could be adopted as a radiation quality descriptor for carbon-ion treatment planning. The novel method proposed herein not only might be helpful for rapid physical verification of biological-oriented ion-beam treatment plans with the development of experimental nanodosimetry, but also makes the direct comparison of ion-beam treatment plans in different institutions possible. Thus, our proposed method might be potentially developed to be a new strategy for carbon-ion treatment planning and improve patient safety for carbon-ion radiotherapy.

Knockout of the VPS22 component of the ESCRT-II complex in rice (Oryza sativa L.) causes chalky endosperm and early seedling lethality.

In both yeast and mammals, the major constituent of the endosomal sorting complex required for transport-II (ESCRT-II) is the VPS22/EAP30 protein, which plays an important role in ubiquitin-mediated degradation of membrane proteins through the multivesicular body pathway. However, the functions of ESCRT-II subunits in plants are largely unknown. In this work, we report the genetic analysis and phenotypic characterization of mutants in OsVPS22 gene, which encodes a functional VPS22 homolog in rice. On the basis of a collection of T-DNA lines, we identified a T-DNA insertion mutant, which showed abnormal segregation ratios; we then found that the T-DNA insertion is located within the sixth intron of the OsVPS22 gene. Compared with the wild type, this vps22 mutant exhibited seedling lethality and severe reduction in shoot and root growth. In addition, the vps22 mutant had a chalky endosperm in the grain. In summary, our data suggest that OsVPS22 may be required for seedling viability and grain filling in rice, thus providing a valuable resource for further exploration of the functions of the ESCRTing machinery in plants.

Time-dependent quantum wave packet study of the Ar+H2+→ArH(+)+H reaction on a new ab initio potential energy surface for the ground electronic state (1(2)A').

A new global potential energy surface for the ground electronic state (1(2)A') of the Ar+H2(+)→ArH(+)+H reaction has been constructed by multi-reference configuration interaction method with Davidson correction and a basis set of aug-cc-pVQZ. Using 6080 ab initio single-point energies of all the regions for the dynamics, a many-body expansion function form has been used to fit these points. The quantum reactive scattering dynamics calculations taking into account the Coriolis coupling (CC) were carried out on the new potential energy surface over a range of collision energies (0.03-1.0 eV). The reaction probabilities and integral cross sections for the title reaction were calculated. The significance of including the CC quantum scattering calculation has been revealed by the comparison between the CC and the centrifugal sudden approximation calculation. The calculated cross section is in agreement with the experimental result at collision energy 1.0 eV.

Selection of breathing phase number in 4D scanned proton treatment planning optimization for lung tumors.

The standard four-dimensional (4D) treatment planning includes all breathing states in the optimization process, which is time-consuming. This work was aimed to optimize the number of intermediate phases needed for 4D proton treatment planning optimization to reduce the computational cost. Five 4D optimization strategies adopting different numbers of intermediate states and one three-dimensional (3D) optimization plan were studied for fifteen lung cancer patients treated with scanned protons, optimizing on all ten phases (4D_10), two extreme phases (4D_2), six phases during the exhalation stage (4D_6EX), six phases during the inhalation stage (4D_6IN), two extreme phases plus an intermediate state (4D_3) and average computed tomography image (3D), respectively. The 4D dose evaluation was conducted on all the ten phases, considering the interplay effect. The resulting doses accumulated on the reference phase were computed and compared. Compared to the 4D optimization plans, the 3D optimization plan performed inferiorly in target coverage, but superiorly in organ at risks (OARs) sparing. For the 4D optimization, all the five 4D plans showed similar performance in OARs protection. However, the 4D_6EX and 4D_6IN strategies out-performed the 4D_2 and 4D_3 plans in dose homogeneity. The computing times of the 4D_2, 4D_3, 4D_6EX and 4D_6IN approaches decreased to 32%, 41%, 66% and 67% of the 4D_10 method, respectively. Thus, our study suggested that the use of all phases during inhalation or exhalation stage might be a feasible approach substituting for the full phase strategy to reduce the calculation load while guaranteeing the plan quality for scanned proton therapy.

Deep learning-based internal gross target volume definition in 4D CT images of lung cancer patients.

BACKGROUND: Contouring of internal gross target volume (iGTV) is an essential part of treatment planning in radiotherapy to mitigate the impact of intra-fractional target motion. However, it is usually time-consuming and easily subjected to intra-observer and inter-observer variability. So far, few studies have been explored to directly predict iGTV by deep learning technique, because the iGTV contains not only the gross target volume (GTV) but also the motion information of the GTV. PURPOSE: This work was an exploratory study to present a deep learning-based framework to segment iGTV rapidly and accurately in 4D CT images for lung cancers. METHODS: Five models, including 3D UNet, mmUNet with point-wise add merging approach (mmUNet-add), mmUNet with concatenate fusion strategy (mmUNet-cat), gruUNet with point-wise add fusion approach (gruUNet-add), and gruUNet with concatenate method (gruUNet-cat), were adopted for iGTV segmentation. All the models originated from the 3D UNet network, with multi-channel multi-path and convolutional gated recurrent unit (GRU) added in the mmUNet and gruUNet networks, respectively. Seventy patients with lung cancers were collected and 55 cases were randomly selected as the training set, and 15 cases as the testing set. In addition, the segmentation results of the five models were compared with the ground truths qualitatively and quantitatively. RESULTS: In terms of Dice Similarity Coefficient (DSC), the proposed four networks (mmUNet-add, mmUNet-cat, gruUNet-add, and gruUNet-cat) increased the DSC score of 3D UNet from 0.6945 to 0.7342, 0.7253, 0.7405, and 0.7365, respectively. However, the differences were not statistically significant (p > 0.05). After a simple post-processing to remove the small isolated connected regions, the mean 95th percentile Hausdorff distances (HD_95s) of the 3D UNet, mmUNet-add, mmUNet-cat, gruUNet-add, and gruUNet-cat networks were 19.70, 15.75, 15.84, 15.61, and 15.83 mm, respectively, corresponding to 25.35, 25.96, 25.11, 28.23, and 24.47 mm before the post-processing. With regard to runtime, significant elapsed time growths (about 70s and 230s) were observed both in the mmUNet and gruUNet architectures due to the increasing parameters. But the mmUNet structure showed less growth. CONCLUSION: Our study demonstrated the ability of the deep learning technique to predict iGTVs directly. With the introduction of multi-channel multi-path and convolutional GRU, the segmentation accuracy was improved under certain conditions with a reduced segmentation efficiency and a further research topic when the 3D UNet network would lead to poor performance is elicited. Less efficiency degradation was observed in the mmUNet structure. Besides, the element-wise add fusing strategy was favorable to increase DSC, whereas HD_95 benefited from the concentrate merging approach. Nevertheless, the segmentation accuracy by deep learning still remains to be improved.

Insights into the selectivity of a brain-penetrant CDK4/6 vs CDK1/2 inhibitor for glioblastoma used in multiple replica molecular dynamics simulations.

Cyclin dependent kinases (CDKs) play an important role in cell cycle regulation and their dysfunction is associated with many cancers. That is why CDKs have been attractive targets for the treatment of cancer. Glioblastoma is a cancer caused by the aberrant expression of CDK4/6, so exploring the mechanism of the selection of CDK4/6 toward inhibitors relative to the other family members CDK1/2 is essential. In this work, multiple replica molecular dynamics (MRMD) simulations, principal component analysis (PCA), free energy landscapes (FELs), molecular mechanics Poisson-Boltzmann/Generalized Born surface area (MM-PB/GBSA) and other methods were integrated to decipher the selectively binding mechanism of the inhibitor N1J to CDK4/6 and CDK1/2. Molecular electrostatic potential (MESP) analysis provides an explanation for the N1J selectivity. Residue-based free energy decomposition reveals that most of the hot residues are located at the same location of CDKs proteins, but the different types of residues in different proteins cause changes in binding energy, which is considered as a potential developmental direction to improve the selectivity of inhibitors to CDK4/6. These results provide insights into the source of inhibitor and CDK4/6 selectivity for the future development of more selective inhibitors.Communicated by Ramaswamy H. Sarma.

Deep learning-based fast denoising of Monte Carlo dose calculation in carbon ion radiotherapy.

BACKGROUND: Plan verification is one of the important steps of quality assurance (QA) in carbon ion radiotherapy. Conventional methods of plan verification are based on phantom measurement, which is labor-intensive and time-consuming. Although the plan verification method based on Monte Carlo (MC) simulation provides a more accurate modeling of the physics, it is also time-consuming when simulating with a large number of particles. Therefore, how to ensure the accuracy of simulation results while reducing simulation time is the current difficulty and focus. PURPOSE: The purpose of this work was to evaluate the feasibility of using deep learning-based MC denoising method to accelerate carbon-ion radiotherapy plan verification. METHODS: Three models, including CycleGAN, 3DUNet and GhostUNet with Ghost module, were used to denoise the 1 × 106 carbon ions-based MC dose distribution to the accuracy of 1 × 108 carbon ions-based dose distribution. The CycleGAN's generator, 3DUNet and GhostUNet were all derived from the 3DUNet network. A total of 59 cases including 29 patients with head-and-neck cancers and 30 patients with lung cancers were collected, and 48 cases were randomly selected as the training set of the CycleGAN network and six cases as the test set. For the 3DUNet and GhostUNet models, the numbers of training set, validation set, and test set were 47, 6, and 6, respectively. Finally, the three models were evaluated qualitatively and quantitatively using RMSE and three-dimensional gamma analysis (3 mm, 3%). RESULTS: The three end-to-end trained models could be used for denoising the 1 × 106 carbon ions-based dose distribution, and their generalization was proved. The GhostUNet obtained the lowest RMSE value of 0.075, indicating the smallest difference between its denoised and 1 × 108 carbon ions-based dose distributions. The average gamma passing rate (GPR) between the GhostUNet denoising-based versus 1 × 108 carbon ions-based dose distributions was 99.1%, higher than that of the CycleGAN at 94.3% and the 3DUNet at 96.2%. Among the three models, the GhostUNet model had the fewest parameters (4.27 million) and the shortest training time (99 s per epoch) but achieved the best denoising results. CONCLUSION: The end-to-end deep network GhostUNet outperforms the CycleGAN, 3DUNet models in denoising MC dose distributions for carbon ion radiotherapy. The network requires less than 5 s to denoise a sample of MC simulation with few particles to obtain a qualitative and quantitative result comparable to the dose distribution simulated by MC with relatively large number particles, offering a significant reduction in computation time.

Ferroptosis: a novel regulated cell death participating in cellular stress response, radiotherapy, and immunotherapy.

BACKGROUND: Ferroptosis is a regulated cell death mode triggered by iron-dependent toxic membrane lipid peroxidation. As a novel cell death modality that is morphologically and mechanistically different from other forms of cell death, such as apoptosis and necrosis, ferroptosis has attracted extensive attention due to its association with various diseases. Evidence on ferroptosis as a potential therapeutic strategy has accumulated with the rapid growth of research on targeting ferroptosis for tumor suppression in recent years. METHODS: We summarize the currently known characteristics and major regulatory mechanisms of ferroptosis and present the role of ferroptosis in cellular stress responses, including ER stress and autophagy. Furthermore, we elucidate the potential applications of ferroptosis in radiotherapy and immunotherapy, which will be beneficial in exploring new strategies for clinical tumor treatment. RESULT AND CONCLUSION: Based on specific biomarkers and precise patient-specific assessment, targeting ferroptosis has great potential to be translated into practical new approaches for clinical cancer therapy, significantly contributing to the prevention, diagnosis, prognosis, and treatment of cancer.

Towards an uncertainty evaluation model for marine track considering local sampling cloud and sample information.

To evaluate the practical ability of crews during the navigation of an inward-port single ship, a track evaluation model was developed on a planar forward normal cloud chart under sample information based on the forward normal and the backward normal cloud generator. Since the track sampling cloud may be too divergent, a track belt division method based on the contributions of normal cloud drops was proposed. Combining the track evaluation model with the track belt division method, a comprehensive track evaluation scheme of the local sampling cloud based on sampling information was established. The results of an example of M.V. DAQING 257 unloaded into Dalian Port demonstrated the effectiveness of the model and showed its consistency with expert evaluation results based on subjective information. The proposed uncertainty evaluation model provides a new approach for intelligent evaluation under sample information.

Free Energy Profiles Relating With Conformational Transition of the Switch Domains Induced by G12 Mutations in GTP-Bound KRAS.

Mutations of G12 in KRAS have been involved in different cancers. Multiple replica-Gaussian accelerated molecular dynamics (MR-GaMD) simulations are applied to investigate conformational changes of the switch domains caused by G12C, G12D and G12R. Free energy landscapes suggest that G12C, G12D and G12R induce more energetic states compared to the GTP-bound WT KRAS and make the conformations of the switch domains more disordered, which disturbs bindings of KRAS to effectors. Dynamics analyses based on MR-GaMD trajectory show that G12C, G12D and G12R not only change structural flexibility of the switch domains but also affect their motion behavior, indicating that these three mutations can be used to tune the activity of KRAS. The analyses of interaction networks verify that the instability in interactions of the GTP with the switch SⅠ plays an important role in the high disorder states of the switch domain. This work is expected to provide useful information for deeply understanding the function of KRAS.