Enhancing the carbon content of coal gangue for composting through sludge amendment: A feasibility study

Cocomposting coal gangue and sludge eliminates the challenge of utilizing coal gangue. However, there is limited understanding about the feasibility of cocomposting sludge and coal gangue, as well as the composting indicators, functional microorganisms, and safety risks involved. Therefore, this study evaluated the feasibility of enhancing carbon composting in coal gangue by incorporating sludge along with sawdust as a conditioner. Three laboratory-scale reactors were designed and labeled as T1 (20 % coal gangue, 60 % sludge, and 20 % sawdust), T2 (40 % coal gangue, 40 % sludge, and 20 % sawdust), and T3 (60 % coal gangue, 20 % sludge, and 20 % sawdust). Seed germination and plant growth assessments were conducted to ensure compost stability and assess phytotoxicity to cabbage (Brassica rapa chinensis L.) in terms of growth and biomass. The results indicated that the temperature, pH, EC and ammonia nitrogen of all three reactor conditions met the requirements for product decomposition. Composting was successfully achieved when the sludge proportion was 20 % (T3). However, when the sludge proportion was markedly high (T1), the harmlessness of the compost was reduced. The germination indices of T1, T2, and T3 reached 95 %, 122 %, and 119 % at maturity, respectively. This confirmed that the harmless cycle, which involved promoting condensation and aromatization, enhancing decay, and reducing composting time, was shorter in T2 and T3 than in T1. Coal gangue can also serve as a beneficial habitat for microorganisms, promoting an increase in their population and activity. Potting experiments in sandy soil revealed that the mechanism of action of compost products in soil included not only the enhancement of soil nutrients but also the improvement of soil texture. The results of this study suggest that using coal gangue as a raw material for composting is an efficient and environmentally friendly approach for producing organic fertilizers.

A comprehensive review of advances in hepatocyte microencapsulation: selecting materials and preserving cell viability

Liver failure represents a critical medical condition with a traditionally grim prognosis, where treatment options have been notably limited. Historically, liver transplantation has stood as the sole definitive cure, yet the stark disparity between the limited availability of liver donations and the high demand for such organs has significantly hampered its feasibility. This discrepancy has necessitated the exploration of hepatocyte transplantation as a temporary, supportive intervention. In light of this, our review delves into the burgeoning field of hepatocyte transplantation, with a focus on the latest advancements in maintaining hepatocyte function, co-microencapsulation techniques, xenogeneic hepatocyte transplantation, and the selection of materials for microencapsulation. Our examination of hepatocyte microencapsulation research highlights that, to date, most studies have been conducted in vitro or using liver failure mouse models, with a notable paucity of experiments on larger mammals. The functionality of microencapsulated hepatocytes is primarily inferred through indirect measures such as urea and albumin production and the rate of ammonia clearance. Furthermore, research on the mechanisms underlying hepatocyte co-microencapsulation remains limited, and the practicality of xenogeneic hepatocyte transplantation requires further validation. The potential of hepatocyte microencapsulation extends beyond the current scope of application, suggesting a promising horizon for liver failure treatment modalities. Innovations in encapsulation materials and techniques aim to enhance cell viability and function, indicating a need for comprehensive studies that bridge the gap between small-scale laboratory success and clinical applicability. Moreover, the integration of bioengineering and regenerative medicine offers novel pathways to refine hepatocyte transplantation, potentially overcoming the challenges of immune rejection and ensuring the long-term functionality of transplanted cells. In conclusion, while hepatocyte microencapsulation and transplantation herald a new era in liver failure therapy, significant strides must be made to translate these experimental approaches into viable clinical solutions. Future research should aim to expand the experimental models to include larger mammals, thereby providing a clearer understanding of the clinical potential of these therapies. Additionally, a deeper exploration into the mechanisms of cell survival and function within microcapsules, alongside the development of innovative encapsulation materials, will be critical in advancing the field and offering new hope to patients with liver failure.

Constant-pH Simulations with the Polarizable Atomic Multipole AMOEBA Force Field

Accurately predicting protein behavior across diverse pH environments remains a significant challenge in biomolecular simulations. Existing constant-pH molecular dynamics (CpHMD) algorithms are limited to fixed-charge force fields, hindering their application to biomolecular systems described by permanent atomic multipoles or induced dipoles. This work overcomes these limitations by introducing the first polarizable CpHMD algorithm in the context of the Atomic Multipole Optimized Energetics for Biomolecular Applications (AMOEBA) force field. Additionally, our implementation in the open-source Force Field X (FFX) software has the unique ability to handle titration state changes for crystalline systems including flexible support for all 230 space groups. The evaluation of constant-pH molecular dynamics (CpHMD) with the AMOEBA force field was performed on 11 crystalline peptide systems that span the titrating amino acids (Asp, Glu, His, Lys, and Cys). Titration states were correctly predicted for 15 out of the 16 amino acids present in the 11 systems, including for the coordination of Zn2+ by cysteines. The lone exception was for a HIS-ALA peptide where CpHMD predicted both neutral histidine tautomers to be equally populated, whereas the experimental model did not consider multiple conformers and diffraction data are unavailable for rerefinement. This work demonstrates the promise polarizable CpHMD simulations for pKa predictions, the study of biochemical mechanisms such as the catalytic triad of proteases, and for improved protein-ligand binding affinity accuracy in the context of pharmaceutical lead optimization.

Oxygen tracer diffusion in amorphous hafnia films for resistive memory

The oxygen diffusion rate in hafnia (HfO2)-based resistive memory plays a pivotal role in enabling nonvolatile data retention. However, the information retention times obtained in HfO2 resistive memory devices are many times higher than the expected values obtained from oxygen diffusion measurements in HfO2 materials. In this study, we resolve this discrepancy by conducting oxygen isotope tracer diffusion measurements in amorphous hafnia (a-HfO2) thin films. Our results show that the oxygen tracer diffusion in amorphous HfO2 films is orders of magnitude lower than that of previous measurements on monoclinic hafnia (m-HfO2) pellets. Moreover, oxygen tracer diffusion is much lower in denser a-HfO2 films deposited by atomic layer deposition (ALD) than in less dense a-HfO2 films deposited by sputtering. The ALD films yield similar oxygen diffusion times as experimentally measured device retention times, reconciling this discrepancy between oxygen diffusion and retention time measurements. More broadly, our work shows how processing conditions can be used to control oxygen transport characteristics in amorphous materials without long-range crystal order.

Incorporating Neural Networks into the AMOEBA Polarizable Force Field

Neural network potentials (NNPs) offer significant promise to bridge the gap between the accuracy of quantum mechanics and the efficiency of molecular mechanics in molecular simulation. Most NNPs rely on the locality assumption that ensures the model's transferability and scalability and thus lack the treatment of long-range interactions, which are essential for molecular systems in the condensed phase. Here we present an integrated hybrid model, AMOEBA+NN, which combines the AMOEBA potential for the short- and long-range noncovalent atomic interactions and an NNP to capture the remaining local covalent contributions. The AMOEBA+NN model was trained on the conformational energy of the ANI-1x data set and tested on several external data sets ranging from small molecules to tetrapeptides. The hybrid model demonstrated substantial improvements over the baseline models in term of accuracy as the molecule size increased, suggesting its potential as a next-generation approach for chemically accurate molecular simulations.

RUNX3 inhibits KSHV lytic replication by binding to the viral genome and repressing transcription

Kaposi's sarcoma-associated herpesvirus (KSHV) belongs to the gamma herpesvirus family, which can cause human malignancies including Kaposi sarcoma, primary effusion lymphoma, and multicentric Castleman's diseases. KSHV typically maintains a persistent latent infection within the host. However, after exposure to intracellular or extracellular stimuli, KSHV lytic replication can be reactivated. The reactivation process of KSHV triggers the innate immune response to limit viral replication. Here, we found that the transcriptional regulator RUNX3 is transcriptionally upregulated by the NF-κB signaling pathway in KSHV-infected SLK cells and B cells during KSHV reactivation. Notably, knockdown of RUNX3 significantly promotes viral lytic replication as well as the gene transcription of KSHV. Consistent with this finding, overexpression of RUNX3 impairs viral lytic replication. Mechanistically, RUNX3 binds to the KSHV genome and limits viral replication through transcriptional repression, which is related to its DNA- and ATP-binding ability. However, KSHV has also evolved corresponding strategies to antagonize this inhibition by using the viral protein RTA to target RUNX3 for ubiquitination and proteasomal degradation. Altogether, our study suggests that RUNX3, a novel host-restriction factor of KSHV that represses the transcription of viral genes, may serve as a potential target to restrict KSHV transmission and disease development.IMPORTANCEThe reactivation of Kaposi's sarcoma-associated herpesvirus (KSHV) from latent infection to lytic replication is important for persistent viral infection and tumorigenicity. However, reactivation is a complex event, and the regulatory mechanisms of this process are not fully elucidated. Our study revealed that the host RUNX3 is upregulated by the NF-κB signaling pathway during KSHV reactivation, which can repress the transcription of KSHV genes. At the late stage of lytic replication, KSHV utilizes a mechanism involving RTA to degrade RUNX3, thus evading host inhibition. This finding helps elucidate the regulatory mechanism of the KSHV life cycle and may provide new clues for the development of therapeutic strategies for KSHV-associated diseases.

Gga-let-7a-3p inhibits the proliferation and differentiation of chicken intramuscular preadipocytes

1. Intramuscular fat (IMF) is a key parameter for chicken meat quality. IMF deposition is driven by genetic, nutritional and management factors, with genetics being the determining factor. Previous whole transcriptome sequencing revealed that microRNA gga-let-7a-3p was related to lipid metabolism in breast muscle. This study further investigated the potential role of gga-let-7a-3p in IMF deposition.2. The mimic and inhibitor of gga-let-7a-3p were individually transfected into chicken intramuscular preadipocytes. Subsequently, the proliferation and differentiation states of the cells were detected. Transcriptome sequencing was performed on cells transfected with gga-let-7a-3p mimic.3. The results indicated that gga-let-7a-3p suppressed the mRNA levels of proliferation and differentiation-related genes, as well as the protein levels. EdU and Oil Red O assays revealed that gga-let-7a-3p restrained preadipocyte proliferation and differentiation. In addition, a total of 333 up-regulated genes and 807 down-regulated genes were identified in cells transfected with gga-let-7a-3p mimic. Using Kyoto Encyclopaedia of Genes and Genomes (KEGG) and Gene Ontology (GO) analysis, differential genes were found to be enriched in processes such as the peroxisome proliferator activated receptor (PPAR) signalling pathway and oxidative phosphorylation.4. The study demonstrated that gga-let-7a-3p inhibits the proliferation and differentiation of chicken intramuscular preadipocytes, which provides new understanding to further unravel the function of gga-let-7a-3p.

Unleashing the Potential of 1,3-Diketone Analogues as Selective LH2 Inhibitors

Lysyl hydroxylase 2 (LH2) catalyzes the formation of highly stable hydroxylysine aldehyde-derived collagen cross-links (HLCCs), thus promoting lung cancer metastasis through its capacity to modulate specific types of collagen cross-links within the tumor stroma. Using 1 and 2 from our previous high-throughput screening (HTS) as lead probes, we prepared a series of 1,3-diketone analogues, 1-18, and identified 12 and 13 that inhibit LH2 with IC50's of approximately 300 and 500 nM, respectively. Compounds 12 and 13 demonstrate selectivity for LH2 over LH1 and LH3. Quantum mechanics/molecular mechanics (QM/MM) modeling indicates that the selectivity of 12 and 13 may stem from noncovalent interactions like hydrogen bonding between the morpholine/piperazine rings with the LH2-specific Arg661. Treatment of 344SQ WT cells with 13 resulted in a dose-dependent reduction in their migration potential, whereas the compound did not impede the migration of the same cell line with an LH2 knockout (LH2KO).

Negative regulation of angiogenesis and the MAPK pathway may be a shared biological pathway between IS and epilepsy

Ischemia stroke and epilepsy are two neurological diseases that have significant patient and societal burden, with similar symptoms of neurological deficits. However, the underlying mechanism of their co-morbidity are still unclear. In this study, we performed a combined analysis of six gene expression profiles (GSE58294, GSE22255, GSE143272, GSE88723, GSE163654, and GSE174574) to reveal the common mechanisms of IS and epilepsy. In the mouse datasets, 74 genes were co-upregulated and 7 genes were co-downregulated in the stroke and epilepsy groups. Further analysis revealed that the co-expressed differentially expressed genes (DEGs) were involved in negative regulation of angiogenesis and the MAPK signaling pathway, and this was verified by Gene Set Enrichment Analysis of human datasets and single cell RNA sequence of middle cerebral artery occlusion mice. In addition, combining DEGs of human and mouse, PTGS2, TMCC3, KCNJ2, and GADD45B were identified as cross species conserved hub genes. Meanwhile, molecular docking results revealed that trichostatin A and valproic acid may be potential therapeutic drugs. In conclusion, to our best knowledge, this study conducted the first comorbidity analysis of epilepsy and ischemic stroke to identify the potential common pathogenic mechanisms and drugs. The findings may provide an important reference for the further studies on post-stroke epilepsy.

Impact of Learning Burnout on Mobile Phone Dependence Among Adolescents in Western China Under the "Double Reduction" Policy: The Mediating Role of Social Support

Purpose

This paper investigated the prevalence of mobile phone dependence (MPD) and its associated with learning burnout under the "double reduction" policy among adolescents in Guizhou Province in western China. In addition, the influence of the mediating mechanism of social support on this relationship was investigated.

Methods

The sample was collected from 16,216 adolescents in West China's Guizhou province, from December 2021 to January 2022 via multistage stratified random sampling. The Self-rating Questionnaire for Adolescent Problematic Mobile Phone Use (SQAPMPU) was used to assess the MPD, the Adolescent Student Burnout Scale (ASBI) was used to assess the learning burnout, and the Social Support Scale (SSS) was used to assess the social support. A hierarchical linear regression model was used to analyze the relationship between MPD, learning burnout, and social support. The mediating effect of social support between MPD and learning burnout was analyzed by structural equation model.

Results

Prevalence of MPD was 26.4% among adolescents in Guizhou province in western China. After adjusting for confounding variables like demographics, multiple linear regression model has revealed that learning burnout positively predicted MPD and social support negatively predicted MPD. The structural equation model showed that 10.9% of the effect was explained by the mediating effect of social support.

Conclusion

These findings could inform service delivery and policy formulation to reduce learning and avoid MPD in adolescents.

A generalized Kirkwood implicit solvent for the polarizable AMOEBA protein model

Computational simulation of biomolecules can provide important insights into protein design, protein-ligand binding interactions, and ab initio biomolecular folding, among other applications. Accurate treatment of the solvent environment is essential in such applications, but the use of explicit solvents can add considerable cost. Implicit treatment of solvent effects using a dielectric continuum model is an attractive alternative to explicit solvation since it is able to describe solvation effects without the inclusion of solvent degrees of freedom. Previously, we described the development and parameterization of implicit solvent models for small molecules. Here, we extend the parameterization of the generalized Kirkwood (GK) implicit solvent model for use with biomolecules described by the AMOEBA force field via the addition of corrections to the calculation of effective radii that account for interstitial spaces that arise within biomolecules. These include element-specific pairwise descreening scale factors, a short-range neck contribution to describe the solvent-excluded space between pairs of nearby atoms, and finally tanh-based rescaling of the overall descreening integral. We then apply the AMOEBA/GK implicit solvent to a set of ten proteins and achieve an average coordinate root mean square deviation for the experimental structures of 2.0 Å across 500 ns simulations. Overall, the continued development of implicit solvent models will help facilitate the simulation of biomolecules on mechanistically relevant timescales.

Molecular Cage Reports on Its Contents: Spectroscopic Signatures of Cryo-Cooled K+- and Ba2+-Benzocryptand Complexes

UV photofragment spectroscopy and IR-UV double resonance methods are used to determine the structure and spectroscopic responses of a three-dimensional [2.2.2]-benzocryptand cage to the incorporation of a single K⁺ or Ba²⁺ imbedded inside it (labeled as K⁺-BzCrypt, Ba²⁺-BzCrypt). We studied the isolated ion-cryptand complex under cryo-cooled conditions, brought into the gas phase by nano-electrospray ionization. Incorporation of a phenyl ring in place of the central ethyl group in one of the three N-CH₂-CH₂-O-CH₂-CH₂-O-CH₂-CH₂-N chains provides a UV chromophore whose S₀-S₁ transition we probe. K⁺-BzCrypt and Ba²⁺-BzCrypt have their S₀-S₁ origin transitions at 35,925 and 36,446 cm^-1, respectively, blue-shifted by 174 and 695 cm^-1 from that of 1,2-dimethoxybenzene. These origins are used to excite a single conformation of each complex selectively and record their IR spectra using IR-UV dip spectroscopy. The alkyl CH stretch region (2800-3000 cm^-1) is surprisingly sensitive to the presence and nature of the encapsulated ion. We carried out an exhaustive conformational search of cage conformations for K⁺-BzCrypt and Ba²⁺-BzCrypt, identifying two conformations (A and B) that lie below all others in energy. We extend our local mode anharmonic model of the CH stretch region to these strongly bound ion-cage complexes to predict conformation-specific alkyl CH stretch spectra, obtaining quantitative agreement with experiment for conformer A, the gas-phase global minimum. The large electrostatic effect of the charge on the O- and N-lone pairs affects the local mode frequencies of the CH₂ groups adjacent to these atoms. The localized CH₂ scissors modes are pushed up in frequency by the adjacent O/N-atoms so that their overtones have little effect on the alkyl CH stretch region. However, the localized CH₂ wags are nearly degenerate and strongly coupled to one another, producing an array of delocalized wag normal modes, whose highest frequency members reach up above 1400 cm^-1. As such, their overtones mix significantly with the CH stretch modes, most notably involving the CH₂ symmetric stretch fundamentals of the central ethyl groups in the all-alkyl chains and the CH stretches adjacent to the N-atoms and antiperiplanar to the nitrogen lone pair.

AMOEBA Force Field Predicts Accurate Hydrogen Bond Counts of Nitriles in SNase by Revealing Water-Protein Interaction in Vibrational Absorption Frequencies

Precisely quantifying the magnitude and direction of electric fields in proteins has long been an outstanding challenge in understanding biological functions. Nitrile vibrational Stark effect probes have been shown to be minimally disruptive to the protein structure and can be better direct reporters of local electrostatic field in the native state of a protein than other measures such as pKa shifts of titratable residues. However, interpretations of the connection between measured vibrational energy and electric field rely on the accurate molecular understanding of interactions of the nitrile group and its environment, particularly from hydrogen bonding. In this work, we compared the extent of hydrogen bonding calculated in two common force fields, the fixed charge force field Amber03 and polarizable force field AMOEBA, at 10 locations of cyanocysteine (CNC) in staphylococcal nuclease (SNase) against the experimental nitrile absorption frequency in terms of full width at half-maximum (FWHM) and frequency temperature line slope (FTLS). We observed that the number of hydrogen bonds correlated well in AMOEBA trajectories with respect to both the FWHM (r = 0.88) and the FTLS (r = -0.85), whereas the correlation of Amber03 trajectories was less reliable because the Amber03 force field predicted more hydrogen bonds in some mutants. Moreover, we demonstrated that contributions from the interactions between CNC and nearby water molecules were significant in AMOEBA trajectories but were not predicted by Amber03. We conclude that although the nitrile absorption peak shape could be qualitatively predicted by the fixed charge Amber03 force field, the detailed electrostatic environment measured by the nitrile probe in terms of the extent of hydrogen bonding could only be accurately observed in the AMOEBA trajectories, where the permanent dipole, quadrupole, and dipole-induced-dipole polarizable interactions were all taken into account. The significance of this finding to the goal of accurately predicting electric fields in complex biomolecular environments is discussed.

Integrative analysis of single-cell and bulk RNA sequencing unveils the senescence landscape in ischemic stroke

Ischemic stroke (IS) is a fatal neurological disease that occurs when the blood flow to the brain is disrupted, leading to brain tissue damage and functional impairment. Cellular senescence, a vital characteristic of aging, is associated with a poor prognosis for IS. This study explores the potential role of cellular senescence in the pathological process following IS by analyzing transcriptome data from multiple datasets (GSE163654, GSE16561, GSE119121, and GSE174574). By using bioinformatics methods, we identified hub-senescence-related genes such as ANGPTL4, CCL3, CCL7, CXCL16, and TNF and verified them using quantitative reverse transcription polymerase chain reaction. Further analysis of single-cell RNA sequencing data suggests that MG4 microglial is highly correlated with cellular senescence in MCAO, and might play a crucial role in the pathological process after IS. Additionally, we identified retinoic acid as a potential drug for improving the prognosis of IS. This comprehensive investigation of cellular senescence in various brain tissues and peripheral blood cell types provides valuable insights into the underlying mechanisms of the pathology of IS and identifies potential therapeutic targets for improving patient outcomes.

Scalable hybrid deep neural networks/polarizable potentials biomolecular simulations including long-range effects

Deep-HP is a scalable extension of the Tinker-HP multi-GPU molecular dynamics (MD) package enabling the use of Pytorch/TensorFlow Deep Neural Network (DNN) models. Deep-HP increases DNNs' MD capabilities by orders of magnitude offering access to ns simulations for 100k-atom biosystems while offering the possibility of coupling DNNs to any classical (FFs) and many-body polarizable (PFFs) force fields. It allows therefore the introduction of the ANI-2X/AMOEBA hybrid polarizable potential designed for ligand binding studies where solvent-solvent and solvent-solute interactions are computed with the AMOEBA PFF while solute-solute ones are computed by the ANI-2X DNN. ANI-2X/AMOEBA explicitly includes AMOEBA's physical long-range interactions via an efficient Particle Mesh Ewald implementation while preserving ANI-2X's solute short-range quantum mechanical accuracy. The DNN/PFF partition can be user-defined allowing for hybrid simulations to include key ingredients of biosimulation such as polarizable solvents, polarizable counter ions, etc.… ANI-2X/AMOEBA is accelerated using a multiple-timestep strategy focusing on the model's contributions to low-frequency modes of nuclear forces. It primarily evaluates AMOEBA forces while including ANI-2X ones only via correction-steps resulting in an order of magnitude acceleration over standard Velocity Verlet integration. Simulating more than 10 μs, we compute charged/uncharged ligand solvation free energies in 4 solvents, and absolute binding free energies of host-guest complexes from SAMPL challenges. ANI-2X/AMOEBA average errors are discussed in terms of statistical uncertainty and appear in the range of chemical accuracy compared to experiment. The availability of the Deep-HP computational platform opens the path towards large-scale hybrid DNN simulations, at force-field cost, in biophysics and drug discovery.

The sustainable antihypertensive and target organ damage protective effect of transcranial focused ultrasound stimulation in spontaneously hypertensive rats

OBJECTIVE

In this study, we aimed to investigate the sustainable antihypertensive effects and protection against target organ damage caused by low-intensity focused ultrasound (LIFU) stimulation and the underlying mechanism in spontaneously hypertensive rats (SHRs) model.

METHODS AND RESULTS

SHRs were treated with ultrasound stimulation of the ventrolateral periaqueductal gray (VlPAG) for 20 min every day for 2 months. Systolic blood pressure (SBP) was compared among normotensive Wistar-Kyoto rats, SHR control group, SHR Sham group, and SHR LIFU stimulation group. Cardiac ultrasound imaging and hematoxylin-eosin and Masson staining of the heart and kidney were performed to assess target organ damage. The c-fos immunofluorescence analysis and plasma levels of angiotensin II, aldosterone, hydrocortisone, and endothelin-1 were measured to investigate the neurohumoral and organ systems involved. We found that SBP was reduced from 172 ± 4.2 mmHg to 141 ± 2.1 mmHg after 1 month of LIFU stimulation, P  < 0.01. The next month of treatment can maintain the rat's blood pressure at 146 ± 4.2 mmHg at the end of the experiment. LIFU stimulation reverses left ventricular hypertrophy and improves heart and kidney function. Furthermore, LIFU stimulation enhanced the neural activity from the VLPAG to the caudal ventrolateral medulla and reduced the plasma levels of ANGII and Aldo.

CONCLUSION

We concluded that LIFU stimulation has a sustainable antihypertensive effect and protects against target organ damage by activating antihypertensive neural pathways from VLPAG to the caudal ventrolateral medulla and further inhibiting the renin-angiotensin system (RAS) activity, thereby supporting a novel and noninvasive alternative therapy to treat hypertension.

Research on quality evaluation of innovation and entrepreneurship education for college students based on random forest algorithm and logistic regression model

The quality evaluation of innovation and entrepreneurship (I&E) in the education sector is achieving worldwide attention as empowering nations with high quality talents is quintessential for economic progress. China, a pioneer in the world market in almost all sectors have transformed its educational policies and incorporated entrepreneurial skills as a part of their education models to further catalyst the country's economic progress. This research focuses on building a novel hybrid Machine Learning (ML) model by integrating two powerful algorithms namely Random Forest (RF) and Logistic Regression (LR) to assess the intensity of the I&E in education from the data acquired from 25 leading Higher Educational Institution's (HEI) in different provinces. The major contributions to the work are, (1) construction of quality index for each topic of interest using individual RF, (2) ranking the indicators based on the quality index to assess the strength and weaknesses, (3) and finally use the LR algorithm study the quality of each indicator. The efficacy of the proposed hybrid model is validated using the benchmark classification metrics to assess its learning and prediction performance in evaluating the quality of I&E education. The result of the research portrays that the universities have now started to integrate entrepreneurship skills as a part of the curriculum, which is evident from the better ranking of the topic curriculum development which is followed by the enrichment of skills. This comprehensive research will help the institutions to identify the potential areas of growth to boost the economic development and improve the skill set necessary for I&E education among college students.

Identification of a novel acylthiourea-based potent broad-spectrum inhibitor for enterovirus 3D polymerase in vitro and in vivo

Enterovirus infections have become a serious public health threat to young children, leading to hand-foot-and-mouth disease and more severe nervous system diseases. Due to the lack of licensed anti enterovirus drugs, we reported herein that a Tenovin-1 analog, acylthiourea-based 4-(tert-butyl)-N-((4-(4-(tert-butyl)benzamido)phenyl)carbamothioyl) benzamide (AcTU), displayed low nanomolar anti-EV-A71 activity with an EC₅₀ of 1.0 nM in RD cells. Moreover, AcTU exhibited nanomolar to picomolar inhibitory activity against a series of enteroviruses including EV-D68, CV-A21, CV-A16 and CV-B1 (EC₅₀ = 0.75-17.15 nM). Mechanistic studies indicated that AcTU inhibited enterovirus proliferation by targeting 3D polymerase. In addition, AcTU displayed moderate pharmacokinetic properties in rats (F = 7.4%, T_1/2 = 3.26 h), and in vivo protection studies demonstrated that AcTU orally administered at 0.6 mg/kg/d was highly protective against lethal EV-A71 challenge in mice, potentially reducing mortality from 100% to 20% as well as alleviating symptoms. These results suggested that AcTU could be a potent clinical candidate for the treatment of enterovirus infections.

Development of AMOEBA Polarizable Force Field for Rare-Earth La3+ Interaction with Bioinspired Ligands

Rare-earth metals (REMs) are crucial for many important industries, such as power generation and storage, in addition to cancer treatment and medical imaging. One promising new REM refinement approach involves mimicking the highly selective and efficient binding of REMs observed in relatively recently discovered proteins. However, realizing any such bioinspired approach requires an understanding of the biological recognition mechanisms. Here, we developed a new classical polarizable force field based on the AMOEBA framework for modeling a lanthanum ion (La3+) interacting with water, acetate, and acetamide, which have been found to coordinate the ion in proteins. The parameters were derived by comparing to high-level ab initio quantum mechanical (QM) calculations that include relativistic effects. The AMOEBA model, with advanced atomic multipoles and electronic polarization, is successful in capturing both the QM distance-dependent La3+-ligand interaction energies and experimental hydration free energy. A new scheme for pairwise polarization damping (POLPAIR) was developed to describe the polarization energy in La3+ interactions with both charged and neutral ligands. Simulations of La3+ in water showed water coordination numbers and ion-water distances consistent with previous experimental and theoretical findings. Water residence time analysis revealed both fast and slow kinetics in water exchange around the ion. This new model will allow investigation of fully solvated lanthanum ion-protein systems using GPU-accelerated dynamics simulations to gain insights on binding selectivity, which may be applied to the design of synthetic analogues.

Computational Investigation of a Series of Small Molecules as Potential Compounds for Lysyl Hydroxylase-2 (LH2) Inhibition

The catalytic function of lysyl hydroxylase-2 (LH2), a member of the Fe(II)/αKG-dependent oxygenase superfamily, is to catalyze the hydroxylation of lysine to hydroxylysine in collagen, resulting in stable hydroxylysine aldehyde-derived collagen cross-links (HLCCs). Reports show that high amounts of LH2 lead to the accumulation of HLCCs, causing fibrosis and specific types of cancer metastasis. Some members of the Fe(II)/αKG-dependent family have also been reported to have intramolecular O2 tunnels, which aid in transporting one of the required cosubstrates into the active site. While LH2 can be a promising target to combat these diseases, efficacious inhibitors are still lacking. We have used computational simulations to investigate a series of 44 small molecules as lead compounds for LH2 inhibition. Tunneling analyses indicate the existence of several intramolecular tunnels. The lengths of the calculated O2-transporting tunnels in holoenzymes are relatively longer than those in the apoenzyme, suggesting that the ligands may affect the enzyme's structure and possibly block (at least partially) the tunnels. The sequence alignment analysis between LH enzymes from different organisms shows that all of the amino acid residues with the highest occurrence rate in the oxygen tunnels are conserved. Our results suggest that the enolate form of diketone compounds establishes stronger interactions with the Fe(II) in the active site. Branching the enolate compounds with functional groups such as phenyl and pyridinyl enhances the interaction with various residues around the active site. Our results provide information about possible leads for further LH2 inhibition design and development.

Molecular Dynamics and Free Energy Calculations of Dicyclohexano-18-crown-6 Diastereoisomers with Sm2+, Eu2+, Dy2+, Yb2+, Cf2+, and Three Halide Salts in Tetrahydrofuran and Acetonitrile Using the AMOEBA Force Field

With the continual development of lanthanides (Ln) in current technological devices, an efficient separation process is needed that can recover greater amounts of these rare elements. Dicyclohexano-18-crown-6 (DCH18C6) is a crown ether that may be a promising candidate for Ln separation, but additional research is required. As such, molecular dynamics (MD) simulations have been performed on four divalent lanthanide halide salts (Sm2+, Eu2+, Dy2+, and Yb2+) and one divalent actinide halide salt (Cf2+) bound to three diastereoisomers of DCH18C6. Dy2+, Yb2+, Cf2+, DCH18C6, and tetrahydrofuran (THF) solvent were parameterized for the AMOEBA polarizable force field for the first time, whereas existing parameters for Sm2+ and Eu2+ were utilized from our previous efforts. A coordination number (CN) of six for Ln2+/An2+-O solvated in THF indicated that the cations interacted almost entirely with the oxygens of the polyether ring. A CN of one for Ln2+/An2+-N solvated in acetonitrile for systems containing iodide suggested that the N atom of acetonitrile was competitive with I- for cation interactions. Fluctuation between five and six CNs for Dy2+ and Yb2+ suggested that although the cations remained in the polyether ring, the size of the ring may not be an ideal fit as these cations possess comparatively smaller ionic radii. Gibbs binding free energies of Sm2+ in all DCH18C6 diastereoisomers solvated in THF were calculated. The binding free energy of the cis-syn-cis diastereoisomer was the most favorable, followed by cis-anti-cis, and then trans-anti-trans. Finally, two major types of conformation were observed for each diastereoisomer that were related to the electrostatic interactions and charge density of the cations.

An Efficient Approach to Large-Scale Ab Initio Conformational Energy Profiles of Small Molecules

Accurate conformational energetics of molecules are of great significance to understand maby chemical properties. They are also fundamental for high-quality parameterization of force fields. Traditionally, accurate conformational profiles are obtained with density functional theory (DFT) methods. However, obtaining a reliable energy profile can be time-consuming when the molecular sizes are relatively large or when there are many molecules of interest. Furthermore, incorporation of data-driven deep learning methods into force field development has great requirements for high-quality geometry and energy data. To this end, we compared several possible alternatives to the traditional DFT methods for conformational scans, including the semi-empirical method GFN2-xTB and the neural network potential ANI-2x. It was found that a sequential protocol of geometry optimization with the semi-empirical method and single-point energy calculation with high-level DFT methods can provide satisfactory conformational energy profiles hundreds of times faster in terms of optimization.

Clean Household Energy Consumption and Residents' Well-Being: Empirical Analysis and Mechanism Test

Improving well-being is a critical problem for worldwide social progress. Research on well-being effects of clean household energy consumption is of great significance for the realization of United Nation's Sustainable Development Goals (SDGs). Due to the multifaceted role of clean household energy in enhancing well-being as a sustainable development goal, it has attracted extensive academic attention and research but still needs to be further refined through new perspectives. This paper uses data from the 2018 China Health and Retirement Longitudinal Study to conduct an empirical analysis of clean household energy consumption and residents' well-being using an ordered probit model, the instrumental variables method, a conditional mixed process (CMP) method, and a mechanism analysis model. The results show that (1) household clean energy consumption contributes to residents' well-being, and the results remained significant after selecting "Do you have an electric Bicycle?" as an instrumental variable to address the endogeneity question. (2) According to heterogeneity research, women's wellbeing is impacted by clean energy consumption in the household more than men's. Only in rural locations can clean household energy consumption significantly boost well-being; urban and suburban areas are unaffected. (3) In the mechanism analysis, health condition and depression played a mediating role on the impact of clean household energy consumption on well-being, and social contact played a moderating role on the impact of clean household energy consumption on well-being. The findings and policy recommendations in this paper are suggestive of how we can improve the well-being of residents in low- and middle-income countries and provide reference values for research in related fields around the world.

AMOEBA Force Field Trajectories Improve Predictions of Accurate pKa Values of the GFP Fluorophore: The Importance of Polarizability and Water Interactions

Precisely quantifying the magnitude, direction, and biological functions of electric fields in proteins has long been an outstanding challenge in the field. The most widely implemented experimental method to measure such electric fields at a particular residue in a protein has been through changes in pKa of titratable residues. While many computational strategies exist to predict these values, it has been difficult to do this accurately or connect predicted results to key structural or mechanistic features of the molecule. Here, we used experimentally determined pKa values of the fluorophore in superfolder green fluorescent protein (GFP) with amino acid mutations made at position Thr 203 to evaluate the pKa prediction ability of molecular dynamics (MD) simulations using a polarizable force field, AMOEBA. Structure ensembles from AMOEBA were used to calculate pKa values of the GFP fluorophore. The calculated pKa values were then compared to trajectories using a conventional fixed charge force field (Amber03 ff). We found that the position of water molecules included in the pKa calculation had opposite effects on the pKa values between the trajectories from AMOEBA and Amber03 force fields. In AMOEBA trajectories, the inclusion of water molecules within 35 Å of the fluorophore decreased the difference between the predicted and experimental values, resulting in calculated pKa values that were within an average of 0.8 pKa unit from the experimental results. On the other hand, in Amber03 trajectories, including water molecules that were more than 5 Å from the fluorophore increased the differences between the calculated and experimental pKa values. The inaccuracy of pKa predictions determined from Amber03 trajectories was caused by a significant stabilization of the deprotonated chromophore's free energy compared to the result in AMOEBA. We rationalize the cutoffs for explicit water molecules when calculating pKa to better predict the electrostatic environment surrounding the fluorophore buried in GFP. We discuss how the results from this work will assist the prospective prediction of pKa values or other electrostatic effects in a wide variety of folded proteins.

Molecular Dynamics Simulations of Protein RNA Complexes by Using an Advanced Electrostatic Model

Protein-RNA interactions are integral to the biological functions of RNA. It is well recognized that molecular dynamics (MD) simulations of protein-RNA complexes are more challenging than those of each component. The difficulty arises from the strong electrostatic interactions and the delicate balance between various types of physical forces at the interface. Previously, MD simulations of protein-RNA complexes have predominantly employed fixed-charge force fields. Although force field modifications have been developed to address problems identified in the simulations, some protein-RNA structures are still hard to reproduce by simulations. Here, we present MD simulations of two representative protein-RNA complexes using the AMOEBA polarizable force field. The van der Waals parameters were refined to reproduce accurate quantum-mechanical data of base-base and base-amino acid interactions. It was found that the refined parameters produced a more stable hydrogen-bond network in the interface. One of the complexes remained stable during the short simulations, whereas it could quickly break down in previous simulations using fixed-charge force fields. There was reversible breaking and formation of hydrogen bonds that are observed in the crystal structure, which may indicate the difference in solution and crystal structures. While further improvement and validation of the force fields are still needed, this work demonstrates that polarizable force fields are promising for the study of protein-RNA complexes.

A novel risk score model based on fourteen chromatin regulators-based genes for predicting overall survival of patients with lower-grade gliomas

Background: Low grade gliomas(LGGs) present vexatious management issues for neurosurgeons. Chromatin regulators (CRs) are emerging as a focus of tumor research due to their pivotal role in tumorigenesis and progression. Hence, the goal of the current work was to unveil the function and value of CRs in patients with LGGs.

Methods: RNA-Sequencing and corresponding clinical data were extracted from The Cancer Genome Atlas (TCGA) and the Chinese Glioma Genome Atlas (CGGA) database. A single-cell RNA-seq dataset was sourced from the Gene Expression Omnibus (GEO) database. Altogether 870 CRs were retrieved from the published articles in top academic journals. The least absolute shrinkage and selection operator (LASSO) algorithm and Cox regression analysis were applied to construct the prognostic risk model. Patients were then assigned into high- and low-risk groups based on the median risk score. The Kaplan–Meier (K-M) survival curve and receiver operating characteristic curve (ROC) were performed to assess the prognostic value. Sequentially, functional enrichment, tumor immune microenvironment, tumor mutation burden, drug prediction, single cell analysis and so on were analyzed to further explore the value of CR-based signature. Finally, the expression of signature genes were validated by immunohistochemistry (IHC) and quantitative real-time PCR (qRT-PCR).

Results: We successfully constructed and validated a 14 CRs-based model for predicting the prognosis of patients with LGGs. Moreover, we also found 14 CRs-based model was an independent prognostic factor. Functional analysis revealed that the differentially expressed genes were mainly enriched in tumor and immune related pathways. Subsequently, our research uncovered that LGGs patients with higher risk scores exhibited a higher TMB and were less likely to be responsive to immunotherapy. Meanwhile, the results of drug analysis offered several potential drug candidates. Furthermore, tSNE plots highlighting the magnitude of expression of the genes of interest in the cells from the scRNA-seq assay. Ultimately, transcription expression of six representative signature genes at the mRNA level was consistent with their protein expression changes.

Conclusion: Our findings provided a reliable biomarker for predicting the prognosis, which is expected to offer new insight into LGGs management and would hopefully become a promising target for future research.

Automation of AMOEBA polarizable force field for small molecules: Poltype 2

A next-generation protocol (Poltype 2) has been developed which automatically generates AMOEBA polarizable force field parameters for small molecules. Both features and computational efficiency have been drastically improved. Notable advances include improved database transferability using SMILES, robust torsion fitting, non-aromatic ring torsion parameterization, coupled torsion-torsion parameterization, Van der Waals parameter refinement using ab initio dimer data and an intelligent fragmentation scheme that produces parameters with dramatically reduced ab initio computational cost. Additional improvements include better local frame assignment for atomic multipoles, automated formal charge assignment, Zwitterion detection, smart memory resource defaults, parallelized fragment job submission, incorporation of Psi4 quantum package, ab initio error handling, ionization state enumeration, hydration free energy prediction and binding free energy prediction. For validation, we have applied Poltype 2 to ~1000 FDA approved drug molecules from DrugBank. The ab initio molecular dipole moments and electrostatic potential values were compared with Poltype 2 derived AMOEBA counterparts. Parameters were further substantiated by calculating hydration free energy (HFE) on 40 small organic molecules and were compared with experimental data, resulting in an RMSE error of 0.59 kcal/mol. The torsion database has expanded to include 3543 fragments derived from FDA approved drugs. Poltype 2 provides a convenient utility for applications including binding free energy prediction for computational drug discovery. Further improvement will focus on automated parameter refinement by experimental liquid properties, expansion of the Van der Waals parameter database and automated parametrization of modified bio-fragments such as amino and nucleic acids.

Computational study on the binding of Mango-II RNA aptamer and fluorogen using the polarizable force field AMOEBA

Fluorescent light-up aptamers (FLAPs) are well-performed biosensors for cellular imaging and the detection of different targets of interest, including RNA, non-nucleic acid molecules, metal ions, and so on. They could be easily designed and emit a strong fluorescence signal once bound to specified fluorogens. Recently, one unique aptamer called Mango-II has been discovered to possess a strong affinity and excellent fluorescent properties with fluorogens TO1-Biotin and TO3-Biotin. To explore the binding mechanisms, computational simulations have been performed to obtain structural and thermodynamic information about FLAPs at atomic resolution. AMOEBA polarizable force field, with the capability of handling the highly charged and flexible RNA system, was utilized for the simulation of Mango-II with TO1-Biotin and TO3-Biotin in this work. The calculated binding free energy using published crystal structures is in excellent agreement with the experimental values. Given the challenges in modeling complex RNA dynamics, our work demonstrates that MD simulation with a polarizable force field is valuable for understanding aptamer-fluorogen binding and potentially designing new aptamers or fluorogens with better performance.

[The working principle and prototype construction of the Chinese vestibular prosthesis]

Objective:To complete the working principle design and prototype construction of the Chinese multichannel vestibular prosthesis （CMVP） with independent intellectual property rights, and verify its working performance, so as to lay the foundation for the clinical promotion and application of CMVP. Methods:On the basis of previous research, the working principle of CMVP was constructed based on the information encoding principle of vestibular nervous system, and the circuit was designed according to the principle. Then, appropriate electronic components and software systems were selected to construct a CMVP prototype according to the design. Finally, the input and output characteristics of the CMVP prototype were verified through the performance test. Results:In the present study, a block diagram of the working principle of the CMVP was successfully designed and drawn, and the working principle was explained in detail according to the block diagram. Further, the circuit diagram of the CMVP was designed and drawn based on the working principle, then the selected electronic components and software systems were combined one by one to complete the construction of a prototype. Finally, the performance test for the prototype was completed, which showed that all stimulus electrodes of the prototype could output biphasic pulse current, and the frequency of biphasic pulse current was modulated by the spatial rotation data input sensed by a motion sensor. Conclusion:The working principle and circuit design of the CMVP are reasonable; the CMVP prototype in China has been successfully constructed; the spatial rotation motion sensing input and the modulated pulse current output are stable and reliable.

Bilateral Dysfunction of Otolith Pathway in Patients With Unilateral Idiopathic BPPV Detected by ACS-VEMPs

Objective

To observe the functional status of the otolith pathway in patients with unilateral idiopathic benign paroxysmal positional vertigo (BPPV) by combining air-conducted sound elicited cervical vestibular-evoked myogenic potential (ACS-cVEMP) and ocular vestibular-evoked myogenic potential (ACS-oVEMP).

Methods

One hundred and eighty patients with BPPV were recruited for conventional cVEMP and oVEMP tests. The abnormal rates of VEMPs were compared between BPPV patients and control participants.

Results

The abnormal rates of cVEMP and oVEMP in BPPV patients were 46.7% (84/180) and 57.2% (103/180) in affected ears, respectively, and 45.0% (81/180) and 56.7% (102/180) in unaffected ears, respectively; both were significantly higher than the abnormal rates of cVEMP and oVEMP in normal control ears. Compared with normal subjects, the cVEMP response rate was lower in affected and unaffected ears in BPPV patients. The abnormal rates of cVEMP and oVEMP were 48.1% (76/158) and 57.6% (91/158) in patients with posterior semicircular canal BPPV, and 36.4% (8/22) and 54.5% (12/22) in lateral semicircular canal BPPV. There was no significant difference in VEMP abnormalities between posterior semicircular canal BPPV and lateral semicircular canal BPPV.

Conclusion

The prevalence of abnormal cVEMPs and oVEMPs in both affected and unaffected ears of patients with BPPV was significantly higher than that observed in the control group. The pathological mechanism of unilateral idiopathic BPPV may be associated with bilateral degeneration of otolith pathways.

Hydration of divalent lanthanides, Sm2+ and Eu2+ : A molecular dynamics study with polarizable AMOEBA force field

The chemistry of divalent lanthanides, Ln2+ , is a growing sub-field of heavy element chemistry owing to new synthetic approaches. However, some theoretical aspects of these unusual cations are currently underdeveloped, especially as they relate to their dynamic properties in solution. In this work, we address the hydration of two of the classical Ln2+ cations, Sm2+ and Eu2+ , using atomic multipole optimized energetic for biomolecular applications (AMOEBA) force fields. These cations have not been parameterized to date with AMOEBA, and few studies are available because of their instability with respect to oxidation in aqueous media. Coordination numbers (CN's) of 8.2 and 8.1 respectively for Sm2+ and Eu2+ , and 8.8 for both Sm3+ and Eu3+ have been obtained and are in good agreement with the few available AIMD and X-ray absorption fine structures studies. The decreased CN of Ln2+ compared with Ln3+ arises from progressive water exchange events that indicates the gradual stabilization of 8-coordinate structures with respect to 9-coordinate geometries. Moreover, the effects of the chloride counter anions on the coordination of Ln2+ cations have been studied at different chloride concentrations in this work. Lastly, water exchange times of Ln2+ cations have been calculated to provide a comprehensive understanding of the behavior of Eu2+ and Sm2+ in aqueous chloride media.

O-235 ACLY entry into the nucleus to bind with EP300 and regulate the genomic DNA and ribosomal DNA transcription during early embryonic development

Study question

ATP citrate lyase (ACLY) is one of the key enzymes to produce acetyl coenzyme A. What is its role in early embryonic development?

Summary answer

ACLY affected by activation of T447, S451 and S455 sites into nucleus and band with EP300 to affect genomic DNA and ribosomal DNA transcription.

What is known already

ACLY is one of the key enzymes to produce acetyl coenzyme A, which participates in fatty acid metabolism and regulates acetylation. Studies have found that ACLY knockout cells can survive, while the cell proliferation capacity will be weakened. ACLY is highly expressed in the cytoplasm of Hela cells colon cancer cells and other somatic cells. When AKT phosphorylates S455 site, ACLY enters the nucleus and participates in DNA damage repair. Its localization and role in embryos have not been studied. The only thing we know is homozygous mice could not be detected at E8.5 days after implantation.

Study design, size, duration

Twenty-four human day 6 developmentally arrested embryos were used to detect ACLY localization. Dozens of ACLY heterozygous female mice were used to test the time points of ACLY homozygous lethality after mating with heterozygous male mice. After zygotes were taken from wild-type C57BL6 mice, microinjection technology was used to explore the effects of knocking down ACLY, changing ACLY location or mutating ACLY. HeLa cells were used to explore at the molecular level.

Participants/materials, setting, methods

In this experiment, we used human abandoned or donated embryos, gene knockout mice, wild-type C57BL6 mice and Hela cell.

The experimental techniques used include in vitro embryo culture, EU staining, plasmid transfection, immunoprecipitation, western blot analysis, in vitro transcription, microinjection, immunofluorescence, RT-PCR and RNA-seq.

Main results and the role of chance

Firstly, our study found its special nuclear localization in early embryos for the first time and proved that its localization was related to embryonic development potential.

Secondly, the specific period of ACLY homozygous death was identified by ACLY heterozygous mice. And by knocking down ACLY or adding inhibitor SB204990, It was found that ACLY played an important role in cell proliferation, blastocyst formation and lineage differentiation of early embryo development.

Thirdly, we proved that ACLY nucleation played an important role in the development of early embryos. and the activated form of T447, S451 and S455 three sites led to ACLY nucleation.

Finally, we found that ACLY affected the level of histone acetylation and UBF acetylation. And then genomic DNA and ribosomal DNA transcription level were decreased. We further explored how the acetyl group produced by ACLY regulated transcription. The interaction between ACLY and EP300 was confirmed by immunofluorescence and COIP experiments. After knocking down EP300 in embryos or inhibiting its function with inhibitor c646, the results are consistent with knocking down ACLY. The above experimental data prove that ACLY affects histone acetylation level and regulates gene opening through the interaction with EP300.

Limitations, reasons for caution

We need to further explore the factors that promote ACLY into the nucleus, so as to add it to in vitro culture to enhance the embryonic development potential.

Wider implications of the findings

The study provides a new mechanism for insulin or lysophosphatidic acid to promote embryo development, and helps to improve the quality of embryo development in the process of embryo culture in vitro.

Trial registration number

not applicable

Development of the Quantum-Inspired SIBFA Many-Body Polarizable Force Field: Enabling Condensed-Phase Molecular Dynamics Simulations

We present the extension of the Sum of Interactions Between Fragments Ab initio Computed (SIBFA) many-body polarizable force field to condensed-phase molecular dynamics (MD) simulations. The quantum-inspired SIBFA procedure is grounded on simplified integrals obtained from localized molecular orbital theory and achieves full separability of its intermolecular potential. It embodies long-range multipolar electrostatics (up to quadrupole) coupled to a short-range penetration correction (up to charge-quadrupole), exchange repulsion, many-body polarization, many-body charge transfer/delocalization, exchange dispersion, and dispersion (up to C10). This enables the reproduction of all energy contributions of ab initio symmetry-adapted perturbation theory (SAPT(DFT)) gas-phase reference computations. The SIBFA approach has been integrated within the Tinker-HP massively parallel MD package. To do so, all SIBFA energy gradients have been derived and the approach has been extended to enable periodic boundary conditions simulations using smooth particle mesh Ewald. This novel implementation also notably includes a computationally tractable simplification of the many-body charge transfer/delocalization contribution. As a proof of concept, we perform a first computational experiment defining a water model fitted on a limited set of SAPT(DFT) data. SIBFA is shown to enable a satisfactory reproduction of both gas-phase energetic contributions and condensed-phase properties highlighting the importance of its physically motivated functional form.

High Order Ab Initio Valence Force Field with Chemical Pattern Based Parameter Assignment

Bonded (or valence) interactions, which directly determine the local structures of the molecules, are fundamental parts of molecular mechanics force fields (FFs). Most popular classical FFs adopt the simple harmonic models for bond stretching and angle bending and ignore cross-coupling effects among the valence terms. This may lead to less accurate vibrational properties and configurations in molecular dynamics (MD) simulations. AMOEBA models utilize an MM3(MM4)-style bonded interaction model, in which the vibrational anharmonicity, the coupling effects among different energy terms, and the out-of-plane bending for sp2-hybridized atoms are considered. In this work, we report the development of bonded interaction parameters for a wide range of chemistry based on quantum mechanics (QM). About 270 atomic types defined by SMARTS strings were used to model the valence interactions. Our results indicate that the resulting valence parameters produce accurate vibrational frequencies (RMSD from QM is less than ~36.6 cm^-1) over a large set of molecules with diverse functional groups (445 molecules). By contrast, the harmonic models usually give an RMS error greater than 60 cm^-1. Meanwhile, this model accurately reflects the potential energy surface of the out-of-plane bending. Our model can generally be applied to the AMOEBA family and any MM3(MM4)-based molecular mechanics FFs.

Inhibition of FAM83D displays antitumor effects in glioblastoma via down-regulation of the AKT/Wnt/β-catenin pathway

Up-regulation of family with sequence similarity 83 member D (FAM83D) has been acknowledged as a vital contributor for the carcinogenesis of numerous cancers. The relevance of FAM83D in glioblastoma (GBM), however, is not well understood. This current work aimed to determine the possible roles and mechanisms of FAM83D in GBM. By analyzing The Cancer Genome Atlas (TCGA) data, we found dramatic increases in FAM83D expression in GBM tissue. We also observed elevated levels of FAM83D in the clinical specimens of GBM. In vitro data showed that silencing FAM83D resulted in remarkable antitumor effects via inhibiting the proliferation, invasion and epithelial-mesenchymal transition of GBM cells. Moreover, the knockdown of FAM83D improved sensitivity to the chemotherapy drug temozolomide. In-depth mechanism research revealed that the silencing of FAM83D strikingly decreased the phosphorylation levels of AKT and glycogen synthase kinase-3β, and prohibited activation of the Wnt/β-catenin pathway. The suppression of AKT abolished FAM83D-mediated activation of the Wnt/β-catenin pathway. The re-expression of β-catenin reversed FAM83D-silencing-induced antitumor effects in GBM cells. In addition, GBM cells with FAM83D silencing exhibited reduced tumorigenic potential in vivo. Overall, the data from this work show that the inhibition of FAM83D displays antitumor effects in GBM via down-regulation of the AKT/Wnt/β-catenin pathway and propose FAM83D as a new therapeutic target for GBM.

B81 Bone Vibrator-Induced Vestibular-Evoked Myogenic Potentials: Normal Values and the Effect of Age

Objective

To define the normal values and examine the influence of aging on B81 bone vibrator-induced cervical vestibular-evoked myogenic potentials (B81-cVEMPs) and ocular vestibular-evoked myogenic potentials (B81-oVEMPs).

Methods

Seventy healthy subjects, divided into seven groups according to their ages, were enrolled in this study. The 4–9-, 10–19-, 20–29-, 30–39-, 40–49-, 50–59-, and 60–70-year-old participants were divided into groups I–VII, respectively. B81-cVEMP and B81-oVEMP were recorded in each group.

Results

The B81-cVEMP response rates for groups I–VII were 100, 100, 100, 100, 95, 95, and 75%, respectively, with significant differences only between groups I–VI and group VII (p = 0.047, p &lt; 0.05). The B81-oVEMP response rates for groups I–VII were 100, 100, 100, 100, 70, 65, and 40%, respectively, with significant differences only between groups I–IV and groups V–VII (p = 0.020, p = 0.008, p = 0.000; p &lt; 0.05). The threshold, P13, and N23 latencies of B81-cVEMP positively correlated with age (r = 0.756, p = 0.000; r = 0.357, p = 0.003; r = 0.316, p = 0.009; p &lt; 0.05). The raw amplitudes and corrected amplitudes negatively correlated with age (r = −0.641, p = 0.000; r = −0.609, p = 0.000, p &lt; 0.05). For B81-oVEMP, the corrected amplitudes negatively correlated with age (r = −0.638, p = 0.000, p&lt;0.05), but the threshold and N10 latency positively correlated with age (r = 0.768, p = 0.000; r = 0.334, p = 0.009, p &lt; 0.05). Moreover, the interaural asymmetry ratio did not significantly correlate with age for B81-cVEMP and B81-oVEMP.

Conclusion

As age increased, the B81-cVEMP response rate decreased, the thresholds increased, P13 and N23 latencies were prolonged, and the raw amplitude and corrected amplitude decreased. The B81-oVEMP response rate and corrected amplitude decreased, the thresholds increased, and N10 latency was prolonged with age. These changes are probably due to the occurrence of morphological and functional changes in the vestibular system with aging. Therefore, we suggest establishing different reference values according to different age groups when evaluating the VEMP results in patients with vestibular diseases.

Trapping Ca+ inside a molecular cavity: computational study of the potential energy surfaces for Ca+-[n]cycloparaphenylene, n = 5-12

Ion trap quantum computing utilizes electronic states of atomic ions such as Ca⁺ to encode information on to a qubit. To explore the fundamental properties of Ca⁺ inside molecular cavities, we describe here a computational study of Ca⁺ bound inside neutral [n]-cycloparaphenylenes (n = 5-12), often referred to as "nanohoops". This ab initio study characterizes optimized structures, harmonic vibrational frequencies, potential energy surfaces, and ion molecular orbital distortion as functions of increasing nanohoop size. The results of this work provide a first step in guiding experimental studies of the spectroscopy of these ion-molecular cavity complexes.

[Clinical application of three-dimensional printed preformed titanium mesh combined with free latissimus dorsi muscle flap in the treatment of squamous cell carcinoma with skull defect in the vertex]

Objective: To explore the clinical effects of three-dimensional printed preformed titanium mesh combined with latissimus dorsi muscle flap free transplantation in the treatment of wounds with skull defect after radical surgery of squamous cell carcinoma in the vertex. Methods: A retrospective observational study was conducted. From January 2010 to December 2019, 5 patients with squamous cell carcinoma in the vertex accompanied with skull invasion who met the inclusion criteria were admitted to the Department of Burns and Plastic Surgery of the Second Affiliated Hospital of Air Force Medical University, including four males and one female, aged 50 to 65 years. The original lesion areas ranged from 5 cm×4 cm to 15 cm×8 cm. The titanium mesh was prefabricated via three-dimensional technic based on the result the scope of skull resection predicted with computerized tomography three-dimensional reconstruction before surgery. During the first stage, the soft tissue defect area of scalp (8 cm×7 cm to 18 cm×11 cm) after tumor enlargement resection was repaired with the preformed titanium mesh, and the titanium mesh was covered with latissimus dorsi muscle flap, with area of 10 cm×9 cm to 20 cm×13 cm. The thoracodorsal artery/vein was anastomosed with the superficial temporal artery/vein on one side. The muscle ends in the donor site were sutured together or performed with transfixion, and then the skin on the back were covered back to the donor site. On the 10th day after the first-stage surgery, the second-stage surgery was performed. The thin intermediate thickness skin graft was taken from the anterolateral thigh to cover the latissimus dorsi muscle flap. The duration and intraoperative blood loss of first-stage surgery were recorded. The postoperative muscle flap survival after the first-stage surgery and skin graft survival after the second-stage surgery was observed. The occurrence of complications, head appearance, and recurrence of tumor were followed up. Results: The average first-stage surgery duration of patients was 12.1 h, and the intraoperative blood loss was not more than 1 200 mL. The muscle flaps in the first-stage surgery and the skin grafts in the second-stage surgery all survived well. During the follow-up of 6-18 months, no complications such as exposure of titanium mesh or infection occurred, with good shape in the recipient sites in the vertex, and no recurrence of tumor. Conclusions: Three-dimensional printed preformed titanium mesh combined with latissimus dorsi muscle flap free transplantation and intermediate thickness skin graft cover is an effective and reliable method for repairing the wound with skull defect after extended resection of squamous cell carcinoma in the vertex. This method can cover the wound effectively as well as promote both recipient and donor sites to obtain good function and appearance.

[Research advances on application of botulinum toxin type A in scar prevention and treatment]

The wound healing time, tension of wound edge, proliferation of fibroblast, and extracellular matrix deposition are the important factors of scar formation, and botulinum toxin type A can regulate the above. Prevention and treatment of scar with botulinum toxin type A is one of the hot topics of clinical research in recent years. This paper briefly reviews researches by scholars at home and abroad on the mechanism, clinical application, complications, and adverse effects of botulinum toxin type A in scar prevention and treatment.

Long non-coding RNA AGAP2-AS1 promotes cell proliferation and invasion through regulating miR-193a-3p/LOXL4 axis in laryngeal squamous cell carcinoma

Laryngeal squamous cell carcinoma (LSCC) is an aggressive malignancy with highly mortality rate. Long non-coding RNA (lncRNA) AGAP2-AS1 is an identified oncogene in several types of cancers. However, the role of AGAP2-AS1 in LSCC remains unclear. The expression levels of AGAP2-AS1 in LSCC tissues and cell lines were measured using qRT-PCR. AGAP2-AS1 was knocked down in LSCC cells through transfection with siRNA-AGAP2-AS1. Cell proliferation and invasion were detected using MTT and transwell assays. Dual-luciferase reporter gene assay was performed to confirm the interaction with AGAP2-AS1 and downstream genes. Our results showed that AGAP2-AS1 expression was remarkably increased in human LSCC tissues and cell lines. Knockdown of AGAP2-AS1 significantly inhibited the proliferation and invasion of LSCC cells. In addition, AGAP2-AS1 sponged miR-193a-3p and regulated its expression in LSCC cells. Inhibition of miR-193a-3p reversed the effects of AGAP2-AS1 knockdown on LSCC cells. Furthermore, Lysyl oxidase-like 4 (LOXL4) was a target gene of miR-193a-3p and the role of miR-193a-3p was mediated by LOXL4. In conclusion, these findings suggest that knockdown of AGAP2-AS1 inhibited the proliferation and invasion of LSCC cells through regulating the miR-193a-3p/LOXL4 axis.

Computationally driven discovery of SARS-CoV-2 Mpro inhibitors: from design to experimental validation

We report a fast-track computationally driven discovery of new SARS-CoV-2 main protease (M^pro) inhibitors whose potency ranges from mM for the initial non-covalent ligands to sub-μM for the final covalent compound (IC₅₀ = 830 ± 50 nM). The project extensively relied on high-resolution all-atom molecular dynamics simulations and absolute binding free energy calculations performed using the polarizable AMOEBA force field. The study is complemented by extensive adaptive sampling simulations that are used to rationalize the different ligand binding poses through the explicit reconstruction of the ligand–protein conformation space. Machine learning predictions are also performed to predict selected compound properties. While simulations extensively use high performance computing to strongly reduce the time-to-solution, they were systematically coupled to nuclear magnetic resonance experiments to drive synthesis and for in vitro characterization of compounds. Such a study highlights the power of in silico strategies that rely on structure-based approaches for drug design and allows the protein conformational multiplicity problem to be addressed. The proposed fluorinated tetrahydroquinolines open routes for further optimization of M^pro inhibitors towards low nM affinities.

The dominant binding mode of the QUB-00006-Int-07 main protease inhibitor during absolute binding free energy simulations.

Recent progress in general force fields of small molecules

Recent advances in computational hardware and free energy algorithms enable a broader application of molecular simulation of binding interactions between receptors and small-molecule ligands. The underlying molecular mechanics force fields (FFs) for small molecules have also achieved advancements in accuracy, user-friendliness, and speed during the past several years (2018-2020). Besides the expansion of chemical space coverage of ligand-like molecules among major popular classical additive FFs and polarizable FFs, new charge models have been proposed for better accuracy and transferability, new chemical perception of avoiding predefined atom types have been applied, and new automated parameterization toolkits, including machine learning approaches, have been developed for users' convenience.

Atomic Polarizabilities for Interactive Dipole Induction Models

Thole-style mutual induction models for molecular polarization have been adopted by several popular polarizable force fields (FFs) for their simplicity and transferability. The atomic polarizability parameters of these models are typically derived by fitting to ab initio or/and experimental molecular polarizabilities. In this work, we improve upon Thole polarizability parameters by employing both high-level quantum mechanics molecular polarizabilities and electrostatic potential (ESP) responses on three-dimensional grids. Our results indicate that the two approaches to derive atomic polarizability parameters are both effective, while the ESP approaches can also capture the polarization for the atoms with lone pair electrons. The resulting polarizability parameters have been validated on a set of over 7200 molecules covering the most common elements found in organic molecules (C, H, O, N, P, S, F, Cl, Br, and I). These parameters have also been tested on the experimentally measured molecular polarizabilities of 422 molecules. The final set of parameters derived in this work show notable improvement over the current AMOEBA set. The result is a highly transferable, expanded set of atomic polarizabilities defined by the local chemical environment in the form of SMARTS patterns. These parameters can be used directly in molecular mechanics polarizable potential energy functions such as AMOEBA, AMOEBA+, and other Thole-style models.

Non-invasive ultrasound modulation of solitary tract nucleus exerts a sustainable antihypertensive effect in spontaneously hypertensive rats

GOAL

We applied the method of non-invasive ultrasound (US) neuromodulation to regulate blood pressure (BP) by stimulating the solitary tract nucleus (NTS) of spontaneously hypertensive rats (SHRs).

METHODS

The rats were exposed to US stimulation for 20 mins every day for two months. Morphology and function of the hypertensive target organs (heart and kidney) were then examined by echocardiography and immunohistochemical staining. C-fos immunofluorescence assays were used to evaluate neuronal activity in the US stimulated areas and to explore related neural pathways. Moreover, the effects of US stimulation on biochemical indicators (ANGII, Aldo, ANF, ET-1, and Cor) in SHRs were detected. In addition, HE, TUNEL, and Nissl staining were performed to evaluate the safety of long-term transcranial US stimulation.

RESULTS

After two months of US stimulation, systolic blood pressure (SBP) decreased from 170±1.1 mmHg to 158±5.3 mmHg, p<0.01. What's more, US stimulation effectively inhibited the pathological process of target organs from both morphological and functional levels. With US stimulation, neuronal activities were also significantly enhanced in the NTS, ventrolateral periaqueductal gray (vlPAG), and the caudal ventrolateral medulla (CVLM) region. And US stimulation did not cause brain tissue damage. Meanwhile, the plasma levels of ANF, ATCH, and Cor content were inhibited.

CONCLUSION

US stimulation of the NTS could significantly lower BP in SHRs.

SIGNIFICANCE

Non-invasive transcranial US stimulation acting on the NTS might be a potential therapeutic intervention due to its efficacy and safety.

[Application effects of free transplantation of autologous omentum in treating maxillofacial infection complicated with complex sinus tract]

Objective: To investigate the application effects of free transplantation of autologous omentum in treating maxillofacial infection complicated with complex sinus tract. Methods: The retrospective observational study method was used. Four patients with maxillofacial infection complicated with complex sinus tract were admitted to Department of Burns and Plastic Surgery of the Second Affiliated Hospital of Air Force Medical University from July 2017 to the December 2019, including 3 males and 1 female aged 36-60 years. Preoperative facial computed tomography (CT) was performed on patients for calculating the volume of sinus tract. During the operation, the sinus tract was thoroughly debrided, and the volume of the tissue defect was about 5 cm×3 cm×2 cm-10 cm×5 cm×3 cm after debridement. The tissue defect area was filled with omentum of 100-300 mL which was cut under laparoscopy. The artery and vein on the right side of the omentum were reserved as the vascular pedicle of the donor area, which were anastomosed with the facial artery and external jugular vein of the recipient area. The survival of omentum, and the occurrences of reinfection and complication were observed after operation, respectively. On the 10th day and in 1 month after the operation, the blood supply of omentum was examined by colored Doppler ultrasound and CT angiography, and the filling of tissue defect area was examined by head and face CT. During follow-up after the operation, the recoveries of face appearance and function and scar hyperplasia in the donor area. Results: The transplanted omentums in 4 patients survived after the operation with no reinfection and complication. On the 10th day and in 1 month after the operation, the transplanted omentums had good blood supply, and the filled area with omentum was in good shape, without formation of dead cavity. During follow-up of 6-10 months after surgery, the appearance and function of face recovered well, and there was no obvious scar hyperplasia in the donor area. Conclusions: After free transplantation of omentum in treating maxillofacial infection with complex sinus tract, the patients have good facial appearance and function, and the application of laparoscopy results in little damage to the patients and quick postoperative recovery.

E2EDNA: Simulation Protocol for DNA Aptamers with Ligands

We present E2EDNA, a simulation protocol and accompanying code for the molecular biophysics and materials science communities. This protocol is both easy to use and sufficiently efficient to simulate single-stranded (ss)DNA and small analyte systems that are important to cellular processes and nanotechnologies such as DNA aptamer-based sensors. Existing computational tools used for aptamer design focus on cost-effective secondary structure prediction and motif analysis in the large data sets produced by SELEX experiments. As a rule, they do not offer flexibility with respect to the choice of the theoretical engine or direct access to the simulation platform. Practical aptamer optimization often requires higher accuracy predictions for only a small subset of sequences suggested, e.g., by SELEX experiments, but in the absence of a streamlined procedure, this task is extremely time and expertise intensive. We address this gap by introducing E2EDNA, a computational framework that accepts a DNA sequence in the FASTA format and the structures of the desired ligands and performs approximate folding followed by a refining step, analyte complexation, and molecular dynamics sampling at the desired level of accuracy. As a case study, we simulate a DNA-UTP (uridine triphosphate) complex in water using the state-of-the-art AMOEBA polarizable force field. The code is available at https://github.com/InfluenceFunctional/E2EDNA.

Imaging sub-diffuse optical properties of cancerous and normal skin tissue using machine learning-aided spatial frequency domain imaging

SIGNIFICANCE

Sub-diffuse optical properties may serve as useful cancer biomarkers, and wide-field heatmaps of these properties could aid physicians in identifying cancerous tissue. Sub-diffuse spatial frequency domain imaging (sd-SFDI) can reveal such wide-field maps, but the current time cost of experimentally validated methods for rendering these heatmaps precludes this technology from potential real-time applications.

AIM

Our study renders heatmaps of sub-diffuse optical properties from experimental sd-SFDI images in real time and reports these properties for cancerous and normal skin tissue subtypes.

APPROACH

A phase function sampling method was used to simulate sd-SFDI spectra over a wide range of optical properties. A machine learning model trained on these simulations and tested on tissue phantoms was used to render sub-diffuse optical property heatmaps from sd-SFDI images of cancerous and normal skin tissue.

RESULTS

The model accurately rendered heatmaps from experimental sd-SFDI images in real time. In addition, heatmaps of a small number of tissue samples are presented to inform hypotheses on sub-diffuse optical property differences across skin tissue subtypes.

CONCLUSION

These results bring the overall process of sd-SFDI a fundamental step closer to real-time speeds and set a foundation for future real-time medical applications of sd-SFDI such as image guided surgery.

Molecular dynamics free energy simulations of ATP:Mg2+ and ADP:Mg2+ using the polarizable force field AMOEBA

ATPases and GTPases are two important classes of protein that play critical roles in energy transduction, cellular signaling, gene regulation and catalysis. These proteins use cofactors such as nucleoside di and tri-phosphates (NTP, NDP) and can detect the difference between NDP and NTP which then induce different protein conformations. Mechanisms that drive proteins into the NTP or NDP conformation may depend on factors such as ligand structure and how Mg²⁺ coordinates with the ligand, amino acids in the pocket and water molecules. Here, we have used the advanced electrostatic and polarizable force field AMOEBA and molecular dynamics free energy simulations (MDFE) to examine the various binding mechanisms of ATP:Mg²⁺ and ADP:Mg²⁺.We compared the ATP:Mg²⁺ binding with previous studies using non-polarizable force fields and experimental data on the binding affinity. It was found that the total free energy of binding for ATP:Mg²⁺ (-7.00 ± 2.13 kcal/mol) is in good agreement with experimental values (-8.6 ± .2 kcal/mol)¹. In addition, parameters for relevant protonation states of ATP, ADP, GTP and GDP have been derived. These parameters will allow for researchers to investigate biochemical phenomena involving NTP's and NDP's with greater accuracy than previous studies involving non-polarizable force fields.

Interfacial Water Many-Body Effects Drive Structural Dynamics and Allosteric Interactions in SARS-CoV-2 Main Protease Dimerization Interface

Following our previous work ( Chem. Sci. 2021, 12, 4889-4907), we study the structural dynamics of the SARS-CoV-2 Main Protease dimerization interface (apo dimer) by means of microsecond adaptive sampling molecular dynamics simulations (50 μs) using the AMOEBA polarizable force field (PFF). This interface is structured by a complex H-bond network that is stable only at physiological pH. Structural correlations analysis between its residues and the catalytic site confirms the presence of a buried allosteric site. However, noticeable differences in allosteric connectivity are observed between PFFs and non-PFFs. Interfacial polarizable water molecules are shown to appear at the heart of this discrepancy because they are connected to the global interface H-bond network and able to adapt their dipole moment (and dynamics) to their diverse local physicochemical microenvironments. The water-interface many-body interactions appear to drive the interface volume fluctuations and to therefore mediate the allosteric interactions with the catalytic cavity.

Thermodynamics of ion binding and occupancy in potassium channels

Potassium channels modulate various cellular functions through efficient and selective conduction of K⁺ ions. The mechanism of ion conduction in potassium channels has recently emerged as a topic of debate. Crystal structures of potassium channels show four K⁺ ions bound to adjacent binding sites in the selectivity filter, while chemical intuition and molecular modeling suggest that the direct ion contacts are unstable. Molecular dynamics (MD) simulations have been instrumental in the study of conduction and gating mechanisms of ion channels. Based on MD simulations, two hypotheses have been proposed, in which the four-ion configuration is an artifact due to either averaged structures or low temperature in crystallographic experiments. The two hypotheses have been supported or challenged by different experiments. Here, MD simulations with polarizable force fields validated by ab initio calculations were used to investigate the ion binding thermodynamics. Contrary to previous beliefs, the four-ion configuration was predicted to be thermodynamically stable after accounting for the complex electrostatic interactions and dielectric screening. Polarization plays a critical role in the thermodynamic stabilities. As a result, the ion conduction likely operates through a simple single-vacancy and water-free mechanism. The simulations explained crystal structures, ion binding experiments and recent controversial mutagenesis experiments. This work provides a clear view of the mechanism underlying the efficient ion conduction and demonstrates the importance of polarization in ion channel simulations.

Polarization shapes the energy landscape of ion conduction in potassium channels.