Novel Insights into the Molecular Genetic Background of Colorectal Tumors

Many of the genes predisposing to highly penetrant colorectal cancer (CRC) syndromes, including hereditary non-polyposis colorectal cancer (MLH1, MSH2, MSH6, PMS2), familial adenomatous polyposis (APC), Peutz-Jeghers syndrome (LKB1), juvenile polyposis (SMAD4, BMPR1A), MYH-associated polyposis (MYH), and Cowden syndrome (PTEN) have already been discovered. Identification of these genes has allowed a more precise classification of the hereditary CRC syndromes and provided a means for predictive genetic testing and surveillance. Some of the genes are also involved in sporadic cancer forms, and therefore the investigation of the rare CRC syndromes has been a breakthrough for general cancer research. Despite the accumulating knowledge on hereditary cancer syndromes, a significant number of familial CRCs remain molecularly unexplained after genetic testing, reflecting the possibility of other predisposing genes or existence of novel syndromes. Moreover, genetic variants conferring low-penetrance risk are still largely unknown. In this study, we examined the role of some new high- and low-penetrance alleles on CRC predisposition. We identified disease causing MYH mutations in a subset (9%) of patients with APC and AXIN2 mutation negative adenomatous polyposis. Due to differences in the pattern of inheritance and clinical manifestation, screening for mutations in MYH is beneficial in view of genetic counselling and surveillance. A novel functionally deficient MYH founder mutation A459D was identified in the Finnish population, and this finding had immediate clinical implications for genetic counselling of at risk families. Many patients with hamartomatous polyposis remain without molecular diagnosis due to atypical phenotypes. We therefore sought to classify 49 patients with unexplained hamartomatous or hyperplastic/mixed polyposis by extensive molecular analyses of PTEN, LKB1, BMPR1A, SMAD4, ENG, BRAF, MYH, and BHD along with revision of polyp histology. Mutations were identified in 11/49 (22%) of the patients. In 6 cases the molecular diagnosis was re-classified guiding surveillance and decisions for prophylactic surgery. Re-evaluation of polyp histology with subsequent more accurate selection of candidate gene analyses is beneficial and can be recommended for patients with unexplained polyposis. Furthermore, germline mutations in ENG underlying juvenile polyposis were described for the first time, characterizing a possible novel genetically defined form of hereditary CRC. Association analyses on two putative low-penetrance alleles, NOD2 3020insC and MDM2 SNP309 were performed in a population-based series of 1042 Finnish CRC patients and in cancer-free controls. In contrast to previous results, NOD2 3020insC did not associate with CRC or age at disease onset in the Finnish population. These data suggest that NOD2 3020insC alone might not be sufficient for CRC predisposition. MDM2 SNP309 was as common in the CRC cohort as in the healthy controls. Interesting trends, however, were observed, which after correction for multiple testing did not reach statistical significance. SNP309 was more common in female CRC patients and a trend towards an earlier age at disease onset was observed in women with SNP309. Subsequent studies have supported this observation and SNP309 could affect gender- or hormone-related tumorigenesis. Finally, a large-scale unbiased effort was designed to characterize the complete mutatome of CRC with microsatellite instability (MSI). Using an approach combining expression microarray and genome database searches, we were able to identify putative MSI target genes. Further characterization of one of the genes suggested that it might play a role also in microsatellite stable CRC and Peutz-Jeghers syndrome pathogenesis.

Mathematical models on the impact of noise and dyadic molecular structures on the properties of a cardiac myocyte

In cardiac myocytes (heart muscle cells), coupling of electric signal known as the action potential to contraction of the heart depends crucially on calcium-induced calcium release (CICR) in a microdomain known as the dyad. During CICR, the peak number of free calcium ions (Ca) present in the dyad is small, typically estimated to be within range 1-100. Since the free Ca ions mediate CICR, noise in Ca signaling due to the small number of free calcium ions influences Excitation-Contraction (EC) coupling gain. Noise in Ca signaling is only one noise type influencing cardiac myocytes, e.g., ion channels playing a central role in action potential propagation are stochastic machines, each of which gates more or less randomly, which produces gating noise present in membrane currents. How various noise sources influence macroscopic properties of a myocyte, how noise is attenuated and taken advantage of are largely open questions. In this thesis, the impact of noise on CICR, EC coupling and, more generally, macroscopic properties of a cardiac myocyte is investigated at multiple levels of detail using mathematical models. Complementarily to the investigation of the impact of noise on CICR, computationally-efficient yet spatially-detailed models of CICR are developed. The results of this thesis show that (1) gating noise due to the high-activity mode of L-type calcium channels playing a major role in CICR may induce early after-depolarizations associated with polymorphic tachycardia, which is a frequent precursor to sudden cardiac death in heart failure patients; (2) an increased level of voltage noise typically increases action potential duration and it skews distribution of action potential durations toward long durations in cardiac myocytes; and that (3) while a small number of Ca ions mediate CICR, Excitation-Contraction coupling is robust against this noise source, partly due to the shape of ryanodine receptor protein structures present in the cardiac dyad.

Mechanisms and molecular regulation of mammalian tooth replacement

In most non-mammalian vertebrates, such as fish and reptiles, teeth are replaced continuously. However, tooth replacement in most mammals, including human, takes place only once and further renewal is apparently inhibited. It is not known how tooth replacement is genetically regulated, and little is known on the physiological mechanism and evolutionary reduction of tooth replacement in mammals. In this study I have attempted to address these questions. In a rare human condition cleidocranial dysplasia, caused by a mutation in a Runt domain transcription factor Runx2, tooth replacement is continued. Runx2 mutant mice were used to investigate the molecular mechanisms of Runx2 function. Microarray analysis from dissected embryonic day 14 Runx2 mutant and wild type dental mesenchymes revealed many downstream targets of Runx2, which were validated using in situ hybridization and tissue culture methods. Wnt signaling inhibitor Dkk1 was identified as a candidate target, and in tissue culture conditions it was shown that Dkk1 is induced by FGF4 and this induction is Runx2 dependent. These experiments demonstrated a connection between Runx2, FGF and Wnt signaling in tooth development and possibly also in tooth replacement. The role of Wnt signaling in tooth replacement was further investigated by using a transgenic mouse model where Wnt signaling mediator β-catenin is continuously stabilized in dental epithelium. This stabilization led to activated Wnt signaling and to the formation of multiple enamel knots. In vitro and transplantation experiments were performed to examine the process of extra tooth formation. We showed that new teeth were continuously generated and that new teeth form from pre-existing teeth. A morphodynamic activator-inhibitor model was used to simulate enamel knot formation. By increasing the intrinsic production rate of the activator (β-catenin), the multiple enamel knot phenotype was reproduced by computer simulations. It was thus concluded that β-catenin acts as an upstream activator of enamel knots, closely linking Wnt signaling to the regulation of tooth renewal. As mice do not normally replace teeth, we used other model animals to investigate the physiological and genetic mechanisms of tooth replacement. Sorex araneus, the common shrew was earlier reported to have non-functional tooth replacement in all antemolar tooth positions. We showed by histological and gene expression studies that there is tooth replacement only in one position, the premolar 4 and that the deciduous tooth is diminished in size and disappears during embryogenesis without becoming functional. The growth rates of deciduous and permanent premolar 4 were measured and it was shown by competence inference that the early initiation of the replacement tooth in relation to the developmental stage of the deciduous tooth led to the inhibition of deciduous tooth morphogenesis. It was concluded that the evolutionary loss of deciduous teeth may involve the early activation of replacement teeth, which in turn suppress their predecessors. Mustela putorius furo, the ferret, has a dentition that resembles that of the human as ferrets have teeth that belong to all four tooth families, and all the antemolar teeth are replaced once. To investigate the replacement mechanism, histological serial sections from different embryonic stages were analyzed. It was noticed that tooth replacement is a process which involves the growth and detachment of the dental lamina from the lingual cervical loop of the deciduous tooth. Detachment of the deciduous tooth leads to a free successional dental lamina, which grows deeper into the mesenchyme, and later buds the replacement tooth. A careful 3D analysis of serial histological sections was performed and it was shown that replacement teeth are initiated from the successional dental lamina and not from the epithelium of the deciduous tooth. The molecular regulation of tooth replacement was studied and it was shown by examination of expression patterns of candidate regulatory genes that BMP/Wnt inhibitor Sostdc1 was strongly expressed in the buccal aspect of the dental lamina, and in the intersection between the detaching deciduous tooth and the successional dental lamina, suggesting a role for Sostdc1 in the process of detachment. Shh was expressed in the enamel knot and in the inner enamel epithelium in both generations of teeth supporting the view that the morphogenesis of both generations of teeth is regulated by similar mechanisms. In summary, histological and molecular studies on different model animals and transgenic mouse models were used to investigate tooth replacement. This thesis work has significantly contributed to the knowledge on the physiological mechanisms and molecular regulation of tooth replacement and its evolutionary suppression in mammals.

Prognostic molecular factors and algorithms in diffuse large B-cell lymphoma

Diffuse large B-cell lymphoma (DLBCL) is the most common of the non-Hodgkin lymphomas. As DLBCL is characterized by heterogeneous clinical and biological features, its prognosis varies. To date, the International Prognostic Index has been the strongest predictor of outcome for DLBCL patients. However, no biological characters of the disease are taken into account. Gene expression profiling studies have identified two major cell-of-origin phenotypes in DLBCL with different prognoses, the favourable germinal centre B-cell-like (GCB) and the unfavourable activated B-cell-like (ABC) phenotypes. However, results of the prognostic impact of the immunohistochemically defined GCB and non-GCB distinction are controversial. Furthermore, since the addition of the CD20 antibody rituximab to chemotherapy has been established as the standard treatment of DLBCL, all molecular markers need to be evaluated in the post-rituximab era. In this study, we aimed to evaluate the predictive value of immunohistochemically defined cell-of-origin classification in DLBCL patients. The GCB and non-GCB phenotypes were defined according to the Hans algorithm (CD10, BCL6 and MUM1/IRF4) among 90 immunochemotherapy- and 104 chemotherapy-treated DLBCL patients. In the chemotherapy group, we observed a significant difference in survival between GCB and non-GCB patients, with a good and a poor prognosis, respectively. However, in the rituximab group, no prognostic value of the GCB phenotype was observed. Likewise, among 29 high-risk de novo DLBCL patients receiving high-dose chemotherapy and autologous stem cell transplantation, the survival of non-GCB patients was improved, but no difference in outcome was seen between GCB and non-GCB subgroups. Since the results suggested that the Hans algorithm was not applicable in immunochemotherapy-treated DLBCL patients, we aimed to further focus on algorithms based on ABC markers. We examined the modified activated B-cell-like algorithm based (MUM1/IRF4 and FOXP1), as well as a previously reported Muris algorithm (BCL2, CD10 and MUM1/IRF4) among 88 DLBCL patients uniformly treated with immunochemotherapy. Both algorithms distinguished the unfavourable ABC-like subgroup with a significantly inferior failure-free survival relative to the GCB-like DLBCL patients. Similarly, the results of the individual predictive molecular markers transcription factor FOXP1 and anti-apoptotic protein BCL2 have been inconsistent and should be assessed in immunochemotherapy-treated DLBCL patients. The markers were evaluated in a cohort of 117 patients treated with rituximab and chemotherapy. FOXP1 expression could not distinguish between patients, with favourable and those with poor outcomes. In contrast, BCL2-negative DLBCL patients had significantly superior survival relative to BCL2-positive patients. Our results indicate that the immunohistochemically defined cell-of-origin classification in DLBCL has a prognostic impact in the immunochemotherapy era, when the identifying algorithms are based on ABC-associated markers. We also propose that BCL2 negativity is predictive of a favourable outcome. Further investigational efforts are, however, warranted to identify the molecular features of DLBCL that could enable individualized cancer therapy in routine patient care.

Development of Novel Assays for Measuring Different Molecular Forms of Prostate Specific Antigen

Prostate cancer (PCa) is the most commonly diagnosed non-skin cancer and second leading cause of cancer-related death of men in developed countries. Measurement of prostate specific antigen (PSA) is a very sensitive method for diagnosing and monitoring of prostate cancer (PCa), but the specificity needs improvement. Measurements of different molecular forms of PSA have been shown to improve differentiation between PCa and benign prostatic diseases. However, accurate measurement of some isoforms has not been achieved in previous assays. The aim of the present study was to develop new assays that reliably measure enzymatically active PSA, PSA-α1-chymotryposin (PSA-ACT) and PSA-α1-protease inhibitor (PSA-API), and to evaluate their diagnostic value. Double-label immunofluorometric assays using a novel monoclonal antibody (MAb) and another antibody to either free PSA (fPSA) or total PSA (tPSA) were developed and used to measure PSA-ACT and fPSA or tPSA at the same time. These assays provide enough sensitivity for measurement of PSA-ACT in sera with low PSA levels. The results obtained confirmed that proportion of PSA-ACT to tPSA (%PSA-ACT) was as useful as proportion of fPSA to tPSA (%fPSA) for discrimination between PCa and benign prostatic hyperplasia (BPH). We developed an immunoassay for detection of PSA-API based on proximity ligation, which improved assay sensitivity 10-fold compared with conventional assays. Our results confirmed previous findings that the PSA-API level is somewhat lower in men with than without PCa, and the combination of %fPSA and proportion of PSA-API to tPSA (%PSA-API) provides diagnostic improvement compared with either method alone. Assays based on this principle should be applicable to other immunoassays in which the nonspecific background is a problem. An immunopeptidometric sandwich assay (IPMA) was developed to measure the enzymatically active PSA. This assay showed high specificity, but sensitivity was not good enough for measurement of PSA concentrations in the gray zone, 2-10 µg/L, in which tPSA does not efficiently differentiate between PCa and BPH. We further developed a solid-phase proximity ligation immunoassay, which provided a 10-fold improvement in sensitivity. This proof of concept study shows that peptides reacting with proteins are potentially useful for sensitive and specific measurement of protein variants for which specific MAbs cannot be obtained.

Neurological manifestations and molecular genetics of acute intermittent porphyria in North Western Russia

Acute intermittent porphyria (AIP, MIM #176000) is an inherited metabolic disease due to a partial deficiency of the third enzyme, hydroxymethylbilane synthase (HMBS, EC: 4.3.1.8), in the haem biosynthesis. Neurological symptoms during an acute attack, which is the major manifestation of AIP, are variable and relatively rare, but may endanger a patient's life. In the present study, 12 Russian and two Finnish AIP patients with severe neurological manifestations during an acute attack were studied prospectively from 1995 to 2006. Autonomic neuropathy manifested as abdominal pain (88%), tachycardia (94%), hypertension (75%) and constipation (88%). The most common neurological sign was acute motor peripheral neuropathy (PNP, 81%) often associated with neuropathic sensory loss (54%) and CNS involvement (85%). Despite heterogeneity of the neurological manifestations in our patients with acute porphyria, the major pattern of PNP associated with abdominal pain, dysautonomia, CNS involvement and mild hepatopathy could be demonstrated. If more strict inclusion criteria for biochemical abnormalities (>10-fold increase in excretion of urinary PBG) are applied, neurological manifestations in an acute attack are probably more homogeneous than described previously, which suggests that some of the neurological patients described previously may not have acute porphyria but rather secondary porphyrinuria. Screening for acute porphyria using urinary PBG is useful in a selected group of neurological patients with acute PNP or encephalopathy and seizures associated with pain and dysautonomia. Clinical manifestations and the outcome of acute attacks were used as a basis for developing a 30-score scale of the severity of an acute attack. This scale can easily be used in clinical practice and to standardise the outcome of an attack. Degree of muscle weakness scored by MRC, prolonged mechanical ventilation, bulbar paralysis, impairment of consciousness and hyponatraemia were important signs of a poor prognosis. Arrhythmia was less important and autonomic dysfunction, severity of pain and mental symptoms did not affect the outcome. The delay in the diagnosis and repeated administrations of precipitating factors were the main cause of proceeding of an acute attack into pareses and severe CNS involvement and a fatal outcome in two patients. Nerve conduction studies and needle EMG were performed in eleven AIP patients during an acute attack and/or in remission. Nine patients had severe PNP and two patients had an acute encephalopathy but no clinically evident PNP. In addition to axonopathy, features suggestive of demyelination could be demonstrated in patients with severe PNP during an acute attack. PNP with a moderate muscle weakness was mainly pure axonal. Sensory involvement was common in acute PNP and could be subclinical. Decreased conduction velocities with normal amplitudes of evoked potentials during acute attacks with no clinically evident PNP indicated subclinical polyneuropathy. Reversible symmetrical lesions comparable with posterior reversible encephalopathy syndrome (PRES) were revealed in two patients' brain CT or MRI during an acute attack. In other five patients brain MRI during or soon after the symptoms was normal. The frequency of reversible brain oedema in AIP is probably under-estimated since it may be short-lasting and often indistinguishable on CT or MRI. In the present study, nine different mutations were identified in the HMBS gene in 11 unrelated Russian AIP patients from North Western Russia and their 32 relatives. AIP was diagnosed in nine symptom-free relatives. The majority of the mutations were family-specific and confirmed allelic heterogeneity also among Russian AIP patients. Three mutations, c.825+5G>C, c.825+3_825+6del and c.770T>C, were novel. Six mutations, c.77G>A (p.R26H), c.517C>T (p.R173W), c.583C>T (p.R195C), c.673C>T (p.R225X), c.739T>C (p.C247R) and c.748G>C (p.E250A), have previously been identified in AIP patients from Western and other Eastern European populations. The effects of novel mutations were studied by amplification and sequencing of the reverse-transcribed total RNA obtained from the patients' lymphoblastoid or fibroblast cell lines. The mutations c.825+5G>C and c.770T>C resulted in varyable amounts of abnormal transcripts, r.822_825del (p.C275fsX2) and [r.770u>c, r.652_771del, r.613_771del (p.L257P, p.G218_L257del, p.I205_L257del)]. All mutations demonstrated low residual activities (0.1-1.3 %) when expressed in COS-1 cells confirming the causality of the mutations and the enzymatic defect of the disease. The clinical outcome, prognosis and correlation between the HMBS genotype and phenotype were studied in 143 Finnish and Russian AIP patients with ten mutations (c.33G>T, c.97delA, InsAlu333, p.R149X, p.R167W, p.R173W, p.R173Q, p.R225G, p.R225X, c.1073delA) and more than six patients in each group. The patients were selected from the pool of 287 Finnish AIP patients presented in a Finnish Porphyria Register (1966-2003) and 23 Russian AIP patients (diagnosed 1995-2003). Patients with the p.R167W and p.R225G mutations showed lower penetrance (19% and 11%) and the recurrence rate (33% and 0%) in comparison to the patients with other mutations (range 36 to 67% and 0 to 66%, respectively), as well as milder biochemical abnormalities [urinary porphobilinogen 47±10 vs. 163±21 mol/L, p<0.001; uroporphyrin 130±40 vs. 942±183 nmol/L, p<0.001] suggesting a milder form of AIP in these patients. Erythrocyte HMBS activity did not correlate with the porphobilinogen excretion in remission or the clinical of the disease. In all AIP severity patients, normal PBG excretion predicted freedom from acute attacks. Urinary PBG excretion together with gender, age at the time of diagnosis and mutation type could predict the likelihood of acute attacks in AIP patients.

Monte Carlo Simulations of Molecular Clusters in Nucleation

A better understanding of the limiting step in a first order phase transition, the nucleation process, is of major importance to a variety of scientific fields ranging from atmospheric sciences to nanotechnology and even to cosmology. This is due to the fact that in most phase transitions the new phase is separated from the mother phase by a free energy barrier. This barrier is crossed in a process called nucleation. Nowadays it is considered that a significant fraction of all atmospheric particles is produced by vapor-to liquid nucleation. In atmospheric sciences, as well as in other scientific fields, the theoretical treatment of nucleation is mostly based on a theory known as the Classical Nucleation Theory. However, the Classical Nucleation Theory is known to have only a limited success in predicting the rate at which vapor-to-liquid nucleation takes place at given conditions. This thesis studies the unary homogeneous vapor-to-liquid nucleation from a statistical mechanics viewpoint. We apply Monte Carlo simulations of molecular clusters to calculate the free energy barrier separating the vapor and liquid phases and compare our results against the laboratory measurements and Classical Nucleation Theory predictions. According to our results, the work of adding a monomer to a cluster in equilibrium vapour is accurately described by the liquid drop model applied by the Classical Nucleation Theory, once the clusters are larger than some threshold size. The threshold cluster sizes contain only a few or some tens of molecules depending on the interaction potential and temperature. However, the error made in modeling the smallest of clusters as liquid drops results in an erroneous absolute value for the cluster work of formation throughout the size range, as predicted by the McGraw-Laaksonen scaling law. By calculating correction factors to Classical Nucleation Theory predictions for the nucleation barriers of argon and water, we show that the corrected predictions produce nucleation rates that are in good comparison with experiments. For the smallest clusters, the deviation between the simulation results and the liquid drop values are accurately modelled by the low order virial coefficients at modest temperatures and vapour densities, or in other words, in the validity range of the non-interacting cluster theory by Frenkel, Band and Bilj. Our results do not indicate a need for a size dependent replacement free energy correction. The results also indicate that Classical Nucleation Theory predicts the size of the critical cluster correctly. We also presents a new method for the calculation of the equilibrium vapour density, surface tension size dependence and planar surface tension directly from cluster simulations. We also show how the size dependence of the cluster surface tension in equimolar surface is a function of virial coefficients, a result confirmed by our cluster simulations.

Monte Carlo simulations of molecular clusters in nucleation

A better understanding of the limiting step in a first order phase transition, the nucleation process, is of major importance to a variety of scientific fields ranging from atmospheric sciences to nanotechnology and even to cosmology. This is due to the fact that in most phase transitions the new phase is separated from the mother phase by a free energy barrier. This barrier is crossed in a process called nucleation. Nowadays it is considered that a significant fraction of all atmospheric particles is produced by vapor-to liquid nucleation. In atmospheric sciences, as well as in other scientific fields, the theoretical treatment of nucleation is mostly based on a theory known as the Classical Nucleation Theory. However, the Classical Nucleation Theory is known to have only a limited success in predicting the rate at which vapor-to-liquid nucleation takes place at given conditions. This thesis studies the unary homogeneous vapor-to-liquid nucleation from a statistical mechanics viewpoint. We apply Monte Carlo simulations of molecular clusters to calculate the free energy barrier separating the vapor and liquid phases and compare our results against the laboratory measurements and Classical Nucleation Theory predictions. According to our results, the work of adding a monomer to a cluster in equilibrium vapour is accurately described by the liquid drop model applied by the Classical Nucleation Theory, once the clusters are larger than some threshold size. The threshold cluster sizes contain only a few or some tens of molecules depending on the interaction potential and temperature. However, the error made in modeling the smallest of clusters as liquid drops results in an erroneous absolute value for the cluster work of formation throughout the size range, as predicted by the McGraw-Laaksonen scaling law. By calculating correction factors to Classical Nucleation Theory predictions for the nucleation barriers of argon and water, we show that the corrected predictions produce nucleation rates that are in good comparison with experiments. For the smallest clusters, the deviation between the simulation results and the liquid drop values are accurately modelled by the low order virial coefficients at modest temperatures and vapour densities, or in other words, in the validity range of the non-interacting cluster theory by Frenkel, Band and Bilj. Our results do not indicate a need for a size dependent replacement free energy correction. The results also indicate that Classical Nucleation Theory predicts the size of the critical cluster correctly. We also presents a new method for the calculation of the equilibrium vapour density, surface tension size dependence and planar surface tension directly from cluster simulations. We also show how the size dependence of the cluster surface tension in equimolar surface is a function of virial coefficients, a result confirmed by our cluster simulations.

Molecular Dynamics Simulations of Homogeneous Nucleation

Nucleation is the first step in a phase transition where small nuclei of the new phase start appearing in the metastable old phase, such as the appearance of small liquid clusters in a supersaturated vapor. Nucleation is important in various industrial and natural processes, including atmospheric new particle formation: between 20 % to 80 % of atmospheric particle concentration is due to nucleation. These atmospheric aerosol particles have a significant effect both on climate and human health. Different simulation methods are often applied when studying things that are difficult or even impossible to measure, or when trying to distinguish between the merits of various theoretical approaches. Such simulation methods include, among others, molecular dynamics and Monte Carlo simulations. In this work molecular dynamics simulations of the homogeneous nucleation of Lennard-Jones argon have been performed. Homogeneous means that the nucleation does not occur on a pre-existing surface. The simulations include runs where the starting configuration is a supersaturated vapor and the nucleation event is observed during the simulation (direct simulations), as well as simulations of a cluster in equilibrium with a surrounding vapor (indirect simulations). The latter type are a necessity when the conditions prevent the occurrence of a nucleation event in a reasonable timeframe in the direct simulations. The effect of various temperature control schemes on the nucleation rate (the rate of appearance of clusters that are equally able to grow to macroscopic sizes and to evaporate) was studied and found to be relatively small. The method to extract the nucleation rate was also found to be of minor importance. The cluster sizes from direct and indirect simulations were used in conjunction with the nucleation theorem to calculate formation free energies for the clusters in the indirect simulations. The results agreed with density functional theory, but were higher than values from Monte Carlo simulations. The formation energies were also used to calculate surface tension for the clusters. The sizes of the clusters in the direct and indirect simulations were compared, showing that the direct simulation clusters have more atoms between the liquid-like core of the cluster and the surrounding vapor. Finally, the performance of various nucleation theories in predicting simulated nucleation rates was investigated, and the results among other things highlighted once again the inadequacy of the classical nucleation theory that is commonly employed in nucleation studies.