Tetrahydrofolate attenuates cognitive impairment after hemorrhagic stroke by promoting hippocampal neurogenesis via PTEN signaling.

Intracerebral hemorrhage (ICH), the most common subtype of hemorrhagic stroke, leads to cognitive impairment and imposes significant psychological burdens on patients. Hippocampal neurogenesis has been shown to play an essential role in cognitive function. Our previous study has shown that tetrahydrofolate (THF) promotes the proliferation of neural stem cells (NSCs). However, the effect of THF on cognition after ICH and the underlying mechanisms remain unclear. Here, we demonstrated that administration of THF could restore cognition after ICH. Using Nestin-GFP mice, we further revealed that THF enhanced the proliferation of hippocampal NSCs and neurogenesis after ICH. Mechanistically, we found that THF could prevent ICH-induced elevated level of PTEN and decreased expressions of phosphorylated AKT and mTOR. Furthermore, conditional deletion of PTEN in NSCs of hippocampus attenuated the inhibitory effect of ICH on the proliferation of NSCs and abnormal neurogenesis. Taken together, these results provide molecular insights into ICH-induced cognitive impairment and suggest translational clinical therapeutic strategy for hemorrhagic stroke.Significance Statement Intracerebral hemorrhage (ICH) has been associated with cognitive dysfunction, yet its underlying mechanism remains elusive. Tetrahydrofolate (THF) has shown potential in promoting the proliferation of neural stem cells (NSCs), but its specific impact on cognitive recovery following ICH is still to be confirmed. Through the utilization of the Nestin-GFP genetic marker to track endogenous NSCs in mice, our study revealed that THF could regulate PTEN pathway to ameliorate cognitive impairment post-ICH by enhancing the proliferation of NSCs and sustaining neurogenesis. These findings contribute to valuable insights into the molecular mechanisms involved and suggest potential clinical applications for enhancing cognitive function recovery after ICH.

A DNA demethylase reduces seed size by decreasing the DNA methylation of AT-rich transposable elements in soybean.

Understanding how to increase soybean yield is crucial for global food security. The genetic and epigenetic factors influencing seed size, a major crop yield determinant, are not fully understood. We explore the role of DNA demethylase GmDMEa in soybean seed size. Our research indicates that GmDMEa negatively correlates with soybean seed size. Using CRISPR-Cas9, we edited GmDMEa in the Dongnong soybean cultivar, known for small seeds. Modified plants had larger seeds and greater yields without altering plant architecture or seed nutrition. GmDMEa preferentially demethylates AT-rich transposable elements, thus activating genes and transcription factors associated with the abscisic acid pathway, which typically decreases seed size. Chromosomal substitution lines confirm that these modifications are inheritable, suggesting a stable epigenetic method to boost seed size in future breeding. Our findings provide insights into epigenetic seed size control and suggest a strategy for improving crop yields through the epigenetic regulation of crucial genes. This work implies that targeted epigenetic modification has practical agricultural applications, potentially enhancing food production without compromising crop quality.

Computational study on the endocrine-disrupting metabolic activation of Benzophenone-3 catalyzed by cytochrome P450 1A1: A QM/MM approach.

Predicting the metabolic activation mechanism and potential hazardous metabolites of environmental endocrine-disruptors is a challenging and significant task in risk assessment. Here the metabolic activation mechanism of benzophenone-3 catalyzed by P450 1A1 was investigated by using Molecular Dynamics, Quantum Mechanics/Molecular Mechanics and Density Functional Theory approaches. Two elementary reactions involved in the metabolic activation of BP-3 with P450 1A1: electrophilic addition and hydrogen abstraction reactions were both discussed. Further conversion reactions of epoxidation products, ketone products and the formaldehyde formation reaction were investigated in the non-enzymatic environment based on previous experimental reports. Binding affinities analysis of benzophenone-3 and its metabolites to sex hormone binding globulin indirectly demonstrates that they all exhibit endocrine-disrupting property. Toxic analysis shows that the eco-toxicity and bioaccumulation values of the benzophenone-3 metabolites are much lower than those of benzophenone-3. However, the metabolites are found to have skin-sensitization effects. The present study provides a deep insight into the biotransformation process of benzophenone-3 catalyzed by P450 1A1 and alerts us to pay attention to the adverse effects of benzophenone-3 and its metabolites in human livers.

Ultrasensitive Graphene Field-Effect Biosensors Based on Ferroelectric Polarization of Lithium Niobate for Breast Cancer Marker Detection.

Herein, we present a novel ultrasensitive graphene field-effect transistor (GFET) biosensor based on lithium niobate (LiNbO3) ferroelectric substrate for the application of breast cancer marker detection. The electrical properties of graphene are varied under the electrostatic field, which is generated through the spontaneous polarization of the ferroelectric substrate. It is demonstrated that the properties of interface between graphene and solution are also altered due to the interaction between the electrostatic field and ions. Compared with the graphene field-effect biosensor based on the conventional Si/SiO2 gate structure, our biosensor achieves a higher sensitivity to 64.7 mV/decade and shows a limit of detection down to 1.7 fM (equivalent to 12 fg·mL-1) on the detection of microRNA21 (a breast cancer marker). This innovative design combining GFETs with ferroelectric substrates holds great promise for developing an ultrahigh-sensitivity biosensing platform based on graphene that enables rapid and early disease diagnosis.

Characterization of VOC emissions and health risk assessment in the plastic manufacturing industry.

Volatile organic compounds (VOCs) significantly contribute to ozone pollution formation, and many VOCs are known to be harmful to human health. Plastic has become an indispensable material in various industries and daily use scenarios, yet the VOC emissions and associated health risks in the plastic manufacturing industry have received limited attention. In this study, we conducted sampling in three typical plastic manufacturing factories to analyze the emission characteristics of VOCs, ozone formation potential (OFP), and health risks for workers. Isopropanol was detected at relatively high concentrations in all three factories, with concentrations in organized emissions reaching 322.3 µg/m3, 344.8 µg/m3, and 22.6 µg/m3, respectively. Alkanes are the most emitted category of VOCs in plastic factories. However, alkenes and oxygenated volatile organic compounds (OVOCs) exhibit higher OFP. In organized emissions of different types of VOCs in the three factories, alkenes and OVOCs contributed 22.8%, 67%, and 37.8% to the OFP, respectively, highlighting the necessity of controlling them. The hazard index (HI) for all three factories was less than 1, indicating a low non-carcinogenic toxic risk; however, there is still a possibility of non-cancerous health risks in two of the factories, and a potential lifetime cancer risk in all of the three factories. For workers with job tenures exceeding 5 years, there may be potential health risks, hence wearing masks with protective capabilities is necessary. This study provides evidence for reducing VOC emissions and improving management measures to ensure the health protection of workers in the plastic manufacturing industry.

Climate-related naturally occurring epimutation and their roles in plant adaptation in A. thaliana.

DNA methylation has been proposed to be an important mechanism that allows plants to respond to their environments sometimes entirely uncoupled from genetic variation. To understand the genetic basis, biological functions and climatic relationships of DNA methylation at a population scale in Arabidopsis thaliana, we performed a genome-wide association analysis with high-quality single nucleotide polymorphisms (SNPs), and found that ~56% on average, especially in the CHH sequence context (71%), of the differentially methylated regions (DMRs) are not tagged by SNPs. Among them, a total of 3235 DMRs are significantly associated with gene expressions and potentially heritable. 655 of the 3235 DMRs are associated with climatic variables, and we experimentally verified one of them, HEI10 (HUMAN ENHANCER OF CELL INVASION NO.10). Such epigenetic loci could be subjected to natural selection thereby affecting plant adaptation, and would be expected to be an indicator of accessions at risk. We therefore incorporated these climate-related DMRs into a gradient forest model, and found that the natural A. thaliana accessions in Southern Europe that may be most at risk under future climate change. Our findings highlight the importance of integrating DNA methylation that is independent of genetic variations, and climatic data to predict plants' vulnerability to future climate change.

A Retrospective Study of Biological Risk Factors Associated with Primary Knee Osteoarthritis and the Development of a Nomogram Model.

Aim: A high percentage of the elderly suffer from knee osteoarthritis (KOA), which imposes a certain economic burden on them and on society as a whole. The purpose of this study is to examine the risk of KOA and to develop a KOA nomogram model that can timely intervene in this disease to decrease patient psychological burdens. Methods: Data was collected from patients with KOA and without KOA at our hospital from February 2021 to February 2023. Initially, a comparison was conducted between the variables, identifying statistical differences between the two groups. Subsequently, the risk of KOA was evaluated using the Least Absolute Shrinkage and Selection Operator method and multivariate logistic regression to determine the most effective predictive index and develop a prediction model. The examination of the disease risk prediction model in KOA includes the corresponding nomogram, which encompasses various potential predictors. The assessment of disease risk entails the application of various metrics, including the consistency index (C index), the area under the curve (AUC) of the receiver operating characteristic curve, the calibration chart, the GiViTi calibration band, and the model for predicting KOA. Furthermore, the potential clinical significance of the model is explored through decision curve analysis (DCA) and clinical influence curve analysis. Results: The study included a total of 582 patients, consisting of 392 patients with KOA and 190 patients without KOA. The nomogram utilized age, haematocrit, platelet count, apolipoprotein a1, potassium, magnesium, hydroxybutyrate dehydrogenase, creatine kinase, and estimated glomerular filtration rate as predictors. The C index, AUC, calibration plot, Giviti calibration band, DCA and clinical influence KOA indicated the ability of nomogram model to differentiate KOA. Conclusion: Using nomogram based on disease risk, high-risk KOA can be identified directly without imaging.

Updating vehicle VOCs emissions characteristics under clean air actions in a tropical city of China.

In the context of clean air actions in China, vehicle emission limits have been continuously tightened, which has facilitated the reduction of volatile organic compounds (VOCs) emissions. However, the characteristics of VOC emissions from vehicles with strict emission limits are poorly understood. This study investigated the VOC emission characteristics from vehicles under the latest standards based on tunnel measurements, and identified future control strategies for vehicle emissions. The results showed that the highest percentage of VOCs from vehicle consisted of alkanes (80.9 %), followed by aromatics (15.8 %) and alkenes (3.1 %). Alkanes had the most significant ozone formation potential due to their high concentrations, in contrast to the aromatics that have been dominant in previous studies. The measured fleet-average VOC emission factor was 71.3 mg·km-1, including tailpipe emissions of 39.6 mg·km-1 and evaporative emissions of 31.7 mg·km-1. The VOC emission factors of the subgroups were obtained. The emission of evaporated VOCs accounted for 44.5 % of the total vehicle VOC emissions, which have increased substantially from previous studies. In addition, the emission characteristics of vehicles that are under the latest emission threshold values have changed significantly, and the mixing ratio of toluene/benzene (T/B) has been updated to 3:1. This study updates the VOCs emission factors of vehicles under clean air actions and highlights the future mitigation policies should focus on reducing evaporative VOC emissions.

Semirational Design Strategy To Enhance the Thermostability and Catalytic Activity of Cytochrome P450 105D7 for the Degradation of the Pharmaceutically Active Compounds: Diclofenac.

Pharmaceutically active compounds are an important category of emerging pollutants, and their biological transformation processes in the environment are crucial for understanding and evaluating the migration, transformation, and environmental fate of emerging pollutants. The cytochrome P450 105 enzyme family has been proven to play an important role in the degradation of exogenous environmental pollutants. However, its thermostability and catalytic activity still need to be improved to better adapt to complex environmental conditions. This work elucidates the key mechanisms and important residues of the degradation reaction through multiple computational strategies, establishes a mutation library, and obtains 21 single-point mutation designs. Experimental verification showed that 16 single mutants had enhanced thermostability, with the R89F and L197Y mutants showing the highest increases in thermostability at 135 and 119% relative to the wild-type enzyme, respectively. Additionally, as a result of the higher specific activity of D390Q, it was selected for combination mutagenesis, ultimately resulting in three combination mutants (R89F/L197Y, R89F/D390Q, and R89F/L197Y/D390Q) with enhanced thermostability and catalytic activity. This study provides a modification approach for constructing efficient enzyme variants through semirational design and can contribute to the development of control technologies for emerging pollutants.

Ultrasensitive 3D Stacked Silicon Nanosheet Field-Effect Transistor Biosensor with Overcoming Debye Shielding Effect for Detection of DNA.

Silicon nanowire field effect (SiNW-FET) biosensors have been successfully used in the detection of nucleic acids, proteins and other molecules owing to their advantages of ultra-high sensitivity, high specificity, and label-free and immediate response. However, the presence of the Debye shielding effect in semiconductor devices severely reduces their detection sensitivity. In this paper, a three-dimensional stacked silicon nanosheet FET (3D-SiNS-FET) biosensor was studied for the high-sensitivity detection of nucleic acids. Based on the mainstream Gate-All-Around (GAA) fenestration process, a three-dimensional stacked structure with an 8 nm cavity spacing was designed and prepared, allowing modification of probe molecules within the stacked cavities. Furthermore, the advantage of the three-dimensional space can realize the upper and lower complementary detection, which can overcome the Debye shielding effect and realize high-sensitivity Point of Care Testing (POCT) at high ionic strength. The experimental results show that the minimum detection limit for 12-base DNA (4 nM) at 1 × PBS is less than 10 zM, and at a high concentration of 1 µM DNA, the sensitivity of the 3D-SiNS-FET is approximately 10 times higher than that of the planar devices. This indicates that our device provides distinct advantages for detection, showing promise for future biosensor applications in clinical settings.

Optimization of CRISPR-Cas9 system in Eustoma grandiflorum.

The optimization of the CRISPR-Cas9 system for enhancing editing efficiency holds significant value in scientific research. In this study, we optimized single guide RNA and Cas9 promoters of the CRISPR-Cas9 vector and established an efficient protoplast isolation and transient transformation system in Eustoma grandiflorum, and we successfully applied the modified CRISPR-Cas9 system to detect editing efficiency of the EgPDS gene. The activity of the EgU6-2 promoter in E. grandiflorum protoplasts was approximately three times higher than that of the GmU6 promoter. This promoter, along with the EgUBQ10 promoter, was applied in the CRISPR-Cas9 cassette, the modified CRISPR-Cas9 vectors that pEgU6-2::sgRNA-2/pEgUBQ10::Cas9-2 editing efficiency was 37.7%, which was 30.3% higher than that of the control, and the types of mutation are base substitutions, small fragment deletions and insertions. Finally we obtained an efficient gene editing vector for E. grandiflorum. This project provides an important technical platform for the study of gene function in E. grandiflorum.

New structure transistors for advanced technology node CMOS ICs.

Over recent decades, advancements in complementary metal-oxide-semiconductor integrated circuits (ICs) have mainly relied on structural innovations in transistors. From planar transistors to the fin field-effect transistor (FinFET) and gate-all-around FET (GAAFET), more gate electrodes have been added to three-dimensional (3D) channels with enhanced control and carrier conductance to provide higher electrostatic integrity and higher operating currents within the same device footprint. Beyond the 1-nm node, Moore's law scaling is no longer expected to be applicable to geometrical shrinkage. Vertical transistor stacking, e.g. in complementary FETs (CFET), 3D stack (3DS) FETs and vertical-channel transistors (VFET), for enhanced density and variable circuit or system design represents a revolutionary scaling approach for sustained IC development. Herein, innovative works on specific structures, key process breakthroughs, shrinking cell sizes and design methodologies for transistor structure research and development are reviewed. Perspectives on future innovations in advanced transistors with new channel materials and operating theories are also discussed.

New solution for predicting the removal efficiency of organic contaminants by UV/ H2O2: From a case study of 1H-benzotriazole.

Rapid prediction of the removal efficiency and energy consumption of organic contaminants under various operating conditions is crucial for advanced oxidation processes (AOPs) in industrial application. In this study, 1H-Benzotriazole (BTZ, CAS: 95-14-7) is selected as a model micropollutant, a validated incorporated Computational Fluid Dynamics (CFD) model is employed to comprehensively investigate the impacts of initial concentrations of H2O2, BTZ and dissolved organic carbon (DOC) (i.e., [DOC]0, [BTZ]0 and [DOC]0), as well as the effective UV lamp power P and volumetric flow rate Qv. Generally, the operation performance depends on [DOC]0 and [BTZ]0 in similar trends, but with quantitatively different ways. The increase in [H2O2]0 and P/Qv can promote •OH generation, leading to the elimination of BTZ. It is worth noting that P/Qv is found to be linearly correlated with the removal order of BTZ (ROBTZ) under specific conditions. Based on this finding, the degradation of other potential organic contaminants with a wide range of rate constants by UV/H2O2 is further investigated. A model for predicting energy consumption for target removal rates of organic pollutants is established from massive simulation data for the first time. Additionally, a handy Matlab app is first developed for convenient application in water treatment. This work proposes a new operable solution for fast predicting operation performance and energy consumption for the removal of organic contaminants in industrial applications of advanced oxidation processes.

Reinforced hydrogel building via formation of alginate-chitosan double network with pH & salt-responsiveness and electric conductivity for soft actuators.

Aiming at green and friendly environmental protection, polyvinyl alcohol/sodium alginate/chitosan (PSCS) double network hydrogel was successfully prepared through diffusing the high molecular weight chitosan into PVA/sodium alginate (PS) hydrogel without any other toxic reagents. The polyanion hydrogels could be significantly enhanced by immersing the polyanion hydrogel in high molecular weight chitosan solution without requiring specific structure. The PSCS hydrogel had a compact and rough surface structure with the smaller porosities and larger crystallization degree compared with polyvinyl alcohol/sodium alginate hydrogels and polyvinyl alcohol/sodium alginate/Ca2+ (PSCa) hydrogels. The PSCS hydrogel possessed excellent hydrolysis resistance, the significant pH-sensitive and salt-sensitive swelling. In addition, the flexibility, Young's modulus and mechanical properties of PSCS hydrogel can be adjusted through the changing the content of sodium alginate. Moreover, PS, PSCa and PSCS had electric conductivity, and PSCS showed twice the conductivity compared to PS hydrogel. Based on differences of swelling ratio, a PSCS bilayer hydrogel was designed and showed excellent pH-driven deformation ability. The PSCS hydrogel is expected to expand the application of hydrogels in conditions involving stimulus response, and might serve as a promising intelligent actuators or soft robots.

PM2.5 pollution modulates the response of ozone formation to VOC emitted from various sources: Insights from machine learning.

Numerous studies have linked ozone (O3) production to its precursors and fine particulate matter (PM2.5), while the complex interaction effects of PM2.5 and volatile organic compounds (VOCs) on O3 remain poorly understood. A systematic approach based on an interpretable machine learning (ML) model was utilized to evaluate the primary driving factors that impact O3 and to elucidate how changes in PM2.5, VOCs from different sources, NOx, and meteorological conditions either promote or inhibit O3 formation through their individual and synergistic effects in a tropical coastal city, Haikou, from 2019 to 2020. The results suggest that under low PM2.5 levels, alongside the linear O3-PM2.5 relationship observed, O3 formation is suppressed by PM2.5 with higher proportions of traffic-derived aerosol. Vehicle VOC emissions contributed maximally to O3 formation at midday, despite the lowest concentration. VOCs from fossil fuel combustion and industry emissions, which have opposing effects on O3, act as inhibitors and promoters by inducing diverse photochemical regimes. As PM2.5 pollution escalates, the impact of these VOCs reverses, becoming more pronounced in shaping O3 variation. Sensitivity analysis reveals that the O3 formation regime is VOC-limited, and effective regional O3 mitigation requires prioritizing substantial VOC reductions to offset enhanced VOC sensitivity induced by the co-reduction in PM2.5, with a focus on industrial and vehicular emissions, and subsequently, fossil fuel combustion once PM2.5 is effectively controlled. This study underscores the potential of the SHAP-based ML approach to decode the intricate O3-NOx-VOCs-PM2.5 interplay, considering both meteorological and atmospheric compositional variations.

New mechanism for the participation of aromatic oxidation products in atmospheric nucleation.

Oxygenated organic molecules (OOMs) are recognized as important precursors for new particle formation (NPF) in the urban atmosphere. The paper theoretically studied the formation of OOMs by styrene oxidation processes initiated by OH radicals, focusing on the OOMs nucleation mechanism. The results found that in the presence of an H2SO4 molecule, lowly oxygenated organic molecules containing a benzene ring (LOMBs) can form stable clusters and grow to the scale of a critical nucleus through pi-pi stacking and OH hydrogen bonding. In addition, LOMBs are more readily generated in a styrene-oxidized system in the presence/absence of NOx than highly oxygenated organic molecules (HOMs). The reaction of OH radicals with other aromatics containing a branched chain on the benzene ring produces LOMBs to varying degrees, with pi-pi stacking playing an essential role. This result suggests that, in the presence of H2SO4 molecules, LOMBs may play a more significant role in promoting nucleation than HOMs. Our findings serve as a pivotal foundation for future investigations into the oxidation and nucleation processes of diverse aromatics in urban environments.

[Ozone Pollution Characteristics and Formation Mechanism in a Typical Tropical Seaside City].

To investigate the pollution characteristics and formation mechanism of ambient air ozone(O3) in a typical tropical seaside city, we conducted an observational experiment on O3 and its precursors at an urban site in Haikou, Hainan Province, from June to October 2019. The O3 pollution characteristics were analyzed comprehensively; the O3-NOx-VOCs sensitivities and key precursors were determined, and the control strategies for O3 pollution were carried out. The results were as follows:1 O3 pollution in Haikou mainly occurred in September and October, with daily maximum 8-h O3 concentrations in the range of 39-190 µg·m-3, and the daily variation in O3 was unimodal, peaking at approximately 14:00. 2 The concentrations of NO2 and VOCs were higher during O3 pollution episodes than their respective mean values in Haikou City. The increased O3 precursor concentrations were an important factor leading to O3 pollution, whereas O3 pollution was also influenced by regional transport, with pollutants mainly transported from the northeastern part of Haikou City. 3 O3-NOx-VOCs sensitivity in Haikou City was in the VOCs and NOx transitional regime, and the most sensitive precursors in various months were different. O3 formation in September was sensitive to anthropogenic VOCs the most; however, in October it was sensitive to NOx. 4 In the future, the reduction ratio of VOCs to NOx should be 1:1-4:1 to control O3 pollution effectively in Haikou.

New insight into biodegradation mechanism of phenylurea herbicides by cytochrome P450 enzymes: Successive N-demethylation mechanism.

Phenylurea herbicides (PUHs) present one of the most important herbicides, which have cause serious effects on ecological environment and humans. Nowadays enzyme strategy shows great advantages in degradation of PUHs. Here density functional theory (DFT), quantitative structure - activity relationship (QSAR) and quantum mechanics/molecular mechanics (QM/MM) approaches are used to investigate the degradation mechanism of PUHs catalyzed by P450 enzymes. Two successive N-demethylation pathways are identified and two hydrogen abstraction (H-abstraction) reaction pathways are identified as the rate-determining step through high-throughput DFT calculations. The Boltzmann-weighted average energy barrier of the second H-abstraction pathway (19.95 kcal/mol) is higher than that of the first H-abstraction pathway (16.80 kcal/mol). Two QSAR models are established to predict the energy barriers of the two H-abstraction pathways based on the quantum chemical descriptors and mordred molecular descriptors. The determination coefficient (R2) values of QSAR models are > 0.9, which reveal that the established QSAR models have great predictive capability. QM/MM calculations indicate that human P450 enzymes are more efficient in degradation of PUHs than crop and weed P450 enzymes. Correlations between energy barriers and key structural/charge parameters are revealed and key parameters that have influence on degradation efficiency of PUHs are identified. This study provides lateral insights into the biodegradation strategy and removal method of PUHs and valuable information for designing or engineering of highly efficient degradation enzymes and genetically modified crops.

Impacts of QM region sizes and conformation numbers on modelling enzyme reactions: a case study of polyethylene terephthalate hydrolase.

A quantum mechanics/molecular mechanics (QM/MM) approach is a broadly used tool in computational enzymology. Treating the QM region with a high-level DFT method is one of the important branches. Here, taking leaf-branch compost cutinase-catalyzed polyethylene terephthalate depolymerization as an example, the convergence behavior of energy barriers as well as key structural and charge features with respect to the size of the QM region (up to 1000 atoms) is systematically investigated. BP86/6-31G(d)//CHARMM and M06-2X/6-311G(d,p)//CHARMM level of theories were applied for geometry optimizations and single-point energy calculations, respectively. Six independent enzyme conformations for all the four catalytic steps (steps (i)-(iv)) were considered. Most of the twenty-four cases show that at least 500 QM atoms are needed while only two rare cases show that ∼100 QM atoms are sufficient for convergence when only a single conformation was considered. This explains why most previous studies showed that 500 or more QM atoms are required while a few others showed that ∼100 QM atoms are sufficient for DFT/MM calculations. More importantly, average energy barriers and key structural/charge features from six conformations show an accelerated convergence than that in a single conformation. For instance, to reach energy barrier convergence (within 2.0 kcal mol-1) for step (ii), only ∼100 QM atoms are required if six conformations are considered while 500 or more QM atoms are needed with a single conformation. The convergence is accelerated to be more rapid if hundreds and thousands of conformations were considered, which aligns with previous findings that only several dozens of QM atoms are required for convergence with semi-empirical QM/MM MD simulations.