Salmonella Typhimurium sequence type (ST) 313 causes invasive nontyphoidal Salmonella (iNTS) disease in sub-Saharan Africa, targeting susceptible HIV+, malarial, or malnourished individuals. An in-depth genomic comparison between the ST313 isolate D23580 and the well-characterized ST19 isolate 4/74 that causes gastroenteritis across the globe revealed extensive synteny. To understand how the 856 nucleotide variations generated phenotypic differences, we devised a large-scale experimental approach that involved the global gene expression analysis of strains D23580 and 4/74 grown in 16 infection-relevant growth conditions. Comparison of transcriptional patterns identified virulence and metabolic genes that were differentially expressed between D23580 versus 4/74, many of which were validated by proteomics. We also uncovered the S. Typhimurium D23580 and 4/74 genes that showed expression differences during infection of murine macrophages. Our comparative transcriptomic data are presented in a new enhanced version of the Salmonella expression compendium, SalComD23580: http://bioinf.gen.tcd.ie/cgi-bin/salcom_v2.pl. We discovered that the ablation of melibiose utilization was caused by three independent SNP mutations in D23580 that are shared across ST313 lineage 2, suggesting that the ability to catabolize this carbon source has been negatively selected during ST313 evolution. The data revealed a novel, to our knowledge, plasmid maintenance system involving a plasmid-encoded CysS cysteinyl-tRNA synthetase, highlighting the power of large-scale comparative multicondition analyses to pinpoint key phenotypic differences between bacterial pathovariants.

We describe here the PRIMO (PRotein Intermediate Model) force field, a physics-based fully transferable additive coarse-grained potential energy function that is compatible with an all-atom force field for multi-scale simulations. The energy function consists of standard molecular dynamics energy terms plus a hydrogen-bonding potential term and is mainly parameterized based on the CHARMM22/CMAP force field in a bottom-up fashion. The solvent is treated implicitly via the generalized Born model. The bonded interactions are either harmonic or distance-based spline interpolated potentials. These potentials are defined on the basis of all-atom molecular dynamics (MD) simulations of dipeptides with the CHARMM22/CMAP force field. The non-bonded parameters are tuned by matching conformational free energies of diverse set of conformations with that of CHARMM all-atom results. PRIMO is designed to provide a correct description of conformational distribution of the backbone (ϕ/ψ) and side chains (χ1) for all amino acids with a CMAP correction term. The CMAP potential in PRIMO is optimized based on the new CHARMM C36 CMAP. The resulting optimized force field has been applied in MD simulations of several proteins of 36–155 amino acids and shown that the root-mean-squared-deviation of the average structure from the corresponding crystallographic structure varies between 1.80 and 4.03 Å. PRIMO is shown to fold several small peptides to their native-like structures from extended conformations. These results suggest the applicability of the PRIMO force field in the study of protein structures in aqueous solution, structure predictions as well as ab initio folding of small peptides.

Bloodstream infections caused by nontyphoidal Salmonella are a major public health concern in Africa, causing ~49,600 deaths every year. The most common Salmonella enterica pathovariant associated with invasive nontyphoidal Salmonella disease is Salmonella Typhimurium sequence type (ST)313. It has been proposed that antimicrobial resistance and genome degradation has contributed to the success of ST313 lineages in Africa, but the evolutionary trajectory of such changes was unclear. Here, to define the evolutionary dynamics of ST313, we sub-sampled from two comprehensive collections of Salmonella isolates from African patients with bloodstream infections, spanning 1966 to 2018. The resulting 680 genome sequences led to the discovery of a pan-susceptible ST313 lineage (ST313 L3), which emerged in Malawi in 2016 and is closely related to ST313 variants that cause gastrointestinal disease in the United Kingdom and Brazil. Genomic analysis revealed degradation events in important virulence genes in ST313 L3, which had not occurred in other ST313 lineages. Despite arising only recently in the clinic, ST313 L3 is a phylogenetic intermediate between ST313 L1 and L2, with a characteristic accessory genome. Our in-depth genotypic and phenotypic characterization identifies the crucial loss-of-function genetic events that occurred during the stepwise evolution of invasive S. Typhimurium across Africa.

The extension of life span driven by 40% caloric restriction (CR) in rodents causes trade-offs in growth, reproduction, and immune defense that make it difficult to identify therapeutically relevant CR-mimetic targets. We report that about 14% CR for 2 years in healthy humans improved thymopoiesis and was correlated with mobilization of intrathymic ectopic lipid. CR-induced transcriptional reprogramming in adipose tissue implicated pathways regulating mitochondrial bioenergetics, anti-inflammatory responses, and longevity. Expression of the gene Pla2g7 encoding platelet activating factor acetyl hydrolase (PLA2G7) is inhibited in humans undergoing CR. Deletion of Pla2g7 in mice showed decreased thymic lipoatrophy, protection against age-related inflammation, lowered NLRP3 inflammasome activation, and improved metabolic health. Therefore, the reduction of PLA2G7 may mediate the immunometabolic effects of CR and could potentially be harnessed to lower inflammation and extend the health span.

Macromolecular crowding is recognized as an important factor influencing folding and conformational dynamics of proteins and nucleic acids. Previous views of crowding have focused on the mostly entropic volume exclusion effect of crowding but recent studies are indicating the importance of enthalpic effects, in particular changes in electrostatic interactions due to crowding. Here, temperature replica exchange molecular dynamics simulations of trp-cage and melittin in the presence of explicit protein crowders are presented to further examine the effect of protein crowders on peptide dynamics. The simulations involve a three-component multiscale modeling scheme where the peptides are represented at atomistic level, the crowder proteins at a coarse-grained level, and the surrounding aqueous solvent as implicit solvent. This scheme optimally balances a physically realistic description for the peptide with computational efficiency. The multiscale simulations were compared with simulations of the same peptides in different dielectric environments with dielectric constants ranging from 5 to 80. It is found that the sampling in the presence of the crowders resembles sampling with reduced dielectric constants between 10 and 40. Furthermore, diverse conformational ensembles are generated in the presence of crowders including partially unfolded states for trp-cage. These findings emphasize the importance of enthalpic interactions over volume exclusion effects in describing the effects of cellular crowding.