Compact sum-of-products form of the molecular electronic Hamiltonian based on canonical polyadic decomposition.

We propose an approach to represent the second-quantized electronic Hamiltonian in a compact sum-of-products (SOP) form. The approach is based on the canonical polyadic decomposition of the original Hamiltonian projected onto the sub-Fock spaces formed by groups of spin-orbitals. The algorithm for obtaining the canonical polyadic form starts from an exact sum-of-products, which is then optimally compactified using an alternating least squares procedure. We discuss the relation of this specific SOP with related forms, namely the Tucker format and the matrix product operator often used in conjunction with matrix product states. We benchmark the method on the electronic dynamics of an excited water molecule, trans-polyenes, and the charge migration in glycine upon inner-valence ionization. The quantum dynamics are performed with the multilayer multiconfiguration time-dependent Hartree method in second quantization representation. Other methods based on tree-tensor Ansätze may profit from this general approach.

Quantum and classical molecular dynamics for H atom scattering from graphene.

This work presents systematic comparisons between classical molecular dynamics (cMD) and quantum dynamics (QD) simulations of 15-dimensional and 75-dimensional models in their description of H atom scattering from graphene. We use an experimentally validated full-dimensional neural network potential energy surface of a hydrogen atom interacting with a large cell of graphene containing 24 carbon atoms. For quantum dynamics simulations, we apply Monte Carlo canonical polyadic decomposition to transform the original potential energy surface (PES) into a sum of products form and use the multi-layer multi-configuration time-dependent Hartree method to simulate the quantum scattering of a hydrogen or deuterium atom with an initial kinetic energy of 1.96 or 0.96 eV and an incident angle of 0°, i.e., perpendicular to the graphene surface. The cMD and QD initial conditions have been carefully chosen in order to be as close as possible. Our results show little differences between cMD and QD simulations when the incident energy of the H atom is equal to 1.96 eV. However, a large difference in sticking probability is observed when the incident energy of the H atom is equal to 0.96 eV, indicating the predominance of quantum effects. To the best of our knowledge, our work provides the first benchmark of quantum against classical simulations for a system of this size with a realistic PES. Additionally, new projectors are implemented in the Heidelberg multi-configuration time-dependent Hartree package for the calculation of the atom scattering energy transfer distribution as a function of outgoing angles.

Inhibition of nonsense-mediated mRNA decay reduces the tumorigenicity of human fibrosarcoma cells.

Nonsense-mediated mRNA decay (NMD) is a eukaryotic RNA decay pathway with roles in cellular stress responses, differentiation, and viral defense. It functions in both quality control and post-transcriptional regulation of gene expression. NMD has also emerged as a modulator of cancer progression, although available evidence supports both a tumor suppressor and a pro-tumorigenic role, depending on the model. To further investigate the role of NMD in cancer, we knocked out the NMD factor SMG7 in the HT1080 human fibrosarcoma cell line, resulting in suppression of NMD function. We then compared the oncogenic properties of the parental cell line, the SMG7-knockout, and a rescue cell line in which we re-introduced both isoforms of SMG7. We also tested the effect of a drug inhibiting the NMD factor SMG1 to distinguish NMD-dependent effects from putative NMD-independent functions of SMG7. Using cell-based assays and a mouse xenograft tumor model, we showed that suppression of NMD function severely compromises the oncogenic phenotype. Molecular pathway analysis revealed that NMD suppression strongly reduces matrix metalloprotease 9 (MMP9) expression and that MMP9 re-expression partially rescues the oncogenic phenotype. Since MMP9 promotes cancer cell migration and invasion, metastasis and angiogenesis, its downregulation may contribute to the reduced tumorigenicity of NMD-suppressed cells. Collectively, our results highlight the potential value of NMD inhibition as a therapeutic approach.

Epitope-engineered human hematopoietic stem cells are shielded from CD123-targeted immunotherapy.

Targeted eradication of transformed or otherwise dysregulated cells using monoclonal antibodies (mAb), antibody-drug conjugates (ADC), T cell engagers (TCE), or chimeric antigen receptor (CAR) cells is very effective for hematologic diseases. Unlike the breakthrough progress achieved for B cell malignancies, there is a pressing need to find suitable antigens for myeloid malignancies. CD123, the interleukin-3 (IL-3) receptor alpha-chain, is highly expressed in various hematological malignancies, including acute myeloid leukemia (AML). However, shared CD123 expression on healthy hematopoietic stem and progenitor cells (HSPCs) bears the risk for myelotoxicity. We demonstrate that epitope-engineered HSPCs were shielded from CD123-targeted immunotherapy but remained functional, while CD123-deficient HSPCs displayed a competitive disadvantage. Transplantation of genome-edited HSPCs could enable tumor-selective targeted immunotherapy while rebuilding a fully functional hematopoietic system. We envision that this approach is broadly applicable to other targets and cells, could render hitherto undruggable targets accessible to immunotherapy, and will allow continued posttransplant therapy, for instance, to treat minimal residual disease (MRD).

Sustained agreement rates in the longitudinal assessment of lupus patients to a Broad Consent for personal data and specimen usage in the RHINEVIT biobank.

Background: Biobanks are essential structures for scientific research. The RHINEVIT biobank is used to recruit biomaterials from rheumatology patients in outpatient care and to conduct clinical research studies (e.g., cohort studies) and basic research. RHINEVIT established Broad Consents (BC) to allow extensive and relevant usage of data and biospecimens without the need for specific project restrictions. For quality assurance, we compared the consent rate of individual items of the BC versions in patients with systemic lupus erythematosus (SLE) in the longitudinal study. Methods: BCs were used for biomaterial donation. Informed consent data from RHINEVIT were analyzed. Due to the content restructuring of the BC items due to changes from the templates of the working group of the Medical Ethics Commissions in the Federal Republic of Germany and GDPR requirements, content mapping of the items was performed for the analysis. Results: From September 2015 to March 2022, 291 SLE outpatients donated biomaterials. In 119 patients, the BC was renewed at least once in a subsequent biomaterial donation. Three biomaterial donations were obtained from 21 patients and four from six patients using the respective BC. However, one consent was later revoked. Consent to the BC topics showed consistently high rates of agreement (range 97.5%-100%), with only some patients disagreeing with individual topics. This remained stable over time (median 526 days [Q1 400, Q3 844]). None of the patients disagreed with a certain topic in two consecutive visits. Conclusion: Modifications to the BC did not result in any relevant changes in the approval rates for SLE patients. RHINEVIT's BC is successfully used for the quality-assured handling of excellently annotated biomaterial. The long-term use of these highly valuable biospecimens for unrestricted research, also in an international context, remains assured.

Simultaneous bilateral TKA in the context of fast track surgery - Do patients meet discharge criteria as anticipated?

BACKGROUND: The efficacy and safety of enhanced recovery after surgery (ERAS) protocols for patients undergoing total knee arthroplasty (TKA) have been generally proven. Previous studies investigating patients undergoing simultaneous bilateral TKA (SBTKA) focused on complications, mortality, and pain and did not examine patients' functional limitations. Therefore, the aim of this study was to investigate to what extent patients undergoing SBTKA are able to meet functional discharge criteria originally designed for their counterparts undergoing unilateral TKA (UTKA) in an ERAS setting. MATERIALS AND METHODS: All patients who received primary SBTKA between June 2015 and December 2018 were included in this retrospective analysis. For comparison, UTKA patients were matched 1:1 to SBTKA patients using Propensity Score Matching based on age, gender, and BMI. The times to achieving the rehabilitation checkpoints of walking 150 m, walking a flight of stairs, and 90° knee flexion were evaluated. RESULTS: 63 (SBTKA group) and 64 (UTKA group) patients were included. Due to the Propensity-Score-Matching there were no differences regarding age, gender, and BMI. The mean length of stay (LOS) was 9.1 days in the SBTKA and 7.6 days in the UTKA group (p = 0.003). On average, it took SBTKA patients 5.4 days to achieve an uninterrupted walking distance of at least 150 m, while it took UTKA patients 4.1 days (p < 0.001). Mean time to walking a flight of stairs was 6.3 days for SBTKA patients and 4.7 days for UTKA patients (p < 0.001). 90° flexion was achieved after 4.1 days by SBTKA patients and 3.5 days by UTKA patients (p = 0.241). CONCLUSION: The vast majority of SBTKA patients were able to achieve functional discharge criteria within their inpatient stay when allowed about 30% extra time. Therefore, functional discharge criteria in ERAS protocols designed for UTKA can be considered appropriate for SBTKA patients. LEVEL OF EVIDENCE: Therapeutic Level III.

Direct determination of glyphosate, aminomethylphosphonic acid and glufosinate in food samples with ion chromatography coupled to electrospray ionization tandem mass spectrometry.

A fast and reliable method for the direct determination of the herbicide glyphosate, its major degradation product aminomethylphosphonic acid (AMPA) and glufosinate is presented for a variety of food matrices. The Quick Polar Pesticides in food of Plant Origin method (QuPPe-PO-Method) was used for extraction without further preconcentration or clean-up steps involving e.g. solid phase extraction (SPE). The method makes use of a commercially available high performance liquid chromatograph coupled to a tandem mass spectrometer with electrospray ionization (LC-ESI-MS/MS) - as present in many laboratories - equipped with an ion chromatography (IC)-column using an MS-compatible eluent made of 0.8% formic acid in water. Due to the absence of time-consuming clean-up procedures, strong matrix effects (ME) of up to 91% for AMPA in grapefruit can be observed, when comparing its sensitivity to that obtained for solvent-based standards. The limits of detection (LODs) were determined for the sample matrices apple, mushrooms, grapefruit, linseed, red lentils and wheat and they were found to be in the range of 0.09 to 0.8, 0.04 to 1 and 0.2 to 2 µg/kg for glyphosate, AMPA and glufosinate, respectively. For the same matrices the validation was carried out according to SANTE guidelines for different commodity groups by spiking them up prior to extraction to concentrations ranging from 10 to 400 µg/kg for matrices with high water content and from 10 to 800 µg/kg for matrices with low water content. When using solvent-based calibration under the use of isotopically labelled internal standards (ILIS) the recoveries were found to range from 84% to 120% and the relative standard deviations (RSD) range between 1% and 19% for glyphosate, AMPA and glufosinate at all fortification levels for all matrices investigated. Accordingly, the method was successfully introduced in our laboratory with limits of quantification (LOQs) of 10 µg/kg for glyphosate, AMPA and glufosinate in samples from SANTE commodity groups 1, 2, 4a and 5. The reliability and robustness of the method are demonstrated by showing a recovery control chart obtained for glyphosate in randomly selected samples from different commodity groups. Therefore, the samples were spiked up with 10 µg/kg of glyphosate during routine analysis, whereby all recoveries were found to be in the range between 70 and 120%.

State-resolved infrared spectrum of the protonated water dimer: revisiting the characteristic proton transfer doublet peak.

The infrared (IR) spectra of protonated water clusters encode precise information on the dynamics and structure of the hydrated proton. However, the strong anharmonic coupling and quantum effects of these elusive species remain puzzling up to the present day. Here, we report unequivocal evidence that the interplay between the proton transfer and the water wagging motions in the protonated water dimer (Zundel ion) giving rise to the characteristic doublet peak is both more complex and more sensitive to subtle energetic changes than previously thought. In particular, hitherto overlooked low-intensity satellite peaks in the experimental spectrum are now unveiled and mechanistically assigned. Our findings rely on the comparison of IR spectra obtained using two highly accurate potential energy surfaces in conjunction with highly accurate state-resolved quantum simulations. We demonstrate that these high-accuracy simulations are important for providing definite assignments of the complex IR signals of fluxional molecules.

The coupling of the hydrated proton to its first solvation shell.

The Zundel ([Formula: see text]) and Eigen ([Formula: see text]) cations play an important role as intermediate structures for proton transfer processes in liquid water. In the gas phase they exhibit radically different infrared (IR) spectra. The question arises: is there a least common denominator structure that explains the IR spectra of both, the Zundel and Eigen cations, and hence of the solvated proton? Full dimensional quantum simulations of these protonated cations demonstrate that two dynamical water molecules and an excess proton constitute this fundamental subunit. Embedded in the static environment of the parent Eigen cation, this subunit reproduces the positions and broadenings of its main excess-proton bands. In isolation, its spectrum reverts to the well-known Zundel ion. Hence, the dynamics of this subunit polarized by an environment suffice to explain the spectral signatures and anharmonic couplings of the solvated proton in its first solvation shell.

Recursive Editing improves homology-directed repair through retargeting of undesired outcomes.

CRISPR-Cas induced homology-directed repair (HDR) enables the installation of a broad range of precise genomic modifications from an exogenous donor template. However, applications of HDR in human cells are often hampered by poor efficiency, stemming from a preference for error-prone end joining pathways that yield short insertions and deletions. Here, we describe Recursive Editing, an HDR improvement strategy that selectively retargets undesired indel outcomes to create additional opportunities to produce the desired HDR allele. We introduce a software tool, named REtarget, that enables the rational design of Recursive Editing experiments. Using REtarget-designed guide RNAs in single editing reactions, Recursive Editing can simultaneously boost HDR efficiencies and reduce undesired indels. We also harness REtarget to generate databases for particularly effective Recursive Editing sites across the genome, to endogenously tag proteins, and to target pathogenic mutations. Recursive Editing constitutes an easy-to-use approach without potentially deleterious cell manipulations and little added experimental burden.

Representation of Diabatic Potential Energy Matrices for Multiconfiguration Time-Dependent Hartree Treatments of High-Dimensional Nonadiabatic Photodissociation Dynamics.

Conventional quantum mechanical characterization of photodissociation dynamics is restricted by steep scaling laws with respect to the dimensionality of the system. In this work, we examine the applicability of the multi-configurational time-dependent Hartree (MCTDH) method in treating nonadiabatic photodissociation dynamics in two prototypical systems, taking advantage of its favorable scaling laws. To conform to the sum-of-product form, elements of the ab initio diabatic potential energy matrix (DPEM) are re-expressed using the recently proposed Monte Carlo canonical polyadic decomposition method, with enforcement of proper symmetry. The MCTDH absorption spectra and product branching ratios are shown to compare well with those calculated using conventional grid-based methods, demonstrating its promise for treating high-dimensional nonadiabatic photodissociation problems.

High-Dimensional Quantum Dynamics Study on Excitation-Specific Surface Scattering Including Lattice Effects of a Five-Atom Surface Cell.

In this work, high-dimensional (21D) quantum dynamics calculations on the mode-specific surface scattering of a carbon monoxide molecule on a copper(100) surface with lattice effects of a five-atom surface cell are performed through the multilayer multiconfiguration time-dependent Hartree (ML-MCTDH) method. We employ a surface model in which five surface atoms near the impact site are treated as fully flexible quantum particles, while all other more distant atoms are kept at fixed locations. To efficiently perform the 21D ML-MCTDH wave packet propagation, the potential energy surface is transferred to a canonical polyadic decomposition form with the aid of a Monte Carlo-based method. Excitation-specific sticking probabilities of CO on Cu(100) are computed, and lattice effects caused by the flexible surface atoms are demonstrated by comparison with sticking probabilities computed for a rigid surface. The dependence of the sticking probability of the initial state of the system is studied, and it is found that the sticking probability is reduced when the surface atom on the impact site is initially vibrationally excited.

PnB Designer: a web application to design prime and base editor guide RNAs for animals and plants.

BACKGROUND: The rapid expansion of the CRISPR toolbox through tagging effector domains to either enzymatically inactive Cas9 (dCas9) or Cas9 nickase (nCas9) has led to several promising new gene editing strategies. Recent additions include CRISPR cytosine or adenine base editors (CBEs and ABEs) and the CRISPR prime editors (PEs), in which a deaminase or reverse transcriptase are fused to nCas9, respectively. These tools hold great promise to model and correct disease-causing mutations in animal and plant models. But so far, no widely-available tools exist to automate the design of both BE and PE reagents. RESULTS: We developed PnB Designer, a web-based application for the design of pegRNAs for PEs and guide RNAs for BEs. PnB Designer makes it easy to design targeting guide RNAs for single or multiple targets on a variant or reference genome from organisms spanning multiple kingdoms. With PnB Designer, we designed pegRNAs to model all known disease causing mutations available in ClinVar. Additionally, PnB Designer can be used to design guide RNAs to install or revert a SNV, scanning the genome with one CBE and seven different ABE PAM variants and returning the best BE to use. PnB Designer is publicly accessible at http://fgcz-shiny.uzh.ch/PnBDesigner/ CONCLUSION: With PnB Designer we created a user-friendly design tool for CRISPR PE and BE reagents, which should simplify choosing editing strategy and avoiding design complications.

Controlled Cycling and Quiescence Enables Efficient HDR in Engraftment-Enriched Adult Hematopoietic Stem and Progenitor Cells.

Genome editing often takes the form of either error-prone sequence disruption by non-homologous end joining (NHEJ) or sequence replacement by homology-directed repair (HDR). Although NHEJ is generally effective, HDR is often difficult in primary cells. Here, we use a combination of immunophenotyping, next-generation sequencing, and single-cell RNA sequencing to investigate and reprogram genome editing outcomes in subpopulations of adult hematopoietic stem and progenitor cells. We find that although quiescent stem-enriched cells mostly use NHEJ, non-quiescent cells with the same immunophenotype use both NHEJ and HDR. Inducing quiescence before editing results in a loss of HDR in all cell subtypes. We develop a strategy of controlled cycling and quiescence that yields a 6-fold increase in the HDR/NHEJ ratio in quiescent stem cells ex vivo and in vivo. Our results highlight the tension between editing and cellular physiology and suggest strategies to manipulate quiescent cells for research and therapeutic genome editing.

Integrated genomic analysis reveals mutated ELF3 as a potential gallbladder cancer vaccine candidate.

Gallbladder cancer (GBC) is an aggressive gastrointestinal malignancy with no approved targeted therapy. Here, we analyze exomes (n = 160), transcriptomes (n = 115), and low pass whole genomes (n = 146) from 167 gallbladder cancers (GBCs) from patients in Korea, India and Chile. In addition, we also sequence samples from 39 GBC high-risk patients and detect evidence of early cancer-related genomic lesions. Among the several significantly mutated genes not previously linked to GBC are ETS domain genes ELF3 and EHF, CTNNB1, APC, NSD1, KAT8, STK11 and NFE2L2. A majority of ELF3 alterations are frame-shift mutations that result in several cancer-specific neoantigens that activate T-cells indicating that they are cancer vaccine candidates. In addition, we identify recurrent alterations in KEAP1/NFE2L2 and WNT pathway in GBC. Taken together, these define multiple targetable therapeutic interventions opportunities for GBC treatment and management.

ATF4 Regulates MYB to Increase γ-Globin in Response to Loss of ß-Globin.

ß-Hemoglobinopathies can trigger rapid production of red blood cells in a process known as stress erythropoiesis. Cellular stress prompts differentiating erythroid precursors to express high levels of fetal γ-globin. However, the mechanisms underlying γ-globin production during cellular stress are still poorly defined. Here, we use CRISPR-Cas genome editing to model the stress caused by reduced levels of adult ß-globin. We find that decreased ß-globin is sufficient to induce robust re-expression of γ-globin, and RNA sequencing (RNA-seq) of differentiating isogenic erythroid precursors implicates ATF4 as a causal regulator of this response. ATF4 binds within the HBS1L-MYB intergenic enhancer and regulates expression of MYB, a known γ-globin regulator. Overall, the reduction of ATF4 upon ß-globin knockout decreases the levels of MYB and BCL11A. Identification of ATF4 as a key regulator of globin compensation adds mechanistic insight to the poorly understood phenomenon of stress-induced globin compensation and could inform strategies to treat hemoglobinopathies.

Effect of progesterone on Candida albicans biofilm formation under acidic conditions: A transcriptomic analysis.

Vulvovaginal candidiasis (VVC) caused by Candida albicans is a common disease worldwide. A very important C. albicans virulence factor is its ability to form biofilms on epithelium and/or on intrauterine devices promoting VVC. It has been shown that VVC has a hormonal dependency and that progesterone affects virulence traits of C. albicans cells. To understand how the acidic environment (pH 4) and progesterone (either alone and in combination) modulate C. albicans response during formation of biofilm, a transcriptomic analysis was performed together with characterization of the biofilm properties. Compared to planktonic cells, acidic biofilm-cells exhibited major changes in their transcriptome, including modifications in the expression of 286 genes that were not previously associated with biofilm formation in C. albicans. The vast majority of the genes up-regulated in the acidic biofilm cells (including those uniquely identified in our study) are known targets of Sfl1, and consistently, Sfl1 deletion is herein shown to impair the formation of acidic biofilms (pH 4). Under the acidic conditions used, the presence of progesterone reduced C. albicans biofilm biomass and structural cohesion. Transcriptomic analysis of biofilms developed in the presence of progesterone led to the identification of 65 down-regulated genes including, among others, the regulator Tec1 and several of its target genes, suggesting that the function of this transcription factor is inhibited by the presence of the hormone. Additionally, progesterone reduced the susceptibility of biofilm cells to fluconazole, consistent with an up-regulation of efflux pumps. Overall, the results of this study show that progesterone modulates C. albicans biofilm formation and genomic expression under acidic conditions, which may have implications for C. albicans pathogenicity in the vaginal environment.

Transforming high-dimensional potential energy surfaces into a canonical polyadic decomposition using Monte Carlo methods.

A Monte Carlo method is proposed for transforming high-dimensional potential energy surfaces evaluated on discrete grid points into a sum-of-products form, more precisely into a Canonical Polyadic Decomposition form. To this end, a modified existing ansatz based on the alternating least squares method is used, in which numerically exact integrals are replaced with Monte Carlo integrals. This largely reduces the numerical cost by avoiding the evaluation of the potential on all grid points and allows the treatment of surfaces with many degrees of freedom. Calculations on the 15D potential of the protonated water dimer (Zundel cation) in a sum-of-products form are presented and compared to the results obtained in a previous work [M. Schröder and H.-D. Meyer, J. Chem. Phys. 147, 064105 (2017)], where a sum-of-products form of the potential was obtained in the Tucker format.

The Indian cobra reference genome and transcriptome enables comprehensive identification of venom toxins.

Snakebite envenoming is a serious and neglected tropical disease that kills ~100,000 people annually. High-quality, genome-enabled comprehensive characterization of toxin genes will facilitate development of effective humanized recombinant antivenom. We report a de novo near-chromosomal genome assembly of Naja naja, the Indian cobra, a highly venomous, medically important snake. Our assembly has a scaffold N50 of 223.35 Mb, with 19 scaffolds containing 95% of the genome. Of the 23,248 predicted protein-coding genes, 12,346 venom-gland-expressed genes constitute the 'venom-ome' and this included 139 genes from 33 toxin families. Among the 139 toxin genes were 19 'venom-ome-specific toxins' (VSTs) that showed venom-gland-specific expression, and these probably encode the minimal core venom effector proteins. Synthetic venom reconstituted through recombinant VST expression will aid in the rapid development of safe and effective synthetic antivenom. Additionally, our genome could serve as a reference for snake genomes, support evolutionary studies and enable venom-driven drug discovery.

Correlating Genotype and Phenotype in the Asexual Yeast Candida orthopsilosis Implicates ZCF29 in Sensitivity to Caffeine.

Candida orthopsilosis is diploid asexual yeast that causes human disease. Most C. orthopsilosis isolates arose from at least four separate hybridizations between related, but not identical, parents. Here, we used population genomics data to correlate genotypic and phenotypic variation in 28 C. orthopsilosis isolates. We used cosine similarity scores to identify 65 variants with potential high-impact (deleterious effects) that correlated with specific phenotypes. Of these, 19 were Single Nucleotide Polymorphisms (SNPs) that changed stop or start codons, or splice sites. One variant resulted in a premature stop codon in both alleles of the gene ZCF29 in C. orthopsilosis isolate 185, which correlated with sensitivity to nystatin and caffeine. We used CRISPR-Cas9 editing to introduce this polymorphism into two resistant C. orthopsilosis isolates. Introducing the stop codon resulted in sensitivity to caffeine and to ketoconazole, but not to nystatin. Our analysis shows that it is possible to associate genomic variants with phenotype in asexual Candida species, but that only a small amount of genomic variation can be easily explored.