Abstract Recent work has introduced a new fluorescent voltage sensor, ASAP1, which can monitor rapid trains of action potentials in cultured neurons. This indicator is based on the Gallus gallus voltage-sensitive phosphatase with the phosphatase domain removed and a circularly permuted GFP placed in the S3-S4 linker. However, many of the biophysical details of this indicator remain unknown. In this work, we study the biophysical properties of ASAP1. Using the cut-open voltage clamp technique, we have simultaneously recorded fluorescence signals and gating currents from Xenopus laevis oocytes expressing ASAP1. Gating charge movement and fluorescence kinetics track closely with each other, although ASAP1 gating currents are significantly faster than those of Ciona intestinalis voltage-sensitive phosphatase. Altering the residue before the first gating charge removes a split in the ASAP1 QV curve, but preserves the accelerated kinetics that allow for the faithful tracking of action potentials in neurons.

Abstract Integrins are bidirectional, allosteric transmembrane receptors that play a central role in hemostasis and arterial thrombosis. Using cryo-electron microscopy, multireference single-particle reconstruction methods, and statistics-based computational fitting approaches, we determined three-dimensional structures of human integrin αIIbβ3 embedded in a lipid bilayer (nanodiscs) while bound to domains of the cytosolic regulator talin and to extracellular ligands. We also determined the conformations of integrin in solution by itself to localize the membrane and the talin-binding site. To our knowledge, our data provide unprecedented three-dimensional information about the conformational states of intact, full-length integrin within membrane bilayers under near-physiological conditions and in the presence of cytosolic activators and extracellular ligands. We show that αIIbβ3 integrins exist in a conformational equilibrium clustered around four main states. These conformations range from a compact bent nodule to two partially extended intermediate conformers and finally to a fully upright state. In the presence of nanodiscs and the two ligands, the equilibrium is significantly shifted toward the upright conformation. In this conformation, the receptor extends ∼20 nm upward from the membrane. There are no observable contacts between the two subunits other than those in the headpiece near the ligand-binding pocket, and the α- and β-subunits are well separated with their cytoplasmic tails ∼8 nm apart. Our results indicate that extension of the ectodomain is possible without separating the legs or extending the hybrid domain, and that the ligand-binding pocket is not occluded by the membrane in any conformations of the equilibrium. Further, they suggest that integrin activation may be influenced by equilibrium shifts.

Abstract Prediction of gene expression levels from regulatory sequences is one of the major challenges of genomic biology today. A particularly promising approach to this problem is that taken by thermodynamics-based models that interpret an enhancer sequence in a given cellular context specified by transcription factor concentration levels and predict precise expression levels driven by that enhancer. Such models have so far not accounted for the effect of chromatin accessibility on interactions between transcription factor and DNA and consequently on gene-expression levels. Here, we extend a thermodynamics-based model of gene expression, called GEMSTAT (Gene Expression Modeling Based on Statistical Thermodynamics), to incorporate chromatin accessibility data and quantify its effect on accuracy of expression prediction. In the new model, called GEMSTAT-A, accessibility at a binding site is assumed to affect the transcription factor’s binding strength at the site, whereas all other aspects are identical to the GEMSTAT model. We show that this modification results in significantly better fits in a data set of over 30 enhancers regulating spatial expression patterns in the blastoderm-stage Drosophila embryo. It is important to note that the improved fits result not from an overall elevated accessibility in active enhancers but from the variation of accessibility levels within an enhancer. With whole-genome DNA accessibility measurements becoming increasingly popular, our work demonstrates how such data may be useful for sequence-to-expression models. It also calls for future advances in modeling accessibility levels from sequence and the transregulatory context, so as to predict accurately the effect of cis and trans perturbations on gene expression.

Abstract Denaturant-induced unfolding of helical membrane proteins provides insights into their mechanism of folding and domain organization, which take place in the chemically heterogeneous, anisotropic environment of a lipid membrane. Rhomboid proteases are intramembrane proteases that play key roles in various diseases. Crystal structures have revealed a compact helical bundle with a buried active site, which requires conformational changes for the cleavage of transmembrane substrates. A dimeric form of the rhomboid protease has been shown to be important for activity. In this study, we examine the mechanism of refolding for two distinct rhomboids to gain insight into their secondary structure-activity relationships. Although helicity is largely abolished in the unfolded states of both proteins, unfolding is completely reversible for HiGlpG but only partially reversible for PsAarA. Refolding of both proteins results in reassociation of the dimer, with a 90% regain of catalytic activity for HiGlpG but only a 70% regain for PsAarA. For both proteins, a broad, gradual transition from the native, folded state to the denatured, partly unfolded state was revealed with the aid of circular dichroism spectroscopy as a function of denaturant concentration, thus arguing against a classical two-state model as found for many globular soluble proteins. Thermal denaturation has irreversible destabilizing effects on both proteins, yet reveals important functional details regarding substrate accessibility to the buried active site. This concerted biophysical and functional analysis demonstrates that HiGlpG, with a simple six-transmembrane-segment organization, is more robust than PsAarA, which has seven predicted transmembrane segments, thus rendering HiGlpG amenable to in vitro studies of membrane-protein folding.

Abstract Glucosylceramide (GlcCer), one of the simplest glycosphingolipids, plays key roles in physiology and pathophysiology. It has been suggested that GlcCer modulates cellular events by forming specialized domains. In this study, we investigated the interplay between GlcCer and cholesterol (Chol), an important lipid involved in the formation of liquid-ordered (lo) phases. Using fluorescence microscopy and spectroscopy, and dynamic and electrophoretic light scattering, we characterized the interaction between these lipids in different pH environments. A quantitative description of the phase behavior of the ternary unsaturated phospholipid/Chol/GlcCer mixture is presented. The results demonstrate coexistence between lo and liquid-disordered (ld) phases. However, the extent of lo/ld phase separation is sparse, mainly due to the ability of GlcCer to segregate into tightly packed gel domains. As a result, the phase diagram of these mixtures is characterized by an extensive three-phase coexistence region of fluid (ld-phospholipid enriched)/lo (Chol enriched)/gel (GlcCer enriched). Moreover, the results show that upon acidification, GlcCer solubility in the lo phase is increased, leading to a larger lo/ld coexistence region. Quantitative analyses allowed us to determine the differences in the composition of the phases at neutral and acidic pH. These results predict the impact of GlcCer on domain formation and membrane organization in complex biological membranes, and provide a background for unraveling the relationship between the biophysical properties of GlcCer and its biological action.

Abstract Cu+-specific P-type ATPase membrane protein transporters regulate cellular copper levels. The lack of crystal structures in Cu+-binding states has limited our understanding of how ion entry and binding are achieved. Here, we characterize the molecular basis of Cu+ entry using molecular-dynamics simulations, structural modeling, and in vitro and in vivo functional assays. Protein structural rearrangements resulting in the exposure of positive charges to bulk solvent rather than to lipid phosphates indicate a direct molecular role of the putative docking platform in Cu+ delivery. Mutational analyses and simulations in the presence and absence of Cu+ predict that the ion-entry path involves two ion-binding sites: one transient Met148-Cys382 site and one intramembranous site formed by trigonal coordination to Cys384, Asn689, and Met717. The results reconcile earlier biochemical and x-ray absorption data and provide a molecular understanding of ion entry in Cu+-transporting P-type ATPases.

Abstract We describe the ways in which Galaxy, a web-based reproducible research platform, can be used for web-based sharing of complex computational models. Galaxy allows users to seamlessly customize and run simulations on cloud computing resources, a concept we refer to as Models and Simulations as a Service (MaSS). To illustrate this application of Galaxy, we have developed a tool suite for simulating a high spatial-resolution model of the cardiac Ca2+ spark that requires supercomputing resources for execution. We also present tools for simulating models encoded in the SBML and CellML model description languages, thus demonstrating how Galaxy’s reproducible research features can be leveraged by existing technologies. Finally, we demonstrate how the Galaxy workflow editor can be used to compose integrative models from constituent submodules. This work represents an important novel approach, to our knowledge, to making computational simulations more accessible to the broader scientific community.

Abstract Tightly clustered inositol trisphosphate receptors (IP3Rs) control localized Ca2+ liberation from the endoplasmic reticulum to generate repetitive Ca2+ puffs. Distributions of the interpuff interval (IPI), i.e., the waiting time between successive puffs, are found to be well characterized by a probability density function involving only two parameters, λ and ξ, which represent the basal rate of puff generation and the recovery rate from refractoriness, respectively. However, how the two parameters depend on the kinetic parameters of single IP3Rs in a cluster is still unclear. In this article, using a stochastic puff model and a single-channel data-based IP3R model, we establish the dependencies of λ and ξ on two important IP3R model parameters, IP3 concentration ([IP3]) and the recovery rate from Ca2+ inhibition (rlow). By varying [IP3] and rlow in physiologically plausible ranges, we find that the ξ-λ plane is comprised of only two disjoint regions, a biologically impermissible region and a region where each parameter set (ξ, λ) can be caused by using two different combinations of [IP3] and rlow. The two combinations utilize very different mechanisms to maintain the same IPI distribution, and the mechanistic difference provides a way of identifying IP3R kinetic parameters by observing properties of the IPI.

Abstract Changing the helical propensity of a polypeptide sequence might be expected to affect the conformational properties of the denatured state of a protein. To test this hypothesis, alanines at positions 83 and 87 near the center of helix 3 of cytochrome c′ from Rhodopseudomonas palustris were mutated to serine to decrease the stability of this helix. A set of 13 single histidine variants in the A83S/A87S background were prepared to permit assessment of the conformational properties of the denatured state using histidine-loop formation in 3 M guanidine hydrochloride. The data are compared with previous histidine-heme loop formation data for wild-type cytochrome c′. As expected, destabilization of helix 3 decreases the global stabilities of the histidine variants in the A83S/A87S background relative to the wild-type background. Loop stability versus loop size data yields a scaling exponent of 2.1 ± 0.2, similar to the value of 2.3 ± 0.2 obtained for wild-type cytochrome c′. However, the stabilities of all histidine-heme loops, which contain the helix 3 sequence segment, are increased in the A83S/A87S background compared to the wild-type background. Rate constants for histidine-heme loop breakage are similar for the wild-type and A83S/A87S variants. However, for histidine-heme loops that contain the helix 3 sequence segment, the rate constants for loop formation increase in the A83S/A87S background compared to the wild-type background. Thus, residual helical structure appears to stiffen the polypeptide chain slowing loop formation in the denatured state. The implications of these results for protein folding mechanisms are discussed.

Abstract Biological mixtures such as the cytosol may consist of thousands of distinct components. There is now a substantial body of evidence showing that, under physiological conditions, intracellular mixtures can phase separate into spatially distinct regions with differing compositions. In this article we present numerical evidence indicating that such spontaneous compartmentalization exploits general features of the phase diagram of a multicomponent biomolecular mixture. In particular, we show that demixed domains are likely to segregate when the variance in the intermolecular interaction strengths exceeds a well-defined threshold. Multiple distinct phases are likely to become stable under very similar conditions, which can then be tuned to achieve multiphase coexistence. As a result, only minor adjustments to the composition of the cytosol or the strengths of the intermolecular interactions are needed to regulate the formation of different domains with specific compositions, implying that phase separation is a robust mechanism for creating spatial organization. We further predict that this functionality is only weakly affected by increasing the number of components in the system. Our model therefore suggests that, for purely physico-chemical reasons, biological mixtures are naturally poised to undergo a small number of demixing phase transitions.

Abstract Recent work has introduced a new fluorescent voltage sensor, ASAP1, which can monitor rapid trains of action potentials in cultured neurons. This indicator is based on the Gallus gallus voltage-sensitive phosphatase with the phosphatase domain removed and a circularly permuted GFP placed in the S3-S4 linker. However, many of the biophysical details of this indicator remain unknown. In this work, we study the biophysical properties of ASAP1. Using the cut-open voltage clamp technique, we have simultaneously recorded fluorescence signals and gating currents from Xenopus laevis oocytes expressing ASAP1. Gating charge movement and fluorescence kinetics track closely with each other, although ASAP1 gating currents are significantly faster than those of Ciona intestinalis voltage-sensitive phosphatase. Altering the residue before the first gating charge removes a split in the ASAP1 QV curve, but preserves the accelerated kinetics that allow for the faithful tracking of action potentials in neurons.

Abstract Integrins are bidirectional, allosteric transmembrane receptors that play a central role in hemostasis and arterial thrombosis. Using cryo-electron microscopy, multireference single-particle reconstruction methods, and statistics-based computational fitting approaches, we determined three-dimensional structures of human integrin αIIbβ3 embedded in a lipid bilayer (nanodiscs) while bound to domains of the cytosolic regulator talin and to extracellular ligands. We also determined the conformations of integrin in solution by itself to localize the membrane and the talin-binding site. To our knowledge, our data provide unprecedented three-dimensional information about the conformational states of intact, full-length integrin within membrane bilayers under near-physiological conditions and in the presence of cytosolic activators and extracellular ligands. We show that αIIbβ3 integrins exist in a conformational equilibrium clustered around four main states. These conformations range from a compact bent nodule to two partially extended intermediate conformers and finally to a fully upright state. In the presence of nanodiscs and the two ligands, the equilibrium is significantly shifted toward the upright conformation. In this conformation, the receptor extends ∼20 nm upward from the membrane. There are no observable contacts between the two subunits other than those in the headpiece near the ligand-binding pocket, and the α- and β-subunits are well separated with their cytoplasmic tails ∼8 nm apart. Our results indicate that extension of the ectodomain is possible without separating the legs or extending the hybrid domain, and that the ligand-binding pocket is not occluded by the membrane in any conformations of the equilibrium. Further, they suggest that integrin activation may be influenced by equilibrium shifts.

Abstract Prediction of gene expression levels from regulatory sequences is one of the major challenges of genomic biology today. A particularly promising approach to this problem is that taken by thermodynamics-based models that interpret an enhancer sequence in a given cellular context specified by transcription factor concentration levels and predict precise expression levels driven by that enhancer. Such models have so far not accounted for the effect of chromatin accessibility on interactions between transcription factor and DNA and consequently on gene-expression levels. Here, we extend a thermodynamics-based model of gene expression, called GEMSTAT (Gene Expression Modeling Based on Statistical Thermodynamics), to incorporate chromatin accessibility data and quantify its effect on accuracy of expression prediction. In the new model, called GEMSTAT-A, accessibility at a binding site is assumed to affect the transcription factor’s binding strength at the site, whereas all other aspects are identical to the GEMSTAT model. We show that this modification results in significantly better fits in a data set of over 30 enhancers regulating spatial expression patterns in the blastoderm-stage Drosophila embryo. It is important to note that the improved fits result not from an overall elevated accessibility in active enhancers but from the variation of accessibility levels within an enhancer. With whole-genome DNA accessibility measurements becoming increasingly popular, our work demonstrates how such data may be useful for sequence-to-expression models. It also calls for future advances in modeling accessibility levels from sequence and the transregulatory context, so as to predict accurately the effect of cis and trans perturbations on gene expression.

Abstract Denaturant-induced unfolding of helical membrane proteins provides insights into their mechanism of folding and domain organization, which take place in the chemically heterogeneous, anisotropic environment of a lipid membrane. Rhomboid proteases are intramembrane proteases that play key roles in various diseases. Crystal structures have revealed a compact helical bundle with a buried active site, which requires conformational changes for the cleavage of transmembrane substrates. A dimeric form of the rhomboid protease has been shown to be important for activity. In this study, we examine the mechanism of refolding for two distinct rhomboids to gain insight into their secondary structure-activity relationships. Although helicity is largely abolished in the unfolded states of both proteins, unfolding is completely reversible for HiGlpG but only partially reversible for PsAarA. Refolding of both proteins results in reassociation of the dimer, with a 90% regain of catalytic activity for HiGlpG but only a 70% regain for PsAarA. For both proteins, a broad, gradual transition from the native, folded state to the denatured, partly unfolded state was revealed with the aid of circular dichroism spectroscopy as a function of denaturant concentration, thus arguing against a classical two-state model as found for many globular soluble proteins. Thermal denaturation has irreversible destabilizing effects on both proteins, yet reveals important functional details regarding substrate accessibility to the buried active site. This concerted biophysical and functional analysis demonstrates that HiGlpG, with a simple six-transmembrane-segment organization, is more robust than PsAarA, which has seven predicted transmembrane segments, thus rendering HiGlpG amenable to in vitro studies of membrane-protein folding.

Abstract Glucosylceramide (GlcCer), one of the simplest glycosphingolipids, plays key roles in physiology and pathophysiology. It has been suggested that GlcCer modulates cellular events by forming specialized domains. In this study, we investigated the interplay between GlcCer and cholesterol (Chol), an important lipid involved in the formation of liquid-ordered (lo) phases. Using fluorescence microscopy and spectroscopy, and dynamic and electrophoretic light scattering, we characterized the interaction between these lipids in different pH environments. A quantitative description of the phase behavior of the ternary unsaturated phospholipid/Chol/GlcCer mixture is presented. The results demonstrate coexistence between lo and liquid-disordered (ld) phases. However, the extent of lo/ld phase separation is sparse, mainly due to the ability of GlcCer to segregate into tightly packed gel domains. As a result, the phase diagram of these mixtures is characterized by an extensive three-phase coexistence region of fluid (ld-phospholipid enriched)/lo (Chol enriched)/gel (GlcCer enriched). Moreover, the results show that upon acidification, GlcCer solubility in the lo phase is increased, leading to a larger lo/ld coexistence region. Quantitative analyses allowed us to determine the differences in the composition of the phases at neutral and acidic pH. These results predict the impact of GlcCer on domain formation and membrane organization in complex biological membranes, and provide a background for unraveling the relationship between the biophysical properties of GlcCer and its biological action.

Abstract Cu+-specific P-type ATPase membrane protein transporters regulate cellular copper levels. The lack of crystal structures in Cu+-binding states has limited our understanding of how ion entry and binding are achieved. Here, we characterize the molecular basis of Cu+ entry using molecular-dynamics simulations, structural modeling, and in vitro and in vivo functional assays. Protein structural rearrangements resulting in the exposure of positive charges to bulk solvent rather than to lipid phosphates indicate a direct molecular role of the putative docking platform in Cu+ delivery. Mutational analyses and simulations in the presence and absence of Cu+ predict that the ion-entry path involves two ion-binding sites: one transient Met148-Cys382 site and one intramembranous site formed by trigonal coordination to Cys384, Asn689, and Met717. The results reconcile earlier biochemical and x-ray absorption data and provide a molecular understanding of ion entry in Cu+-transporting P-type ATPases.

Abstract We describe the ways in which Galaxy, a web-based reproducible research platform, can be used for web-based sharing of complex computational models. Galaxy allows users to seamlessly customize and run simulations on cloud computing resources, a concept we refer to as Models and Simulations as a Service (MaSS). To illustrate this application of Galaxy, we have developed a tool suite for simulating a high spatial-resolution model of the cardiac Ca2+ spark that requires supercomputing resources for execution. We also present tools for simulating models encoded in the SBML and CellML model description languages, thus demonstrating how Galaxy’s reproducible research features can be leveraged by existing technologies. Finally, we demonstrate how the Galaxy workflow editor can be used to compose integrative models from constituent submodules. This work represents an important novel approach, to our knowledge, to making computational simulations more accessible to the broader scientific community.

Abstract Tightly clustered inositol trisphosphate receptors (IP3Rs) control localized Ca2+ liberation from the endoplasmic reticulum to generate repetitive Ca2+ puffs. Distributions of the interpuff interval (IPI), i.e., the waiting time between successive puffs, are found to be well characterized by a probability density function involving only two parameters, λ and ξ, which represent the basal rate of puff generation and the recovery rate from refractoriness, respectively. However, how the two parameters depend on the kinetic parameters of single IP3Rs in a cluster is still unclear. In this article, using a stochastic puff model and a single-channel data-based IP3R model, we establish the dependencies of λ and ξ on two important IP3R model parameters, IP3 concentration ([IP3]) and the recovery rate from Ca2+ inhibition (rlow). By varying [IP3] and rlow in physiologically plausible ranges, we find that the ξ-λ plane is comprised of only two disjoint regions, a biologically impermissible region and a region where each parameter set (ξ, λ) can be caused by using two different combinations of [IP3] and rlow. The two combinations utilize very different mechanisms to maintain the same IPI distribution, and the mechanistic difference provides a way of identifying IP3R kinetic parameters by observing properties of the IPI.

Abstract Changing the helical propensity of a polypeptide sequence might be expected to affect the conformational properties of the denatured state of a protein. To test this hypothesis, alanines at positions 83 and 87 near the center of helix 3 of cytochrome c′ from Rhodopseudomonas palustris were mutated to serine to decrease the stability of this helix. A set of 13 single histidine variants in the A83S/A87S background were prepared to permit assessment of the conformational properties of the denatured state using histidine-loop formation in 3 M guanidine hydrochloride. The data are compared with previous histidine-heme loop formation data for wild-type cytochrome c′. As expected, destabilization of helix 3 decreases the global stabilities of the histidine variants in the A83S/A87S background relative to the wild-type background. Loop stability versus loop size data yields a scaling exponent of 2.1 ± 0.2, similar to the value of 2.3 ± 0.2 obtained for wild-type cytochrome c′. However, the stabilities of all histidine-heme loops, which contain the helix 3 sequence segment, are increased in the A83S/A87S background compared to the wild-type background. Rate constants for histidine-heme loop breakage are similar for the wild-type and A83S/A87S variants. However, for histidine-heme loops that contain the helix 3 sequence segment, the rate constants for loop formation increase in the A83S/A87S background compared to the wild-type background. Thus, residual helical structure appears to stiffen the polypeptide chain slowing loop formation in the denatured state. The implications of these results for protein folding mechanisms are discussed.

Abstract Biological mixtures such as the cytosol may consist of thousands of distinct components. There is now a substantial body of evidence showing that, under physiological conditions, intracellular mixtures can phase separate into spatially distinct regions with differing compositions. In this article we present numerical evidence indicating that such spontaneous compartmentalization exploits general features of the phase diagram of a multicomponent biomolecular mixture. In particular, we show that demixed domains are likely to segregate when the variance in the intermolecular interaction strengths exceeds a well-defined threshold. Multiple distinct phases are likely to become stable under very similar conditions, which can then be tuned to achieve multiphase coexistence. As a result, only minor adjustments to the composition of the cytosol or the strengths of the intermolecular interactions are needed to regulate the formation of different domains with specific compositions, implying that phase separation is a robust mechanism for creating spatial organization. We further predict that this functionality is only weakly affected by increasing the number of components in the system. Our model therefore suggests that, for purely physico-chemical reasons, biological mixtures are naturally poised to undergo a small number of demixing phase transitions.