Interatomic potentials for fusion reactor material simulations

Fusion energy is a clean and safe solution for the intricate question of how to produce non-polluting and sustainable energy for the constantly growing population. The fusion process does not result in any harmful waste or green-house gases, since small amounts of helium is the only bi-product that is produced when using the hydrogen isotopes deuterium and tritium as fuel. Moreover, deuterium is abundant in seawater and tritium can be bred from lithium, a common metal in the Earth's crust, rendering the fuel reservoirs practically bottomless. Due to its enormous mass, the Sun has been able to utilize fusion as its main energy source ever since it was born. But here on Earth, we must find other means to achieve the same. Inertial fusion involving powerful lasers and thermonuclear fusion employing extreme temperatures are examples of successful methods. However, these have yet to produce more energy than they consume. In thermonuclear fusion, the fuel is held inside a tokamak, which is a doughnut-shaped chamber with strong magnets wrapped around it. Once the fuel is heated up, it is controlled with the help of these magnets, since the required temperatures (over 100 million degrees C) will separate the electrons from the nuclei, forming a plasma. Once the fusion reactions occur, excess binding energy is released as energetic neutrons, which are absorbed in water in order to produce steam that runs turbines. Keeping the power losses from the plasma low, thus allowing for a high number of reactions, is a challenge. Another challenge is related to the reactor materials, since the confinement of the plasma particles is not perfect, resulting in particle bombardment of the reactor walls and structures. Material erosion and activation as well as plasma contamination are expected. Adding to this, the high energy neutrons will cause radiation damage in the materials, causing, for instance, swelling and embrittlement. In this thesis, the behaviour of a material situated in a fusion reactor was studied using molecular dynamics simulations. Simulations of processes in the next generation fusion reactor ITER include the reactor materials beryllium, carbon and tungsten as well as the plasma hydrogen isotopes. This means that interaction models, {\it i.e. interatomic potentials}, for this complicated quaternary system are needed. The task of finding such potentials is nonetheless nearly at its end, since models for the beryllium-carbon-hydrogen interactions were constructed in this thesis and as a continuation of that work, a beryllium-tungsten model is under development. These potentials are combinable with the earlier tungsten-carbon-hydrogen ones. The potentials were used to explain the chemical sputtering of beryllium due to deuterium plasma exposure. During experiments, a large fraction of the sputtered beryllium atoms were observed to be released as BeD molecules, and the simulations identified the swift chemical sputtering mechanism, previously not believed to be important in metals, as the underlying mechanism. Radiation damage in the reactor structural materials vanadium, iron and iron chromium, as well as in the wall material tungsten and the mixed alloy tungsten carbide, was also studied in this thesis. Interatomic potentials for vanadium, tungsten and iron were modified to be better suited for simulating collision cascades that are formed during particle irradiation, and the potential features affecting the resulting primary damage were identified. Including the often neglected electronic effects in the simulations was also shown to have an impact on the damage. With proper tuning of the electron-phonon interaction strength, experimentally measured quantities related to ion-beam mixing in iron could be reproduced. The damage in tungsten carbide alloys showed elemental asymmetry, as the major part of the damage consisted of carbon defects. On the other hand, modelling the damage in the iron chromium alloy, essentially representing steel, showed that small additions of chromium do not noticeably affect the primary damage in iron. Since a complete assessment of the response of a material in a future full-scale fusion reactor is not achievable using only experimental techniques, molecular dynamics simulations are of vital help. This thesis has not only provided insight into complicated reactor processes and improved current methods, but also offered tools for further simulations. It is therefore an important step towards making fusion energy more than a future goal.

Dynamic nature of proteins - interpretation of residual dipolar couplings

Protein conformations and dynamics can be studied by nuclear magnetic resonance spectroscopy using dilute liquid crystalline samples. This work clarifies the interpretation of residual dipolar coupling data yielded by the experiments. It was discovered that unfolded proteins without any additional structure beyond that of a mere polypeptide chain exhibit residual dipolar couplings. Also, it was found that molecular dynamics induce fluctuations in the molecular alignment and doing so affect residual dipolar couplings. The finding clarified the origins of low order parameter values observed earlier. The work required the development of new analytical and computational methods for the prediction of intrinsic residual dipolar coupling profiles for unfolded proteins. The presented characteristic chain model is able to reproduce the general trend of experimental residual dipolar couplings for denatured proteins. The details of experimental residual dipolar coupling profiles are beyond the analytical model, but improvements are proposed to achieve greater accuracy. A computational method for rapid prediction of unfolded protein residual dipolar couplings was also developed. Protein dynamics were shown to modulate the effective molecular alignment in a dilute liquid crystalline medium. The effects were investigated from experimental and molecular dynamics generated conformational ensembles of folded proteins. It was noted that dynamics induced alignment is significant especially for the interpretation of molecular dynamics in small, globular proteins. A method of correction was presented. Residual dipolar couplings offer an attractive possibility for the direct observation of protein conformational preferences and dynamics. The presented models and methods of analysis provide significant advances in the interpretation of residual dipolar coupling data from proteins.

The GAF Domain of the cGMP-Binding, cGMP-Specific Phosphodiesterase (PDE5) Is a Sensor and a Sink for cGMP

We describe here a novel sensor for cGMP based on the GAF domain of the cGMP-binding, cGMP-specific phosphodiesterase 5 (PDE5) using bioluminescence resonance energy transfer (BRET). The wild type GAFa domain, capable of binding cGMP with high affinity, and a mutant (GAFaF163A) unable to bind cGMP were cloned as fusions between GFP and Rluc for BRET2 assays. BRET2 ratios of the wild type GAFa fusion protein, but not GAFaF163A, increased in the presence of cGMP but not cAMP. Higher basal BRET2 ratios were observed in cells expressing the wild type GAFa domain than in cells expressing GAFaF163A. This was correlated with elevated basal intracellular levels of cGMP, indicating that the GAF domain could act as a sink for cGMP. The tandem GAF domains in full length PDE5 could also sequester cGMP when the catalytic activity of PDE5 was inhibited. Therefore, these results describe a cGMP sensor utilizing BRET2 technology and experimentally demonstrate the reservoir of cGMP that can be present in cells that express cGMP-binding GAF domain-containing proteins. PDE5 is the target for the anti-impotence drug sildenafil citrate; therefore, this GAF-BRET2 sensor could be used for the identification of novel compounds that inhibit cGMP binding to the GAF domain, thereby regulating PDE5 catalytic activity.

Proteomic profiling of forskolin-induced differentiated BeWo cells: an in-vitro model of cytotrophoblast differentiation

Placental trophoblastic differentiation is characterized by the fusion of monolayer cytotrophoblasts into syncytiotrophoblasts. During this process of differentiation, several morphological and biochemical changes are known to occur, and this model has been employed to investigate the changes that occur at the gene and protein level during differentiation. Using the sensitive technique of proteomics [two-dimensional gel electrophoresis (2DGE)], changes in protein profile were evaluated in the control and forskolin-induced differentiated cells of trophoblastic choriocarcinoma BeWo cell line. Several proteins were differentially expressed in control and differentiated cells. Four major proteins were up-regulated as assessed by silver staining, and were further characterized as c-h-ras p 21 (phosphorylated), retinoblastoma susceptibility protein I and integrase interactor protein 1. These proteins are known to play an important role in growth arrest of cells, and thus may play a role in initiating the process of differentiation.

Translational fusion of two beta-subunits of human chorionic gonadotropin results in production of a novel antagonist of

The strategy of translationally fusing the alpha-and beta-subunits of human chorionic gonadotropin (hCG) into a single-chain molecule has been used to produce novel analogs of hCG. Previously we reported expression of a biologically active singlechain analog hCG alpha beta expressed using Pichia expression system. Using the same expression system, another analog, in which the alpha-subunit was replaced with the second beta-subunit, was expressed (hCG beta beta) and purified. hCG beta beta could bind to LH receptor with an affinity three times lower than that of hCG but failed to elicit any response. However, it could inhibit response to the hormone in vitro in a dose- dependent manner. Furthermore, it inhibited response to hCG in vivo indicating the antagonistic nature of the analog. However, it was unable inhibit human FSH binding or response to human FSH, indicating the specificity of the effect. Characterization of hCG alpha beta and hCG beta beta using immunological tools showed alterations in the conformation of some of the epitopes, whereas others were unaltered. Unlike hCG, hCG beta beta interacts with two LH receptor molecules. These studies demonstrate that the presence of the second beta-subunit in the single-chain molecule generated a structure that can be recognized by the receptor. However, due to the absence of alpha-subunit, the molecule is unable to elicit response. The strategy of fusing two beta-subunits of glycoprotein hormones can be used to produce antagonists of these hormones.

Fusion of Multisensor Data:Review and Comparative Analysis

Image fusion is a formal framework which is expressed as means and tools for the alliance of multisensor, multitemporal, and multiresolution data. Multisource data vary in spectral, spatial and temporal resolutions necessitating advanced analytical or numerical techniques for enhanced interpretation capabilities. This paper reviews seven pixel based image fusion techniques - intensity-hue-saturation, brovey, high pass filter (HPF), high pass modulation (HPM), principal component analysis, fourier transform and correspondence analysis.Validation of these techniques on IKONOS data (Panchromatic band at I m spatial resolution and Multispectral 4 bands at 4 in spatial resolution) reveal that HPF and HPM methods synthesises the images closest to those the corresponding multisensors would observe at the high resolution level.

Multiscale study on hydrogen mobility in metallic fusion divertor material

For achieving efficient fusion energy production, the plasma-facing wall materials of the fusion reactor should ensure long time operation. In the next step fusion device, ITER, the first wall region facing the highest heat and particle load, i.e. the divertor area, will mainly consist of tiles based on tungsten. During the reactor operation, the tungsten material is slowly but inevitably saturated with tritium. Tritium is the relatively short-lived hydrogen isotope used in the fusion reaction. The amount of tritium retained in the wall materials should be minimized and its recycling back to the plasma must be unrestrained, otherwise it cannot be used for fueling the plasma. A very expensive and thus economically not viable solution is to replace the first walls quite often. A better solution is to heat the walls to temperatures where tritium is released. Unfortunately, the exact mechanisms of hydrogen release in tungsten are not known. In this thesis both experimental and computational methods have been used for studying the release and retention of hydrogen in tungsten. The experimental work consists of hydrogen implantations into pure polycrystalline tungsten, the determination of the hydrogen concentrations using ion beam analyses (IBA) and monitoring the out-diffused hydrogen gas with thermodesorption spectrometry (TDS) as the tungsten samples are heated at elevated temperatures. Combining IBA methods with TDS, the retained amount of hydrogen is obtained as well as the temperatures needed for the hydrogen release. With computational methods the hydrogen-defect interactions and implantation-induced irradiation damage can be examined at the atomic level. The method of multiscale modelling combines the results obtained from computational methodologies applicable at different length and time scales. Electron density functional theory calculations were used for determining the energetics of the elementary processes of hydrogen in tungsten, such as diffusivity and trapping to vacancies and surfaces. Results from the energetics of pure tungsten defects were used in the development of an classical bond-order potential for describing the tungsten defects to be used in molecular dynamics simulations. The developed potential was utilized in determination of the defect clustering and annihilation properties. These results were further employed in binary collision and rate theory calculations to determine the evolution of large defect clusters that trap hydrogen in the course of implantation. The computational results for the defect and trapped hydrogen concentrations were successfully compared with the experimental results. With the aforedescribed multiscale analysis the experimental results within this thesis and found in the literature were explained both quantitatively and qualitatively.

Computer simulation of H and He effects in fusion reactor materials

Fusion power is an appealing source of clean and abundant energy. The radiation resistance of reactor materials is one of the greatest obstacles on the path towards commercial fusion power. These materials are subject to a harsh radiation environment, and cannot fail mechanically or contaminate the fusion plasma. Moreover, for a power plant to be economically viable, the reactor materials must withstand long operation times, with little maintenance. The fusion reactor materials will contain hydrogen and helium, due to deposition from the plasma and nuclear reactions because of energetic neutron irradiation. The first wall divertor materials, carbon and tungsten in existing and planned test reactors, will be subject to intense bombardment of low energy deuterium and helium, which erodes and modifies the surface. All reactor materials, including the structural steel, will suffer irradiation of high energy neutrons, causing displacement cascade damage. Molecular dynamics simulation is a valuable tool for studying irradiation phenomena, such as surface bombardment and the onset of primary damage due to displacement cascades. The governing mechanisms are on the atomic level, and hence not easily studied experimentally. In order to model materials, interatomic potentials are needed to describe the interaction between the atoms. In this thesis, new interatomic potentials were developed for the tungsten-carbon-hydrogen system and for iron-helium and chromium-helium. Thus, the study of previously inaccessible systems was made possible, in particular the effect of H and He on radiation damage. The potentials were based on experimental and ab initio data from the literature, as well as density-functional theory calculations performed in this work. As a model for ferritic steel, iron-chromium with 10% Cr was studied. The difference between Fe and FeCr was shown to be negligible for threshold displacement energies. The properties of small He and He-vacancy clusters in Fe and FeCr were also investigated. The clusters were found to be more mobile and dissociate more rapidly than previously assumed, and the effect of Cr was small. The primary damage formed by displacement cascades was found to be heavily influenced by the presence of He, both in FeCr and W. Many important issues with fusion reactor materials remain poorly understood, and will require a huge effort by the international community. The development of potential models for new materials and the simulations performed in this thesis reveal many interesting features, but also serve as a platform for further studies.

Rotavirus nonstructural proteins: a structural perspective

Rotavirus is a major cause of acute infantile diarrhoea worldwide. The virus genome consists of 11 segments of double-stranded RNA that codesfor six structural proteins (VP1-6) and six non-structural proteins(NSP1-6). NSPs are proteins expressed from the virus genome in the infected cell, but are not incorporated into the mature virus article. NSPs play an essential role in virus replication, morphogenesis and pathogenesis, and most of them exhibit multifunctional properties. Structure-function analysis of the NSPs is essential for understanding the molecular mechanisms by which the virus circumvents host innate immune responses, inhibits cellular protein synthesis, hijacks the protein synthetic machinery for its own propagation and manifests the disease process. Because of their essential roles in virus biology, NSPs represent potential targets for the development of antiviral agents. Determination of the three-dimensional structure of NSPs has been hindered due to low-level expression and aggregation. To date, the complete three-dimensional structure of only NSP2 has been determined. The structures of the N- and C-terminal domains of NSP3 and the diarrhoea-inducing domain of NSP4 have also been determined. This review primarily covers the structural and biological functions of the NSPs whose three-dimensional structural aspects have been fully or partially understood, but provides a brief account of other NSPs and the structural features of the mature virion as determined by electron cryomicroscopy.

Characterization of hydrophobin proteins at interfaces and in solutions using x rays

Hydrophobins are a group of particularly surface active proteins. The surface activity is demonstrated in the ready adsorption of hydrophobins to hydrophobic/hydrophilic interfaces such as the air/water interface. Adsorbed hydrophobins self-assemble into ordered films, lower the surface tension of water, and stabilize air bubbles and foams. Hydrophobin proteins originate from filamentous fungi. In the fungi the adsorbed hydrophobin films enable the growth of fungal aerial structures, form protective coatings and mediate the attachment of fungi to solid surfaces. This thesis focuses on hydrophobins HFBI, HFBII, and HFBIII from a rot fungus Trichoderma reesei. The self-assembled hydrophobin films were studied both at the air/water interface and on a solid substrate. In particular, using grazing-incidence x-ray diffraction and reflectivity, it was possible to characterize the hydrophobin films directly at the air/water interface. The in situ experiments yielded information on the arrangement of the protein molecules in the films. All the T. reesei hydrophobins were shown to self-assemble into highly crystalline, hexagonally ordered rafts. The thicknesses of these two-dimensional protein crystals were below 30 Å. Similar films were also obtained on silicon substrates. The adsorption of the proteins is likely to be driven by the hydrophobic effect, but the self-assembly into ordered films involves also specific protein-protein interactions. The protein-protein interactions lead to differences in the arrangement of the molecules in the HFBI, HFBII, and HFBIII protein films, as seen in the grazing-incidence x-ray diffraction data. The protein-protein interactions were further probed in solution using small-angle x-ray scattering. Both HFBI and HFBII were shown to form mainly tetramers in aqueous solution. By modifying the solution conditions and thereby the interactions, it was shown that the association was due to the hydrophobic effect. The stable tetrameric assemblies could tolerate heating and changes in pH. The stability of the structure facilitates the persistence of these secreted proteins in the soil.

Efficient content-based image retrieval based on multi-feature fusion

The increased availability of image capturing devices has enabled collections of digital images to rapidly expand in both size and diversity. This has created a constantly growing need for efficient and effective image browsing, searching, and retrieval tools. Pseudo-relevance feedback (PRF) has proven to be an effective mechanism for improving retrieval accuracy. An original, simple yet effective rank-based PRF mechanism (RB-PRF) that takes into account the initial rank order of each image to improve retrieval accuracy is proposed. This RB-PRF mechanism innovates by making use of binary image signatures to improve retrieval precision by promoting images similar to highly ranked images and demoting images similar to lower ranked images. Empirical evaluations based on standard benchmarks, namely Wang, Oliva & Torralba, and Corel datasets demonstrate the effectiveness of the proposed RB-PRF mechanism in image retrieval.

Conformations of disulfide bridges in proteins

The conformational characteristics of disulfide bridges in proteins have been analyzed using a dataset of 22 protein structures, available at a resolution of 2.0 Å, containing a total of 72 disulfide crosslinks. The parameters used in the analysis include (φ, Ψ) values at Cys residues, bridge dihedral angles χss, χ1i, χ1j, χ2i and χ2j the distances Cαi-Cαj and Cβi-Cβj between the Cα and Cβ atoms of Cys(i) and Cys(j). Eight families of bridge conformations with three or more occurrences have been identified on the basis of these stereochemical parameters. The most populated family corresponds to the "left handed spiral" identified earlier by Richardson ((1981) Adv. Protein Chem. 34, 167–330). Disulfide bridging across antiparallel extended strands is observed in α-lytic protease, crambin, and β-trypsin and this structure is shown to be very similar to those obtained in small cystine peptides. Solvent accessible surface area calculations show that the overwhelming majority of disulfide bridges are inaccessible to solvent.

Selection of intrusion detection system threshold bounds for effective sensor fusion

The motivation behind the fusion of Intrusion Detection Systems was the realization that with the increasing traffic and increasing complexity of attacks, none of the present day stand-alone Intrusion Detection Systems can meet the high demand for a very high detection rate and an extremely low false positive rate. Multi-sensor fusion can be used to meet these requirements by a refinement of the combined response of different Intrusion Detection Systems. In this paper, we show the design technique of sensor fusion to best utilize the useful response from multiple sensors by an appropriate adjustment of the fusion threshold. The threshold is generally chosen according to the past experiences or by an expert system. In this paper, we show that the choice of the threshold bounds according to the Chebyshev inequality principle performs better. This approach also helps to solve the problem of scalability and has the advantage of failsafe capability. This paper theoretically models the fusion of Intrusion Detection Systems for the purpose of proving the improvement in performance, supplemented with the empirical evaluation. The combination of complementary sensors is shown to detect more attacks than the individual components. Since the individual sensors chosen detect sufficiently different attacks, their result can be merged for improved performance. The combination is done in different ways like (i) taking all the alarms from each system and avoiding duplications, (ii) taking alarms from each system by fixing threshold bounds, and (iii) rule-based fusion with a priori knowledge of the individual sensor performance. A number of evaluation metrics are used, and the results indicate that there is an overall enhancement in the performance of the combined detector using sensor fusion incorporating the threshold bounds and significantly better performance using simple rule-based fusion.

Peptide mimics for structural features in proteins. Crystal structures of three heptapeptide helices with a C-terminal 6-->1 hydrogen bond.

The crystal structure determination of three heptapeptides containing alpha-aminoisobutyryl (Aib) residues as a means of helix stabilization provides a high-resolution characterization of 6-->1 hydrogen-bonded conformations, reminiscent of helix-terminating structural features in proteins. The crystal parameters for the three peptides, Boc-Val-Aib-X-Aib-Ala-Aib-Y-OMe, where X and Y are Phe, Leu (I), Leu, Phe (II) and Leu, Leu (III) are: (I) space group P1, Z = 1, a = 9.903 A, b = 10.709 A, c = 11.969 A, alpha = 102.94 degrees, beta = 103.41 degrees, gamma = 92.72 degrees, R = 4.55%; (II) space group P21, Z = 2, a = 10.052 A, b = 17.653 A, c = 13.510 A, beta = 108.45 degrees, R = 4.49%; (III) space group P1, Z = 2 (two independent molecules IIIa and IIIb in the asymmetric unit), a = 10.833 A, b = 13.850 A, c = 16.928 A, alpha = 99.77 degrees, beta = 105.90 degrees, gamma = 90.64 degrees, R = 8.54%. In all cases the helices form 3(10)/alpha-helical (or 3(10)helical) structures, with helical columns formed by head-to-tail hydrogen bonding. The helices assemble in an all-parallel motif in crystals I and III and in an antiparallel motif in II. In the four crystallographically characterized molecules, I, II, IIIa and IIIb, Aib(6) adopts a left-handed helical (hL) conformation with positive phi, psi values, resulting in 6-->1 hydrogen-bond formation between Aib(2) CO and Leu(7)/Phe(7) NH groups. In addition a 4-->1 hydrogen bond is seen between Aib(3) CO and Aib(6) NH groups. This pattern of hydrogen bonding is often observed at the C-terminus of helices proteins, with the terminal pi-type turn being formed by four residues adopting the hRhRhRhL conformation.