Mammalian ribosomal and chaperone protein RPS3A counteracts α-synuclein aggregation and toxicity in a yeast model system.

Accumulation of aggregated forms of αSyn (α-synuclein) into Lewy bodies is a known hallmark associated with neuronal cell death in Parkinson's disease. When expressed in the yeast Saccharomyces cerevisiae, αSyn interacts with the plasma membrane, forms inclusions and causes a concentration-dependent growth defect. We have used a yeast mutant, cog6Δ, which is particularly sensitive to moderate αSyn expression, for screening a mouse brain-specific cDNA library in order to identify mammalian proteins that counteract αSyn toxicity. The mouse ribosomal and chaperone protein RPS3A was identified as a suppressor of αSyn [WT (wild-type) and A53T] toxicity in yeast. We demonstrated that the 50 N-terminal amino acids are essential for this function. The yeast homologues of RPS3A were not effective in suppressing the αSyn-induced growth defect, illustrating the potential of our screening system to identify modifiers that would be missed using yeast gene overexpression as the first screening step. Co-expression of mouse RPS3A delayed the formation of αSyn-GFP inclusions in the yeast cells. The results of the present study suggest that the recently identified extraribosomal chaperonin function of RPS3A also acts on the neurodegeneration-related protein αSyn and reveal a new avenue for identifying promising candidate mammalian proteins involved in αSyn functioning.

In vitro phosphorylation does not influence the aggregation kinetics of WT α-synuclein in contrast to its phosphorylation mutants.

The aggregation of alpha-synuclein (α-SYN) into fibrils is characteristic for several neurodegenerative diseases, including Parkinson's disease (PD). Ninety percent of α-SYN deposited in Lewy Bodies, a pathological hallmark of PD, is phosphorylated on serine129. α-SYN can also be phosphorylated on tyrosine125, which is believed to regulate the membrane binding capacity and thus possibly its normal function. A better understanding of the effect of phosphorylation on the aggregation of α-SYN might shed light on its role in the pathogenesis of PD. In this study we compare the aggregation properties of WT α-SYN with the phospho-dead and phospho-mimic mutants S129A, S129D, Y125F and Y125E and in vitro phosphorylated α-SYN using turbidity, thioflavin T and circular dichroism measurements as well as transmission electron microscopy. We show that the mutants S129A and S129D behave similarly compared to wild type (WT) α-SYN, while the mutants Y125F and Y125E fibrillate significantly slower, although all mutants form fibrillar structures similar to the WT protein. In contrast, in vitro phosphorylation of α-SYN on either S129 or Y125 does not significantly affect the fibrillization kinetics. Moreover, FK506 binding proteins (FKBPs), enzymes with peptidyl-prolyl cis-trans isomerase activity, still accelerate the aggregation of phosphorylated α-SYN in vitro, as was shown previously for WT α-SYN. In conclusion, our results illustrate that phosphorylation mutants can display different aggregation properties compared to the more biologically relevant phosphorylated form of α-SYN.

The characterization of the nuclear dynamics of syntenin-2, a PIP2 binding PDZ protein.

Cellular signaling is largely dependent on the presence, that is, assembly/disassembly, of supramolecular complexes. Postsynaptic density protein, Discs-large, Zona occludens (PDZ) domains play important roles in the assembly of various signaling complexes. Syntenin-2 (S2) is a PDZ protein that interacts with nuclear phosphatidylinositol 4,5-bisphosphate (PIP2 ). Although nuclear lipids emerge as key players in nuclear processes, the global significance of nuclear phosphoinositide-protein interactions is still poorly understood. Those phosphoinositide-protein interactions that have been studied in detail appear to have profound physiological effects. To our knowledge none of these were investigated by dynamic studies such as Fluorescence Correlation Spectroscopy (FCS), Fluorescence Cross-Correlation Spectroscopy (FCCS), or Fluorescence Recovery After Photobleaching (FRAP). Although the exact function of S2 is unknown, siRNA experiments suggest that this PDZ protein plays a role in the organization of nuclear PIP2 , cell division, and cell survival. As a consequence of its PIP2 interaction, its reported self-association in a yeast two-hybrid study and its speculated interaction with many, yet unidentified, proteins one can hypothesize that S2 plays an important role in cell signaling. Therefore, we studied the dynamics of S2 using FCS, FCCS, and FRAP, utilizing an active truncated form deleted for the first 94 amino acids (S2-ΔN). We showed that S2-ΔN self-associates and is distributed in three groups. One immobile group, one slow diffusing group, which interacts with the nuclear environment and one fast diffusing group, which is not incorporated in high molecular weight complexes. In addition, our FCS and FRAP measurements on S2-ΔN mutants affected in their PIP2 binding showed that PIP2 plays an important role in the distribution of S2-ΔN among these groups, and favors the enrichment of S2-ΔN in the slow diffusing and immobile group. This work indicates that S2 relies on nuclear PIP2 to interact with practically immobile structures, possibly chromatin.

The transcriptional co-activator LEDGF/p75 displays a dynamic scan-and-lock mechanism for chromatin tethering.

Nearly all cellular and disease related functions of the transcriptional co-activator lens epithelium-derived growth factor (LEDGF/p75) involve tethering of interaction partners to chromatin via its conserved integrase binding domain (IBD), but little is known about the mechanism of in vivo chromatin binding and tethering. In this work we studied LEDGF/p75 in real-time in living HeLa cells combining different quantitative fluorescence techniques: spot fluorescence recovery after photobleaching (sFRAP) and half-nucleus fluorescence recovery after photobleaching (hnFRAP), continuous photobleaching, fluorescence correlation spectroscopy (FCS) and an improved FCS method to study diffusion dependence of chromatin binding, tunable focus FCS. LEDGF/p75 moves about in nuclei of living cells in a chromatin hopping/scanning mode typical for transcription factors. The PWWP domain of LEDGF/p75 is necessary, but not sufficient for in vivo chromatin binding. After interaction with HIV-1 integrase via its IBD, a general protein-protein interaction motif, kinetics of LEDGF/p75 shift to 75-fold larger affinity for chromatin. The PWWP is crucial for locking the complex on chromatin. We propose a scan-and-lock model for LEDGF/p75, unifying paradoxical notions of transcriptional co-activation and lentiviral integration targeting.

Tryptophan side chain conformers monitored by NMR and time-resolved fluorescence spectroscopies.

We have inserted a tryptophan (F77W) in the core of the regulatory domain of cardiac troponin C (cNTnC), and previously determined the structure of this mutant with and without the cosolvent trifluoroethanol (TFE). Interestingly, the orientations of the indole side chain of the Trp are in opposite directions in the two structures (Julien et al., Protein Sci 2009; 18:1165-1174). Fluorescence decay experiments for single Trp-containing proteins often show several lifetimes, which have been interpreted as reflecting conformational heterogeneity of the Trp side chain resulting from different rotamers. To test this interpretation, we monitored the effect of TFE on wild type, F77W and F77W-V82A calcium-saturated cNTnC using 2D (13)C-HSQC NMR and time-correlated single photon counting fluorescence spectroscopies. The time dependence of the Trp fluorescence decay was fit with three lifetimes. Addition of TFE caused a gradual, but limited decrease of the lifetimes due to dynamic quenching. For F77W cNTnC, the amplitude fractions of the lifetimes also changed upon addition of TFE-the long lifetime increased from 13 to 29%, while the middle lifetime decreased from 63 to 50% and the short lifetime remained relatively unchanged. For F77W-V82A cNTnC, comparable NMR changes are observed, confirming the switch in rotamer conformation, but only much smaller changes in fluorescence decay parameters were detected. These data indicate that the balance between the rotamer states can be changed without changing the lifetime amplitude fractions appreciably, and suggest that the environment(s) of the indole ring, responsible for the different lifetimes, can result from factors other than the intrinsic rotamer state of the tryptophan.

Inhibition of FK506 binding proteins reduces alpha-synuclein aggregation and Parkinson's disease-like pathology.

alpha-Synuclein (alpha-SYN) is a key player in the pathogenesis of Parkinson's disease (PD). In pathological conditions, the protein is present in a fibrillar, aggregated form inside cytoplasmic inclusions called Lewy bodies. Members of the FK506 binding protein (FKBP) family are peptidyl-prolyl isomerases that were shown recently to accelerate the aggregation of alpha-SYN in vitro. We now established a neuronal cell culture model for synucleinopathy based on oxidative stress-induced alpha-SYN aggregation and apoptosis. Using high-content analysis, we examined the role of FKBPs in aggregation and apoptotic cell death. FK506, a specific inhibitor of this family of proteins, inhibited alpha-SYN aggregation and neuronal cell death in this synucleinopathy model dose dependently. Knockdown of FKBP12 or FKBP52 reduced the number of alpha-SYN aggregates and protected against cell death, whereas overexpression of FKBP12 or FKBP52 accelerated both aggregation of alpha-SYN and cell death. Thus, FK506 likely targets FKBP members in the cell culture model. Furthermore, oral administration of FK506 after viral vector-mediated overexpression of alpha-SYN in adult mouse brain significantly reduced alpha-SYN aggregate formation and neuronal cell death. Our data explain previously described neuroregenerative and neuroprotective effects of immunophilin ligands and validate FKBPs as a novel drug target for the causative treatment of PD.

SAP-like domain in nucleolar spindle associated protein mediates mitotic chromosome loading as well as interphase chromatin interaction.

Nucleolar spindle associated protein (NuSAP) is a microtubule-stabilizing protein that localizes to chromosome arms and chromosome-proximal microtubules during mitosis and to the nucleus, with enrichment in the nucleoli, during interphase. The critical function of NuSAP is underscored by the finding that its depletion in HeLa cells results in various mitotic defects. Moreover, NuSAP is found overexpressed in multiple cancers and its expression levels often correlate with the aggressiveness of cancer. Due to its localization on chromosome arms and combination of microtubule-stabilizing and DNA-binding properties, NuSAP takes a special place within the extensive group of spindle assembly factors. In this study, we identify a SAP-like domain that shows DNA binding in vitro with a preference for dsDNA. Deletion of the SAP-like domain abolishes chromosome arm binding of NuSAP during mitosis, but is not sufficient to abrogate its chromosome-proximal localization after anaphase onset. Fluorescence recovery after photobleaching experiments revealed the highly dynamic nature of this NuSAP-chromatin interaction during mitosis. In interphase cells, NuSAP also interacts with chromatin through its SAP-like domain, as evident from its enrichment on dense chromatin regions and intranuclear mobility, measured by fluorescence correlation spectroscopy. The obtained results are in agreement with a model where NuSAP dynamically stabilizes newly formed microtubules on mitotic chromosomes to enhance chromosome positioning without immobilizing these microtubules. Interphase NuSAP-chromatin interaction suggests additional functions for NuSAP, as recently identified for other nuclear spindle assembly factors with a role in gene expression or DNA damage response.

Early aggregation steps in alpha-synuclein as measured by FCS and FRET: evidence for a contagious conformational change.

The kinetics of aggregation of alpha-synuclein are usually studied by turbidity or Thio-T fluorescence. Here we follow the disappearance of monomers and the formation of early oligomers using fluorescence correlation spectroscopy. Alexa488-labeled A140C-synuclein was used as a fluorescent probe in trace amounts in the presence of excess unlabeled alpha-synuclein. Repeated short measurements produce a distribution of diffusion coefficients. Initially, a sharp peak is obtained corresponding to monomers, followed by a distinct transient population and the gradual formation of broader-sized distributions of higher oligomers. The kinetics of aggregation can be followed by the decreasing number of fast-diffusing species. Both the disappearance of fast-diffusing species and the appearance of turbidity can be fitted to the Finke-Watzky equation, but the apparent rate constants obtained are different. This reflects the fact that the disappearance of fast species occurs largely during the lag phase of turbidity development, due to the limited sensitivity of turbidity to the early aggregation process. The nucleation of the early oligomers is concentration-dependent and accompanied by a conformational change that precedes beta-structure formation, and can be visualized using fluorescence resonance energy transfer between the donor-labeled N-terminus and the acceptor-labeled cysteine in the mutant A140C.

The conformation and the aggregation kinetics of α-synuclein depend on the proline residues in its C-terminal region.

The neuronal protein α-synuclein (α-syn) plays a central role in Parkinson's disease (PD). The pathological features of PD are the loss of dopaminergic neurons in the substantia nigra pars compacta and the presence of Lewy bodies. The C-terminal domain of α-syn is characterized by the presence of 15 acidic amino acids and all five proline residues of the protein (P108, P117, P120, P128, and P138). The aggregation of this natively unfolded protein is accelerated in vitro by FK506 binding proteins (FKBPs) showing peptidyl-prolyl cis-trans isomerase activity. These proteins catalyze the cis-trans conformational change of the X-Pro peptide bond, often a rate-limiting step in protein folding. The acceleration of the folding of α-syn by FKBPs may accelerate disease-associated aggregation. To further elucidate the role of the proline residues in the conformation and aggregation of α-syn, we constructed several mutants of α-syn in which one or more proline residues are mutated to alanine via site-directed mutagenesis. For this purpose, we produced and purified His-WT α-syn, a recombinant α-syn with a polyhistidine tag (six His residues) and a linker, and a number of Pro-to-Ala mutants. The aggregation kinetics of these mutants and His-WT α-syn were studied by turbidity, thioflavin T fluorescence, and CD measurements. We can conclude that mutation of the proline residues to alanine accelerates the aggregation kinetics of α-syn while all proline mutants formed fibrils similar to His-WT α-syn, as visualized via transmission electron microscopy. We also demonstrate that the accelerating effect of hFKBP12 is abolished via removal of the proline residues from the C-terminus. Finally, we show that the mutant of His α-syn with all five proline residues mutated to alanine is more structured (more α-helix) than His-WT α-syn, indicating the role of the Pro residues as potential helix breakers in the inhibitory conformation of the C-terminus.

Tryptophan conformations associated with partial unfolding in ribonuclease T1.

The origin of the biexponential fluorescence decay of Trp in ribonuclease T1 under mildly destabilizing conditions, such as increased pH or temperature, or the presence of detergent, is still not understood. We have performed two extended replica-exchange molecular dynamics simulations to obtain a detailed representation of the native state at two protonation states corresponding to a high and low pH. At high pH, the appearance of partially unfolded states is evident. We found that this pH-induced destabilization originates from increased global repulsion as well as reduced local favorable electrostatic interactions and reduced H-bonding strength of His(27), His(40), and His(92). At high pH, alternative tryptophan rotamers appear and are linked to a distorted environment of the tryptophan, which also acts as a separate source of ground-state heterogeneity. The total population of these alternative conformations agrees well with the amplitude of the experimentally observed secondary fluorescence lifetime.

The high-affinity peptidoglycan binding domain of Pseudomonas phage endolysin KZ144.

The binding affinity of the N-terminal peptidoglycan binding domain of endolysin KZ144 (PBD(KZ)), originating from Pseudomonas aeruginosa bacteriophage varphiKZ, has been examined using a fusion protein of PBD(KZ) and green fluorescent protein (PBD(KZ)-GFP). A fluorescence recovery after photobleaching analysis of bound PBD(KZ)-GFP molecules showed less than 10% fluorescence recovery in the bleached area within 15 min. Surface plasmon resonance analysis confirmed this apparent high binding affinity revealing an equilibrium affinity constant of 2.95 x 10(7)M(-1) for the PBD(KZ)-peptidoglycan interaction. This unique domain, which binds to the peptidoglycan of all tested Gram-negative species, was harnessed to improve the specific activity of the peptidoglycan hydrolase domain KMV36C. The chimeric peptidoglycan hydrolase (PBD(KZ)-KMV36C) exhibits a threefold higher specific activity than the native catalytic domain (KMV36C). These results demonstrate that the modular assembly of functional domains is a rational approach to improve the specific activity of endolysins from phages infecting Gram-negatives.

Virus evolution reveals an exclusive role for LEDGF/p75 in chromosomal tethering of HIV.

Retroviruses by definition insert their viral genome into the host cell chromosome. Although the key player of retroviral integration is viral integrase, a role for cellular cofactors has been proposed. Lentiviral integrases use the cellular protein LEDGF/p75 to tether the preintegration complex to the chromosome, although the existence of alternative host proteins substituting for the function of LEDGF/p75 in integration has been proposed. Truncation mutants of LEDGF/p75 lacking the chromosome attachment site strongly inhibit HIV replication by competition for the interaction with integrase. In an attempt to select HIV strains that can overcome the inhibition, we now have used T-cell lines that stably express a C-terminal fragment of LEDGF/p75. Despite resistance development, the affinity of integrase for LEDGF/p75 is reduced and replication kinetics in human primary T cells is impaired. Detection of the integrase mutations A128T and E170G at key positions in the LEDGF/p75-integrase interface provides in vivo evidence for previously reported crystallographic data. Moreover, the complementary inhibition by LEDGF/p75 knockdown and mutagenesis at the integrase-LEDGF/p75 interface points to the incapability of HIV to circumvent LEDGF/p75 function during proviral integration. Altogether, the data provide a striking example of the power of viral molecular evolution. The results underline the importance of the LEDGF/p75 HIV-1 interplay as target for innovative antiviral therapy. Moreover, the role of LEDGF/p75 in targeting integration will stimulate research on strategies to direct gene therapy vectors into safe landing sites.

Biliverdin reductase is a transporter of haem into the nucleus and is essential for regulation of HO-1 gene expression by haematin.

hBVR (human biliverdin reductase) is an enzyme that reduces biliverdin (the product of haem oxygenases HO-1 and HO-2 activity) to the antioxidant bilirubin. It also functions as a kinase and as a transcription factor in the MAPK (mitogen-activated protein kinase) signalling cascade. Fluorescence correlation spectroscopy was used to investigate the mobility of hBVR in living cells and its function in the nuclear transport of haematin for induction of HO-1. In transiently transfected HeLa cells only kinase-competent hBVR translocates to the nucleus. A reduced mobility in the nucleus of haematin-treated cells suggests formation of an hBVR-haematin complex and its further association with large nuclear components. The binding of haematin is specific, with the formation of a 1:1 molar complex, and the C-terminal 7-residue fragment KYCCSRK(296) of hBVR contributes to the binding. The following data suggest formation of dynamic complexes of hBVR-haematin with chromatin: (i) the reduction of hBVR mobility in the presence of haematin is greater in heterochromatic regions than in euchromatic domains and (ii) hBVR mobility is not retarded by haematin in nuclear lysates that contain only soluble factors. Moreover, hBVR kinase activity is stimulated in the presence of double-stranded DNA fragments corresponding to HO-1 antioxidant and HREs (hypoxia response elements), as well as by haematin. Experiments with nuclear localization, export signal mutants and si-hBVR [siRNA (small interfering RNA) specific to hBVR] indicate that nuclear localization of hBVR is required for induction of HO-1 by haematin. Because gene regulation is energy-dependent and haematin regulates gene expression, our data suggest that hBVR functions as an essential component of the regulatory mechanisms for haem-responsive transcriptional activation.

Dark states in monomeric red fluorescent proteins studied by fluorescence correlation and single molecule spectroscopy.

Monomeric red fluorescent proteins (mRFPs) have become indispensable tools for studying protein dynamics, interactions and functions in the cellular environment. Their emission spectrum can be well separated from other fluorescent proteins, and their monomeric structure preserves the natural function of fusion proteins. However, previous photophysical studies of some RFPs have shown the presence of light-induced dark states that can complicate the interpretation of cellular experiments. In this article, we extend these studies to mRFP1, mCherry, and mStrawberry by means of fluorescence correlation spectroscopy and prove that this light-driven intensity flickering also occurs in these proteins. Furthermore, we show that the flickering in these proteins is pH-dependent. Single molecule spectroscopy revealed reversible transitions from a bright to a dark state in several timescales, even up to seconds. Time-resolved fluorescence spectroscopy showed multiexponential decays, consistent with a "loose" conformation. We offer a structural basis for the fluorescence flickering using known crystal structures and point out that the environment of Glu-215 is critical for the pH dependence of the flickering in RFPs. We apply dual-color fluorescence correlation spectroscopy inside live cells to prove that this flickering can seriously hamper cellular measurements if the timescales of the flickering and diffusion are not well separated.

An unusual red-edge excitation and time-dependent Stokes shift in the single tryptophan mutant protein DD-carboxypeptidase from Streptomyces: the role of dynamics and tryptophan rotamers.

The fluorescence emission of the single tryptophan (W233) of the mutant protein DD-carboxypeptidase from streptomyces is characterized by a red-edge excitation shift (REES), i.e., the phenomenon that the wavelength of maximum emission depends on the excitation wavelength. This phenomenon is an indication for a strongly reduced dynamic environment of the single tryptophan, which has a very low accessibility to the solvent. The REES shows, however, an unusual temperature and time dependence. This, together with the fluorescence lifetime analysis, showing three resolvable lifetimes, can be explained by the presence of three rotameric states that can be identified using the Dead-End Elimination method. The three individual lifetimes increase with increasing emission wavelength, indicating the presence of restricted protein dynamics within the rotameric states. This is confirmed by time-resolved anisotropy measurements that show dynamics within the rotamers but not among the rotamers. The global picture is that of a protein with a single buried tryptophan showing strongly restricted dynamics within three distinct rotameric states with different emission spectra and an anisotropic environment.

How do rotameric conformations influence the time-resolved fluorescence of tryptophan in proteins? A perspective based on molecular modeling and quantum chemistry.

We discuss the dynamics of tryptophan rotamers in the context of the non-exponential fluorescence decay in proteins. The central question is: how does the ground-state conformational heterogeneity influence the time evolution of tryptophan fluorescence? This problem is examined here from the theoretical perspective. Three methods at different levels of theory, and with different scopes and computational requirements are reviewed. The Dead-end elimination method is limited to side-chain dynamics and provides an efficient way to detect the stable tryptophan rotamers in a protein. Its application to the study of heterogeneous emission characteristics is illustrated. Molecular dynamics is aimed at the full phase space of the macromolecule in solution, but must rely on classical force fields and laws of evolution. We examine to what extent the molecular mechanics paradigm yields sufficiently accurate thermodynamic results, and what are the possible kinetic implications. Finally Quantum Chemistry is the only theoretical method that allows a direct assessment of the excited states. It is necessarily restricted to small molecular systems, and thus must be used in a hybrid combination with classical methods and electrostatic models. So far understanding of the emitting state has greatly progressed as a result of these calculations, but the actual treatment of the photophysical decay processes at the quantum level has not yet really started.

FK506 binding protein 12 differentially accelerates fibril formation of wild type alpha-synuclein and its clinical mutants A30P or A53T.

Aggregation of alpha-synuclein (alpha-SYN) plays a key role in Parkinson's disease. We have previously shown that aggregation of alpha-SYN in vitro is accelerated by addition of FK506 binding proteins (FKBP) and that this effect can be counteracted by FK506, a specific inhibitor of these enzymes. In this paper, we investigated in detail the effect of FKBP12 on early aggregation and on fibril formation of wild-type, A53T and A30P alpha-SYN. FKBP12 has a much smaller effect on the fibril formation of these two clinical mutants alpha-SYN. Using an inactive enzyme, we were able to discriminate between catalytic and non-catalytic effects that differentially influence the two processes. A model explaining non-linear concentration dependencies is proposed.

Prolyl oligopeptidase stimulates the aggregation of alpha-synuclein.

Despite its thorough enzymological and biochemical characterization the exact function of prolyl oligopeptidase (PO, E.C. 3.4.21.26) remains unclear. The positive effect of PO inhibitors on learning and memory in animal models for amnesia, enzyme activity measurements in patient samples and (neuro)peptide degradation studies link the enzyme with neurodegenerative disorders. The brain protein alpha-synuclein currently attracts much attention because of its proposed role in the pathology of Parkinson's disease. A fundamental question concerns how the essentially disordered protein is transformed into the highly organized fibrils that are found in Lewy bodies, the hallmarks of Parkinson's disease. Using gel electrophoresis and MALDI TOF/TOF mass spectrometry we investigated the possibility of alpha-synuclein as a PO substrate. We found that in vitro incubation of the protein with PO did not result in truncation of full-length alpha-synuclein. Surprisingly, however, we found an acceleration of the aggregation process of alpha-synuclein using turbidity measurements that was reversed by specific inhibitors of PO enzymatic activity. If PO displays this activity also in vivo, PO inhibitors might have an effect on neurodegenerative disorders through a decrease in the aggregation of alpha-synuclein.

Exploration of the activation pathway of Deltaalpha-Chymotrypsin with molecular dynamics simulations and correlation with kinetic experiments.

Correlating the experimentally observed kinetics of protein conformational changes with theoretical predictions is a formidable and challenging task, due to the multitude of degrees of freedom (>5,000) in a protein and the huge gap between the timescale of the kinetic event of interest (ms) and the typical timescale of computer simulations (ns). In this study we show that using the targeted molecular dynamics (TMD) method it is possible to simulate conformational changes of the ms time range and to correlate multiple simulations of single pathways with ensemble experiments on both the structural and energetic basis. As a model system we chose to study the conformational change of rat-Deltaalpha-chymotrypsin from its inactive to its active conformation. This activation process has been analyzed previously by experimental and theoretical methods, i.e. fluorescence stopped-flow spectroscopy (FSF), molecular dynamics (MD) and TMD. Inspired by the results of these studies on the wild type (WT) enzyme, several mutants were constructed to alter the conformational pathway and studied by FSF measurements. In the present work WT and mutant N18G were subjected to multiple MD and subsequent TMD simulations. We report the existence of two main activation pathways, a feature of chymotrypsin activation that has been hitherto unknown. A method to correlate the energetics of the different pathways calculated by TMD and the kinetic parameters observed by experimental methods such as FSF is presented. Our work is relevant for experimental single molecule studies of enzymes in general.