Assessing the Reproducibility of Clustering of Molecular Dynamics Conformations on Self-Organizing Maps

The goal of most clustering algorithms is to find the optimal number of clusters (i.e. fewest number of clusters). However, analysis of molecular conformations of biological macromolecules obtained from computer simulations may benefit from a larger array of clusters. The Self-Organizing Map (SOM) clustering method has the advantage of generating large numbers of clusters, but often gives ambiguous results. In this work, SOMs have been shown to be reproducible when the same conformational dataset is independently clustered multiple times (~100), with the help of the Cramérs V-index (C_v). The ability of C_v to determine which SOMs are reproduced is generalizable across different SOM source codes. The conformational ensembles produced from MD (molecular dynamics) and REMD (replica exchange molecular dynamics) simulations of the penta peptide Met-enkephalin (MET) and the 34 amino acid protein human Parathyroid Hormone (hPTH) were used to evaluate SOM reproducibility. The training length for the SOM has a huge impact on the reproducibility. Analysis of MET conformational data definitively determined that toroidal SOMs cluster data better than bordered maps due to the fact that toroidal maps do not have an edge effect. For the source code from MATLAB, it was determined that the learning rate function should be LINEAR with an initial learning rate factor of 0.05 and the SOM should be trained by a sequential algorithm. The trained SOMs can be used as a supervised classification for another dataset. The toroidal 10×10 hexagonal SOMs produced from the MATLAB program for hPTH conformational data produced three sets of reproducible clusters (27%, 15%, and 13% of 100 independent runs) which find similar partitionings to those of smaller 6×6 SOMs. The χ^2 values produced as part of the C_v calculation were used to locate clusters with identical conformational memberships on independently trained SOMs, even those with different dimensions. The χ^2 values could relate the different SOM partitionings to each other.

Relationship between physical activity and resting secretory immunity in children

This study examined relationships among physical activity, body fat and salivary immonoglobulin A (sIgA) levels in adolescent children of Southern Ontario. Gender differences on these factors were also assessed. Sixty-one grade-five students (10-1 lyrs), males (n=29) and females (n=31), who had not received a flu vaccination in the past 12 months, participated in the study. They were assessed for: aerobic power (20-m shuttle run), relative body fat (bioelectrical impedance analysis), sIgA, sIgA/albumin ratio, and salivary Cortisol. Each subject completed the Habitual Activity Estimation Scale and the Participation Questioimaire. Students wore a pedometer for 48h to estimate their average total distance traveled per day. The results show 40% of the children were over 25% body fat and 50% of them spend less than five hours per day in any physical activities. Salivary IgA was not related to salivary Cortisol, physical activity, fitness level or body fat in this age group. There were no gender differences in sIgA and Cortisol levels. Boys had a significantly higher aerobic power and daily distance traveled, but reported similar organized and fi-ee time activity participation levels as the girls. The test-retest reproducibility for salivary Cortisol was 0.663 (p<0.01), while long term sIgA and sIgA/albumin ratio reproducibility was non-significant for repeated measurements taken after six weeks. It was found that salivary IgA has not been shovm to be a stable measure in children, in contrast to the results found in the literatiu-e that tested adults and the relationship with physical activity, fitness level and body fat.

Investigation into the determination of arsenic by direct current plasma atomic emission spectrometer

Improvements have been made on the currently available hydride generator system manufactured by SpectraMetrics Incorporated, because the system was found to be unsatisfactory with respect to the following: 1. the drying agent, anhydrous calcium chloride, 2. the special sample tube, 3. the direction of argon flow through the Buchner funnel when it came to dealing with real sample, that is, with reference only to aqueous extracts of soil samples. Changes that were made on the system included the replacement of anhydrous calcium chloride with anhydrous calcium sulphate and the replacement of the special sample tube with a modified one made from silica. Re-directing the flow of argon through the top of the Buchner funnel appeared to make the system compatible with aqueous extracts of soil samples. The interferences from 1000 ~g/mL of nickel(II) , cobalt(II), iron(III), copper(II) have been eliminated with the aid of 1.4 M hydrochloric acid and 1% (weight/volume) L-cystine. Greater than 90% recovery of 0.3 ~g/mL arsenic signal was achieved in each case. Furthermore, 103% of arsenic signal was accomplished in the presence of 1000 ~g/mL cadmium with 5 M Hel. tVhen each of the interferents was present in solution at 1000 ppm, a recovery of 85% was achieved by using 5 M hydrochloric acid and 3% (weight/volume) L-cystine. Without L-cystine and when 1.4 M hydrochloric acid was used, the recoveries were 0% (Ni), 0% (Co), 88% (Fe), 15% (Cu), 18% (Cd). Similarly, a solution containing 1000 ppm of each interferent gave a zero percent recovery of arsenic. The reduction of trivalent and pentavalent arsenic at a pH less than one has also been investigated and shown to be quantitative if peak areas are measured. The reproducibility determination of a 0.3 Vg/mL standard arsenic solution by hydride generation shows a relative standard deviation of 3.4%. The detection limits with and without Porapak Q have been found to be 0.6 ng/mL and 1.0 ng/mL, respectively.

Investigations into the determinations of hydride-forming elements

Analytical methods for the determination of trace amounts of germanium, tin and arsenic were established using hydride generation coupled with direct current plasma atomic emission spectrometry. A continuous gas flowing batch system for the hydride generation was investigated and was applied to the determination of germanium(Ge), tin(Sn), antimony(Sb) and lead(Pb) (Preliminary results suggest that it is also applicable to arsenic)As) ). With this system, the reproducibility of signals was improved and the determination was speeded up, compared with the conventional batch type hydride generation system. Each determination was complete within one minute. Interferences from a number of transition metal ions, especially from Pd(II), Pt(IV), Ni(II), Cu(II), Co(II), and Fe(II, III), have proven to be very serious under normal conditions, in the determination of germanium, tin, and arsenic. These interference effects were eliminated or significantly reduced in the presence of L-cystine or L-cysteine. Thus, a 10-1000 fold excess of Ni(II), Cu(II), Co(II), Fe(II), Pt(IV), Pd(II), etc. can be tolerated without interference, In the presence of L-cystine or L-cysteine, compared with absence of interference reducing agent. The methods for the determination of Ge, Sn, and As were examined by the analyses of standard reference materials. Interference effects from the sample matrix, for example, in transition metal-rich samples, copper, iron and steel, were eliminated by L-cystine (for As and Sn) and by LI cysteine (for Ge). The analysis of a number of standard reference materials gave excellent results of As and Sn contents in agreement with the certified values, showing there was no systematic interference. The detection limits for both germanium and tin were 20 pg ml- I . Preliminary studies were carried out for the determination of antimony and lead. Antimony was found to react with NaBH4, remaInIng from the previous determinations, giving an analytical signal. A reversed injection manner, i.e., injection of the NaBH4 solution prior to the analyte solution was used to avoid uncertainty caused by residual NaBH4 present and to ensure that an excess of NaB H4 was available. A solution of 0.4% L-cysteine was found to reduce the interference from selected transition metal ions, Co(II), Cu(II), Ni(II) and Pt(IV). Hydrochloric acid - hydrogen peroxide, nitric acid - ammonium persulphate, and potassium dichromate malic acid reaction systems for lead hydride generation were compared. The potassium dichromate - malic acid reaction medium proved to be the best with respect to reproducibility and minimal interference. Cu(II), Ni(II), and Fe(II) caused strong interference In lead determinations, which was not reduced by L-cysteine or Lcystine. Sodium citrate, ascorbic acid, dithizone, thiosemicarbazide and penicillamine reduced interferences to some extent. Further interference reduction studies were carried out uSIng a number of amino acids, glycine, alanine, valine, leucine and histidine, as possible interference reducing agents in the determination of germanium. From glycine, alanine, valine to leucine, the interference reduction effect in germanium determinations decreased. Histidine II was found to be very promising In the reduction of interference. In fact, histidine proved more efficient than L-cystine and L-cysteine In the reduction of interference from Ni(II) in the determination of germanium. Signal enhancement by easily ionized elements (EIEs), usually regarded as an interference effect in analysis by DCP-AES, was studied and successfully applied to advantage in improving the sensitivity and detection limit in the determination of As, Ge, Sn, Sb, and Pb. The effect of alkali and alkaline-earth elements on the determination of the above five hydride forming elements was studied. With the appropriate EIE, a signal enhancement of 40-115% was achieved. Linear calibration and good reproducibility were also obtained in the presence of EIEs. III

The determination of sinigrin in Brassica, and investigations into the use of allyl-isothiocyanate as a nematicide

The goal of this thesis was to study factors related to the development of Brassica juncea as a sustainable nematicide. Brassica juncea is characterized by the glycoside (glucosinolate) sinigrin. Various methods were developed for the determination of sinigrin in Brassica juncea tissue extracts. Sinigrin concentrations in plant tissues at various stages of growth were monitored. Sinigrin enzymatically breaks down into allylisothiocyanate (AITC). AITC is unstable in aqueous solution and degradation was studied in water and in soil. Finally, the toxicity of AITC against the root-lesion nematode (Pratylenchus penetrans) was determined. A method was developed to extract sinigrin from whole Brassica j uncea tissues. The optimal time of extraction wi th boiling phosphate buffer (0.7mM, pH=6.38) and methanol/water (70:30 v/v) solutions were both 25 minutes. Methanol/water extracted 13% greater amount of sinigrin than phosphate buffer solution. Degradation of sinigrin in boiling phosphate buffer solution (0.13%/minute) was similar to the loss of sinigrin during the extraction procedure. The loss of sinigrin from boiling methanol/water was estimated to be O.Ol%/minute. Brassica juncea extract clean up was accomplished by an ion-pair solid phase extraction (SPE) method. The recovery of sinigrin was 92.6% and coextractive impurities were not detected in the cleaned up extract. Several high performance liquid chromatography (HPLC) methods were developed for the determination of sinigrin. All the developed methods employed an isocratic mobile phase system wi th a low concentration of phosphate buffer solution, ammonium acetate solution or an ion-pair reagent solution. A step gradient system was also developed. The method involved preconditioning the analytical column with phosphate buffer solution and then switching the mobile phase to 100% water after sample injection.Sinigrin and benzyl-glucosinolate were both studied by HPLC particle beam negative chemical ionization mass spectrometry (HPLCPB- NCI-MS). Comparison of the mass spectra revealed the presence of fragments arising from the ~hioglucose moiety and glucosinolate side-chain. Variation in the slnlgrin concentration within Brassica juncea plants was studied (Domo and Cutlass cuItivars). The sinigrin concentration in the top three leaves was studied during growth of each cultivar. For Cutlass, the minimum (200~100~g/g) and maximum (1300~200~g/g) concentrations were observed at the third and seventh week after planting, respectively. For Domo, the minimum (190~70~g/g) and maximum (1100~400~g/g) concentrations were observed at the fourth and eighth week after planting, respectively. The highest sinigrin concentration was observed in flower tissues 2050±90~g/g and 2300±100~g/g for Cutlass and Domo cultivars, respectively. Physical properties of AITC were studied. The solubility of AITC in water was determined to be approximately 1290~g/ml at 24°C. An HPLC method was developed for the separation of degradation compounds from aqueous AITC sample solutions. Some of the degradation compounds identified have not been reported in the literature: allyl-thiourea, allyl-thiocyanate and diallyl-sulfide. In water, AITC degradation to' diallyl-thiourea was favored at basic pH (9.07) and degradation to diallyl-sulfide was favored at acidic pH (4 . 97). It wap necessary to amend the aqueous AITC sample solution with acetonitrile ?efore injection into the HPLC system. The acetonitrile amendment considerably improved AITC recovery and the reproducibility of the results. The half-life of aqueous AITC degradation at room temperature did not follow first-order kinetics. Beginning with a 1084~g/ml solution, the half-life was 633 hours. Wi th an ini tial AITC concentration of 335~g/ml the half-life was 865 hours. At 35°C the half-life AITC was 76+4 hours essentially independent of the iiisolution pH over the range of pH=4.97 to 9.07 (1000~g/ml). AITC degradation was also studied in soil at 35°C; after 24 hours approximately 75% of the initial AITC addition was unrecoverable by water extraction. The ECso of aqueous AITC against the root-lesion nematode (Pratylenchus penetrans) was determined to be approximately 20~g/ml at one hour exposure of the nematode to the test solution. The toxicological study was also performed with a myrosinase treated Brassica juncea extract. Myrosinase treatment of the Brassica juncea extract gave nearly quantitative conversion of sinigrin into AITC. The myrosinase treated extract was of the same efficacy as an aqueous AITC solution of equivalent concentration. The work of this thesis was focused upon understanding parameters relevant to the development of Brassica juncea as a sustainable nematicide. The broad range of experiments were undertaken in support of a research priority at Agriculture and Agri-Food Canada.

The measured variation of the Debye-Waller factor of aluminum from 295K to 815K by using the energy dispersive x-ray diffraction technique

The Energy Dispersive X-ray Diffraction System at Brock University has been used to measure the intensities of the diffraction lines of aluminum powder sample as a function of temperature. At first, intensity measurements at high temperature were not reproducible. After some modifications have been made, we were able to measure the intensities of the diffraction lines to 815K, with good accuracy and reproducibility. Therefore the changes of the Debye-Waller factor from room temperature up to 815K for aluminum were determined with precision. Our results are in good agreement with those previously published.

Transmission measurements using a new far-infrared set-up

A system comprised of a Martin-Puplett type polarizing interferometer and a Helium-3 cryostat was developed to study the transmission of materials in the very-far-infrared region of the spectrum. This region is of significant interest due to the low-energy excitations which many materials exhibit. The experimental transmission spectrum contains information concerning the optical properties of the material. The set-up of this system is described in detail along with the adaptations and improvements which have been made to the system to ensure the best results. Transmission experiments carried out with this new set-up for two different varieties of materials: superconducting thin films of lead and biological proteins, are discussed. Several thin films of lead deposited on fused silica quartz substrates were studied. From the ratio of the transmission in the superconducting state to that in the normal state the superconducting energy gap was determined to be approximately 25 cm-1 which corresponds to 2~/kBTc rv 5 in agreement with literature data. Furthermore, in agreement with theoretical predictions, the maximum in the transmission ratio was observed to increase as the film thickness was increased. These results provide verification of the system's ability to accurately measure the optical properties of thin low-Tc superconducting films. Transmission measurements were carried out on double deionized water, and a variety of different concentrations by weight of the globular protein, Bovine Serum Albumin, in the sol, gel and crystalline forms. The results of the water study agree well with literature values and thus further illustrate the reproducibility of the system. The results of the protein experiments, although preliminary, indicate that as the concentration increases the samples become more transparent. Some weak structure in the frequency dependent absorption coefficient, which is more prominent in crystalline samples, may be due to low frequency vibrations of the protein molecules.

Preliminary characterization of VmTPT2, a putative monoterpene indole alkaloid (MIA) transporter from Vinca minor L. (Apocynaceae)

The various steps of monoterpene indole alkaloid (MIA) biosynthesis are known to occur in specialized cell types and subcellular compartments. Numerous MIAs display powerful biological activities that have led to their use as pharmaceutical treatments for cancer, hypertension and malaria. Many of these compounds accumulate on the leaf surface of medicinally important Apocynaceae plants, which led to the recent discovery and characterization of an ABC transporter (CrTPT2) that was shown to mobilize catharanthine from its site of biosynthesis in epidermal cells to the leaf surface of Catharanthus roseus. Bioinformatic analysis of transcriptomes from several geographically distant MIA-producing species led to the identification of proteins with high amino acid sequence identity to CrTPT2. Molecular cloning of a similar transporter (VmTPT2) from Vinca minor was carried out and expressed in a yeast heterologous system for transport experiments and functional characterization. In planta studies involved transcript expression analysis of the early MIA biosynthetic gene VmTDC and putative transporter VmTPT2, and alkaloid profile analyses. RT-qPCR results showed that VmTPT2 expression increased 15-fold between the first two leaf pairs, and high levels were maintained across older leaves. The alkaloid accumulation profile on leaf surfaces matched that of VmTPT2 expression, especially for the MIAs vincadifformine and vincamine. Gene expression and alkaloid profile analyses suggest that the functional protein may act as a similar transporter to CrTPT2. However, although VmTPT2 had 88.4% identity at the amino acid level to CrTPT2, it displayed an altered expression pattern in planta across developing leaves, and functional characterization using a previously developed yeast heterologous system was unsuccessful due to difficulties with reproducibility of transport assays.