An amino-domino model described by a cross-peptide-bond Ramachandran plot defines amino acid pairs as local structural units.

Protein structure, both at the global and local level, dictates function. Proteins fold from chains of amino acids, forming secondary structures, α-helices and ß-strands, that, at least for globular proteins, subsequently fold into a three-dimensional structure. Here, we show that a Ramachandran-type plot focusing on the two dihedral angles separated by the peptide bond, and entirely contained within an amino acid pair, defines a local structural unit. We further demonstrate the usefulness of this cross-peptide-bond Ramachandran plot by showing that it captures ß-turn conformations in coil regions, that traditional Ramachandran plot outliers fall into occupied regions of our plot, and that thermophilic proteins prefer specific amino acid pair conformations. Further, we demonstrate experimentally that the effect of a point mutation on backbone conformation and protein stability depends on the amino acid pair context, i.e., the identity of the adjacent amino acid, in a manner predictable by our method.

Meeting the unmet needs of clinicians from AI systems showcased for cardiology with deep-learning-based ECG analysis.

Despite their great promise, artificial intelligence (AI) systems have yet to become ubiquitous in the daily practice of medicine largely due to several crucial unmet needs of healthcare practitioners. These include lack of explanations in clinically meaningful terms, handling the presence of unknown medical conditions, and transparency regarding the system's limitations, both in terms of statistical performance as well as recognizing situations for which the system's predictions are irrelevant. We articulate these unmet clinical needs as machine-learning (ML) problems and systematically address them with cutting-edge ML techniques. We focus on electrocardiogram (ECG) analysis as an example domain in which AI has great potential and tackle two challenging tasks: the detection of a heterogeneous mix of known and unknown arrhythmias from ECG and the identification of underlying cardio-pathology from segments annotated as normal sinus rhythm recorded in patients with an intermittent arrhythmia. We validate our methods by simulating a screening for arrhythmias in a large-scale population while adhering to statistical significance requirements. Specifically, our system 1) visualizes the relative importance of each part of an ECG segment for the final model decision; 2) upholds specified statistical constraints on its out-of-sample performance and provides uncertainty estimation for its predictions; 3) handles inputs containing unknown rhythm types; and 4) handles data from unseen patients while also flagging cases in which the model's outputs are not usable for a specific patient. This work represents a significant step toward overcoming the limitations currently impeding the integration of AI into clinical practice in cardiology and medicine in general.

A dataset of alternately located segments in protein crystal structures.

Protein Data Bank (PDB) files list the relative spatial location of atoms in a protein structure as the final output of the process of fitting and refining to experimentally determined electron density measurements. Where experimental evidence exists for multiple conformations, atoms are modelled in alternate locations. Programs reading PDB files commonly ignore these alternate conformations by default leaving users oblivious to the presence of alternate conformations in the structures they analyze. This has led to underappreciation of their prevalence, under characterisation of their features and limited the accessibility to this high-resolution data representing structural ensembles. We have trawled PDB files to extract structural features of residues with alternately located atoms. The output includes the distance between alternate conformations and identifies the location of these segments within the protein chain and in proximity of all other atoms within a defined radius. This dataset should be of use in efforts to predict multiple structures from a single sequence and support studies investigating protein flexibility and the association with protein function.

Codon-specific Ramachandran plots show amino acid backbone conformation depends on identity of the translated codon.

Synonymous codons translate into chemically identical amino acids. Once considered inconsequential to the formation of the protein product, there is evidence to suggest that codon usage affects co-translational protein folding and the final structure of the expressed protein. Here we develop a method for computing and comparing codon-specific Ramachandran plots and demonstrate that the backbone dihedral angle distributions of some synonymous codons are distinguishable with statistical significance for some secondary structures. This shows that there exists a dependence between codon identity and backbone torsion of the translated amino acid. Although these findings cannot pinpoint the causal direction of this dependence, we discuss the vast biological implications should coding be shown to directly shape protein conformation and demonstrate the usefulness of this method as a tool for probing associations between codon usage and protein structure. Finally, we urge for the inclusion of exact genetic information into structural databases.

Machine learning approaches demonstrate that protein structures carry information about their genetic coding.

Synonymous codons translate into the same amino acid. Although the identity of synonymous codons is often considered inconsequential to the final protein structure, there is mounting evidence for an association between the two. Our study examined this association using regression and classification models, finding that codon sequences predict protein backbone dihedral angles with a lower error than amino acid sequences, and that models trained with true dihedral angles have better classification of synonymous codons given structural information than models trained with random dihedral angles. Using this classification approach, we investigated local codon-codon dependencies and tested whether synonymous codon identity can be predicted more accurately from codon context than amino acid context alone, and most specifically which codon context position carries the most predictive power.

Opening the Schrödinger Box: Short- and Long-Range Mammalian Heart Rate Variability.

BACKGROUND: The interactions between the autonomic nervous system (ANS), intrinsic systems (e.g., endocrine), and internal pacemaker mechanisms govern short (milliseconds-seconds)- and long (seconds-minutes)-range heart rate variability (HRV). However, there is a debate regarding the identity of the mechanism underlying HRV on each time scale. We aim to design a general method that accurately differentiates between the relative contribution of the ANS and pacemaker mechanisms to HRV in various mammals, without the need for drug perturbations or organ isolation. Additionally, we aim to explore the universality of the relative contribution of the ANS and pacemaker system of different mammals. METHODS: This work explored short- and long-range HRVs using published ECG data from dogs, rabbits, and mice. To isolate the effects of ANS on HRV, ECG segments recorded before and after ANS-blockade were compared. RESULTS: Differentiation of the ANS from extrinsic and intrinsic pacemaker mechanisms was successfully achieved. In dogs, the internal pacemaker mechanisms were the main contributors to long-range and the ANS to short-range HRV. In rabbits and mice, the ANS and the internal pacemaker mechanisms affected both time scales, and anesthesia changed the relative contribution of the pacemaker mechanism to short- and long-range HRVs. In mice, the extrinsic mechanisms affected long-range HRV, while their effect was negligible in rabbits. CONCLUSION: We offer a novel approach to determine the relative contributions of ANS and extrinsic and intrinsic pacemaker mechanisms to HRV and highlight the importance of selecting mammalian research models with HRV mechanisms representative of the target species of interest.

Digital oximetry biomarkers for assessing respiratory function: standards of measurement, physiological interpretation, and clinical use.

Pulse oximetry is routinely used to non-invasively monitor oxygen saturation levels. A low oxygen level in the blood means low oxygen in the tissues, which can ultimately lead to organ failure. Yet, contrary to heart rate variability measures, a field which has seen the development of stable standards and advanced toolboxes and software, no such standards and open tools exist for continuous oxygen saturation time series variability analysis. The primary objective of this research was to identify, implement and validate key digital oximetry biomarkers (OBMs) for the purpose of creating a standard and associated reference toolbox for continuous oximetry time series analysis. We review the sleep medicine literature to identify clinically relevant OBMs. We implement these biomarkers and demonstrate their clinical value within the context of obstructive sleep apnea (OSA) diagnosis on a total of n = 3806 individual polysomnography recordings totaling 26,686 h of continuous data. A total of 44 digital oximetry biomarkers were implemented. Reference ranges for each biomarker are provided for individuals with mild, moderate, and severe OSA and for non-OSA recordings. Linear regression analysis between biomarkers and the apnea hypopnea index (AHI) showed a high correlation, which reached [Formula: see text]. The resulting python OBM toolbox, denoted "pobm", was contributed to the open software PhysioZoo ( physiozoo.org ). Studying the variability of the continuous oxygen saturation time series using pbom may provide information on the underlying physiological control systems and enhance our understanding of the manifestations and etiology of diseases, with emphasis on respiratory diseases.

Signatures of the autonomic nervous system and the heart's pacemaker cells in canine electrocardiograms and their applications to humans.

Heart rate and heart rate variability (HRV) are mainly determined by the autonomic nervous system (ANS), which interacts with receptors on the sinoatrial node (SAN; the heart's primary pacemaker), and by the "coupled-clock" system within the SAN cells. HRV changes are associated with cardiac diseases. However, the relative contributions of the ANS and SAN to HRV are not clear, impeding effective treatment. To discern the SAN's contribution, we performed HRV analysis on canine electrocardiograms containing basal and ANS-blockade segments. We also analyzed human electrocardiograms of atrial fibrillation and heart failure patients, as well as healthy aged subjects. Finally, we used a mathematical model to simulate HRV under decreased "coupled-clock" regulation. We found that (a) in canines, the SAN and ANS contribute mainly to long- and short-term HRV, respectively; (b) there is evidence suggesting a similar relative SAN contribution in humans; (c) SAN features can be calculated from beat-intervals obtained in-vivo, without intervention; (d) ANS contribution can be modeled by sines embedded in white noise; (e) HRV changes associated with cardiac diseases and aging can be interpreted as deterioration of both SAN and ANS; and (f) SAN clock-coupling can be estimated from changes in HRV. This may enable future non-invasive diagnostic applications.

A Universal Scaling Relation for Defining Power Spectral Bands in Mammalian Heart Rate Variability Analysis.

Background: Power spectral density (PSD) analysis of the heartbeat intervals in the three main frequency bands [very low frequency (VLF), low frequency (LF), and high frequency (HF)] provides a quantitative non-invasive tool for assessing the function of the cardiovascular control system. In humans, these frequency bands were standardized following years of empirical evidence. However, no quantitative approach has justified the frequency cutoffs of these bands and how they might be adapted to other mammals. Defining mammal-specific frequency bands is necessary if the PSD analysis of the HR is to be used as a proxy for measuring the autonomic nervous system activity in animal models. Methods: We first describe the distribution of prominent frequency peaks found in the normalized PSD of mammalian data using a Gaussian mixture model while assuming three components corresponding to the traditional VLF, LF and HF bands. We trained the algorithm on a database of human electrocardiogram recordings (n = 18) and validated it on databases of dogs (n = 17) and mice (n = 8). Finally, we tested it to predict the bands for rabbits (n = 4) for the first time. Results: Double-logarithmic analysis demonstrates a scaling law between the GMM-identified cutoff frequencies and the typical heart rate (HRm): fVLF-LF = 0.0037⋅ HRm0.58 , fLF-HF = 0.0017⋅ HRm1.01 and fHFup = 0.0128⋅ HRm0.86 . We found that the band cutoff frequencies and Gaussian mean scale with a power law of 1/4 or 1/8 of the typical body mass (BMm), thus revealing allometric power laws. Conclusion: Our automated data-driven approach allowed us to define the frequency bands in PSD analysis of beat-to-beat time series from different mammals. The scaling law between the band frequency cutoffs and the HRm can be used to approximate the PSD bands in other mammals.

PhysioZoo: A Novel Open Access Platform for Heart Rate Variability Analysis of Mammalian Electrocardiographic Data.

Background: The time variation between consecutive heartbeats is commonly referred to as heart rate variability (HRV). Loss of complexity in HRV has been documented in several cardiovascular diseases and has been associated with an increase in morbidity and mortality. However, the mechanisms that control HRV are not well understood. Animal experiments are the key to investigating this question. However, to date, there are no standard open source tools for HRV analysis of mammalian electrocardiogram (ECG) data and no centralized public databases for researchers to access. Methods: We created an open source software solution specifically designed for HRV analysis from ECG data of multiple mammals, including humans. We also created a set of public databases of mammalian ECG signals (dog, rabbit and mouse) with manually corrected R-peaks (>170,000 annotations) and signal quality annotations. The platform (software and databases) is called PhysioZoo. Results: PhysioZoo makes it possible to load ECG data and perform very accurate R-peak detection (F 1 > 98%). It also allows the user to manually correct the R-peak locations and annotate low signal quality of the underlying ECG. PhysioZoo implements state of the art HRV measures adapted for different mammals (dogs, rabbits, and mice) and allows easy export of all computed measures together with standard data representation figures. PhysioZoo provides databases and standard ranges for all HRV measures computed on healthy, conscious humans, dogs, rabbits, and mice at rest. Study of these measures across different mammals can provide new insights into the complexity of heart rate dynamics across species. Conclusion: PhysioZoo enables the standardization and reproducibility of HRV analysis in mammalian models through its open source code, freely available software, and open access databases. PhysioZoo will support and enable new investigations in mammalian HRV research. The source code and software are available on www.physiozoo.com.