Hybrid Modelling of Titan's Interaction with the Magnetosphere of Saturn

In this dissertation we study the interaction between Saturn's moon Titan and the magnetospheric plasma and magnetic field. The method of research is a three-dimensional computer simulation model, that is used to simulate this interaction. The simulation model used is a hybrid model. Hybrid models enable individual tracking or tracing of ions and also take into account the particle motion in the propagation of the electromagnetic fields. The hybrid model has been developed at the Finnish Meteorological Institute. This thesis gives a general description of the effects that the solar wind has on Earth and other planets of our solar system. Planetary satellites can also have similar interactions with the solar wind but also with the plasma flows of planetary magnetospheres. Titan is clearly the largest among the satellites of Saturn and also the only known satellite with a dense atmosphere. It is the atmosphere that makes Titan's plasma interaction with the magnetosphere of Saturn so unique. Nevertheless, comparisons with the plasma interactions of other solar system bodies are valuable. Detecting charged plasma particles requires in situ measurements obtainable through scientific spacecraft. The Cassini mission has been one of the most remarkable international efforts in space science. Since 2004 the measurements and images obtained from instruments onboard the Cassini spacecraft have increased the scientific knowledge of Saturn as well as its satellites and magnetosphere in a way no one was probably able to predict. The current level of science on Titan is practically unthinkable without the Cassini mission. Many of the observations by Cassini instrument teams have influenced this research both the direct measurements of Titan as well as observations of its plasma environment. The theoretical principles of the hybrid modelling approach are presented in connection to the broader context of plasma simulations. The developed hybrid model is described in detail: e.g. the way the equations of the hybrid model are solved is shown explicitly. Several simulation techniques, such as the grid structure and various boundary conditions, are discussed in detail as well. The testing and monitoring of simulation runs is presented as an essential routine when running sophisticated and complex models. Several significant improvements of the model, that are in preparation, are also discussed. A main part of this dissertation are four scientific articles based on the results of the Titan model. The Titan model developed during the course of the Ph.D. research has been shown to be an important tool to understand Titan's plasma interaction. One reason for this is that the structures of the magnetic field around Titan are very much three-dimensional. The simulation results give a general picture of the magnetic fields in the vicinity of Titan. The magnetic fine structure of Titan's wake as seen in the simulations seems connected to Alfvén waves an important wave mode in space plasmas. The particle escape from Titan is also a major part of these studies. Our simulations show a bending or turning of Titan's ionotail that we have shown to be a direct result of the basic principles in plasma physics. Furthermore, the ion flux from the magnetosphere of Saturn into Titan's upper atmosphere has been studied. The modelled ion flux has asymmetries that would likely have a large impact in the heating in different parts of Titan's upper atmosphere.

Modular design of synthetic protein mimics. Crystal structure of two seven-residue helical peptide segments linked by .epsilon.-aminocaproic acid

Two seven-residue helical segments, Val-Ala-Leu-Aib-Val-Ala-Leu, were linked synthetically with an epsilon-aminocaproic acid (Acp) linker with the intention of making a stable antiparallel helix-helix motif. The crystal structure of the linked peptide Boc-Val-Ala-Leu-Aib-Val-Ala-Leu-Acp-Val-Ala-Leu-Aib-Val-Ala-Leu-OMe (1) shows the two helices displaced laterally from each other by the linker, but the linker has not folded the molecule into a close-packed antiparallel conformation. Two strong intermolecular NH...O = C hydrogen bonds are formed between the top of the lower helix of one molecule and the bottom of the upper helix in a laterally adjacent molecule to give the appearance of an extended single helix. The composite peptide with Boc and OMe end groups, C76H137N15O18.H2O, crystallize in space group P2(1) with a = 8.802 (1) angstrom, b = 20.409 (4) angstrom, c = 26.315 (3) angstrom, and beta = 90.72 (1)degrees; overall agreement R = 7.86% for 5030 observed reflections (\F(o)\ > 3-sigma(F)); resolution = 0.93 angstrom. Limited evidence for a more compact conformation in solution consistent with an antiparallel helix arrangement is obtained by comparison of the HPLC retention times and CD spectra of peptide 1 with well-characterized continuous helices of similar length and sequence.

Stereochemical analysis of the antigenic tip of the V3 loop peptide of HIV-1 gp120

A novel multiple turn conformation has been observed for a segment GPGRAFY in the crystal structure of a complex of HIV-1 gp120 V3 loop peptide with the Fab fragment of a neutralizing antibody [Ghiara ct al. (1994) Science 264, 82-85]. A structural motif has been defined for the peptide segment, employing idealized backbone conformations characterized by ranges of virtual C-alpha torsion angles and bond angles. A search of 122 high-resolution protein crystal structures has permitted identification of 24 examples of similar structural motifs. Two major conformational families have been identified, which differ primarily in the conformation at residue 3. The observed conformation at residue 3 in family 1 is left-handed helical (alpha(L)) and that in family 2 is right-handed helical (alpha(R)). Of the 10 examples in family 1, 9 examples have Gly residues at position 3. Of the 12 examples in family 2, 7 examples have Asn/Asp at position 3. Computer modeling of the V3 loop tip sequence using the two backbone conformational families as starting points leads to minimum-energy conformations in which antigenically important side-chains occupy similar spatial arrangements. This stereochemical analysis of the V3 loop tip sequence suggests a rational basis for the design of synthetic analog peptides for use as viral antagonists or synthetic antigens. (C) Munksgaard 1995.

Crystal structure of dimeric FabZ of Plasmodium falciparum reveals conformational switching to active hexamers by peptide flips

The crystal structure of beta-hydroxyacyl acyl carrier protein dehydratase of Plasmodium falciparum (PfFabZ) has been determined at a resolution of 2.4 angstrom. PfFabZ has been found to exist as a homodimer (d-PfFabZ) in the crystals of the present study in contrast to the reported hexameric form (h-PfFabZ) which is a trimer of dimers crystallized in a different condition. The catalytic sites of this enzyme are located in deep narrow tunnel-shaped pockets formed at the dimer interface. A histidine residue from one subunit of the dimer and a glutamate residue from the other subunit lining the tunnel form the catalytic dyad in the reported crystal structures. While the position of glutamate remains unaltered in the crystal structure of d-PffabZ compared to that in b-PfFabZ, the histidine residue takes up an entirely different conformation and moves away from the tunnel leading to a His-Phe cis-trans peptide flip at the histidine residue. In addition, a loop in the vicinity has been observed to undergo a similar flip at a Tyr-Pro peptide bond. These alterations not only prevent the formation of a hexamer but also distort the active site geometry resulting in a dimeric form of FabZ that is incapable of substrate binding. The dimeric state and an altered catalytic site architecture make d-PfFabZ distinctly different from the FabZ structures described so far. Dynamic light scattering and size exclusion chromatographic studies clearly indicate a pH-related switching of the dimers to active hexamers. (c) 2006 Federation of European Biochemical Societies. Published by Elsevier B.V. All rights reserv.

Structure determination and biochemical studies on Bacillus stearothermophilus E53Q serine hydroxymethyltransferase and its complexes provide insights on function and enzyme memory

Serine hydroxymethyltransferase (SHMT) belongs to the alpha-family of pyridoxal 5'-phosphate-dependent enzymes and catalyzes the reversible conversion of L-Ser and etrahydrofolate to Gly and 5,10-methylene tetrahydrofolate. 5,10-Methylene tetrahydrofolate serves as a source of one-carbon fragment in many biological processes. SHMT also catalyzes the tetrahydrofolate-independent conversion of L-allo-Thr to Gly and acetaldehyde. The crystal structure of Bacillus stearothermophilus SHMT (bsSHMT) suggested that E53 interacts with the substrate, L-Ser and etrahydrofolate. To elucidate the role of E53, it was mutated to Q and structural and biochemical studies were carried out with the mutant enzyme. The internal aldimine structure of E53QbsSHMT was similar to that of the except for significant changes at Q53, Y60 and Y61. The wild-type enzyme, carboxyl of Gly and side chain of L-Ser were in two conformations in the respective external aldimine structures. The mutant enzyme was completely inactive for tetrahydrofolate-depen dent cleavage of L-Ser, whereas there was a 1.5-fold increase in the rate of tetrahydrofolate-independent reaction with L-allo-Thr. The results obtained from these studies suggest that E53 plays an essential role in tetrahydrofolate/5-formyl tetrahydrofolate binding and in the proper positioning of C beta of L-Ser for direct attack by N5 of tetrahydrofolate. Most interestingly, the structure of the complex obtained by cocrystallization of E53QbsSHMT with Gly and 5-formyl tetrahydrofolate revealed the gem-diamine form of pyridoxal 5'-phosphate bound to Gly and active site Lys. However, density for 5-formyl tetrahydrofolate was not observed. Gly carboxylate was in a single conformation, whereas pyridoxal 5'-phosphate had two distinct conformations. The differences between the structures of this complex and Gly external aldimine suggest that the changes induced by initial binding of 5-formyl tetrahydrofolate are retained even though 5-formyl tetrahydrofolate is absent in the final structure. Spectral studies carried out with this mutant enzyme also suggest that 5-formyl tetrahydrofolate binds to the E53QbsSHMT-Gly complex forming a quinonoid intermediate and falls off within 4 h of dialysis, leaving behind the mutant enzyme in the gemdiamine form. This is the first report to provide direct evidence for enzyme memory based on the crystal structure of enzyme complexes.

Suppression of spinel formation to induce reversible thermal behavior in the layered double hydroxides (LDHs) of Co with Al, Fe, Ga, and In

The layered double hydroxides (LDHs) of Co with trivalent cations decompose irreversibly to yield oxides with the spinel structure. Spinel formation is aided by the oxidation of Co(II) to Co(III) in the ambient atmosphere. When the decomposition is carried out under N-2, the oxidation of Co(II) is suppressed, and the resulting oxide has the rock salt structure. Thus, the Co-Al-CO32-/Cl- LDHs yield oxides of the type Co1- Al-x(2x/3)rectangle O-x/3, which are highly metastable, given the large defect concentration. This defect oxide rapidly reverts back to the original hydroxide on soaking in a Na2CO3 solution. Interlayer NO3- anions, on the other hand, decompose generating a highly oxidizing atmosphere, whereby the Co-Al-NO3- LDH decomposes to form the spinel phase even in a N-2 atmosphere. The oxide with the defect rock salt structure formed by the thermal decomposition of the Co-Fe-CO32- LDH under N2, on soaking in a Na2CO3 solution, follows a different kinetic pathway and undergoes a solution transformation into the inverse spinel Co(Co, Fe)(2)O-4. Fe3+ has a low octahedral crystal field stabilization energy and therefore prefers the tetrahedral coordination offered by the structure of the inverse spinel rather than the octahedral coordination of the parent LDH. Similar considerations do not hold in the case of Ga- and In-containing LDHs, given the considerable barriers to the diffusion of M3+ (M=Ga, In) from octahedral to tetrahedral sites owing to their large size. Consequently, the In-containing oxide residue reverts back to the parent hydroxide, whereas this reconstruction is partial in the case of the Ga-containing oxide. These studies show that the reversible thermal behavior offers a competing kinetic pathway to spinel formation. Suppression of the latter induces the reversible behavior in an LDH that otherwise decomposes irreversibly to the spinel.

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.

The structure of thiostrepton

Much of the chemical structure of thiostrepton, a sulphur containing metabolic product of the microorganism Streptomyces azureus, has been determined by X-ray crystallographic techniques.

Detection of a periodic structure hidden in random background: the role of signal amplitude in the matched filter detection method

The matched filter method for detecting a periodic structure on a surface hidden behind randomness is known to detect up to (r(0)/Lambda) gt;= 0.11, where r(0) is the coherence length of light on scattering from the rough part and 3 is the wavelength of the periodic part of the surface-the above limit being much lower than what is allowed by conventional detection methods. The primary goal of this technique is the detection and characterization of the periodic structure hidden behind randomness without the use of any complicated experimental or computational procedures. This paper examines this detection procedure for various values of the amplitude a of the periodic part beginning from a = 0 to small finite values of a. We thus address the importance of the following quantities: `(a)lambda) `, which scales the amplitude of the periodic part with the wavelength of light, and (r(0))Lambda),in determining the detectability of the intensity peaks.

Extended self-similarity works for the Burgers equation and why

Extended self-similarity (ESS), a procedure that remarkably extends the range of scaling for structure functions in Navier-Stokes turbulence and thus allows improved determination of intermittency exponents, has never been fully explained. We show that ESS applies to Burgers turbulence at high Reynolds numbers and we give the theoretical explanation of the numerically observed improved scaling at both the IR and UV end, in total a gain of about three quarters of a decade: there is a reduction of subdominant contributions to scaling when going from the standard structure function representation to the ESS representation. We conjecture that a similar situation holds for three-dimensional incompressible turbulence and suggest ways of capturing subdominant contributions to scaling.

Mechanism of crystal structure transformations. Part 3. - Factors affecting the anatase-rutile transformation

The particle size and crystallite size of anatase increase markedly in the region of the crystal structure transformation. The unit cell of anatase seems to expand prior to the transformation to rutile. This expansion has been attributed to a displacive transformation of the type defined by Buerger. Smaller particle size and larger surface area seem to favour the transformation. The kinetics of the transformation of anatase prepared by the hydrolysis of titanium sulphate have been studied at different temperatures and are found to be considerably different from the kinetics of the transformation of pure anatase. The transformation becomes immeasurably slow below ∼695 ± 10°C compared to ∼610°C for pure anatase. An induction period is observed in the transformation of anatase obtained from sulphate hydrolysis and the duration decreases with increase in temperature. The activation energy is ∼120 kcal/mole, a value higher than that for the pure anatase-rutile transformation. The results have been interpreted in terms of the relative rates of nucleation and propagation processes. The activation energy for the nucleation process seems to be much larger than for the propagation process. The kinetics of the transformation of anatase samples doped with different amounts of sulphate ion impurity have also been studied and the transformation is found to be progressively decelerated with increase in the impurity concentration. The energy of activation for the transformation appears to increase progressively with increase in impurity concentration.

Structure of Ce1-xSnxO2 and its relation to oxygen storage property from first-principles analysis

CeO2-SnO2 solid solution has been reported to possess high oxygen storage/release property which possibly originates from local structural distortion. We have performed first-principles based density functional calculations of Ce1-xSnxO2 structure (x=0, 0.25, 0.5, 1) to understand its structural stability in fluorite in comparison to rutile structure of the other end-member SnO2, and studied the local structural distortion induced by the dopant Sn ion. Analysis of relative energies of fluorite and rutile phases of CeO2, SnO2, and Ce1-xSnxO2 indicates that fluorite structure is the most stable for Ce1-xSnxO2 solid solution. An analysis of local structural distortions reflected in phonon dispersion show that SnO2 in fluorite structure is highly unstable while CeO2 in rutile structure is only weakly unstable. Thus, Sn in Ce1-xSnxO2-fluorite structure is associated with high local structural distortion whereas Ce in Ce1-xSnxO2-rutile structure, if formed, will show only marginal local distortion. Determination of M-O (M=Ce or Sn) bond lengths and analysis of Born effective charges for the optimized structure of Ce1-xSnxO2 show that local coordination of these cations changes from ideal eightfold coordination expected of fluorite lattice to 4+4 coordination, leading to generation of long and short Ce-O and Sn-O bonds in the doped structure. Bond valence analyses for all ions show the presence of oxygen with bond valence similar to 1.84. These weakly bonded oxygen ions are relevant for enhanced oxygen storage/release properties observed in Ce1-xSnxO2 solid solution. (C) 2010 American Institute of Physics.

Structure-function relationships of icosahedral plant

X-ray diffraction studies on single crystals of a few viruses have led to the elucidation of their three dimensional structure at near atomic resolution. Both the tertiary structure of the coat protein subunit and the quaternary organization of the icosahedral capsid in these viruses are remarkably similar. These studies have led to a critical re-examination of the structural principles in the architecture of isometric viruses and suggestions of alternative mechanisms of assembly. Apart from their role in the assembly of the virus particle, the coat proteins of certian viruses have been shown to inhibit the replication of the cognate RNA leading to cross-protection. The coat protein amino acid sequence and the genomic sequence of several spherical plant RNA viruses have been determined in the last decade. Experimental data on the mechanisms of uncoating, gene expression and replication of several classes of viruses have also become available. The function of the non-structural proteins of some viruses have been determined. This rapid progress has provided a wealth of information on several key steps in the life cycle of RNA viruses. The function of the viral coat protein, capsid architecture, assembly and disassembly and replication of isometric RNA plant viruses are discussed in the light of this accumulated knowledge.

Protein sequence design on the basis of topology optimization techniques -Using continuous modeling of discrete amino acid types

The notion of optimization is inherent in protein design. A long linear chain of twenty types of amino acid residues are known to fold to a 3-D conformation that minimizes the combined inter-residue energy interactions. There are two distinct protein design problems, viz. predicting the folded structure from a given sequence of amino acid monomers (folding problem) and determining a sequence for a given folded structure (inverse folding problem). These two problems have much similarity to engineering structural analysis and structural optimization problems respectively. In the folding problem, a protein chain with a given sequence folds to a conformation, called a native state, which has a unique global minimum energy value when compared to all other unfolded conformations. This involves a search in the conformation space. This is somewhat akin to the principle of minimum potential energy that determines the deformed static equilibrium configuration of an elastic structure of given topology, shape, and size that is subjected to certain boundary conditions. In the inverse-folding problem, one has to design a sequence with some objectives (having a specific feature of the folded structure, docking with another protein, etc.) and constraints (sequence being fixed in some portion, a particular composition of amino acid types, etc.) while obtaining a sequence that would fold to the desired conformation satisfying the criteria of folding. This requires a search in the sequence space. This is similar to structural optimization in the design-variable space wherein a certain feature of structural response is optimized subject to some constraints while satisfying the governing static or dynamic equilibrium equations. Based on this similarity, in this work we apply the topology optimization methods to protein design, discuss modeling issues and present some initial results.

Ideal Structure of the Silver Code

The Silver code has captured a lot of attention in the recent past,because of its nice structure and fast decodability. In their recent paper, Hollanti et al. show that the Silver code forms a subset of the natural order of a particular cyclic division algebra (CDA). In this paper, the algebraic structure of this subset is characterized. It is shown that the Silver code is not an ideal in the natural order but a right ideal generated by two elements in a particular order of this CDA. The exact minimum determinant of the normalized Silver code is computed using the ideal structure of the code. The construction of Silver code is then extended to CDAs over other number fields.