A structural disulfide of yeast protein-disulfide isomerase destabilizes the active site disulfide of the N-terminal thioredoxin domain

Protein-disulfide isomerase (PDI) is an essential catalyst of disulfide formation and isomerization in the eukaryotic endoplasmic reticulum. PDI has two active sites at either end of the molecule, each containing two cysteines that facilitate thiol-disulfide exchange. In addition to its four catalytic cysteines, PDI possesses two non-active site cysteines whose location and separation distance varies by organism. In higher eukaryotes, the non-active site cysteines are located in the C-terminal half of the protein sequence and are separated by 30 amino acids. In contrast, the internal cysteines of PDI from lower eukaryotes are located near the N-terminal active site and are much closer together in sequence. The function of these cysteines and the significance of their unique location in yeast PDI have been unclear. Previous data (Xiao, R., Wilkinson, B., Solovyov, A., Winther, J. R., Holmgren, A., Lundstrom-Ljung, J., and Gilbert, H. F. (2004) J. Biol. Chem. 279, 49780-49786) suggest that the internal cysteines exist as a disulfide in the endoplasmic reticulum of Saccharomyces cerevisiae. By coupling mass spectrometry with a gel-shift technique that allows us to measure the redox potentials of the PDI active sites in the presence and absence of the non-active site cysteines, we find that the non-active site cysteines form a disulfide that is stable even in a very reducing environment and demonstrate that this disulfide exists to destabilize the N-terminal active site disulfide, making it a better oxidant by 18-fold. Consistent with this finding, we show that mutating the non-active site cysteines to alanines disrupts both the oxidase and isomerase activities of PDI in vitro.     

Catalysis in glycine N-methyltransferase: testing the electrostatic stabilization and compression hypothesis

Glycine N-methyltransferase (GNMT) is an S-adenosyl-l-methionine dependent enzyme that catalyzes glycine transformation to sarcosine. Here, we present a hybrid quantum mechanics/molecular mechanics (QM/MM) computational study of the reaction compared to the counterpart process in water. The process takes place through an SN2 mechanism in both media with a transition state in which the transferring methyl group is placed in between the donor (SAM) and the acceptor (the amine group of glycine). Comparative analysis of structural, electrostatic, and electronic characteristics of the in-solution and enzymatic transition states allows us to get a deeper insight into the origins of the enzyme's catalytic power. We found that the enzyme is able to stabilize the substrate in its more active basic form by means of a positively charged residue (Arg175) placed in the active site. However, the maximum stabilization is attained for the transition state. In this case, the enzyme is able to form stronger hydrogen bonds with the positively charged amine group. Finally, we show that in agreement with previous computational studies on other methyltransferases, there is no computational evidence for the compression hypothesis, as was formulated by Schowen (Hegazi, M. F., Borchardt, R. T., and Schowen, R. L. (1979) J. Am. Chem. Soc. 101, 4359-4365).     

Bayesian mapping of QTL in outbred F2 families allowing inference about whether F0 grandparents are homozygous or heterozygous at QTL

In this paper, we propose a new Bayesian method for QTL analysis in outbred F2 families based on Markov chain Monte Carlo (MCMC) estimation allowing inference about whether each of F0 founders (grandparents) is homozygous or heterozygous at QTL. This, in turn, allows us to select a model accurately explaining observations of phenotypes for F2 individuals. The proposed method performs the fitting a statistical model of the two possible QTL states in each F0 grandparent, that is, homozygous and heterozygous at QTL, and gives a posterior distribution for the QTL states in each F0 grandparent. We confine ourselves to the discrimination of two QTL states, homozygous or heterozygous, for each of the F0 grandparents without taking into consideration whether common alleles are shared by F0 grandparents. The statistical model includes allelic effects and dominance effects for each QTL. The number of parameters representing allelic effects and dominance effects is therefore changed depending on the QTL states. A Reversible Jump MCMC technique is used for transition between the models of different dimensions. The effectiveness of the proposed method was investigated using simulation experiments. It was practicable to estimate the QTL states of F0 grandparents as well as the number, the locations and the effects of QTL segregating in an outbred F2 family.     

Study of Circular Dichroism Modes Through Decomposition of Planar Nanostructures

We extend the plasmon hybridization method from a single nanoparticle to a complex planar nanostructure, decomposing the complex nanostructure into fundamental nanoparticle building blocks. Using gammadion nanostructure as an example, we validated the theory by comparing the field profile in the gammadion’s arms under the influence of an incident circularly polarized wave. This allows us to address the origin of the plasmonics modes in the circular dichroism (CD) spectrum. The use of this hybridization method provides a simple and intuitive explanation on how conductive and inductive coupling may result from complex planar nanostructures, allowing us to study its optical properties. Using our approach, top down hybridization studies can be applied to other complex planar structures to gain further insight on the origin of the CD modes and enhance ultrasensitive sensing of chiral micro and macro molecules.     

A novel galactosyl-binding lectin from the plasma of the blood clam,Anadara granosa (L) and a study of its combining site

Control Analysis has been carried out in the first steps of a rat liver glycolytic system. Attention has been focused on the effect of several glucose concentrations on the control, particularly regarding the role of glucokinase. From kinetic studies of the whole metabolic system we have obtained information on the flux variation under different glucose concentrations. This information together with the kinetics of glucokinase has allowed us to calculate Flux Control and Elasticity Coefficients for glucokinase and the Response Coefficient of the system with respect to glucose. The changes in of the value of Flux Control Coefficients demonstrates that in conditions of low glucose concentration, glucokinase is the main enzyme in controlling the flux through the pathway, but at high glucose concentration the control moves to phosphofructokinase. Next, we have compared our results with those obtained with the shortening and titration method, previously described (Torres, N.V., Mateo, F., Mélendez-Hevia, E. and Kacser, H., (1986) Biochem. J. 234, 169–174; Torres, N.V. and Meléndez-Hevia, E. 1991. Molec. Cell. Biochem. 101, 1–10). Furthermore, from knowledge of the enzyme kinetics of the system we have been able to build a model of the pathway that allows us computer similation of its behavior and calculation of the Flux Control Coefficient profile at different glucose concentrations. By the three methods the results correlate, supporting the use of the pathway substrate as external modulator of the metabolic system as a tool for practical application of Control Analysis.     

Plant metallothionein domains: functional insight into physiological metal binding and protein folding

Plant metallothioneins (MTs) differ from animal MTs by a peculiar sequence organization consisting of two short cysteine-rich terminal domains linked by a long cysteine-devoid spacer. The role of the plant MT domains in the protein structure and functionality is largely unknown. Here, we investigate the separate domain contribution to the in vivo binding of Zn and Cu and to confer metal tolerance to CUP1-null yeast cells of a plant type 2 MT (QsMT). For this purpose, we obtained three recombinant peptides that, respectively, correspond to the single N-terminal (N25) and C-terminal (C18) cysteine-rich domains of QsMT, and a chimera in which the spacer is replaced with a four-glycine bridge (N25-C18). The metal-peptide preparations recovered from Zn- or Cu-enriched cultures were characterized by ESI-MS, ICP-OES and CD and UV-vis spectroscopy and data compared to full length QsMT. Results are consistent with QsMT giving rise to homometallic Zn- or Cu-MT complexes according to a hairpin model in which the two Cys-rich domains interact to form a cluster. In this model the spacer region does not contribute to the metal coordination. However, our data from Zn-QsMT (but not from Cu-QsMT) support a fold of the spacer involving some interaction with the metal core. On the other hand, results from functional complementation assays in endogenous MT-defective yeast cells suggest that the spacer region may play a role in Cu-QsMT stability or subcellular localization. As a whole, our results provide the first insight into the structure/function relationship of plant MTs using the analysis of the separate domain abilities to bind physiological metals.     

Synapse fits neuron: joint reduction by model inversion

In this paper, we introduce a novel simplification method for dealing with physical systems that can be thought to consist of two subsystems connected in series, such as a neuron and a synapse. The aim of our method is to help find a simple, yet convincing model of the full cascade-connected system, assuming that a satisfactory model of one of the subsystems, e.g., the neuron, is already given. Our method allows us to validate a candidate model of the full cascade against data at a finer scale. In our main example, we apply our method to part of the squid’s giant fiber system. We first postulate a simple, hypothetical model of cell-to-cell signaling based on the squid’s escape response. Then, given a FitzHugh-type neuron model, we derive the verifiable model of the squid giant synapse that this hypothesis implies. We show that the derived synapse model accurately reproduces synaptic recordings, hence lending support to the postulated, simple model of cell-to-cell signaling, which thus, in turn, can be used as a basic building block for network models.     

Structural and spectroscopic analysis of the F393H mutant of flavocytochrome P450 BM3.

In the preceding paper in this issue [Ost, T. W. B., Miles, C. S., Munro, A. W., Murdoch, J., Reid, G. A., and Chapman, S. K. (2001) Biochemistry 40, 13421-13429], we have established that the primary role of the phylogenetically conserved phenylalanine in flavocytochrome P450 BM3 (F393) is to control the thermodynamic properties of the heme iron, so as to optimize electron-transfer both to the iron (from the flavin redox partner) and onto molecular oxygen. In this paper, we report a detailed study of the F393H mutant enzyme, designed to probe the structural, spectroscopic, and metabolic profile of the enzyme in an attempt to identify the factors responsible for causing the changes. The heme domain structure of the F393H mutant has been solved to 2.0 A resolution and demonstrates that the histidine replaces the phenylalanine in almost exactly the same conformation. A solvent water molecule is hydrogen bonded to the histidine, but there appears to be little other gross alteration in the environment of the heme. The F393H mutant displays an identical ferric EPR spectrum to wild-type, implying that the degree of splitting of the iron d orbitals is unaffected by the substitution, however, the overall energy of the d-orbitals have changed relative to each other. Magnetic CD studies show that the near-IR transition, diagnostic of heme ligation state, is red-shifted by 40 nm in F393H relative to wild-type P450 BM3, probably reflecting alteration in the strength of the iron-cysteinate bond. Studies of the catalytic turnover of fatty acid (myristate) confirms NADPH oxidation is tightly coupled to fatty acid oxidation in F393H, with a product profile very similar to wild-type. The results indicate that gross conformational changes do not account for the perturbations in the electronic features of the P450 BM3 heme system and that the structural environment on the proximal side of the P450 heme must be conformationally conserved in order to optimize catalytic function.     

The tricarboxylic acid cycle in Dictyostelium discoideum. I. Formulation of alternative kinetic representations.

Enzyme systems within living cells have recently been shown to be highly ordered structures that violate classic assumptions of the Michaelis-Menten formalism, which originally was developed for the characterization of isolated reactions in vitro. This evidence suggests that a thorough examination of alternative kinetic formalisms for integrated biochemical systems is in order. The purpose of this series of papers is to assess the utility of an alternative power-law formalism by carrying out a detailed comparative analysis of a relatively large, representative system--the tricarboxylic acid cycle of Dictyostelium discoideum. This system was chosen because considerable experimental information already has been synthesized into a detailed kinetic model of the intact system. In this first paper, we set the stage for subsequent analysis within the framework of the power-law formalism: we review the underlying theory, emphasizing recent developments, formulate the model in terms that are convenient for the analysis to follow, and develop the system representation in both the Michaelis-Menten and power-law forms. In the second paper (Shiraishi, F., and Savageau, M. A. (1992) J. Biol. Chem. 267, 22919-22925), these alternative representations are shown to be internally consistent and locally equivalent. The third paper (Shiraishi, F., and Savageau, M. A. (1992) J. Biol. Chem. 267, 22926-22933) provides a complete analysis of the steady state behavior and also treats the dynamic behavior of the model.     

Disulfide chromophore and its optical activity

The compounds I-IV derived from α-D-cyclodextrin moiety by bridging and/or interconnecting with various patterns of disulfide bonds were chosen as models for the spectroscopic study of conformation of the disulfide bridge. The energy gap between the disulfide and cyclodextrin's electronic transitions allows us to investigate absorption and electronic circular dichroism spectra without disturbing spectral overlaps with amides or aromatic amino acids in peptides or proteins. Raman optical activity (ROA) spectra were measured and the bands due to S-S and C-S stretching motion identified. Comparison with the quantum mechanical calculations of simple models indicates that sense of disulfide twist follows sign of the measured S-S ROA band.     

Cooperativity of complex salt bridges

The energetic contribution of complex salt bridges, in which one charged residue (anchor residue) forms salt bridges with two or more residues simultaneously, has been suggested to have importance for protein stability. Detailed analysis of the net energetics of complex salt bridge formation using double- and triple-mutant cycle analysis revealed conflicting results. In two cases, it was shown that complex salt bridge formation is cooperative, i.e., the net strength of the complex salt bridge is more than the sum of the energies of individual pairs. In one case, it was reported that complex salt bridge formation is anti-cooperative. To resolve these different findings, we performed analysis of the geometries of salt bridges in a representative set of structures from the PDB and found that over 87% of all complex salt bridges anchored by Arg/Lys have a geometry such that the angle formed by their Calpha atoms, Theta, is <90 degrees . This preferred geometry is observed in the two reported instances when the energetics of complex salt bridge formation is cooperative, while in the reported anti-cooperative complex salt bridge, Theta is close to 160 degrees . Based on these observations, we hypothesized that complex salt bridges are cooperative for Theta < 90 degrees and anti-cooperative for 90 degrees < Theta < 180 degrees . To provide a further experimental test for this hypothesis, we engineered a complex salt bridge with Theta = 150 degrees into a model protein, the activation domain of human procarboxypeptidase A2 (ADA2h). Experimentally derived stabilities of the ADA2h variants allowed us to show that the complex salt bridge in ADA2h is anti-cooperative.     

The tricarboxylic acid cycle in Dictyostelium discoideum. II. Evaluation of model consistency and robustness.

When kinetic models of complex biochemical systems are reconstructed from knowledge of the component reactions that have been characterized in vitro, or when values must be assumed for some of the parameters, errors are invariably encountered, and, as a consequence, the resulting model is frequently internally inconsistent. The simplest and most basic manifestations of such logical inconsistency are the failure of the model to exhibit a steady state or to yield a steady state that is in agreement with the actual steady state of the integrated system, or to yield a steady state that is dynamically stable. Models that are consistent may nonetheless be lacking in robustness, which is manifested as a pathological sensitivity to small changes in the values of their parameters. In this paper, we examine the current model of the tricarboxylic acid cycle in Dictyostelium discoideum (see Shiraishi, F., and Savageau, M. A. (1992) J. Biol. Chem. 267, 22912-22918) with regard to these basic indicators of model quality. This may be viewed as a preliminary analysis; the object is to determine whether or not the model is reasonable and worthy of a more refined analysis and, if not, to diagnose the areas in need of modification before further analysis is undertaken. The results demonstrate that the current model of the tricarboxylic acid cycle is self-consistent and possesses a steady state that is in agreement with experimental evidence. However, the results also suggest that this model is not very robust. The high sensitivities of parameters influencing pyruvate metabolism indicate that the experimental characterization of these reactions might be fruitfully re-examined. These high sensitivities lead us to predict that this model of the tricarboxylic acid cycle should be accurate only over a very narrow range in variation of the independent variables. This is verified by the results presented in the following paper (Shiraishi, F., and Savageau, M. A. (1992) J. Biol. Chem. 267, 22926-22933).     

Estimating Summer Nutrient Concentrations in Northeastern Lakes from SPARROW Load Predictions and Modeled Lake Depth and Volume

Global nutrient cycles have been altered by the use of fossil fuels and fertilizers resulting in increases in nutrient loads to aquatic systems. In the United States, excess nutrients have been repeatedly reported as the primary cause of lake water quality impairments. Setting nutrient criteria that are protective of a lakes ecological condition is one common solution; however, the data required to do this are not always easily available. A useful solution for this is to combine available field data (i.e., The United States Environmental Protection Agency (USEPA) National Lake Assessment (NLA)) with average annual nutrient load models (i.e., USGS SPARROW model) to estimate summer concentrations across a large number of lakes. In this paper we use this combined approach and compare the observed total nitrogen (TN) and total phosphorus (TN) concentrations in Northeastern lakes from the 2007 National Lake Assessment to those predicted by the Northeast SPARROW model. We successfully adjusted the SPARROW predictions to the NLA observations with the use of Vollenweider equations, simple input-output models that predict nutrient concentrations in lakes based on nutrient loads and hydraulic residence time. This allows us to better predict summer concentrations of TN and TP in Northeastern lakes and ponds. On average we improved our predicted concentrations of TN and TP with Vollenweider models by 18.7% for nitrogen and 19.0% for phosphorus. These improved predictions are being used in other studies to model ecosystem services (e.g., aesthetics) and dis-services (e.g. cyanobacterial blooms) for ~18,000 lakes in the Northeastern United States.     

Determination of Fundamental Morphological Parameters of Supported Nanoparticle Ensembles

A new model to extract important morphological parameters of noble metal nanoparticle ensembles with a broad size and shape distribution is presented. The technique is based on a rigorous simulation of the inhomogeneously broadened extinction profiles of nanoparticle ensembles. As input data, only experimentally accessible parameters, such as the amount of deposited material, the nanoparticle number density, and the relative size distribution of the nanoparticles, are used. The model can be applied to oblate nanoparticles, which exhibit a strong correlation between their shape and size, e.g., to supported nanoparticles generated, for example, by deposition of atoms and subsequent nucleation or by gas phase deposition. Both methods are standard preparation techniques to generate well-defined nanoparticle ensembles under ultra high vacuum conditions. We apply our model to gold and silver nanoparticles on sapphire and TiO₂ supports and obtain a perfect agreement between the calculated and experimental data. More importantly, we could extract the functional dependence between the axial ratio and the radius of the nanoparticles within the ensemble and, therewith, the most probable axial ratio in the ensemble. In addition, the extinction spectrum of a nanoparticle ensemble irradiated with nanosecond pulsed laser light during growth has been successfully modeled. This demonstrates, that the model is able to describe shape changes of resonantly heated nanoparticles within the ensemble. By using the coverage as a free parameter, we could calculate from the extinction spectrum the average particle radius as well as the amount of desorbed atoms after irradiation with laser light. In summary, the model allows a fast, easy, but extensive morphological characterization of nanoparticle ensembles that exhibit a broad size and shape distribution.     

Voltage clamping with a single microelectrode.

A technique is described which allows neurons to be voltage clamped with a single microelectrode, and the advantages of this circuit with respect to conventional bridge techniques are discussed. In this circuit, the single microelectrode is rapidly switched from a current passing to a recording mode. The circuitry consists of: (1) an electronic switch; (2) a high impedance, ultralow input capacity amplifier; (3) a sample-and-hold module; (4) conventional voltage clamping circuitry. The closed electronic switch allows current to flow through the electrode. The switch then opens, and the electrode is in a recording mode. The low input capacity of the preamplifier allows the artifact from the current pulse to rapidly abate, after which time the circuit samples the membrane potential. This cycle is repeated at rates up to 10 kHz. The voltage clamping amplifier senses the output of the sample-and-hold module and adjusts the current pulse amplitude to maintain the desired membrane potential. The system was evaluated in Aplysia neurons by inserting two microelectrodes into a cell. One electrode was used to clamp the cell and the other to independently monitor membrane potential at a remote location in the soma.     

EVOLUTION AND MANIPULATION OF PARASITOID EGG LOAD

In proovigenic parasitoids such as Leptopilina boulardi, the female emerges with a limited egg load and no further eggs are produced during its adult life. A female thus runs the risk of exhausting this limited supply of eggs before the end of her life. Given that the production of an egg is costly, what is the evolutionarily stable egg load at emergence? This question has attracted a lot of attention in the last decade. Here, we analyze a model that allows us to track both the evolution and the population dynamics of a solitary, proovigenic parasitoid. First, we show how host–parasitoid dynamics feedbacks on the evolution of parasitoid egg load. Second, we use this model to consider the situation in which the parasitoid can be infected by a virus that manipulates the oviposition behavior of the females. In particular, we model the effect of the LbFV virus in L. boulardi, a virus that is known to enhance its horizontal transmission by increasing superparasitism (i.e., the laying of eggs in a host already parasitized). Specifically, we model (1) the effect of the virus on parasitoid egg load strategies , and (2) the evolution of egg load manipulation by the virus. This analysis yields two alternative, yet not mutually exclusive, adaptive explanations for the observation that females infected by the virus harbor higher egg loads than uninfected females. Infected females could either respond plastically to the infection status, or be manipulated by the virus. Further experimental work is required to distinguish between these two hypotheses. In a broader context, we present a general theoretical framework that allows us to study the epidemiology, the evolution, the coevolution, and the evolution of manipulation of various reproductive strategies of parasitoids.     

A fuzzy expert system design for analysis of body sounds and design of an unique electronic stethoscope (development of HILSA kit)

In this paper we have developed a fuzzy expert system (FES) for different sounds produced by different organs in the human body. We have also constructed a unique electronic stethoscope. The human body sounds produced by different organs like heart, lungs and intestine were analyzed. The doctor provided the data and relation between variables chosen for each organ sound. Using this information a rule base for fuzzy expert system was built. Such FES helps the medical doctor in arriving at appropriate decision in different difficult clinical situations. The examination of body sounds was done using conventional stethoscope (CS) and electronic stethoscope (ES), which was uniquely designed for this study. We have found that unique stethoscope developed by us is far superior to conventional stethoscope by its overall performance.     

Characterization of the Shank family of synaptic proteins. Multiple genes, alternative splicing, and differential expression in brain and development.

Shank1, Shank2, and Shank3 constitute a family of proteins that may function as molecular scaffolds in the postsynaptic density (PSD). Shank directly interacts with GKAP and Homer, thus potentially bridging the N-methyl-D-aspartate receptor-PSD-95-GKAP complex and the mGluR-Homer complex in synapses (Naisbitt, S., Kim, E., Tu, J. C. , Xiao, B., Sala, S., Valtschanoff, J., Weinberg, R. J., Worley, P. F., and Sheng, M. (1999) Neuron 23, 569-582; Tu, J. C., Xiao, B., Naisbitt, S., Yuan, J. P., Petralia, R. S., Brakeman, P., Doan, A., Aakalu, V. K., Lanahan, A. A., Sheng, M., and Worley, P. F. (1999) Neuron 23, 583-592). Shank contains multiple domains for protein-protein interaction including ankyrin repeats, an SH3 domain, a PSD-95/Dlg/ZO-1 domain, a sterile alpha motif domain, and a proline-rich region. By characterizing Shank cDNA clones and RT-PCR products, we found that there are four sites for alternative splicing in Shank1 and another four sites in Shank2, some of which result in deletion of specific domains of the Shank protein. In addition, the expression of the splice variants is differentially regulated in different regions of rat brain during development. Immunoblot analysis of Shank proteins in rat brain using five different Shank antibodies reveals marked heterogeneity in size (120-240 kDa) and differential spatiotemporal expression. Shank1 immunoreactivity is concentrated at excitatory synaptic sites in adult brain, and the punctate staining of Shank1 is seen in developing rat brains as early as postnatal day 7. These results suggest that alternative splicing in the Shank family may be a mechanism that regulates the molecular structure of Shank and the spectrum of Shank-interacting proteins in the PSDs of adult and developing brain.     

Induction of phosphoribomutase in Bacillus cereus growing on nucleosides

In this paper we show that phosphoribomutase is induced in Bacillus cereus by the same metabolizable purine and pyrimidine ribonucleosides previously shown to induce the purine nucleoside phosphorylase (Tozzi, M.G., Sgarrella, F. and Ipata, P.L. (1981) Biochim. Biophys. Acta 678, 460-466). The mutase allows ribose 1-phosphate formed from nucleosides to be utilized by the cell through the pentose cycle, upon transformation to ribose 5-phosphate. The equilibrium constant of the mutase reaction is towards ribose-5-phosphate formation. The coordinate induction of the two enzymes completes the picture of the molecular events leading to the utilization of the sugar moiety of purine nucleotides and nucleosides as an energy source (Mura. U., Sgarrella, F. and Ipata, P.L. (1978) J. Biol. Chem. 253, 7905-7909).