Robust probabilistic superposition and comparison of protein structures

<p>Abstract</p> <p>Background</p> <p>Protein structure comparison is a central issue in structural bioinformatics. The standard dissimilarity measure for protein structures is the root mean square deviation (RMSD) of representative atom positions such as α-carbons. To evaluate the RMSD the structures under comparison must be superimposed optimally so as to minimize the RMSD. How to evaluate optimal fits becomes a matter of debate, if the structures contain regions which differ largely - a situation encountered in NMR ensembles and proteins undergoing large-scale conformational transitions.</p> <p>Results</p> <p>We present a probabilistic method for robust superposition and comparison of protein structures. Our method aims to identify the largest structurally invariant core. To do so, we model non-rigid displacements in protein structures with outlier-tolerant probability distributions. These distributions exhibit heavier tails than the Gaussian distribution underlying standard RMSD minimization and thus accommodate highly divergent structural regions. The drawback is that under a heavy-tailed model analytical expressions for the optimal superposition no longer exist. To circumvent this problem we work with a scale mixture representation, which implies a weighted RMSD. We develop two iterative procedures, an Expectation Maximization algorithm and a Gibbs sampler, to estimate the local weights, the optimal superposition, and the parameters of the heavy-tailed distribution. Applications demonstrate that heavy-tailed models capture differences between structures undergoing substantial conformational changes and can be used to assess the precision of NMR structures. By comparing Bayes factors we can automatically choose the most adequate model. Therefore our method is parameter-free.</p> <p>Conclusions</p> <p>Heavy-tailed distributions are well-suited to describe large-scale conformational differences in protein structures. A scale mixture representation facilitates the fitting of these distributions and enables outlier-tolerant superposition.</p

Atrial Heterogeneity Generates Re-entrant Substrate during Atrial Fibrillation and Anti-arrhythmic Drug Action: Mechanistic Insights from Canine Atrial Models

Anti-arrhythmic drug therapy is a frontline treatment for atrial fibrillation (AF), but its success rates are highly variable. This is due to incomplete understanding of the mechanisms of action of specific drugs on the atrial substrate at different stages of AF progression. We aimed to elucidate the role of cellular, tissue and organ level atrial heterogeneities in the generation of a re-entrant substrate during AF progression, and their modulation by the acute action of selected anti-arrhythmic drugs. To explore the complex cell-to-organ mechanisms, a detailed biophysical models of the entire 3D canine atria was developed. The model incorporated atrial geometry and fibre orientation from high-resolution micro-computed tomography, region-specific atrial cell electrophysiology and the effects of progressive AF-induced remodelling. The actions of multi-channel class III anti-arrhythmic agents vernakalant and amiodarone were introduced in the model by inhibiting appropriate ionic channel currents according to experimentally reported concentration-response relationships. AF was initiated by applied ectopic pacing in the pulmonary veins, which led to the generation of localized sustained re-entrant waves (rotors), followed by progressive wave breakdown and rotor multiplication in both atria. The simulated AF scenarios were in agreement with observations in canine models and patients. The 3D atrial simulations revealed that a re-entrant substrate was typically provided by tissue regions of high heterogeneity of action potential duration (APD). Amiodarone increased atrial APD and reduced APD heterogeneity and was more effective in terminating AF than vernakalant, which increased both APD and APD dispersion. In summary, the initiation and sustenance of rotors in AF is linked to atrial APD heterogeneity and APD reduction due to progressive remodelling. Our results suggest that anti-arrhythmic strategies that increase atrial APD without increasing its dispersion are effective in terminating AF

Understanding biomolecular motion, recognition, and allostery by use of conformational ensembles

We review the role conformational ensembles can play in the analysis of biomolecular dynamics, molecular recognition, and allostery. We introduce currently available methods for generating ensembles of biomolecules and illustrate their application with relevant examples from the literature. We show how, for binding, conformational ensembles provide a way of distinguishing the competing models of induced fit and conformational selection. For allostery we review the classic models and show how conformational ensembles can play a role in unravelling the intricate pathways of communication that enable allostery to occur. Finally, we discuss the limitations of conformational ensembles and highlight some potential applications for the future

The Changing Landscape for Stroke\ua0Prevention in AF: Findings From the GLORIA-AF Registry Phase 2

Background GLORIA-AF (Global Registry on Long-Term Oral Antithrombotic Treatment in Patients with Atrial Fibrillation) is a prospective, global registry program describing antithrombotic treatment patterns in patients with newly diagnosed nonvalvular atrial fibrillation at risk of stroke. Phase 2 began when dabigatran, the first non\u2013vitamin K antagonist oral anticoagulant (NOAC), became available. Objectives This study sought to describe phase 2 baseline data and compare these with the pre-NOAC era collected during phase&nbsp;1. Methods During phase 2, 15,641 consenting patients were enrolled (November 2011 to December 2014); 15,092 were eligible. This pre-specified cross-sectional analysis describes eligible patients\u2019 baseline characteristics. Atrial fibrillation&nbsp;disease characteristics, medical outcomes, and concomitant diseases and medications were collected. Data were analyzed using descriptive statistics. Results Of the total patients, 45.5% were female; median age was 71 (interquartile range: 64, 78) years. Patients were from Europe (47.1%), North America (22.5%), Asia (20.3%), Latin America (6.0%), and the Middle East/Africa (4.0%). Most had high stroke risk (CHA2DS2-VASc [Congestive heart failure, Hypertension, Age&nbsp; 6575 years, Diabetes mellitus, previous Stroke, Vascular disease, Age 65 to 74 years, Sex category] score&nbsp; 652; 86.1%); 13.9% had moderate risk (CHA2DS2-VASc&nbsp;= 1). Overall, 79.9% received oral anticoagulants, of whom 47.6% received NOAC and 32.3% vitamin K antagonists (VKA); 12.1% received antiplatelet agents; 7.8% received no antithrombotic treatment. For comparison, the proportion of phase 1 patients (of N&nbsp;= 1,063 all eligible) prescribed VKA was 32.8%, acetylsalicylic acid 41.7%, and no therapy 20.2%. In Europe in phase 2, treatment with NOAC was more common than VKA (52.3% and 37.8%, respectively); 6.0% of patients received antiplatelet treatment; and 3.8% received no antithrombotic treatment. In North America, 52.1%, 26.2%, and 14.0% of patients received NOAC, VKA, and antiplatelet drugs, respectively; 7.5% received no antithrombotic treatment. NOAC use was less common in Asia (27.7%), where 27.5% of patients received VKA, 25.0% antiplatelet drugs, and 19.8% no antithrombotic treatment. Conclusions The baseline data from GLORIA-AF phase 2 demonstrate that in newly diagnosed nonvalvular atrial fibrillation patients, NOAC have been highly adopted into practice, becoming more frequently prescribed than VKA in&nbsp;Europe and North America. Worldwide, however, a large proportion of patients remain undertreated, particularly in&nbsp;Asia&nbsp;and North America. (Global Registry on Long-Term Oral Antithrombotic Treatment in Patients With Atrial Fibrillation [GLORIA-AF]; NCT01468701

Structure and Dynamics of Biological Systems: Integration of Neutron Scattering with Computer Simulation

The combination of molecular dynamics simulation and neutron scattering techniques has emerged as a highly synergistic approach to elucidate the atomistic details of the structure, dynamics and functions of biological systems. Simulation models can be tested by calculating neutron scattering structure factors and comparing the results directly with experiments. If the scattering profiles agree the simulations can be used to provide a detailed decomposition and interpretation of the experiments, and if not, the models can be rationally adjusted. Comparison with neutron experiment can be made at the level of the scattering functions or, less directly, of structural and dynamical quantities derived from them. Here, we examine the combination of simulation and experiment in the interpretation of SANS and inelastic scattering experiments on the structure and dynamics of proteins and other biopolymers

Efficient Translation of Pelargonium line pattern virus RNAs Relies on a TED-Like 3 '-Translational Enhancer that Communicates with the Corresponding 5 '-Region through a Long-Distance RNA-RNA Interaction

[EN] Cap-independent translational enhancers (CITEs) have been identified at the 3'-terminal regions of distinct plant positive-strand RNA viruses belonging to families Tombusviridae and Luteoviridae. On the bases of their structural and/or functional requirements, at least six classes of CITEs have been defined whose distribution does not correlate with taxonomy. The so-called TED class has been relatively under-studied and its functionality only confirmed in the case of Satellite tobacco necrosis virus, a parasitic subviral agent. The 3' untranslated region of the monopartite genome of Pelargonium line pattern virus (PLPV), the recommended type member of a tentative new genus (Pelarspovirus) in the family Tombusviridae, was predicted to contain a TED-like CITE. Similar CITEs can be anticipated in some other related viruses though none has been experimentally verified. Here, in the first place, we have performed a reassessment of the structure of the putative PLPV-TED through in silico predictions and in vitro SHAPE analysis with the full-length PLPV genome, which has indicated that the presumed TED element is larger than previously proposed. The extended conformation of the TED is strongly supported by the pattern of natural sequence variation, thus providing comparative structural evidence in support of the structural data obtained by in silico and in vitro approaches. Next, we have obtained experimental evidence demonstrating the in vivo activity of the PLPV-TED in the genomic (g) RNA, and also in the subgenomic (sg) RNA that the virus produces to express 3'-proximal genes. Besides other structural features, the results have highlighted the key role of long-distance kissing-loop interactions between the 3'-CITE and 5'-proximal hairpins for gRNA and sgRNA translation. Bioassays of CITE mutants have confirmed the importance of the identified 5'-3' RNA communication for viral infectivity and, moreover, have underlined the strong evolutionary constraints that may operate on genome stretches with both regulatory and coding functions.This work was supported by grants BFU2009-11699 and BFU2012-36095 from the Ministerio de Investigacion, Ciencia e Innovacion (MICINN, Spain, www.micinn.es) and the Ministerio de Economia y Competitividad (MINECO, Spain, http://www.mineco.gob.es), respectively, and ACOMP/2012/100 from the Generalitat Valenciana (http://www.gva.es) (to C.H.). MBP and LR were the recipients of a predoctoral and postdoctoral (Juan de la Cierva program) contract, respectively, from MICINN, and MPC was the recipient of a predoctoral contract from MINECO.Blanco Pérez, M.; Pérez Cañamás, M.; Ruiz, L.; Hernandez Fort, C. (2016). Efficient Translation of Pelargonium line pattern virus RNAs Relies on a TED-Like 3 '-Translational Enhancer that Communicates with the Corresponding 5 '-Region through a Long-Distance RNA-RNA Interaction. PLoS ONE. 11(4):1-24. https://doi.org/10.1371/journal.pone.0152593S12411

Quaternions as a tool for the analysis of molecular systems

International audienceQuaternions are generalized complex numbers and represent rotations in space as ordinary complex numbers represent rotations in a plane. In the context of molecular dynamics (MD) simulations they have been ‘rediscovered’ for the integration of the rotational equation of motion of rigid molecules since they allow one to write down these equations in a singularity-free form. In this paper applications to the analysis of molecular systems are described

Quaternions as a tool for the analysis of molecular systems

Quaternions are generalized complex numbers and represent rotations in space as ordinary complex numbers represent rotations in a plane. In the context of molecular dynamics (MD) simulations they have been ‘rediscovered’ for the integration of the rotational equation of motion of rigid molecules since they allow one to write down these equations in a singularity-free form. In this paper applications to the analysis of molecular systems are described

Liquid-like and solid-like motions in proteins

Recent analyses of molecular dynamics simulations of hydrated C-phycocyanin suggest that the internal single-particle dynamics of this protein can be decomposed into four almost decoupled motion types: (1) diffusion of residues ("beads") in an effective harmonic potential, (2) corresponding vibrations in a local potential well, (3) purely rotational rigid side-chain diffusion, and (4) residue deformations. Each residue bead is represented by the corresponding C[α] carbon atom on the main chain. The effective harmonic residue potential can be imagined as the envelope of many local wells which are separated by small energy barriers. The residue friction matrix is assumed diagonal and the individual friction constants can be related to the density of the surrounding atoms. In this article we show that our model can be applied to lysozyme in solution as well, the only difference being that the side-chain deformations are more important and seem to be strongly correlated with the side-chain rotations. Comparing the simulated coherent scattering function of C-phycocyanin to a neutron spin-echo spectrum we show that our model can also describe collective motions in proteins at the residue level