Direct interaction of Gβγ with a C-terminal Gβγ-binding domain of the Ca2+ channel α1 subunit is responsible for channel inhibition by G protein-coupled receptors

Several classes of voltage-gated Ca2+ channels (VGCCs) are inhibited by G proteins activated by receptors for neurotransmitters and neuromodulatory peptides. Evidence has accumulated to indicate that for non-L-type Ca2+ channels the executing arm of the activated G protein is its βγ dimer (Gβγ). We report below the existence of two Gβγ-binding sites on the A-, B-, and E-type α1 subunits that form non-L-type Ca2+ channels. One, reported previously, is in loop 1 connecting transmembrane domains I and II. The second is located approximately in the middle of the ca. 600-aa-long C-terminal tails. Both Gβγ-binding regions also bind the Ca2+ channel β subunit (CCβ), which, when overexpressed, interferes with inhibition by activated G proteins. Replacement in α1E of loop 1 with that of the G protein-insensitive and Gβγ-binding-negative loop 1 of α1C did not abolish inhibition by G proteins, but the exchange of the α1E C terminus with that of α1C did. This and properties of α1E C-terminal truncations indicated that the Gβγ-binding site mediating the inhibition of Ca2+ channel activity is the one in the C terminus. Binding of Gβγ to this site was inhibited by an α1-binding domain of CCβ, thus providing an explanation for the functional antagonism existing between CCβ and G protein inhibition. The data do not support proposals that Gβγ inhibits α1 function by interacting with the site located in the loop I–II linker. These results define the molecular mechanism by which presynaptic G protein-coupled receptors inhibit neurotransmission.

Folding and activity of circularly permuted forms of a polytopic membrane protein

The transmembrane subunit of the Glc transporter (IICBGlc), which mediates uptake and concomitant phosphorylation of glucose, spans the membrane eight times. Variants of IICBGlc with the native N and C termini joined and new N and C termini in the periplasmic and cytoplasmic surface loops were expressed in Escherichia coli. In vivo transport/in vitro phosphotransferase activities of the circularly permuted variants with the termini in the periplasmic loops 1 to 4 were 35/58, 32/37, 0/3, and 0/0% of wild type, respectively. The activities of the variants with the termini in the cytoplasmic loops 1 to 3 were 0/25, 0/4 and 24/70, respectively. Fusion of alkaline phosphatase to the periplasmic C termini stabilized membrane integration and increased uptake and/or phosphorylation activities. These results suggest that internal signal anchor and stop transfer sequences can function as N-terminal signal sequences in a circularly permuted α-helical bundle protein and that the orientation of transmembrane segments is determined by the amino acid sequence and not by the sequential appearance during translation. Of the four IICBGlc variants with new termini in periplasmic loops, only the one with the discontinuity in loop 4 is inactive. The sequences of loop 4 and of the adjacent TM7 and TM8 are conserved in all phosphoenolpyruvate-dependent carbohydrate:phosphotransferase system transporters of the glucose family.

Alternate protein frameworks for molecular recognition.

In an effort to determine whether proteins with structures other than the immunoglobulin fold can be used to mimic the ligand binding properties of antibodies, we generated a library from the four-helix bundle protein cytochrome b562 in which the two loops were randomized. Panning of this library against the bovine serum albumin (BSA) conjugate of N-methyl-p-nitrobenzylamine derivative 1 by phage display methods yielded cytochromes in which residues Trp-20, Arg-21, and Ser-22 in loop A and Arg-83 and Trp-84 in loop B were conserved. The individual mutants, which fold into native-like structure, bind selectively to the BSA-1 conjugate with micromolar dissociation constants (Kd), in comparison to a monoclonal antibody that binds selectively to 1 with a Kd of 290 nM. These and other antibody-like receptors may prove useful as therapeutic agents or as reagents for both intra- and extracellular studies.

Risco de mercado na comercialização de energia elétrica: uma análise estruturada com foco no ambiente de contratação livre - ACL.

A motivação para o desenvolvimento desse trabalho surge em um momento em que se verifica uma participação cada vez mais significativa das fontes energéticas renováveis não convencionais no País. Não obstante, o cenário de evolução evidencia que o arcabouço regulatório e as regras de mercado não acompanharam as especificidades inerentes à exploração dessas fontes. Assim, para que se mantenha adequado ritmo de inserção na matriz energética, devem ser buscadas opções para que fontes alternativas sejam cada vez mais competitivas na atual configuração do mercado energético. A contribuição dessa pesquisa, portanto, centra-se na análise dos riscos de mercado incorridos por esses geradores de fontes intermitentes de energia ao comercializarem energia no ambiente de contratação livre. Nessa perspectiva, a Dissertação foi desenvolvida abordando tipos de geração de energia e suas características técnicas e econômicas, legislação do setor elétrico, regras de comercialização, balanço energético do sistema, formação de preços no mercado de curto prazo e precificação de contratos no ACL, diferença de preços entre submercados, requisitos de flexibilidade e sazonalidade nos contratos de venda a consumidores livres e seu impacto na precificação de contratos, identificação de comportamento energético complementar para mitigação de riscos de mercado entre fontes renováveis e rebatimento na formulação de mecanismo de hedge, análise de portfólio de contratos e estratégia ótima de contratação de energia para agentes geradores atuando no ACL. Como resposta ao desafio de equacionar o impasse surgido na comercialização de fontes de produção sazonal, propõe-se um modelo para definir estratégias de contratação para agentes geradores e comercializadores a partir da complementação energética entre diferentes tipos de fontes, de forma a maximizar os ganhos de comercialização para um risco estabelecido. Busca-se a composição ideal dessas fontes na carteira de um comercializador para minimizar o risco de exposição à volatilidade dos preços do mercado de curto prazo. Isso é possível em virtude das compensações energéticas feitas entre as diferentes fontes em um portfólio combinado, mitigando a receita em risco decorrente das variações que existem nos preços de curto prazo e na produção energética. De forma complementar, estruturou-se um modelo de negócio no qual uma empresa detentora de ativos de geração hidrelétrica compra os direitos de produção de uma eólica e/ou biomassa para incorporar ao seu portfólio e vender como contrato por quantidade. Determinou-se o volume de energia a ser comprado de cada fonte, o preço, a estratégia mais indicada de contratação e a mitigação de fatores de risco contemplados nos contratos de venda, buscando maximizar os ganhos de comercialização condicionada a critérios de risco pré-fixados.

The influence of magnetic field geometry on magnetars X-ray spectra

Nowadays, the analysis of the X-ray spectra of magnetically powered neutron stars or magnetars is one of the most valuable tools to gain insight into the physical processes occurring in their interiors and magnetospheres. In particular, the magnetospheric plasma leaves a strong imprint on the observed X-ray spectrum by means of Compton up-scattering of the thermal radiation coming from the star surface. Motivated by the increased quality of the observational data, much theoretical work has been devoted to develop Monte Carlo (MC) codes that incorporate the effects of resonant Compton scattering (RCS) in the modeling of radiative transfer of photons through the magnetosphere. The two key ingredients in this simulations are the kinetic plasma properties and the magnetic field (MF) configuration. The MF geometry is expected to be complex, but up to now only mathematically simple solutions (self-similar solutions) have been employed. In this work, we discuss the effects of new, more realistic, MF geometries on synthetic spectra. We use new force-free solutions [14] in a previously developed MC code [9] to assess the influence of MF geometry on the emerging spectra. Our main result is that the shape of the final spectrum is mostly sensitive to uncertain parameters of the magnetospheric plasma, but the MF geometry plays an important role on the angle-dependence of the spectra.

Constraining storm-scale forecasts of deep convective initiation with surface weather observations

Thesis (Ph.D.)--University of Washington, 2016-06

Simulating Nursing: Rhetoric, Materiality, and Disciplinary Learning

Thesis (Ph.D.)--University of Washington, 2016-06

Binomial leap methods for simulating stochastic chemical kinetics

This paper discusses efficient simulation methods for stochastic chemical kinetics. Based on the tau-leap and midpoint tau-leap methods of Gillespie [D. T. Gillespie, J. Chem. Phys. 115, 1716 (2001)], binomial random variables are used in these leap methods rather than Poisson random variables. The motivation for this approach is to improve the efficiency of the Poisson leap methods by using larger stepsizes. Unlike Poisson random variables whose range of sample values is from zero to infinity, binomial random variables have a finite range of sample values. This probabilistic property has been used to restrict possible reaction numbers and to avoid negative molecular numbers in stochastic simulations when larger stepsize is used. In this approach a binomial random variable is defined for a single reaction channel in order to keep the reaction number of this channel below the numbers of molecules that undergo this reaction channel. A sampling technique is also designed for the total reaction number of a reactant species that undergoes two or more reaction channels. Samples for the total reaction number are not greater than the molecular number of this species. In addition, probability properties of the binomial random variables provide stepsize conditions for restricting reaction numbers in a chosen time interval. These stepsize conditions are important properties of robust leap control strategies. Numerical results indicate that the proposed binomial leap methods can be applied to a wide range of chemical reaction systems with very good accuracy and significant improvement on efficiency over existing approaches. (C) 2004 American Institute of Physics.

Effect of granule morphology on breakage behaviour during compression

In the area of dry particle breakage, Discrete Element Method (DEM) simulations have been widely used to analyse the sensitivity of various physical parameters to the behaviour of agglomerates during breakage. This paper looks at the effect of agglomerate shape and structure on the mechanisms and extent of breakage of dry agglomerates under compressive load using DEM simulations. In the simulations, a spherical-shaped agglomerate produced within the DEM code is compared with an irregularly shaped agglomerate, whose structure is that of an actual granule that was characterised with X-ray microtomography (muCT). Both agglomerates have identical particle size distribution, coordination number and surface energy values, with only the agglomerate shape and structure differing between the two. The work here details the breakage behaviour with a number of traditional DEM output parameters (i.e., contact/cluster distributions) with showing vastly different behaviour between the two agglomerates. (C) 2004 Elsevier B.V. All rights reserved.

Properties of the stochastic Gross-Pitaevskii equation: finite temperature Ehrenfest relations and the optimal plane wave representation

We present Ehrenfest relations for the high temperature stochastic Gross-Pitaevskii equation description of a trapped Bose gas, including the effect of growth noise and the energy cutoff. A condition for neglecting the cutoff terms in the Ehrenfest relations is found which is more stringent than the usual validity condition of the truncated Wigner or classical field method-that all modes are highly occupied. The condition requires a small overlap of the nonlinear interaction term with the lowest energy single particle state of the noncondensate band, and gives a means to constrain dynamical artefacts arising from the energy cutoff in numerical simulations. We apply the formalism to two simple test problems: (i) simulation of the Kohn mode oscillation for a trapped Bose gas at zero temperature, and (ii) computing the equilibrium properties of a finite temperature Bose gas within the classical field method. The examples indicate ways to control the effects of the cutoff, and that there is an optimal choice of plane wave basis for a given cutoff energy. This basis gives the best reproduction of the single particle spectrum, the condensate fraction and the position and momentum densities.

Numerical modeling of 11.1 T MRI of a human head using MoM/FDTD methods

In a recent study, severe distortions in the proton images of an excised, fixed, human brain in an 11.1 Tesla/40 cm MR instrument have been observed, and the effect modeled on phantom images using a finite difference time domain (FDTD) model. in the present study, we extend these simulations to that of a complete human head, employing a hybrid FDTD and method of moments (MoM) approach, which provides a validated method for simulating biological samples in coil structures. The effect of fixative on the image distortions is explored. importantly, temperature distributions within the head are also simulated using a bioheat method based on parameters derived from the electromagnetic simulations. The MoM/FDTD simulations confirm that the transverse magnetic field (B,) from a ReCav resonator exhibits good homogeneity in air but strong inhomogeneity when loaded with the head with or without fixative. The fixative serves to increase the distortions, but they are still significant for the in vivo simulations. The simulated signal intensity (SI) distribution within the sample confirm the distortions in the experimental images are caused by the complex interactions of the incident electromagnetic fields with tissue, which is heterogeneous in terms of conductivity and permittivity. The temperature distribution is likewise heterogeneous, raising concerns regarding hot spot generation in the sample that may exceed acceptable levels in future in vivo studies. As human imaging at 11.1 T is some time away, simulations are important in terms of predicting potential safety issues as well as evaluating practical concerns about the quality of images. Simulation on a whole human head at 11.1 T implies the wave behavior presents significant engineering challenges for ultra-high-field (UHF) MRI. Novel strategies will have to be employed in imaging technique and resonator design for UHF MRI to achieve the theoretical signal-to-noise ratio (SNR) improvements it offers over lower field systems. (C) 2005 Wiley Periodicals, Inc.

Parallel implementation of stochastic simulation for large scale cellular processes

Experimental and theoretical studies have shown the importance of stochastic processes in genetic regulatory networks and cellular processes. Cellular networks and genetic circuits often involve small numbers of key proteins such as transcriptional factors and signaling proteins. In recent years stochastic models have been used successfully for studying noise in biological pathways, and stochastic modelling of biological systems has become a very important research field in computational biology. One of the challenge problems in this field is the reduction of the huge computing time in stochastic simulations. Based on the system of the mitogen-activated protein kinase cascade that is activated by epidermal growth factor, this work give a parallel implementation by using OpenMP and parallelism across the simulation. Special attention is paid to the independence of the generated random numbers in parallel computing, that is a key criterion for the success of stochastic simulations. Numerical results indicate that parallel computers can be used as an efficient tool for simulating the dynamics of large-scale genetic regulatory networks and cellular processes

A verifiable solution to the MEG inverse problem

Magnetoencephalography (MEG) is a non-invasive brain imaging technique with the potential for very high temporal and spatial resolution of neuronal activity. The main stumbling block for the technique has been that the estimation of a neuronal current distribution, based on sensor data outside the head, is an inverse problem with an infinity of possible solutions. Many inversion techniques exist, all using different a-priori assumptions in order to reduce the number of possible solutions. Although all techniques can be thoroughly tested in simulation, implicit in the simulations are the experimenter's own assumptions about realistic brain function. To date, the only way to test the validity of inversions based on real MEG data has been through direct surgical validation, or through comparison with invasive primate data. In this work, we constructed a null hypothesis that the reconstruction of neuronal activity contains no information on the distribution of the cortical grey matter. To test this, we repeatedly compared rotated sections of grey matter with a beamformer estimate of neuronal activity to generate a distribution of mutual information values. The significance of the comparison between the un-rotated anatomical information and the electrical estimate was subsequently assessed against this distribution. We found that there was significant (P < 0.05) anatomical information contained in the beamformer images across a number of frequency bands. Based on the limited data presented here, we can say that the assumptions behind the beamformer algorithm are not unreasonable for the visual-motor task investigated.

Investigation of the activation mechanism and structural characterization of the CGRP receptor

The calcitonin gene-related peptide (CGRP) receptor is an unusual G protein-coupled receptor (GPCR) in that it comprises the calcitonin receptor-like receptor (CLR), receptor activity modifying protein 1 (RAMP1) and the receptor component protein (RCP). The RAMP1 has two other homologues – RAMP2 and RAMP3. The endogenous ligand for this receptor is CGRP, a 37 amino acid neuropeptide that act as a vasodilator. This peptide has been implicated in the aetiology of health conditions such as inflammation, Reynaud’s disease and migraine. A clear understanding of the mode of activation of this receptor could be key in developing therapeutic agents for associated health conditions. Although the crystal structure of the N-terminal extracellular domain (ECD) of this receptor (in complex with an antagonist) has been published, the details of receptor-agonist interactions at this domain, and so ultimately the mechanism of receptor activation, are still unclear. Also, the C-terminus of the CLR (in the CGRP receptor), especially around the presumed helix 8 (H8) region, has not been well studied for its role in receptor signalling. This research project investigated these questions. In this study, certain residues making up the putative N-terminal ligand-binding core of the CLR (in the CGRP receptor) were mapped out and found to be crucial for receptor signalling. They included W69 and D70 of the WDG motif in family B GPCRs, as well as Y91, F92, D94 and F95 in loop 2 of CLR N-terminus. Also, F163 at the cytoplasmic end of TM1 and certain residues spanning H8 and associated C-terminal region of CLR were found to be required for CGRP receptor signalling. These residues were investigated by site-directed mutagenesis where they were mutated to alanine (or other residues in specific cases) and the effect of the mutations on receptor pharmacology assessed by evaluating cAMP production, cell surface expression, total cell expression and aCGRP-mediated receptor internalization. Moreover, the N-terminal ECDs of the CLR and RAMPs (RAMP1, RAMP2 and RAMP3) were produced in a yeast host strain (Pichia pastoris) for the purpose of structural interaction study by surface plasmon resonance (SPR). Following expression and purification, these receptor proteins were found to individually retain their secondary structures when analysed by circular dichroism (CD). Results were analysed and interpreted with the knowledge of the secretin family receptor paradigm. The research described in this thesis has produced novel data that contributes to a clearer understanding of CGRP receptor pharmacology. The study on CLR and RAMPs ECDs could be a useful tool in determining novel interacting GPCR partners of RAMPs.

Three-dimensional finite element analysis of the effects of anisotropy on bone mechanical properties measured by nanoindentation

A three-dimensional finite element analysis (FEA) model with elastic-plastic anisotropy was built to investigate the effects of anisotropy on nanoindentation measurements for cortical bone. The FEA model has demonstrated a capability to capture the cortical bone material response under the indentation process. By comparison with the contact area obtained from monitoring the contact profile in FEA simulations, the Oliver-Pharr method was found to underpredict or overpredict the contact area due to the effects of anisotropy. The amount of error (less than 10% for cortical bone) depended on the indentation orientation. The indentation modulus results obtained from FEA simulations at different surface orientations showed a trend similar to experimental results and were also similar to moduli calculated from a mathematical model. The Oliver-Pharr method has been shown to be useful for providing first-order approximations in the analysis of anisotropic mechanical properties of cortical bone, although the indentation modulus is influenced by anisotropy.