A parallel multistate framework for atomistic non-equilibrium reaction dynamics of solutes in strongly interacting organic solvents

We describe a parallel linear-scaling computational framework developed to implement arbitrarily large multi-state empirical valence bond (MS-EVB) calculations within CHARMM. Forces are obtained using the Hellman-Feynmann relationship, giving continuous gradients, and excellent energy conservation. Utilizing multi-dimensional Gaussian coupling elements fit to CCSD(T)-F12 electronic structure theory, we built a 64-state MS-EVB model designed to study the F + CD3CN -> DF + CD2CN reaction in CD3CN solvent. This approach allows us to build a reactive potential energy surface (PES) whose balanced accuracy and efficiency considerably surpass what we could achieve otherwise. We use our PES to run MD simulations, and examine a range of transient observables which follow in the wake of reaction, including transient spectra of the DF vibrational band, time dependent profiles of vibrationally excited DF in CD3CN solvent, and relaxation rates for energy flow from DF into the solvent, all of which agree well with experimental observations. Immediately following deuterium abstraction, the nascent DF is in a non-equilibrium regime in two different respects: (1) it is highly excited, with ~23 kcal mol-1 localized in the stretch; and (2) not yet Hydrogen bonded to the CD3CN solvent, its microsolvation environment is intermediate between the non-interacting gas-phase limit and the solution-phase equilibrium limit. Vibrational relaxation of the nascent DF results in a spectral blue shift, while relaxation of its microsolvation environment results in a red shift. These two competing effects result in a post-reaction relaxation profile distinct from that observed when DF vibration excitation occurs within an equilibrium microsolvation environment. The parallel software framework presented in this paper should be more broadly applicable to a range of complex reactive systems.Comment: 58 pages and 29 Figure

Comparative Fungal Community Analyses Using Metatranscriptomics and Internal Transcribed Spacer Amplicon Sequencing from Norway Spruce

The health, growth, and fitness of boreal forest trees are impacted and improved by their associated microbiomes. Microbial gene expression and functional activity can be assayed with RNA sequencing (RNA-Seq) data from host samples. In contrast, phylogenetic marker gene amplicon sequencing data are used to assess taxonomic composition and community structure of the microbiome. Few studies have considered how much of this structural and taxonomic information is included in transcriptomic data from matched samples. Here, we described fungal communities using both host-derived RNA-Seq and fungal ITS1 DNA amplicon sequencing to compare the outcomes between the methods. We used a panel of root and needle samples from the coniferous tree species Picea abies (Norway spruce) growing in untreated (nutrient-deficient) and nutrient-enriched plots at the Flakaliden forest research site in boreal northern Sweden. We show that the relationship between samples and alpha and beta diversity indicated by the fungal transcriptome is in agreement with that generated by the ITS data, while also identifying a lack of taxonomic overlap due to limitations imposed by current database coverage. Furthermore, we demonstrate how metatranscriptomics data additionally provide biologically informative functional insights. At the community level, there were changes in starch and sucrose metabolism, biosynthesis of amino acids, and pentose and glucuronate interconversions, while processing of organic macromolecules, including aromatic and heterocyclic compounds, was enriched in transcripts assigned to the genus Cortinarius.IMPORTANCE A deeper understanding of microbial communities associated with plants is revealing their importance for plant health and productivity. RNA extracted from plant field samples represents the host and other organisms present. Typically, gene expression studies focus on the plant component or, in a limited number of studies, expression in one or more associated organisms. However, metatranscriptomic data are rarely used for taxonomic profiling, which is currently performed using amplicon approaches. We created an assembly-based, reproducible, and hardware-agnostic workflow to taxonomically and functionally annotate fungal RNA-Seq data obtained from Norway spruce roots, which we compared to matching ITS amplicon sequencing data. While we identified some limitations and caveats, we show that functional, taxonomic, and compositional insights can all be obtained from RNA-Seq data. These findings highlight the potential of metatranscriptomics to advance our understanding of interaction, response, and effect between host plants and their associated microbial communities

FHOD1 and FMNL1 formin proteins in intestinal gastric cancer: correlation with tumor-infiltrating T lymphocytes and molecular subtypes

Background Gastric cancer (GC) is the third most common cause of cancer death. Intestinal type GC is a molecularly diverse disease. Formins control cytoskeletal processes and have been implicated in the progression of many cancers. Their clinical significance in GC remains unclear. Here, we characterize the expression of formin proteins FHOD1 and FMNL1 in intestinal GC tissue samples and investigate their association with clinical parameters, GC molecular subtypes and intratumoral T lymphocytes. Methods The prognostic significance of FHOD1 and FMNL1 mRNA expression was studied with Kaplan-Meier analyses in an online database. The expression of FHOD1 and FMNL1 proteins was characterized in GC cells, and in non-neoplastic and malignant tissues utilizing tumor microarrays of intestinal GC representing different molecular subtypes. FHOD1 and FMNL1 expression was correlated with clinical parameters, molecular features and T lymphocyte infiltration. Immunohistochemical expression of neither formin correlated with survival. Results Kaplan-Meier analysis associated high FHOD1 and FMNL1 mRNA expression with reduced overall survival (OS). Characterization of FHOD1 and FMNL1 in GC cells showed cytoplasmic expression along the actin filaments. Similar pattern was recapitulated in GC tissue samples. Elevated FMNL1 was associated with larger tumor size and higher disease stage. Downregulation of FHOD1 associated with TP53-mutated GC tumors. Tumor cell FHOD1 expression strongly correlated with high numbers of tumor-infiltrating CD8 + lymphocytes. Conclusions FHOD1 and FMNL1 proteins are expressed in the tumor cells of intestinal GC and significantly associate with clinical parameters without direct prognostic significance. FHOD1 correlates with high intratumoral CD8 + T lymphocyte infiltration in this cohort.Peer reviewe

Development of an Industry 4.0 Demonstrator Using Sequence Planner and ROS2

In many modern automation solutions, manual off-line programming is being replaced by online algorithms that dynamically perform tasks based on the state of the environment. Complexities of such systems are pushed even further with collaboration among robots and humans, where intelligent machines and learning algorithms are replacing more traditional automation solutions. This chapter describes the development of an industrial demonstrator using a control infrastructure called Sequence Planner (SP), and presents some lessons learned during development. SP is based on ROS2 and it is designed to aid in handling the increased complexity of these new systems using formal models and online planning algorithms to coordinate the actions of robots and other devices. During development, SP can auto generate ROS nodes and message types as well as support continuous validation and testing. SP is also designed with the aim to handle traditional challenges of automation software development such as safety, reliability and efficiency. In this chapter, it is argued that ROS2 together with SP could be an enabler of intelligent automation for the next industrial revolution

Immune response to SARS-CoV-2 mRNA vaccination in multiple sclerosis patients after rituximab treatment interruption

Peripheral B cell depletion via anti-CD20 treatment is a highly effective disease-modifying treatment for reducing new relapses in multiple sclerosis (MS) patients. A drawback of rituximab (RTX) and other anti-CD20 antibodies is a poor immune response to vaccination. While this can be mitigated by treatment interruption of at least six months prior to vaccination, the timing to resume treatment while maintaining subsequent vaccine responses remains undetermined. Here, we characterized SARS-CoV-2 S-directed antibody and B cell responses throughout three BNT162b2 mRNA vaccine doses in RTX-treated MS patients, with the first two doses given during treatment interruption. We examined B-cell mediated immune responses in blood samples from patients with RTX-treated MS throughout three BNT162b2 vaccine doses, compared to an age- and sex-matched healthy control group. The first vaccine dose was given 1.3 years (median) after the last RTX infusion, the second dose one month after the first, and the third dose four weeks after treatment re-initiation. We analyzed SARS-CoV-2 S-directed antibody levels using enzyme-linked immunosorbent assay (ELISA), and the neutralization capacity of patient serum against SARS-CoV-2 S-pseudotyped lentivirus using luciferase reporter assay. In addition, we assessed switched memory (CD19+CD20+CD27+IgD-), unswitched memory (CD19+CD20+CD27+IgD+), naïve (CD19+CD20+CD27-IgD+), and double negative (DN, CD19+CD20+CD27-IgD-) B cell frequencies, as well as their SARS-CoV-2 S-specific (CoV+) and Decay Accelerating Factor-negative (DAF-) subpopulations, using flow cytometry. After two vaccine doses, S-binding antibody levels and neutralization capacity in SARS-CoV-2-naïve MS patients were comparable to vaccinated healthy controls, albeit with greater variation. Higher antibody response levels and CoV+-DN B cell frequencies after the second vaccine dose were predictive of a boost effect after the third dose, even after re-initiation of rituximab treatment. MS patients also exhibited lower frequencies of DAF- memory B cells, a suggested proxy for germinal centre activity, than control individuals. S-binding antibody levels in RTX-treated MS patients after two vaccine doses could help determine which individuals would need to move up their next vaccine booster dose or postpone their next RTX infusion. Our findings also offer first indications on the potential importance of antigenic stimulation of DN B cells and long-term impairment of germinal centre activity in rituximab-treated MS patients

Sources of variation in cuticular hydrocarbons in the ant formica exsecta

Phenotypic variation arises from interactions between genotype and environment, although how variation is produced and then maintained remains unclear. The discovery of the nest-mate recognition system in Formica exsecta ants has allowed phenotypic variation in chemical profiles to be quantified across a natural population of 83 colonies. We investigated if this variation was correlated or not with intrinsic (genetic relatedness), extrinsic (location, light, temperature) or social (queen number) factors. (Z)-9-Alkenes and n-alkanes showed different patterns of variance: island (location) explained only 0.2% of the variation in (Z)-9-alkenes, but 21¬–29% in n-alkanes, whereas colony of origin explained 96% and 45–49% of the variation in (Z)-9-alkenes and n-alkanes, respectively. By contrast, within-colony variance of (Z)-9-alkenes was 4%, and 23–34% in n-alkanes, supporting the function of the former as recognition cues. (Z)-9-Alkene and n-alkane profiles were correlated with the genetic distance between colonies. Only n-alkane profiles diverged with increasing spatial distance. Sampling year explained a small (5%), but significant, amount of the variation in the (Z)-9-alkenes, but there was no consistent directional trend. Polygynous colonies and populous monogynous colonies were dominated by a rich C23:1 profile. We found no associations between worker size, mound exposure, or humidity, although effect sizes for the latter two factors were considerable. The results support the conjecture that genetic factors are the most likely source of between-colony variation in cuticular hydrocarbons

Enhanced cosmic-ray flux toward zeta Persei inferred from laboratory study of H3+ - e- recombination rate

The H3+ molecular ion plays a fundamental role in interstellar chemistry, as it initiates a network of chemical reactions that produce many interstellar molecules. In dense clouds, the H3+ abundance is understood using a simple chemical model, from which observations of H3+ yield valuable estimates of cloud path length, density, and temperature. On the other hand, observations of diffuse clouds have suggested that H3+ is considerably more abundant than expected from the chemical models. However, diffuse cloud models have been hampered by the uncertain values of three key parameters: the rate of H3+ destruction by electrons, the electron fraction, and the cosmic-ray ionisation rate. Here we report a direct experimental measurement of the H3+ destruction rate under nearly interstellar conditions. We also report the observation of H3+ in a diffuse cloud (towards zeta Persei) where the electron fraction is already known. Taken together, these results allow us to derive the value of the third uncertain model parameter: we find that the cosmic-ray ionisation rate in this sightline is forty times faster than previously assumed. If such a high cosmic-ray flux is indeed ubiquitous in diffuse clouds, the discrepancy between chemical models and the previous observations of H3+ can be resolved.Comment: 6 pages, Nature, in pres

Lymphocyte predominant cells detect Moraxella catarrhalis-derived antigens in nodular lymphocyte-predominant Hodgkin lymphoma.

Nodular lymphocyte-predominant Hodgkin lymphoma (NLPHL) is a rare lymphoma of B-cell origin with frequent expression of functional B-cell receptors (BCRs). Here we report that expression cloning followed by antigen screening identifies DNA-directed RNA polymerase beta' (RpoC) from Moraxella catarrhalis as frequent antigen of BCRs of IgD <sup>+</sup> LP cells. Patients show predominance of HLA-DRB1*04/07 and the IgVH genes encode extraordinarily long CDR3s. High-titer, light-chain-restricted anti-RpoC IgG1/κ-type serum-antibodies are additionally found in these patients. RpoC and MID/hag, a superantigen co-expressed by Moraxella catarrhalis that is known to activate IgD <sup>+</sup> B cells by binding to the Fc domain of IgD, have additive activation effects on the BCR, the NF-κB pathway and the proliferation of IgD <sup>+</sup> DEV cells expressing RpoC-specific BCRs. This suggests an additive antigenic and superantigenic stimulation of B cells with RpoC-specific IgD <sup>+</sup> BCRs under conditions of a permissive MHC-II haplotype as a model of NLPHL lymphomagenesis, implying future treatment strategies

Excitons in a Photosynthetic Light-Harvesting System: A Combined Molecular Dynamics/Quantum Chemistry and Polaron Model Study

The dynamics of pigment-pigment and pigment-protein interactions in light-harvesting complexes is studied with a novel approach which combines molecular dynamics (MD) simulations with quantum chemistry (QC) calculations. The MD simulations of an LH-II complex, solvated and embedded in a lipid bilayer at physiological conditions (with total system size of 87,055 atoms) revealed a pathway of a water molecule into the B800 binding site, as well as increased dimerization within the B850 BChl ring, as compared to the dimerization found for the crystal structure. The fluctuations of pigment (B850 BChl) excitation energies, as a function of time, were determined via ab initio QC calculations based on the geometries that emerged from the MD simulations. From the results of these calculations we constructed a time-dependent Hamiltonian of the B850 exciton system from which we determined the linear absorption spectrum. Finally, a polaron model is introduced to describe quantum mechanically both the excitonic and vibrational (phonon) degrees of freedom. The exciton-phonon coupling that enters into the polaron model, and the corresponding phonon spectral function are derived from the MD/QC simulations. It is demonstrated that, in the framework of the polaron model, the absorption spectrum of the B850 excitons can be calculated from the autocorrelation function of the excitation energies of individual BChls, which is readily available from the combined MD/QC simulations. The obtained result is in good agreement with the experimentally measured absorption spectrum.Comment: REVTeX3.1, 23 pages, 13 (EPS) figures included. A high quality PDF file of the paper is available at http://www.ks.uiuc.edu/Publications/Papers/PDF/DAMJ2001/DAMJ2001.pd