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.

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.

Pollen and grain size records in abyssal sediments of the northwest Pacific Ocean as proxies of Plio-Pleistocene climate change /

A distinctive period of global change occurred during the PUocene between the warm Miocene and subsequent Quaternary cooling. Samples from Ocean Drilling Project Site 11 79 (-5586 mbsl, 41°4'N, 159°57'E), Site 881 (-5765 mbsl, 47°6.133'N, 161°29.490'E) and Site 882 (-3255 mbsl, 50°22'N, 167°36'E) were studied to determine the magnitude and composition ofterrigenous flux to the western mid-latitude North Pacific and its relation to climate change in East Asia since the mid-Pliocene. Dust-sized particles (including pollen), sourced from the arid regions and loess plateaus in East Asia are entrained by prevailing westerly winds and transported to the midlatitude northwest North Pacific Ocean. This is recorded by peaks in the total concentration of pollen and spores, as well as the mean grain size of allochthonous and autochthonous silicate material in abyssal marine sediments. Aridification of the Asian interior due to the phased uplift of the Himalayan-Tibetan Plateau created the modem East Asian Monsoon system dominated by a strengthening of the winter monsoon. The winter monsoon is further enhanced during glacials due to the expansion of desert and steppe environments at the expense ofwoodlands and forests recorded by the composition of palynological assemblages. The late Pliocene-Pleistocene glacials at ODP Sites 1 179, 881, and 882 are characterized by increases in grain size, magnetic susceptibility, pollen and spore concentrations around 3.5-3.3, 2.6-2.4, 1.7-1.6, and 0.9-0.7 Ma (ages based on magnetostratigraphic and biostratigraphic datums). The peaks during these times are relatively rich in pollen taxa derived primarily from steppe and boreal vegetation zones, recording cool, dry climates. The overall size increase of sediment and abundance of terrestrial palynomorphs record enhanced wind strength. The increase in magnitude of pollen and spore concentrations as well as grain size record global cooling and Northern Hemisphere glaciation. The peaks in grain size as well as pollen and spore abundance in marine sediments correlate with the mean grain size of loess in East Asia, consistent with the deflation of unarmoured surfaces during glacials. The transport of limiting nutrients to marine environments enhanced sea surface productivity and increased the rate of sediment accumulation.

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.