The Structure and Functions of Hantavirus Nucleocapsid Protein

Hantaviruses, members of the genus Hantavirus in the Bunyaviridae family, are enveloped single-stranded RNA viruses with tri-segmented genome of negative polarity. In humans, hantaviruses cause two diseases, hemorrhagic fever with renal syndrome (HFRS) and hantavirus pulmonary syndrome (HPS), which vary in severity depending on the causative agent. Each hantavirus is carried by a specific rodent host and is transmitted to humans through excreta of infected rodents. The genome of hantaviruses encodes four structural proteins: the nucleocapsid protein (N), the glycoproteins (Gn and Gc), and the polymerase (L) and also the nonstructural protein (NSs). This thesis deals with the functional characterization of hantavirus N protein with regard to its structure. Structural studies of the N protein have progressed slowly and the crystal structure of the whole protein is still not available, therefore biochemical assays coupled with bioinformatical modeling proved essential for studying N protein structure and functions. Presumably, during RNA encapsidation, the N protein first forms intermediate trimers and then oligomers. First, we investigated the role of N-terminal domain in the N protein oligomerization. The results suggested that the N-terminal region of the N protein forms a coiled-coil, in which two antiparallel alpha helices interact via their hydrophobic seams. Hydrophobic residues L4, I11, L18, L25 and V32 in the first helix and L44, V51, L58 and L65 in the second helix were crucial for stabilizing the structure. The results were consistent with the head-to-head, tail-to-tail model for hantavirus N protein trimerization. We demonstrated that an intact coiled-coil structure of the N terminus is crucial for the oligomerization capacity of the N protein. We also added new details to the head-to-head, tail-to-tail model of trimerization by suggesting that the initial step is based on interaction(s) between intact intra-molecular coiled-coils of the monomers. We further analyzed the importance of charged aa residues located within the coiled-coil for the N protein oligomerization. To predict the interacting surfaces of the monomers we used an upgraded in silico model of the coiled-coil domain that was docked into a trimer. Next the predicted target residues were mutated. The results obtained using the mammalian two-hybrid assay suggested that conserved charged aa residues within the coiled-coil make a substantial contribution to the N protein oligomerization. This contribution probably involves the formation of interacting surfaces of the N monomers and also stabilization of the coiled-coil via intramolecular ionic bridging. We proposed that the tips of the coiled-coils are the first to come into direct contact and thus initiate tight packing of the three monomers into a compact structure. This was in agreement with the previous results showing that an increase in ionic strength abolished the interaction between N protein molecules. We also showed that residues having the strongest effect on the N protein oligomerization are not scattered randomly throughout the coiled-coil 3D model structure, but form clusters. Next we found evidence for the hantaviral N protein interaction with the cytoplasmic tail of the glycoprotein Gn. In order to study this interaction we used the GST pull-down assay in combination with mutagenesis technique. The results demonstrated that intact, properly folded zinc fingers of the Gn protein cytoplasmic tail as well as the middle domain of the N protein (that includes aa residues 80 248 and supposedly carries the RNA-binding domain) are essential for the interaction. Since hantaviruses do not have a matrix protein that mediates the packaging of the viral RNA in other negatve stranded viruses (NSRV), hantaviral RNPs should be involved in a direct interaction with the intraviral domains of the envelope-embedded glycoproteins. By showing the N-Gn interaction we provided the evidence for one of the crucial steps in the virus replication at which RNPs are directed to the site of the virus assembly. Finally we started analysis of the N protein RNA-binding region, which is supposedly located in the middle domain of the N protein molecule. We developed a model for the initial step of RNA-binding by the hantaviral N protein. We hypothesized that the hantaviral N protein possesses two secondary structure elements that initiate the RNA encapsidation. The results suggest that amino acid residues (172-176) presumably act as a hook to catch vRNA and that the positively charged interaction surface (aa residues 144-160) enhances the initial N-RNA interacation. In conclusion, we elucidated new functions of hantavirus N protein. Using in silico modeling we predicted the domain structure of the protein and using experimental techniques showed that each domain is responsible for executing certain function(s). We showed that intact N terminal coiled-coil domain is crucial for oligomerization and charged residues located on its surface form a interaction surface for the N monomers. The middle domain is essential for interaction with the cytoplasmic tail of the Gn protein and RNA binding.

Approximate multipole coefficients of RF ion traps as functions of aperture size

In this study we present approximate analytical expressions for estimating the variation in multipole expansion coefficients as a function of the size of the apertures in the electrodes in axially symmetric (3D) and two-dimensional (2D) ion trap ion traps. Following the approach adopted in our earlier studies which focused on the role of apertures to fields within the traps, here too, the analytical expression we develop is a sum of two terms, A(n,noAperiure), the multipole expansion coefficient for a trap with no apertures and A(n,dueToAperture), the multipole expansion coefficient contributed by the aperture. A(n,noAperture) has been obtained numerically and A(n,dueToAperture) is obtained from the n th derivative of the potential within the trap. The expressions derived have been tested on two 3D geometries and two 2D geometries. These include the quadrupole ion trap (QIT) and the cylindrical ion trap (CIT) for 3D geometries and the linear ion trap (LIT) and the rectilinear ion trap (RIT) for the 2D geometries. Multipole expansion coefficients A(2) to A(12), estimated by our analytical expressions, were compared with the values obtained numerically (using the boundary element method) for aperture sizes varying up to 50% of the trap dimension. In all the plots presented, it is observed that our analytical expression for the variation of multipole expansion coefficients versus aperture size closely follows the trend of the numerical evaluations for the range of aperture sizes considered. The maximum relative percentage errors, which provide an estimate of the deviation of our values from those obtained numerically for each multipole expansion coefficient, are seen to be largely in the range of 10-15%. The leading multipole expansion coefficient, A(2), however, is seen to be estimated very well by our expressions, with most values being within 1% of the numerically determined values, with larger deviations seen for the QIT and the LIT for large aperture sizes. (C) 2010 Elsevier B.V. All rights reserved.

Bound states of two spin-1/2 fermions in a synthetic non-Abelian gauge field

We study the bound states of two spin-1/2 fermions interacting via a contact attraction (characterized by a scattering length) in the singlet channel in three-dimensional space in presence of a uniform non-Abelian gauge field. The configuration of the gauge field that generates a Rashba-type spin-orbit interaction is described by three coupling parameters (lambda(x),lambda(y),lambda(z)). For a generic gauge field configuration, the critical scattering length required for the formation of a bound state is negative, i.e., shifts to the ``BCS side'' of the resonance. Interestingly, we find that there are special high-symmetry configurations (e.g., lambda(x) = lambda(y) = lambda(z)) for which there is a two-body bound state for any scattering length however small and negative. Remarkably, the bound-state wave functions obtained for such configurations have nematic spin structure similar to those found in liquid He-3. Our results show that the BCS-BEC (Bose-Einstein condensation) crossover is drastically affected by the presence of a non-Abelian gauge field. We discuss possible experimental signatures of our findings both at high and low temperatures.

Behavioural and morphological dimorphism of the sexes: an account of two primitively eusocial wasps

In social Hymenoptera, foraging, nest building, brood care and all other colony maintenance functions are carried out by the females while males function solely as reproductives. This asymmetry in social roles of the two sexes has led social insect researchers to focus almost exclusively on the females whereas males have remained relatively neglected. We studied two sympatric, primitively eusocial wasps, Ropalidia marginata and Ropalidia cyathiformis, and compared the morphological and behavioural profiles of males and females. Males of both species are smaller in size and weigh less compared to females. Males of the two species live in the nest for different durations. Borrowing from the ecological literature we use novel methods to compute and compare behavioural diversity and behavioural richness and show that females of both species are behaviourally richer and more diverse than the males.

On an exact populationary model of genetic algorithms

Genetic algorithms (GAs) are search methods that are being employed in a multitude of applications with extremely large search spaces. Recently, there has been considerable interest among GA researchers in understanding and formalizing the working of GAs. In an earlier paper, we have introduced the notion of binomially distributed populations as the central idea behind an exact ''populationary'' model of the large-population dynamics of the GA operators for objective functions called ''functions of unitation.'' In this paper, we extend this populationary model of GA dynamics to a more general class of objective functions called functions of unitation variables. We generalize the notion of a binomially distributed population to a generalized binomially distributed population (GBDP). We show that the effects of selection, crossover, and mutation can be exactly modelled after decomposing the population into GBDPs. Based on this generalized model, we have implemented a GA simulator for functions of two unitation variables-GASIM 2, and the distributions predicted by GASIM 2 match with those obtained from actual GA runs. The generalized populationary model of GA dynamics not only presents a novel and natural way of interpreting the workings of GAs with large populations, but it also provides for an efficient implementation of the model as a GA simulator. (C) Elsevier Science Inc. 1997.

Structure and Mechanistic Insights into Novel Iron-mediated Moonlighting Functions of Human J-protein Cochaperone, Dph4

J-proteins are obligate cochaperones of Hsp70s and stimulate their ATPase activity via the J-domain. Although the functions of J-proteins have been well understood in the context of Hsp70s, their additional co-evolved ``physiological functions'' are still elusive. We report here the solution structure and mechanism of novel iron-mediated functional roles of human Dph4, a type III J-protein playing a vital role in diphthamide biosynthesis and normal development. The NMR structure of Dph4 reveals two domains: a conserved J-domain and a CSL-domain connected via a flexible linker-helix. The linker-helix modulates the conformational flexibility between the two domains, regulating thereby the protein function. Dph4 exhibits a unique ability to bind iron in tetrahedral coordination geometry through cysteines of its CSL-domain. The oxidized Fe-Dph4 shows characteristic UV-visible and electron paramagnetic resonance spectral properties similar to rubredoxins. Iron-bound Dph4 (Fe-Dph4) also undergoes oligomerization, thus potentially functioning as a transient ``iron storage protein,'' thereby regulating the intracellular iron homeostasis. Remarkably, Fe-Dph4 exhibits vital redox and electron carrier activity, which is critical for important metabolic reactions, including diphthamide biosynthesis. Further, we observed that Fe-Dph4 is conformationally better poised to perform Hsp70-dependent functions, thus underlining the significance of iron binding in Dph4. Yeast Jjj3, a functional ortholog of human Dph4 also shows a similar iron-binding property, indicating the conserved nature of iron sequestration across species. Taken together, our findings provide invaluable evidence in favor of additional co-evolved specialized functions of J-proteins, previously not well appreciated.

Junction between surfaces of two topological insulators

We study the properties of a line junction which separates the surfaces of two three-dimensional topological insulators. The velocities of the Dirac electrons on the two surfaces may be unequal and may even have opposite signs. For a time-reversal invariant system, we show that the line junction is characterized by an arbitrary parameter alpha which determines the scattering from the junction. If the surface velocities have the same sign, we show that there can be edge states which propagate along the line junction with a velocity and spin orientation which depend on alpha and the ratio of the velocities. Next, we study what happens if the two surfaces are at an angle phi with respect to each other. We study the scattering and differential conductance through the line junction as functions of phi and alpha. We also find that there are edge states which propagate along the line junction with a velocity and spin orientation which depend on phi. Finally, if the surface velocities have opposite signs, we find that the electrons must transmit into the two-dimensional interface separating the two topological insulators.

Binding of two DNA molecules by type II topoisomerases for decatenation

Topoisomerases (topos) maintain DNA topology and influence DNA transaction processes by catalysing relaxation, supercoiling and decatenation reactions. In the cellular milieu, division of labour between different topos ensures topological homeostasis and control of central processes. In Escherichia coli, DNA gyrase is the principal enzyme that carries out negative supercoiling, while topo IV catalyses decatenation, relaxation and unknotting. DNA gyrase apparently has the daunting task of undertaking both the enzyme functions in mycobacteria, where topo IV is absent. We have shown previously that mycobacterial DNA gyrase is an efficient decatenase. Here, we demonstrate that the strong decatenation property of the enzyme is due to its ability to capture two DNA segments in trans. Topo IV, a strong dedicated decatenase of E. coli, also captures two distinct DNA molecules in a similar manner. In contrast, E. coli DNA gyrase, which is a poor decatenase, does not appear to be able to hold two different DNA molecules in a stable complex. The binding of a second DNA molecule to GyrB/ParE is inhibited by ATP and the non-hydrolysable analogue, AMPPNP, and by the substitution of a prominent positively charged residue in the GyrB N-terminal cavity, suggesting that this binding represents a potential T-segment positioned in the cavity. Thus, after the GyrA/ParC mediated initial DNA capture, GyrB/ParE would bind efficiently to a second DNA in trans to form a T-segment prior to nucleotide binding and closure of the gate during decatenation.

Layerwise decomposition of water dynamics in reverse micelles: A simulation study of two-dimensional infrared spectrum

We present computer simulation study of two-dimensional infrared spectroscopy (2D-IR) of water confined in reverse micelles (RMs) of various sizes. The present study is motivated by the need to understand the altered dynamics of confined water by performing layerwise decomposition of water, with an aim to quantify the relative contributions of different layers water molecules to the calculated 2D-IR spectrum. The 0-1 transition spectra clearly show substantial elongation, due to in-homogeneous broadening and incomplete spectral diffusion, along the diagonal in the surface water layer of different sized RMs. Fitting of the frequency fluctuation correlation functions reveal that the motion of the surface water molecules is sub-diffusive and indicate the constrained nature of their dynamics. This is further supported by two peak nature of the angular analogue of van Hove correlation function. With increasing system size, the water molecules become more diffusive in nature and spectral diffusion almost completes in the central layer of the larger size RMs. Comparisons between experiments and simulations establish the correspondence between the spectral decomposition available in experiments with the spatial decomposition available in simulations. Simulations also allow a quantitative exploration of the relative role of water, sodium ions, and sulfonate head groups in vibrational dephasing. Interestingly, the negative cross correlation between force on oxygen and hydrogen of O-H bond in bulk water significantly decreases in the surface layer of each RM. This negative cross correlation gradually increases in the central water pool with increasing RMs size and this is found to be partly responsible for the faster relaxation rate of water in the central pool. (C) 2013 AIP Publishing LLC.

Structural and functional analysis of two universal stress proteins YdaA and YnaF from Salmonella typhimurium: possible roles in microbial stress tolerance

In many organisms ``Universal Stress Proteins'' CUSPS) are induced in response to a variety of environmental stresses. Here we report the structures of two USPs, YnaF and YdaA from Salmonella typhimurium determined at 1.8 angstrom and 2.4 angstrom resolutions, respectively. YnaF consists of a single USP domain and forms a tetrameric organization stabilized by interactions mediated through chloride ions. YdaA is a larger protein consisting of two tandem USP domains. Two protomers of YdaA associate to form a structure similar to the YnaF tetramer. YdaA showed ATPase activity and an ATP binding motif G-2X-G-9X-G(S/T/N) was found in its C-terminal domain. The residues corresponding to this motif were not conserved in YnaF although YnaF could bind ATP. However, unlike YdaA, YnaF did not hydrolyse ATP in vitro. Disruption of interactions mediated through chloride ions by selected mutations converted YnaF into an ATPase. Residues that might be important for ATP hydrolysis could be identified by comparing the active sites of native and mutant structures. Only the C-terminal domain of YdaA appears to be involved in ATP hydrolysis. The structurally similar N-terminal domain was found to bind a zinc ion near the segment equivalent to the phosphate binding loop of the C-terminal domain. Mass spectrometric analysis showed that YdaA might bind a ligand of approximate molecular weight 800 daltons. Structural comparisons suggest that the ligand, probably related to an intermediate in lipid A biosynthesis, might bind at a site close to the zinc ion. Therefore, the N-terminal domain of YdaA binds zinc and might play a role in lipid metabolism. Thus, USPs appear to perform several distinct functions such as ATP hydrolysis, altering membrane properties and chloride sensing. (C) 2015 Elsevier Inc. All rights reserved.

A perturbation procedure for nonlinear oscillations (The dynamics of two oscillators with weak nonlinear coupling)

This thesis considers in detail the dynamics of two oscillators with weak nonlinear coupling. There are three classes of such problems: non-resonant, where the Poincaré procedure is valid to the order considered; weakly resonant, where the Poincaré procedure breaks down because small divisors appear (but do not affect the O(1) term) and strongly resonant, where small divisors appear and lead to O(1) corrections. A perturbation method based on Cole's two-timing procedure is introduced. It avoids the small divisor problem in a straightforward manner, gives accurate answers which are valid for long times, and appears capable of handling all three types of problems with no change in the basic approach.

One example of each type is studied with the aid of this procedure: for the nonresonant case the answer is equivalent to the Poincaré result; for the weakly resonant case the analytic form of the answer is found to depend (smoothly) on the difference between the initial energies of the two oscillators; for the strongly resonant case we find that the amplitudes of the two oscillators vary slowly with time as elliptic functions of ϵ t, where ϵ is the (small) coupling parameter.

Our results suggest that, as one might expect, the dynamical behavior of such systems varies smoothly with changes in the ratio of the fundamental frequencies of the two oscillators. Thus the pathological behavior of Whittaker's adelphic integrals as the frequency ratio is varied appears to be due to the fact that Whittaker ignored the small divisor problem. The energy sharing properties of these systems appear to depend strongly on the initial conditions, so that the systems not ergodic.

The perturbation procedure appears to be applicable to a wide variety of other problems in addition to those considered here.

Neural and computational representations of decision variables

These studies explore how, where, and when representations of variables critical to decision-making are represented in the brain. In order to produce a decision, humans must first determine the relevant stimuli, actions, and possible outcomes before applying an algorithm that will select an action from those available. When choosing amongst alternative stimuli, the framework of value-based decision-making proposes that values are assigned to the stimuli and that these values are then compared in an abstract “value space” in order to produce a decision. Despite much progress, in particular regarding the pinpointing of ventromedial prefrontal cortex (vmPFC) as a region that encodes the value, many basic questions remain. In Chapter 2, I show that distributed BOLD signaling in vmPFC represents the value of stimuli under consideration in a manner that is independent of the type of stimulus it is. Thus the open question of whether value is represented in abstraction, a key tenet of value-based decision-making, is confirmed. However, I also show that stimulus-dependent value representations are also present in the brain during decision-making and suggest a potential neural pathway for stimulus-to-value transformations that integrates these two results.

More broadly speaking, there is both neural and behavioral evidence that two distinct control systems are at work during action selection. These two systems compose the “goal-directed system”, which selects actions based on an internal model of the environment, and the “habitual” system, which generates responses based on antecedent stimuli only. Computational characterizations of these two systems imply that they have different informational requirements in terms of input stimuli, actions, and possible outcomes. Associative learning theory predicts that the habitual system should utilize stimulus and action information only, while goal-directed behavior requires that outcomes as well as stimuli and actions be processed. In Chapter 3, I test whether areas of the brain hypothesized to be involved in habitual versus goal-directed control represent the corresponding theorized variables.

The question of whether one or both of these neural systems drives Pavlovian conditioning is less well-studied. Chapter 4 describes an experiment in which subjects were scanned while engaged in a Pavlovian task with a simple non-trivial structure. After comparing a variety of model-based and model-free learning algorithms (thought to underpin goal-directed and habitual decision-making, respectively), it was found that subjects’ reaction times were better explained by a model-based system. In addition, neural signaling of precision, a variable based on a representation of a world model, was found in the amygdala. These data indicate that the influence of model-based representations of the environment can extend even to the most basic learning processes.

Knowledge of the state of hidden variables in an environment is required for optimal inference regarding the abstract decision structure of a given environment and therefore can be crucial to decision-making in a wide range of situations. Inferring the state of an abstract variable requires the generation and manipulation of an internal representation of beliefs over the values of the hidden variable. In Chapter 5, I describe behavioral and neural results regarding the learning strategies employed by human subjects in a hierarchical state-estimation task. In particular, a comprehensive model fit and comparison process pointed to the use of "belief thresholding". This implies that subjects tended to eliminate low-probability hypotheses regarding the state of the environment from their internal model and ceased to update the corresponding variables. Thus, in concert with incremental Bayesian learning, humans explicitly manipulate their internal model of the generative process during hierarchical inference consistent with a serial hypothesis testing strategy.

Proposal and validation of methodologies for the categorisation of continuous variables in the development of prediction models

158 p.

Application of neural networks to the ranking of perinatal variables influencing birthweight.

In this paper we compare Multi-Layer Perceptrons (a neural network type) with Multivariate Linear Regression in predicting birthweight from nine perinatal variables which are thought to be related. Results show, that seven of the nine variables, i.e., gestational age, mother's body-mass index (BMI), sex of the baby, mother's height, smoking, parity and gravidity, are related to birthweight. We found no significant relationship between birthweight and each of the two variables, i.e., maternal age and social class.

Fisheries databases for assessing habitat functions of Cape Romain NWR

The US Fish and Wildlife Service Cape Romain National Wildlife Refuge (CRNWR) and the Center for Coastal Environmental Health and Biomolecular Research (CCEHBR) at Charleston are interested in assessing the status of our coastal resources in light of increased coastal development and recreational use. Through an Interagency Agreement (FWS #1448-40181-00-H-001), an ecological characterization was undertaken to describe the status of and potential impacts to resources at CRNWR. This report describes historic fisheries-independent or non-commercial data relevant to CRNWR that can be used to evaluate the role of the Refuge as habitat for nearshore and offshore fish species. The purpose of this document is two-fold, first to give resource managers an understanding of fisheries data that have been collected over the years and, second, to illustrate how these data can be applied to address specific management issues. This report provides an overview of historic fisheries data collected along the southeast coast, as well as basic summaries of that data relevant to CRNWR, indicating how these data can be used to address specific questions of interest to Refuge managers and biologists.