Protein folding: a perspective for biology, medicine and biotechnology

At the present time, protein folding is an extremely active field of research including aspects of biology, chemistry, biochemistry, computer science and physics. The fundamental principles have practical applications in the exploitation of the advances in genome research, in the understanding of different pathologies and in the design of novel proteins with special functions. Although the detailed mechanisms of folding are not completely known, significant advances have been made in the understanding of this complex process through both experimental and theoretical approaches. In this review, the evolution of concepts from Anfinsen's postulate to the "new view" emphasizing the concept of the energy landscape of folding is presented. The main rules of protein folding have been established from in vitro experiments. It has been long accepted that the in vitro refolding process is a good model for understanding the mechanisms by which a nascent polypeptide chain reaches its native conformation in the cellular environment. Indeed, many denatured proteins, even those whose disulfide bridges have been disrupted, are able to refold spontaneously. Although this assumption was challenged by the discovery of molecular chaperones, from the amount of both structural and functional information now available, it has been clearly established that the main rules of protein folding deduced from in vitro experiments are also valid in the cellular environment. This modern view of protein folding permits a better understanding of the aggregation processes that play a role in several pathologies, including those induced by prions and Alzheimer's disease. Drug design and de novo protein design with the aim of creating proteins with novel functions by application of protein folding rules are making significant progress and offer perspectives for practical applications in the development of pharmaceuticals and medical diagnostics.

Determination of serum aluminum, platelet aggregation and lipid peroxidation in hemodialyzed patients

Aluminum (Al3+) overload is frequently associated with lipid peroxidation and neurological disorders. Aluminum accumulation is also reported to be related to renal impairment, anemia and other clinical complications in hemodialysis patients. The aim of the present study was to determine the degree of lipid peroxidation, platelet aggregation and serum aluminum in patients receiving regular hemodialytic treatment. The level of plasma lipid peroxidation was evaluated on the basis of thiobarbituric acid reactive substances (TBARS). Mean platelet peroxidation in patients undergoing hemodialysis was significantly higher than in normal controls (2.7 ± 0.03 vs 1.8 ± 0.06 nmol/l, P<0.05). Platelet aggregation and serum aluminum levels were determined by a turbidimetric method and atomic absorption spectrophotometry, respectively. Serum aluminum was significantly higher in patients than in normal controls (44.5 ± 29 vs 10.8 ± 2.5 µg/l, P<0.05). Human blood platelets were stimulated with collagen (2.2 µg/ml), adenosine diphosphate (6 µM) and epinephrine (6 µM) and showed reduced function with the three agonists utilized. No correlation between aluminum levels and platelet aggregation or between aluminum and peroxidation was observed in hemodialyzed patients.

Role of Tamm-Horsfall protein and uromodulin in calcium oxalate crystallization

One of the defenses against nephrolithiasis is provided by macromolecules that modulate the nucleation, growth, aggregation and retention of crystals in the kidneys. The aim of the present study was to determine the behavior of two of these proteins, Tamm-Horsfall and uromodulin, in calcium oxalate crystallization in vitro. We studied a group of 10 male stone formers who had formed at least one kidney stone composed of calcium oxalate. They were classified as having idiopathic nephrolithiasis and had no well-known metabolic risk factors involved in kidney stone pathogenesis. Ten normal men were used as controls, as was a group consisting of five normal women and another consisting of five pregnant women. Crystallization was induced by a fixed supersaturation of calcium oxalate and measured with a Coulter Counter. All findings were confirmed by light and scanning electron microscopy. The number of particulate material deposited from patients with Tamm-Horsfall protein was higher than that of the controls (P<0.001). However, Tamm-Horsfall protein decreased the particle diameter of the stone formers when analyzed by the mode of the volume distribution curve (P<0.002) (5.64 ± 0.55 µm compared to 11.41 ± 0.48 µm of uromodulin; 15.94 ± 3.93 µm and 12.45 ± 0.97 µm of normal men Tamm-Horsfall protein and uromodulin, respectively; 8.17 ± 1.57 µm and 9.82 ± 0.95 µm of normal women Tamm-Horsfall protein and uromodulin, respectively; 12.17 ± 1.41 µm and 12.99 ± 0.51 µm of pregnant Tamm-Horsfall protein and uromodulin, respectively). Uromodulin produced fewer particles than Tamm-Horsfall protein in all groups. Nonetheless, the total volume of the crystals produced by uromodulin was higher than that produced by Tamm-Horsfall protein. Our results indicate a different effect of Tamm-Horsfall protein and uromodulin. This dual behavior suggests different functions. Tamm-Horsfall protein may act on nucleation and inhibit crystal aggregation, while uromodulin may promote aggregation of calcium oxalate crystals.

Large-eddy simulation of wind flows over complex terrains for wind energy applications

Wind energy has obtained outstanding expectations due to risks of global warming and nuclear energy production plant accidents. Nowadays, wind farms are often constructed in areas of complex terrain. A potential wind farm location must have the site thoroughly surveyed and the wind climatology analyzed before installing any hardware. Therefore, modeling of Atmospheric Boundary Layer (ABL) flows over complex terrains containing, e.g. hills, forest, and lakes is of great interest in wind energy applications, as it can help in locating and optimizing the wind farms. Numerical modeling of wind flows using Computational Fluid Dynamics (CFD) has become a popular technique during the last few decades. Due to the inherent flow variability and large-scale unsteadiness typical in ABL flows in general and especially over complex terrains, the flow can be difficult to be predicted accurately enough by using the Reynolds-Averaged Navier-Stokes equations (RANS). Large- Eddy Simulation (LES) resolves the largest and thus most important turbulent eddies and models only the small-scale motions which are more universal than the large eddies and thus easier to model. Therefore, LES is expected to be more suitable for this kind of simulations although it is computationally more expensive than the RANS approach. With the fast development of computers and open-source CFD software during the recent years, the application of LES toward atmospheric flow is becoming increasingly common nowadays. The aim of the work is to simulate atmospheric flows over realistic and complex terrains by means of LES. Evaluation of potential in-land wind park locations will be the main application for these simulations. Development of the LES methodology to simulate the atmospheric flows over realistic terrains is reported in the thesis. The work also aims at validating the LES methodology at a real scale. In the thesis, LES are carried out for flow problems ranging from basic channel flows to real atmospheric flows over one of the most recent real-life complex terrain problems, the Bolund hill. All the simulations reported in the thesis are carried out using a new OpenFOAM® -based LES solver. The solver uses the 4th order time-accurate Runge-Kutta scheme and a fractional step method. Moreover, development of the LES methodology includes special attention to two boundary conditions: the upstream (inflow) and wall boundary conditions. The upstream boundary condition is generated by using the so-called recycling technique, in which the instantaneous flow properties are sampled on aplane downstream of the inlet and mapped back to the inlet at each time step. This technique develops the upstream boundary-layer flow together with the inflow turbulence without using any precursor simulation and thus within a single computational domain. The roughness of the terrain surface is modeled by implementing a new wall function into OpenFOAM® during the thesis work. Both, the recycling method and the newly implemented wall function, are validated for the channel flows at relatively high Reynolds number before applying them to the atmospheric flow applications. After validating the LES model over simple flows, the simulations are carried out for atmospheric boundary-layer flows over two types of hills: first, two-dimensional wind-tunnel hill profiles and second, the Bolund hill located in Roskilde Fjord, Denmark. For the twodimensional wind-tunnel hills, the study focuses on the overall flow behavior as a function of the hill slope. Moreover, the simulations are repeated using another wall function suitable for smooth surfaces, which already existed in OpenFOAM® , in order to study the sensitivity of the flow to the surface roughness in ABL flows. The simulated results obtained using the two wall functions are compared against the wind-tunnel measurements. It is shown that LES using the implemented wall function produces overall satisfactory results on the turbulent flow over the two-dimensional hills. The prediction of the flow separation and reattachment-length for the steeper hill is closer to the measurements than the other numerical studies reported in the past for the same hill geometry. The field measurement campaign performed over the Bolund hill provides the most recent field-experiment dataset for the mean flow and the turbulence properties. A number of research groups have simulated the wind flows over the Bolund hill. Due to the challenging features of the hill such as the almost vertical hill slope, it is considered as an ideal experimental test case for validating micro-scale CFD models for wind energy applications. In this work, the simulated results obtained for two wind directions are compared against the field measurements. It is shown that the present LES can reproduce the complex turbulent wind flow structures over a complicated terrain such as the Bolund hill. Especially, the present LES results show the best prediction of the turbulent kinetic energy with an average error of 24.1%, which is a 43% smaller than any other model results reported in the past for the Bolund case. Finally, the validated LES methodology is demonstrated to simulate the wind flow over the existing Muukko wind farm located in South-Eastern Finland. The simulation is carried out only for one wind direction and the results on the instantaneous and time-averaged wind speeds are briefly reported. The demonstration case is followed by discussions on the practical aspects of LES for the wind resource assessment over a realistic inland wind farm.

High doses of riboflavin and the elimination of dietary red meat promote the recovery of some motor functions in Parkinson's disease patients

Abnormal riboflavin status in the absence of a dietary deficiency was detected in 31 consecutive outpatients with Parkinson's disease (PD), while the classical determinants of homocysteine levels (B6, folic acid, and B12) were usually within normal limits. In contrast, only 3 of 10 consecutive outpatients with dementia without previous stroke had abnormal riboflavin status. The data for 12 patients who did not complete 6 months of therapy or did not comply with the proposed treatment paradigm were excluded from analysis. Nineteen PD patients (8 males and 11 females, mean age ± SD = 66.2 ± 8.6 years; 3, 3, 2, 5, and 6 patients in Hoehn and Yahr stages I to V) received riboflavin orally (30 mg every 8 h) plus their usual symptomatic medications and all red meat was eliminated from their diet. After 1 month the riboflavin status of the patients was normalized from 106.4 ± 34.9 to 179.2 ± 23 ng/ml (N = 9). Motor capacity was measured by a modification of the scoring system of Hoehn and Yahr, which reports motor capacity as percent. All 19 patients who completed 6 months of treatment showed improved motor capacity during the first three months and most reached a plateau while 5/19 continued to improve in the 3- to 6-month interval. Their average motor capacity increased from 44 to 71% after 6 months, increasing significantly every month compared with their own pretreatment status (P < 0.001, Wilcoxon signed rank test). Discontinuation of riboflavin for several days did not impair motor capacity and yellowish urine was the only side effect observed. The data show that the proposed treatment improves the clinical condition of PD patients. Riboflavin-sensitive mechanisms involved in PD may include glutathione depletion, cumulative mitochondrial DNA mutations, disturbed mitochondrial protein complexes, and abnormal iron metabolism. More studies are required to identify the mechanisms involved.

Rigorous development of safety-critical systems

Nowadays, computer-based systems tend to become more complex and control increasingly critical functions affecting different areas of human activities. Failures of such systems might result in loss of human lives as well as significant damage to the environment. Therefore, their safety needs to be ensured. However, the development of safety-critical systems is not a trivial exercise. Hence, to preclude design faults and guarantee the desired behaviour, different industrial standards prescribe the use of rigorous techniques for development and verification of such systems. The more critical the system is, the more rigorous approach should be undertaken. To ensure safety of a critical computer-based system, satisfaction of the safety requirements imposed on this system should be demonstrated. This task involves a number of activities. In particular, a set of the safety requirements is usually derived by conducting various safety analysis techniques. Strong assurance that the system satisfies the safety requirements can be provided by formal methods, i.e., mathematically-based techniques. At the same time, the evidence that the system under consideration meets the imposed safety requirements might be demonstrated by constructing safety cases. However, the overall safety assurance process of critical computerbased systems remains insufficiently defined due to the following reasons. Firstly, there are semantic differences between safety requirements and formal models. Informally represented safety requirements should be translated into the underlying formal language to enable further veri cation. Secondly, the development of formal models of complex systems can be labour-intensive and time consuming. Thirdly, there are only a few well-defined methods for integration of formal verification results into safety cases. This thesis proposes an integrated approach to the rigorous development and verification of safety-critical systems that (1) facilitates elicitation of safety requirements and their incorporation into formal models, (2) simplifies formal modelling and verification by proposing specification and refinement patterns, and (3) assists in the construction of safety cases from the artefacts generated by formal reasoning. Our chosen formal framework is Event-B. It allows us to tackle the complexity of safety-critical systems as well as to structure safety requirements by applying abstraction and stepwise refinement. The Rodin platform, a tool supporting Event-B, assists in automatic model transformations and proof-based verification of the desired system properties. The proposed approach has been validated by several case studies from different application domains.

Optimization and Measuring Techniques for Collect-and-Place Machines in Printed Circuit Board Industry

This thesis considers optimization problems arising in printed circuit board assembly. Especially, the case in which the electronic components of a single circuit board are placed using a single placement machine is studied. Although there is a large number of different placement machines, the use of collect-and-place -type gantry machines is discussed because of their flexibility and increasing popularity in the industry. Instead of solving the entire control optimization problem of a collect-andplace machine with a single application, the problem is divided into multiple subproblems because of its hard combinatorial nature. This dividing technique is called hierarchical decomposition. All the subproblems of the one PCB - one machine -context are described, classified and reviewed. The derived subproblems are then either solved with exact methods or new heuristic algorithms are developed and applied. The exact methods include, for example, a greedy algorithm and a solution based on dynamic programming. Some of the proposed heuristics contain constructive parts while others utilize local search or are based on frequency calculations. For the heuristics, it is made sure with comprehensive experimental tests that they are applicable and feasible. A number of quality functions will be proposed for evaluation and applied to the subproblems. In the experimental tests, artificially generated data from Markov-models and data from real-world PCB production are used. The thesis consists of an introduction and of five publications where the developed and used solution methods are described in their full detail. For all the problems stated in this thesis, the methods proposed are efficient enough to be used in the PCB assembly production in practice and are readily applicable in the PCB manufacturing industry.

Hygro-elasto-plastic behavior of planar orthotropic material

The mechanical and hygroscopic properties of paper and board are factors affecting the whole lifecycle of a product, including paper/board quality, production, converting, and material and energy savings. The progress of shrinkage profiles, loose edges of web, baggy web causing wrinkling and misregistration in printing are examples of factors affecting runnability and end product quality in the drying section and converting processes, where paper or board is treated as a moving web. The structural properties and internal stresses or plastic strain differences built up during production also cause the end-product defects related to distortion of the shape of the product such as sheet or box. The objective of this work was to construct a model capable of capturing the characteristic behavior of hygroscopic orthotropic material under moisture change, during different external in-plane stretch or stress conditions. Two independent experimental models were constructed: the elasto-plastic material model and the hygroexpansivity-shrinkage model. Both describe the structural properties of the sheet with a fiber orientation probability distribution, and both are functions of the dry solids content and fiber orientation anisotropy index. The anisotropy index, introduced in this work, simplifies the procedure of determining the constitutive parameters of the material model and the hygroexpansion coefficients in different in-plane directions of the orthotropic sheet. The mathematically consistent elasto-plastic material model and the dry solids content dependent hygroexpansivity have been constructed over the entire range from wet to dry. The presented elastoplastic and hygroexpansivity-shrinkage models can be used in an analytical approach to estimate the plastic strain and shrinkage in simple one-dimensional cases. For studies of the combined and more complicated effects of hygro-elasto-plastic behavior, both models were implemented in a finite element program for a numerical solution. The finite element approach also offered possibilities for studying different structural variations of orthotropic planar material, as well as local buckling behavior and internal stress situations of the sheet or web generated by local strain differences. A comparison of the simulation examples presented in this work to results published earlier confirms that the hygro-elasto-plastic model provides at least qualitatively reasonable estimates. The application potential of the hygro-elasto-plastic model is versatile, including several phenomena and defects appearing in the drying, converting and end-use conditions of the paper or board webs and products, or in other corresponding complex planar materials.

The eukaryotic Pso2/Snm1/Artemis proteins and their function as genomic and cellular caretakers

DNA double-strand breaks (DSBs) represent a major threat to the genomic stability of eukaryotic cells. DNA repair mechanisms such as non-homologous end joining (NHEJ) are responsible for the maintenance of eukaryotic genomes. Dysfunction of one or more of the many protein complexes that function in NHEJ can lead to sensitivity to DNA damaging agents, apoptosis, genomic instability, and severe combined immunodeficiency. One protein, Pso2p, was shown to participate in the repair of DSBs induced by DNA inter-strand cross-linking (ICL) agents such as cisplatin, nitrogen mustard or photo-activated bi-functional psoralens. The molecular function of Pso2p in DNA repair is unknown, but yeast and mammalian cell line mutants for PSO2 show the same cellular responses as strains with defects in NHEJ, e.g., sensitivity to ICLs and apoptosis. The Pso2p human homologue Artemis participates in V(D)J recombination. Mutations in Artemis induce a variety of immunological deficiencies, a predisposition to lymphomas, and an increase in chromosomal aberrations. In order to better understand the role of Pso2p in the repair of DSBs generated as repair intermediates of ICLs, an in silico approach was used to characterize the catalytic domain of Pso2p, which led to identification of novel Pso2p homologues in other organisms. Moreover, we found the catalytic core of Pso2p fused to different domains. In plants, a specific ATP-dependent DNA ligase I contains the catalytic core of Pso2p, constituting a new DNA ligase family, which was named LIG6. The possible functions of Pso2p/Artemis/Lig6p in NHEJ and V(D)J recombination and in other cellular metabolic reactions are discussed.

Quantitative evidence for neurofilament heavy subunit aggregation in motor neurons of spinal cords of patients with amyotrophic lateral sclerosis

Amyotrophic lateral sclerosis (ALS), a neurodegenerative disease of unknown etiology, affects motor neurons leading to atrophy of skeletal muscles, paralysis and death. There is evidence for the accumulation of neurofilaments (NF) in motor neurons of the spinal cord in ALS cases. NF are major structural elements of the neuronal cytoskeleton. They play an important role in cell architecture and differentiation and in the determination and maintenance of fiber caliber. They are composed of three different polypeptides: light (NF-L), medium (NF-M) and heavy (NF-H) subunits. In the present study, we performed a morphological and quantitative immunohistochemical analysis to evaluate the accumulation of NF and the presence of each subunit in control and ALS cases. Spinal cords from patients without neurological disease and from ALS patients were obtained at autopsy. In all ALS cases there was a marked loss of motor neurons, besides atrophic neurons and preserved neurons with cytoplasmic inclusions, and extensive gliosis. In control cases, the immunoreaction in the cytoplasm of neurons was weak for phosphorylated NF-H, strong for NF-M and weak for NF-L. In ALS cases, anterior horn neurons showed intense immunoreactivity in focal regions of neuronal perikarya for all subunits, although the difference in the integrated optical density was statistically significant only for NF-H. Furthermore, we also observed dilated axons (spheroids), which were immunopositive for NF-H but negative for NF-M and NF-L. In conclusion, we present qualitative and quantitative evidence of NF-H subunit accumulation in neuronal perikarya and spheroids, which suggests a possible role of this subunit in the pathogenesis of ALS.

Differences in pressure and temperature transitions of proteins and polymer gels

Pressure-driven and temperature-driven transitions of two thermoresponsive polymers, poly(N-isopropylacrylamide) (pNIPAM) and poly(N-vinylisobutyramide) (pNVIBA)), in both a soluble linear polymer form and a cross-linked hydro-gel form, were examined by a dynamic light-scattering method and direct microscopic observation, respectively. Their behavior was compared with that of protein systems. Changes in some characteristic parameters in the time-intensity correlation functions of dynamic light-scattering measurement of aqueous solutions of pNIPAM at various pressures and temperatures showed no essential differences during temperature and pressure scanning and, as a whole, the motions of polymers in aqueous solutions were similar in two types of transitions until chain shrinkage occurred. The gels (cross-linked polymer gels) prepared from the thermoresponsive polymers also showed similar volume transitions responding to the pressure and temperature increase. In temperature transitions, however, gels showed drastic volume shrinkage with loss of transparency, while pressure-induced transition showed a slow recovery of transparency while keeping the size, after first transient drastic volume shrinkage with loss of transparency. At a temperature slightly higher than the transition under atmospheric temperature, so-called reentry of the volume change and recovery of the transparency were observed during the pressure-increasing process, which implies much smaller aggregation or non-aggregated collapsed polymer chains in the gel at higher pressures, indicating a certain mechanistic difference of the dehydration processes induced by temperature and pressure.

Disintegrins: integrin selective ligands which activate integrin-coupled signaling and modulate leukocyte functions

Extracellular matrix proteins and cell adhesion receptors (integrins) play essential roles in the regulation of cell adhesion and migration. Interactions of integrins with the extracellular matrix proteins lead to phosphorylation of several intracellular proteins such as focal adhesion kinase, activating different signaling pathways responsible for the regulation of a variety of cell functions, including cytoskeleton mobilization. Once leukocytes are guided to sites of infection, inflammation, or antigen presentation, integrins can participate in the initiation, maintenance, or termination of the immune and inflammatory responses. The modulation of neutrophil activation through integrin-mediated pathways is important in the homeostatic control of the resolution of inflammatory states. In addition, during recirculation, T lymphocyte movement through distinct microenvironments is mediated by integrins, which are critical for cell cycle, differentiation and gene expression. Disintegrins are a family of low-molecular weight, cysteine-rich peptides first identified in snake venom, usually containing an RGD (Arg-Gly-Asp) motif, which confers the ability to selectively bind to integrins, inhibiting integrin-related functions in different cell systems. In this review we show that, depending on the cell type and the microenvironment, disintegrins are able to antagonize the effects of integrins or to act agonistically by activating integrin-mediated signaling. Disintegrins have proven useful as tools to improve the understanding of the molecular events regulated by integrin signaling in leukocytes and prototypes in order to design therapies able to interfere with integrin-mediated effects.

Differential evolution based classification with pool of distances and aggregation operators

The objective of this thesis is to develop and generalize further the differential evolution based data classification method. For many years, evolutionary algorithms have been successfully applied to many classification tasks. Evolution algorithms are population based, stochastic search algorithms that mimic natural selection and genetics. Differential evolution is an evolutionary algorithm that has gained popularity because of its simplicity and good observed performance. In this thesis a differential evolution classifier with pool of distances is proposed, demonstrated and initially evaluated. The differential evolution classifier is a nearest prototype vector based classifier that applies a global optimization algorithm, differential evolution, to determine the optimal values for all free parameters of the classifier model during the training phase of the classifier. The differential evolution classifier applies the individually optimized distance measure for each new data set to be classified is generalized to cover a pool of distances. Instead of optimizing a single distance measure for the given data set, the selection of the optimal distance measure from a predefined pool of alternative measures is attempted systematically and automatically. Furthermore, instead of only selecting the optimal distance measure from a set of alternatives, an attempt is made to optimize the values of the possible control parameters related with the selected distance measure. Specifically, a pool of alternative distance measures is first created and then the differential evolution algorithm is applied to select the optimal distance measure that yields the highest classification accuracy with the current data. After determining the optimal distance measures for the given data set together with their optimal parameters, all determined distance measures are aggregated to form a single total distance measure. The total distance measure is applied to the final classification decisions. The actual classification process is still based on the nearest prototype vector principle; a sample belongs to the class represented by the nearest prototype vector when measured with the optimized total distance measure. During the training process the differential evolution algorithm determines the optimal class vectors, selects optimal distance metrics, and determines the optimal values for the free parameters of each selected distance measure. The results obtained with the above method confirm that the choice of distance measure is one of the most crucial factors for obtaining higher classification accuracy. The results also demonstrate that it is possible to build a classifier that is able to select the optimal distance measure for the given data set automatically and systematically. After finding optimal distance measures together with optimal parameters from the particular distance measure results are then aggregated to form a total distance, which will be used to form the deviation between the class vectors and samples and thus classify the samples. This thesis also discusses two types of aggregation operators, namely, ordered weighted averaging (OWA) based multi-distances and generalized ordered weighted averaging (GOWA). These aggregation operators were applied in this work to the aggregation of the normalized distance values. The results demonstrate that a proper combination of aggregation operator and weight generation scheme play an important role in obtaining good classification accuracy. The main outcomes of the work are the six new generalized versions of previous method called differential evolution classifier. All these DE classifier demonstrated good results in the classification tasks.