Ensemble reweighting using Cryo-EM particles

Cryo-electron microscopy (cryo-EM) has recently become a premier method for obtaining high-resolution structures of biological macromolecules. However, it is limited to biomolecular samples with low conformational heterogeneity, where all the conformations can be well-sampled at many projection angles. While cryo-EM technically provides single-molecule data for heterogeneous molecules, most existing reconstruction tools cannot extract the full distribution of possible molecular configurations. To overcome these limitations, we build on a prior Bayesian approach and develop an ensemble refinement framework that estimates the ensemble density from a set of cryo-EM particles by reweighting a prior ensemble of conformations, e.g., from molecular dynamics simulations or structure prediction tools. Our work is a general approach to recovering the equilibrium probability density of the biomolecule directly in conformational space from single-molecule data. To validate the framework, we study the extraction of state populations and free energies for a simple toy model and from synthetic cryo-EM images of a simulated protein that explores multiple folded and unfolded conformations

Supplementary guidance: listening to staff: Autumn 2017

Kinases play a critical role in cellular signaling and are dysregulated in a number of diseases, such as cancer, diabetes, and neurodegeneration. Therapeutics targeting kinases currently account for roughly 50% of cancer drug discovery efforts. The ability to explore human kinase biochemistry and biophysics in the laboratory is essential to designing selective inhibitors and studying drug resistance. Bacterial expression systems are superior to insect or mammalian cells in terms of simplicity and cost effectiveness but have historically struggled with human kinase expression. Following the discovery that phosphatase coexpression produced high yields of Src and Abl kinase domains in bacteria, we have generated a library of 52 His-tagged human kinase domain constructs that express above 2 μg/mL of culture in an automated bacterial expression system utilizing phosphatase coexpression (YopH for Tyr kinases and lambda for Ser/Thr kinases). Here, we report a structural bioinformatics approach to identifying kinase domain constructs previously expressed in bacteria and likely to express well in our protocol, experiments demonstrating our simple construct selection strategy selects constructs with good expression yields in a test of 84 potential kinase domain boundaries for Abl, and yields from a high-throughput expression screen of 96 human kinase constructs. Using a fluorescence-based thermostability assay and a fluorescent ATP-competitive inhibitor, we show that the highest-expressing kinases are folded and have well-formed ATP binding sites. We also demonstrate that these constructs can enable characterization of clinical mutations by expressing a panel of 48 Src and 46 Abl mutations. The wild-type kinase construct library is available publicly via Addgene

BHPR research: qualitative1. Complex reasoning determines patients' perception of outcome following foot surgery in rheumatoid arhtritis

Background: Foot surgery is common in patients with RA but research into surgical outcomes is limited and conceptually flawed as current outcome measures lack face validity: to date no one has asked patients what is important to them. This study aimed to determine which factors are important to patients when evaluating the success of foot surgery in RA Methods: Semi structured interviews of RA patients who had undergone foot surgery were conducted and transcribed verbatim. Thematic analysis of interviews was conducted to explore issues that were important to patients. Results: 11 RA patients (9 ♂, mean age 59, dis dur = 22yrs, mean of 3 yrs post op) with mixed experiences of foot surgery were interviewed. Patients interpreted outcome in respect to a multitude of factors, frequently positive change in one aspect contrasted with negative opinions about another. Overall, four major themes emerged. Function: Functional ability & participation in valued activities were very important to patients. Walking ability was a key concern but patients interpreted levels of activity in light of other aspects of their disease, reflecting on change in functional ability more than overall level. Positive feelings of improved mobility were often moderated by negative self perception ("I mean, I still walk like a waddling duck”). Appearance: Appearance was important to almost all patients but perhaps the most complex theme of all. Physical appearance, foot shape, and footwear were closely interlinked, yet patients saw these as distinct separate concepts. Patients need to legitimize these feelings was clear and they frequently entered into a defensive repertoire ("it's not cosmetic surgery; it's something that's more important than that, you know?”). Clinician opinion: Surgeons' post operative evaluation of the procedure was very influential. The impact of this appraisal continued to affect patients' lasting impression irrespective of how the outcome compared to their initial goals ("when he'd done it ... he said that hasn't worked as good as he'd wanted to ... but the pain has gone”). Pain: Whilst pain was important to almost all patients, it appeared to be less important than the other themes. Pain was predominately raised when it influenced other themes, such as function; many still felt the need to legitimize their foot pain in order for health professionals to take it seriously ("in the end I went to my GP because it had happened a few times and I went to an orthopaedic surgeon who was quite dismissive of it, it was like what are you complaining about”). Conclusions: Patients interpret the outcome of foot surgery using a multitude of interrelated factors, particularly functional ability, appearance and surgeons' appraisal of the procedure. While pain was often noted, this appeared less important than other factors in the overall outcome of the surgery. Future research into foot surgery should incorporate the complexity of how patients determine their outcome Disclosure statement: All authors have declared no conflicts of interes

DDR1

Folding@home simulations of Discoidin Domain Receptor Tyrosine Kinase 1 (DDR1

Structural, biochemical and computational studies of TRP channel transmembrane domain modularity

Transient receptor potential (TRP) channels are expressed throughout the central nervous system and have a unique ability to detect a wide range of stimuli including changes in voltage, temperature, pH, lipid environment, small molecule agonists, and mechanical stress. While it is known that TRP channels contain the same six transmembrane helix (S1-S6), tetrameric architecture as voltage-gated channels, the degree to which functional and structural analogies are relevant remains poorly understood. This thesis describes a multidisciplinary approach toward understanding the structure and function of TRP channel transmembrane domains by focusing on the S1-S4 transmembrane helices of the TRPV1. This focus is inspired by the voltage-sensor domain (VSD) of the S1-S4 helices of voltage-gated channels, for which a range of studies show functional and structural independence. While some TRP channels are voltage-sensitive, their S4 helix does not contain the positive string of amino acids of canonical VSDs. However, the S1-S4 helices are functionally significant as the binding site of small molecule ligands in both TRPV1 and TRPM8 (for capsaicin and menthol, respectively). The question of TRP channel transmembrane domain modularity is addressed in this thesis by expression and purification trials as well as radioligand-binding assays. It is demonstrated that the S1-S4 and S1-S6 helices of TRPV1 can be properly inserted, overexpressed, and show signs of stability upon detergent-extraction from Saccharomyces cerevisiae membranes. However the TRPV1 S1-S4 and S1-S6 helices do not show wildtype (WT)-like binding in [3H]-RTX binding assays. These results indicate that the TRPV1 transmembrane domains are likely structural but not functional domains. The S. cerevisiae expression system remained promising for the overexpression of TRP transmembrane domains as well as the production of functional, though not stable upon detergent-extraction, WT TRPV1. This WT TRPV1 was subsequently found to functionally bind both RTX, used in ligand binding assays, as well as the double-knot toxin (DkTx), targeted to the pore domain (the S5-S6 helices). An effect of DkTx on RTX binding affinity demonstrates an allosteric interaction indicative of a possible tighter packing between the two transmembrane domains than is seen in voltage-gated channels containing the canonical VSD. Computational approaches additionally allowed for the investigation of the intramembrane capsaicin binding site in the TRPV1 S1-S4 helices, crucial to the initial motivations of this study. While the literature locates the capsaicin binding site to the TRPV1 S1-S4 helices, a `binding pocket' has yet to be defined, with regards to the orientation of bound capsaicin and its access route to the site via the bilayer. Using molecular dynamics (MD) simulations the preferred location of capsaicin within the bilayer is defined, as well as the elucidiation of capsaicin flip-flop between bilayer leaflets as a key event prior to TRPV1 binding. A transient binding was also observed between a homology model of the TRPV1 S1-S4 helices and capsaicin, possibly encouraging the idea that the S1-S4 helices still contain a partial binding site, though of too low affinity to be observed in the binding experiments performed here.</p

Structural, biochemical and computational studies of TRP channel transmembrane domain modularity

Transient receptor potential (TRP) channels are expressed throughout the central nervous system and have a unique ability to detect a wide range of stimuli including changes in voltage, temperature, pH, lipid environment, small molecule agonists, and mechanical stress. While it is known that TRP channels contain the same six transmembrane helix (S1-S6), tetrameric architecture as voltage-gated channels, the degree to which functional and structural analogies are relevant remains poorly understood. This thesis describes a multidisciplinary approach toward understanding the structure and function of TRP channel transmembrane domains by focusing on the S1-S4 transmembrane helices of the TRPV1. This focus is inspired by the voltage-sensor domain (VSD) of the S1-S4 helices of voltage-gated channels, for which a range of studies show functional and structural independence. While some TRP channels are voltage-sensitive, their S4 helix does not contain the positive string of amino acids of canonical VSDs. However, the S1-S4 helices are functionally significant as the binding site of small molecule ligands in both TRPV1 and TRPM8 (for capsaicin and menthol, respectively). The question of TRP channel transmembrane domain modularity is addressed in this thesis by expression and purification trials as well as radioligand-binding assays. It is demonstrated that the S1-S4 and S1-S6 helices of TRPV1 can be properly inserted, overexpressed, and show signs of stability upon detergent-extraction from Saccharomyces cerevisiae membranes. However the TRPV1 S1-S4 and S1-S6 helices do not show wildtype (WT)-like binding in [3H]-RTX binding assays. These results indicate that the TRPV1 transmembrane domains are likely structural but not functional domains. The S. cerevisiae expression system remained promising for the overexpression of TRP transmembrane domains as well as the production of functional, though not stable upon detergent-extraction, WT TRPV1. This WT TRPV1 was subsequently found to functionally bind both RTX, used in ligand binding assays, as well as the double-knot toxin (DkTx), targeted to the pore domain (the S5-S6 helices). An effect of DkTx on RTX binding affinity demonstrates an allosteric interaction indicative of a possible tighter packing between the two transmembrane domains than is seen in voltage-gated channels containing the canonical VSD. Computational approaches additionally allowed for the investigation of the intramembrane capsaicin binding site in the TRPV1 S1-S4 helices, crucial to the initial motivations of this study. While the literature locates the capsaicin binding site to the TRPV1 S1-S4 helices, a `binding pocket' has yet to be defined, with regards to the orientation of bound capsaicin and its access route to the site via the bilayer. Using molecular dynamics (MD) simulations the preferred location of capsaicin within the bilayer is defined, as well as the elucidiation of capsaicin flip-flop between bilayer leaflets as a key event prior to TRPV1 binding. A transient binding was also observed between a homology model of the TRPV1 S1-S4 helices and capsaicin, possibly encouraging the idea that the S1-S4 helices still contain a partial binding site, though of too low affinity to be observed in the binding experiments performed here

Exploring the Adjugate Matrix Approach to Quaternion Pose Extraction

Quaternions are important for a wide variety of rotation-related problems in computer graphics, machine vision, and robotics. We study the nontrivial geometry of the relationship between quaternions and rotation matrices by exploiting the adjugate matrix of the characteristic equation of a related eigenvalue problem to obtain the manifold of the space of a quaternion eigenvector. We argue that quaternions parameterized by their corresponding rotation matrices cannot be expressed, for example, in machine learning tasks, as single-valued functions: the quaternion solution must instead be treated as a manifold, with different algebraic solutions for each of several single-valued sectors represented by the adjugate matrix. We conclude with novel constructions exploiting the quaternion adjugate variables to revisit several classic pose estimation applications: 2D point-cloud matching, 2D point-cloud-to-projection matching, 3D point-cloud matching, 3D orthographic point-cloud-to-projection matching, and 3D perspective point-cloud-to-projection matching. We find an exact solution to the 3D orthographic least squares pose extraction problem, and apply it successfully also to the perspective pose extraction problem with results that improve on existing methods.Comment: 67 pages, 5 appendices, 9 figure

Extracting thermodynamic properties from van 't Hoff plots with emphasis on temperature-sensing ion channels

Transient receptor potential (TRP) ion channels are among the most well-studied classes of temperature-sensing molecules. Yet, the molecular mechanism and thermodynamic basis for the temperature sensitivity of TRP channels remains to this day poorly understood. One hypothesis is that the temperature-sensing mechanism can simply be described by a difference in heat capacity between the closed and open channel states. While such a two-state model may be simplistic it nonetheless has descriptive value, in the sense that it can be used to to compare overall temperature sensitivity between different channels and mutants. Here, we introduce a mathematical framework based on the two-state model to reliably extract temperature-dependent thermodynamic potentials and heat capacities from measurements of equilibrium constants at different temperatures. Our framework is implemented in an open-source data analysis package that provides a straightforward way to fit both linear and nonlinear van ‘t Hoff plots, thus avoiding some of the previous, potentially erroneous, assumptions when extracting thermodynamic variables from TRP channel electrophysiology data

The Inaugural Flatiron Institute Cryo-EM Heterogeneity Community Challenge

Anonymised submissions to the Flatiron Institute 2023 Heterogeneity Challenge

Microsecond molecular dynamics simulation of kinase domain of the human tyrosine kinase ABL1

This is a one microsecond trajectory of the catalytic domain of the human tyrosine protein kinase ABL1 (UNIPROT:ABL1_HUMAN, residues 242-493). The simulation was generated on Folding@home using OpenMM 6.3.1 with the AMBER99SB-ILDN forcefield, TIP3P water, and neutralizing salt. A Langevin integrator with 2 fs timestep, 300 K temperature, and 1/ps collision rate was used. PME with a nonbonded cutoff of 1.0 nm was used with default tolerance of 5.0e-4. A long-range isotropic dispersion correction was used to correct for the truncation of Lennard-Jones interactions at 1.0 nm. A molecular-scaling Monte Carlo barostat was applied using a 1 atm reference pressure, with MC moves attempted every 50 steps. Coordinates were saved every 250ps, and waters and ions stripped. The trajectory file contains 4000 snapshots. Initial models were generated with Ensembler (see References) using the ABL1 structure from chain B of PDB ID: 2GQG as a template.<div><br></div><div>Manifest:<div>* run17-clone3.jpg - image of initial simulation snapshot</div><div>* run17-clone3.pdb - PDB file of initial simulation snapshot</div><div>* run17-clone3.h5 - MDTraj HDF5 format trajectory</div><div><br></div><div>To convert the trajectory, see the MDTraj `mdconvert` command:</div><div>http://mdtraj.org/latest/mdconvert.html</div><div><br></div><div>For information on installing MDTraj, see:</div><div>http://mdtraj.org/</div></div