FragBuilder:an efficient Python library to setup quantum chemistry calculations on peptides models

We present a powerful Python library to quickly and efficiently generate realistic peptide model structures. The library makes it possible to quickly set up quantum mechanical calculations on model peptide structures. It is possible to manually specify a specific conformation of the peptide. Additionally the library also offers sampling of backbone conformations and side chain rotamer conformations from continuous distributions. The generated peptides can then be geometry optimized by the MMFF94 molecular mechanics force field via convenient functions inside the library. Finally, it is possible to output the resulting structures directly to files in a variety of useful formats, such as XYZ or PDB formats, or directly as input files for a quantum chemistry program. FragBuilder is freely available at https://github.com/jensengroup/fragbuilder/ under the terms of the BSD open source license

Bayesian inference of protein structure from chemical shift data

Protein chemical shifts are routinely used to augment molecular mechanics force fields in protein structure simulations, with weights of the chemical shift restraints determined empirically. These weights, however, might not be an optimal descriptor of a given protein structure and predictive model, and a bias is introduced which might result in incorrect structures. In the inferential structure determination framework, both the unknown structure and the disagreement between experimental and back-calculated data are formulated as a joint probability distribution, thus utilizing the full information content of the data. Here, we present the formulation of such a probability distribution where the error in chemical shift prediction is described by either a Gaussian or Cauchy distribution. The methodology is demonstrated and compared to a set of empirically weighted potentials through Markov chain Monte Carlo simulations of three small proteins (ENHD, Protein G and the SMN Tudor Domain) using the PROFASI force field and the chemical shift predictor CamShift. Using a clustering-criterion for identifying the best structure, together with the addition of a solvent exposure scoring term, the simulations suggests that sampling both the structure and the uncertainties in chemical shift prediction leads more accurate structures compared to conventional methods using empirical determined weights. The Cauchy distribution, using either sampled uncertainties or predetermined weights, did, however, result in overall better convergence to the native fold, suggesting that both types of distribution might be useful in different aspects of the protein structure prediction

LegoNet: Alternating Model Blocks for Medical Image Segmentation

Since the emergence of convolutional neural networks (CNNs), and later vision transformers (ViTs), the common paradigm for model development has always been using a set of identical block types with varying parameters/hyper-parameters. To leverage the benefits of different architectural designs (e.g. CNNs and ViTs), we propose to alternate structurally different types of blocks to generate a new architecture, mimicking how Lego blocks can be assembled together. Using two CNN-based and one SwinViT-based blocks, we investigate three variations to the so-called LegoNet that applies the new concept of block alternation for the segmentation task in medical imaging. We also study a new clinical problem which has not been investigated before, namely the right internal mammary artery (RIMA) and perivascular space segmentation from computed tomography angiography (CTA) which has demonstrated a prognostic value to major cardiovascular outcomes. We compare the model performance against popular CNN and ViT architectures using two large datasets (e.g. achieving 0.749 dice similarity coefficient (DSC) on the larger dataset). We evaluate the performance of the model on three external testing cohorts as well, where an expert clinician made corrections to the model segmented results (DSC>0.90 for the three cohorts). To assess our proposed model for suitability in clinical use, we perform intra- and inter-observer variability analysis. Finally, we investigate a joint self-supervised learning approach to assess its impact on model performance. The code and the pretrained model weights will be available upon acceptance.Comment: 12 pages, 5 figures, 4 table

Inflammatory risk and cardiovascular events in patients without obstructive coronary artery disease: the ORFAN multicentre, longitudinal cohort study

Background: Coronary computed tomography angiography (CCTA) is the first line investigation for chest pain, and it is used to guide revascularisation. However, the widespread adoption of CCTA has revealed a large group of individuals without obstructive coronary artery disease (CAD), with unclear prognosis and management. Measurement of coronary inflammation from CCTA using the perivascular fat attenuation index (FAI) Score could enable cardiovascular risk prediction and guide the management of individuals without obstructive CAD. The Oxford Risk Factors And Non-invasive imaging (ORFAN) study aimed to evaluate the risk profile and event rates among patients undergoing CCTA as part of routine clinical care in the UK National Health Service (NHS); to test the hypothesis that coronary arterial inflammation drives cardiac mortality or major adverse cardiac events (MACE) in patients with or without CAD; and to externally validate the performance of the previously trained artificial intelligence (AI)-Risk prognostic algorithm and the related AI-Risk classification system in a UK population. Methods: This multicentre, longitudinal cohort study included 40 091 consecutive patients undergoing clinically indicated CCTA in eight UK hospitals, who were followed up for MACE (ie, myocardial infarction, new onset heart failure, or cardiac death) for a median of 2·7 years (IQR 1·4–5·3). The prognostic value of FAI Score in the presence and absence of obstructive CAD was evaluated in 3393 consecutive patients from the two hospitals with the longest follow-up (7·7 years [6·4–9·1]). An AI-enhanced cardiac risk prediction algorithm, which integrates FAI Score, coronary plaque metrics, and clinical risk factors, was then evaluated in this population. Findings: In the 2·7 year median follow-up period, patients without obstructive CAD (32 533 [81·1%] of 40 091) accounted for 2857 (66·3%) of the 4307 total MACE and 1118 (63·7%) of the 1754 total cardiac deaths in the whole of Cohort A. Increased FAI Score in all the three coronary arteries had an additive impact on the risk for cardiac mortality (hazard ratio [HR] 29·8 [95% CI 13·9–63·9], p<0·001) or MACE (12·6 [8·5–18·6], p<0·001) comparing three vessels with an FAI Score in the top versus bottom quartile for each artery. FAI Score in any coronary artery predicted cardiac mortality and MACE independently from cardiovascular risk factors and the presence or extent of CAD. The AI-Risk classification was positively associated with cardiac mortality (6·75 [5·17–8·82], p<0·001, for very high risk vs low or medium risk) and MACE (4·68 [3·93–5·57], p<0·001 for very high risk vs low or medium risk). Finally, the AI-Risk model was well calibrated against true events. Interpretation: The FAI Score captures inflammatory risk beyond the current clinical risk stratification and CCTA interpretation, particularly among patients without obstructive CAD. The AI-Risk integrates this information in a prognostic algorithm, which could be used as an alternative to traditional risk factor-based risk calculators