(57)【要約】【課題】広帯域パルスレーザーを用いた２光子励起蛍光観察において同時励起、選択的励起の切り替えを行う。【解決手段】広帯域パルスの位相変調関数として、光源が有する全波長帯域において励起エネルギーが最も大きくなる第１の位相関数φc(ω)と小さくなる第２の位相関数φd(ω)とを組み合わせた関数を用いる。励起したい蛍光分子の２光子吸収断面積が大きな波長帯域１には第１の位相関数を、蛍光を抑制したい分子の２光子吸収断面積が大きな波長帯域２には第２の位相関数を用いることにより選択的励起を行う。さらに、第１の位相関数に０から１の割合で第２の位相関数を加算した位相関数を与え、その割合を波長帯域１と２で各々調整することにより、２種類の蛍光分子から発生する２光子励起蛍光強度を独立に制御する。このとき、波長帯域１，２のいずれかの位相関数に群遅延を与える位相を加えることによって、２種類の蛍光分子を独立に扱うことができる。【選択図】図７
2010014677.tif 000002
【特許請求の範囲】
【請求項１】
広帯域パルスレーザー光源と、
前記パルスレーザー光源から発生されたレーザー光のスペクトル位相を変調する位相変調部と、
前記位相変調部から出射したレーザー光を試料に照射し、試料から発生された２光子励起蛍光強度を測定する蛍光測定部と、
前記位相変調部を制御する制御部とを備え、
前記制御部は、前記パルスレーザー光源の全波長帯域において第二高調波の光強度が最も高くなる第１の位相関数と、前記パルスレーザー光源の全波長帯域において第二高調波の光強度が最も小さくなる第２の位相関数とを保持し、前記第１の位相関数と第２の位相関数を組み合わせて得られる合成位相関数によって前記位相変調部を制御することを特徴とする２光子励起蛍光観察装置。
【請求項２】
請求項１記載の２光子励起蛍光観察装置において、前記試料は少なくとも２種類の蛍光分子を含み、前記制御部は、第１の蛍光分子の２光子吸収断面積が大きな波長帯域に前記第１の位相関数を適用し、他の蛍光分子の２光子吸収断面積が大きな波長帯域に前記第２の位相関数を適用するようにして前記位相変調部を制御し、前記第１の蛍光分子を選択的に２光子励起することを特徴とする２光子励起蛍光観察装置。
【請求項３】
請求項１記載の２光子励起蛍光観察装置において、前記制御部は、所定の境界波長を境に波長を複数の波長帯域に分割し、その波長帯域の一つに前記第１の位相関数を適用し、残りの波長帯域に前記第２の位相関数を適用して生成された合成位相関数によって前記位相変調部を制御し、前記第１の位相関数が適用された波長帯域に大きな２光子吸収断面積を有する蛍光分子を選択的に励起することを特徴とする２光子励起蛍光観察装置。
【請求項４】
請求項１記載の２光子励起蛍光観察装置において、前記試料は少なくとも第１の蛍光分子と第２の蛍光分子を含み、前記制御部は、前記第１の蛍光分子の２光子吸収断面積が大きな第１の波長帯域に対して、前記第１の位相関数に前記第２の位相関数を第１の割合で加算した位相関数を適用し、前記第２の蛍光分子の２光子吸収断面積が大きな第２の波長帯域に対して、前記第１の位相関数に前記第２の位相関数を第２の割合で加算した位相関数を適用し、かつ前記第１の波長帯域に適用する位相関数と第２の波長帯域に適用する位相関数の一方に群遅延時間をもたらす位相を加えて得られた合成位相関数によって前記位相変調部を制御し、前記第１の蛍光分子と前記第２の蛍光分子を同時に２光子励起することを特徴とする２光子励起蛍光観察装置。
【請求項５】
請求項１記載の２光子励起蛍光観察装置において、前記制御部は、所定の境界波長を境に波長を少なくとも第１の波長帯域と第２の波長帯域に分割し、前記第１の波長帯域に対して、前記第１の位相関数に前記第２の位相関数を第１の割合で加算した位相関数を適用し、前記第２の波長帯域に対して、前記第１の位相関数に前記第２の位相関数を第２の割合で加算した位相関数を適用し、かつ前記第１の波長帯域に適用する位相関数と第２の波長帯域に適用する位相関数の一方に群遅延時間をもたらす位相を加えて得られた合成位相関数によって前記位相変調部を制御し、前記第１の波長帯域に大きな２光子吸収断面積を有する第１の蛍光分子と前記第２の波長帯域に大きな２光子吸収断面積を有する第２の蛍光分子を同時に２光子励起することを特徴とする２光子励起蛍光観察装置。
【請求項６】
請求項４又は５記載の２光子励起蛍光観察装置において、前記制御部は、前記第１の割合及び／又は第２の割合を調整することにより前記第１の蛍光分子と第２の蛍光分子から発生する２光子励起蛍光強度を独立に制御することを特徴とする２光子励起蛍光観察装置。
【請求項７】
請求項１〜６のいずれか１項記載の２光子励起蛍光観察装置において、前記位相変調部は、前記パルスレーザー光源から発生されたレーザー光を分光する第１の分光素子と、前記第１の分光素子によって分光されたスペクトル光を平行光にする第１の集光素子と、前記平行光にされたスペクトル光のスペクトル位相を変調する空間光変調器と、前記空間光変調器を通った平行光を集光する第２の集光素子と、前記第２の集光素子の集光点に配置され前記集光された光を光ビームに整形する第２の分光素子とを備えることを特徴とする２光子励起蛍光観察装置。
【請求項８】
請求項７記載の２光子励起蛍光観察装置において、前記空間光変調器は、スペクトル分散方向に並べて配置された複数の位相変調素子を備えることを特徴とする２光子励起蛍光観察装置。
【請求項９】
請求項１記載の２光子励起蛍光観察装置において、前記蛍光測定部に試料に代えて非線形光学結晶を配置し、当該結晶から発生する第二高調波の信号が最も大きくなるように前記制御部によって前記位相変調部を最適制御して前記第１の位相関数を取得し、前記第二高調波の信号が最も小さくなるように前記制御部によって前記位相変調部を最適制御して前記第２の位相関数を取得することを特徴とする２光子励起蛍光観察装置。
【請求項１０】
請求項１記載の２光子励起蛍光観察装置において、前記蛍光測定部の試料位置における光パルスのスペクトル位相を測定し、それを打ち消すような位相を前記第１の位相関数として取得し、前記レーザー光のパルス幅を長くするような位相関数を前記第２の位相関数として取得することを特徴とする２光子励起蛍光観察装置。
【請求項１１】
２光子励起スペクトルが異なる第１の蛍光分子と第２の蛍光分子を含む試料中の前記第１の蛍光分子を選択的に２光子励起して蛍光観察する２光子励起蛍光観察方法において、
広帯域パルスレーザー光源の全波長帯域において第二高調波の光強度が最も高くなる第１の位相関数と、前記第二高調波の光強度が最も小さくなる第２の位相関数とを取得する工程と、
前記広帯域パルスレーザー光源から発生されたレーザー光のスペクトル位相を、前記第１の蛍光分子の２光子吸収断面積が大きな波長帯域には前記第１の位相関数を適用し、前記第２の蛍光分子の２光子吸収断面積が大きな波長帯域には前記第２の位相関数を適用して得られた位相関数で変調する工程と、
前記スペクトル位相が変調されたレーザー光を試料に照射する工程と、
前記試料から発生される蛍光を検出する工程と、
を有することを特徴とする２光子励起蛍光観察方法。
【請求項１２】
２光子励起スペクトルが異なる第１の蛍光分子と第２の蛍光分子を含む試料中の前記第１の蛍光分子と前記第２の蛍光分子を同時に２光子励起して蛍光観察する２光子励起蛍光観察方法において、
広帯域パルスレーザー光源の全波長帯域において第二高調波の光強度が最も高くなる第１の位相関数と、前記第二高調波の光強度が最も小さくなる第２の位相関数とを取得する工程と、
前記広帯域パルスレーザー光源から発生されたレーザー光のスペクトル位相を、前記第１の蛍光分子の２光子吸収断面積が大きな波長帯域には前記第１の位相関数に前記第２の位相関数を第１の割合で加算した位相関数を適用し、前記第２の蛍光分子の２光子吸収断面積が大きな波長帯域には前記第１の位相関数に前記第２の位相関数を第２の割合で加算した位相関数を適用し、かつ前記第１の波長帯域に適用する位相関数と第２の波長帯域に適用する位相関数の一方に群遅延時間をもたらす位相を加えて得られた位相関数で変調する工程と、
前記スペクトル位相が変調されたレーザー光を試料に照射する工程と、
前記試料から発生される蛍光を検出する工程と、
を有することを特徴とする２光子励起蛍光観察方法。
【請求項１３】
請求項１２記載の２光子励起蛍光観察方法において、前記第１の割合及び／又は第２の割合を調整することにより前記第１の蛍光分子と第２の蛍光分子から発生する２光子励起蛍光強度を独立に制御することを特徴とする２光子励起蛍光観察方法。
【請求項１４】
請求項１１〜１３のいずれか１項記載の２光子励起蛍光観察方法において、蛍光測定部に試料に代えて非線形光学結晶を配置し、当該結晶から発生する第二高調波の信号が最も大きくなるように前記レーザー光のスペクトル位相を最適制御して前記第１の位相関数を取得し、前記第二高調波の信号が最も小さくなるように前記レーザー光のスペクトル位相を最適制御して前記第２の位相関数を取得することを特徴とする２光子励起蛍光観察方法。
【請求項１５】
請求項１１〜１３のいずれか１項記載の２光子励起蛍光観察方法において、試料位置における前記レーザー光のスペクトル位相を測定し、それを打ち消すような位相を前記第１の位相関数として取得し、前記レーザー光のパルス幅を長くするような位相関数を前記第２の位相関数として取得することを特徴とする２光子励起蛍光観察方法。
【請求項１６】
請求項１１〜１３のいずれか１項記載の２光子励起蛍光観察方法において、前記第１の位相関数及び第２の位相関数は、記憶媒体から読み出すことによって取得することを特徴とする２光子励起蛍光観察方法。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、２光子励起蛍光観察方法及び装置に関する。
【背景技術】
【０００２】
フェムト秒レーザーを励起光源とした２光子励起蛍光顕微鏡は生細胞内の生命現象可視化に欠かせない技術である。生命現象の解明において重要となっているタンパク質分子間相互作用の可視化においては、マルチカラーイメージング、蛍光相互相関分光法と蛍光共鳴エネルギー移動（FRET）イメージングが必要不可欠である。マルチカラーイメージングは、スペクトル情報を介して回折限界領域内に同時に存在するタンパク質についての情報を与える。蛍光相互相関分光法では、タンパク質の相互作用を定量的に測定可能である。また、FRETイメージングは、タンパク質−タンパク質相互作用及びタンパク質のconformal change（構造変化）についての空間的及び時間的情報を与える。マルチカラーイメージングや蛍光相互相関分光法では複数種類の蛍光分子を同時に２光子励起する必要があり、FRETイメージングでは一種類の蛍光分子のみを選択的に２光子励起する必要がある。
【０００３】
２光子励起蛍光顕微鏡において、特定の蛍光分子を選択的に２光子励起するための広帯域パルスの位相変調法がいくつか報告されている。非特許文献１−３には、正弦波関数の位相変調を用いた選択的励起法が報告されている。また、非特許文献４，５にはバイナリの位相変調を用いた選択的励起法が報告されている。
【０００４】
【特許文献１】米国特許第7105811号明細書
【非特許文献１】D. Meshulach, Y. Silberberg, “Coherent quantum control of two-photon transitions by a femtosecond laser pulse,” Nature 396, 239-242 (1998)
【非特許文献２】K. A. Walowicz, I. Pastirk, V. V. Lozovoy, and M. Dantus, “Multiphoton intrapulse interference. 1. Control of multiphoton processes in condensed phases,” J. Phys. Chem. A 106, 9369-9373 (2002).
【非特許文献３】V. V. Lozovoy, I. Pastirk, K. A. Walowicz, and M. Dantus, “Multiphoton intrapulse interference. II. Control of two- and three-photon laser induced fluorescence with shaped pulses,” J. Chem. Phys. 118, 3187-3196 (2003).
【非特許文献４】M. Comstock, V. V. Lozovoy, I. Pastirk, and M. Dantus, “Multiphoton intrapulse interference 6; binary phase shaping,” Opt. Express 12, 1061-1066 (2004).
【非特許文献５】V. V. Lozovoy, B. Xu, J. C. Shane, and M. Dantus “Selective nonlinear optical excitation with pulses shaped by pseudorandom Galois fields,” Phys. Rev. A 74, 041805(R) (2006).
【発明の開示】
【発明が解決しようとする課題】
【０００５】
マルチカラーイメージングとFRETイメージングには、複数の蛍光分子を選択的あるいは同時に２光子励起する必要がある。選択的励起あるいは同時励起は、それぞれ狭帯域パルスの中心波長を、２種類の蛍光分子における２光子励起スペクトルの比が最大になる波長あるいは等しくなる波長に合わせることによって達成される。複数の蛍光分子によって標識されたサンプルのマルチカラーイメージングにおいては、すべての蛍光分子を同時に励起、あるいはレーザー波長をチューニングして各蛍光分子を順番に励起すればよい。しかし、細胞を複数の蛍光分子で標識する場合、しばしば濃度に大きな差が生じ、ある蛍光分子と他の蛍光分子の２光子励起蛍光強度に大きな差が生じることが問題になる。そのような場合、複数の蛍光分子の２光子励起蛍光強度を個別にかつ任意に制御できると好都合である。また、現状のモードロック・フェムト秒レーザーは波長合わせに数秒以上の時間がかかるため、迅速な制御法が必要になる。
【０００６】
従来の正弦波関数を用いた位相変調の場合、特定の波長λcに励起エネルギーを集中することが可能であるが、λc以外の波長においても励起エネルギーが十分残っているため、選択的励起のコントラストが低い。バイナリ位相変調では、励起エネルギーを集中した波長以外の波長において励起エネルギーが十分低くなるように制御され、正弦波位相変調で問題となったコントラストの低下が抑制される。しかし、０とπの２つの位相のみを用いているため、複数種類の蛍光分子の蛍光強度を独立かつ任意に制御することは困難である。
【０００７】
本発明は、広帯域パルスレーザーを用いた２光子励起蛍光観察において同時励起、選択的励起の切り替えが可能なレーザーパルスの位相変調法を提案するものである。
【課題を解決するための手段】
【０００８】
本発明では、広帯域パルスの位相変調関数として、光源が有する全波長帯域において励起エネルギーが最も大きくなる第１の位相関数と小さくなる第２の位相関数とを組み合わせた関数を用いる。励起したい蛍光分子の２光子吸収断面積が大きな波長帯域１には第１の位相関数を、蛍光を抑制したい分子の２光子吸収断面積が大きな波長帯域２には第２の位相関数を用いることにより選択的励起を行う。この方法により、正弦波位相変調の場合のようなコントラストの低下は抑制される。また、波長帯域１か２の一方に群遅延時間をもたらす位相を与えることにより、波長帯域１と２で構成されるパルスを独立に扱うことを可能とする。さらに、第１の位相関数に０から１の割合で第２の位相関数を加算した位相関数を与え、その割合を波長帯域１と２で各々調整することにより、２種類の蛍光分子から発生する２光子励起蛍光強度を独立かつ任意に制御することが可能となる。また、分割する波長帯域の数を増やすことにより、３種類以上の複数種類の蛍光分子の２光子励起蛍光強度を独立かつ任意に制御することが可能となる。
【０００９】
本発明による２光子励起蛍光観察装置は、広帯域パルスレーザー光源と、パルスレーザー光源から発生されたレーザー光のスペクトル位相を変調する位相変調部と、位相変調部から出射したレーザー光を試料に照射し、試料から発生された２光子励起蛍光強度を測定する蛍光測定部と、位相変調部を制御する制御部とを備え、制御部は、パルスレーザー光源の全波長帯域において第二高調波の光強度が最も高くなる第１の位相関数と、パルスレーザー光源の全波長帯域において第二高調波の光強度が最も小さくなる第２の位相関数とを保持し、第１の位相関数と第２の位相関数を組み合わせて得られる合成位相関数によって位相変調部を制御する。
【００１０】
試料が少なくとも２種類の蛍光分子を含むとき、制御部は、第１の蛍光分子の２光子吸収断面積が大きな波長帯域に第１の位相関数を適用し、他の蛍光分子の２光子吸収断面積が大きな波長帯域に第２の位相関数を適用するようにして位相変調部を制御することで、第１の蛍光分子を選択的に２光子励起することができる。
【００１１】
換言すると、制御部は、所定の境界波長を境に波長を複数の波長帯域に分割し、その波長帯域の一つに第１の位相関数を適用し、残りの波長帯域に第２の位相関数を適用して生成された合成位相関数によって位相変調部を制御することで、第１の位相関数が適用された波長帯域に大きな２光子吸収断面積を有する蛍光分子を選択的に励起することができる。
【００１２】
また、試料が少なくとも第１の蛍光分子と第２の蛍光分子を含むとき、制御部は、第１の蛍光分子の２光子吸収断面積が大きな第１の波長帯域に対して、第１の位相関数に第２の位相関数を第１の割合で加算した位相関数を適用し、第２の蛍光分子の２光子吸収断面積が大きな第２の波長帯域に対して、第１の位相関数に第２の位相関数を第２の割合で加算した位相関数を適用し、かつ第１の波長帯域に適用する位相関数と第２の波長帯域に適用する位相関数の一方に群遅延時間をもたらす位相を加えて得られた合成位相関数によって位相変調部を制御することで、第１の蛍光分子と第２の蛍光分子を同時に２光子励起することができる。
【００１３】
換言すると、制御部は、所定の境界波長を境に波長を少なくとも第１の波長帯域と第２の波長帯域に分割し、第１の波長帯域に対して、第１の位相関数に第２の位相関数を第１の割合で加算した位相関数を適用し、第２の波長帯域に対して、第１の位相関数に第２の位相関数を第２の割合で加算した位相関数を適用し、かつ第１の波長帯域に適用する位相関数と第２の波長帯域に適用する位相関数の一方に群遅延時間をもたらす位相を加えて得られた合成位相関数によって位相変調部を制御することで、第１の波長帯域に大きな２光子吸収断面積を有する第１の蛍光分子と第２の波長帯域に大きな２光子吸収断面積を有する第２の蛍光分子を同時に２光子励起することができる。
【００１４】
このとき、第１の割合及び／又は第２の割合を調整することにより第１の蛍光分子と第２の蛍光分子から発生する２光子励起蛍光強度を独立に制御することができる。
【００１５】
第１の蛍光分子と第２の蛍光分子を含む試料中の第１の蛍光分子を選択的に２光子励起して蛍光観察する本発明の２光子励起蛍光観察方法は、広帯域パルスレーザー光源の全波長帯域において第二高調波の光強度が最も高くなる第１の位相関数と、第二高調波の光強度が最も小さくなる第２の位相関数とを取得する工程と、広帯域パルスレーザー光源から発生されたレーザー光のスペクトル位相を、第１の蛍光分子の２光子吸収断面積が大きな波長帯域には第１の位相関数を適用し、第２の蛍光分子の２光子吸収断面積が大きな波長帯域には第２の位相関数を適用して得られた位相関数で変調する工程と、スペクトル位相が変調されたレーザー光を試料に照射する工程と、試料から発生される蛍光を検出する工程とを有する。
【００１６】
また、第１の蛍光分子と第２の蛍光分子を同時に２光子励起して蛍光観察する本発明の２光子励起蛍光観察方法は、広帯域パルスレーザー光源の全波長帯域において第二高調波の光強度が最も高くなる第１の位相関数と、第二高調波の光強度が最も小さくなる第２の位相関数とを取得する工程と、広帯域パルスレーザー光源から発生されたレーザー光のスペクトル位相を、第１の蛍光分子の２光子吸収断面積が大きな波長帯域には第１の位相関数に第２の位相関数を第１の割合で加算した位相関数を適用し、第２の蛍光分子の２光子吸収断面積が大きな波長帯域には第１の位相関数に第２の位相関数を第２の割合で加算した位相関数を適用し、かつ第１の波長帯域に適用する位相関数と第２の波長帯域に適用する位相関数の一方に群遅延時間をもたらす位相を加えて得られた位相関数で変調する工程と、スペクトル位相が変調されたレーザー光を試料に照射する工程と、試料から発生される蛍光を検出する工程とを有する。
【００１７】
このとき、第１の波長帯域に適用する位相関数と第２の波長帯域に適用する位相関数の一方に群遅延時間をもたらす位相を加える。また、第１の割合及び／又は第２の割合を調整することにより第１の蛍光分子と第２の蛍光分子から発生する２光子励起蛍光強度を独立に制御することができる。
【発明の効果】
【００１８】
本発明によると、選択的励起だけでなく、各々の２光子励起蛍光強度も独立に制御した同時励起が可能となる。これまでマルチカラーイメージングでは、試料作成時において生じる複数種類の蛍光分子濃度差による発光強度の差が問題となっていたが、この問題も本発明により解決できる。
【発明を実施するための最良の形態】
【００１９】
以下、図面を参照して本発明の実施の形態を説明する。
【００２０】
狭帯域パルスに比較して、広帯域パルスは、選択的励起と複数の蛍光分子の２光子励起蛍光強度を独立に制御した同時励起とを迅速かつ簡易に切り替えできる可能性を有する。中間共鳴準位がないと仮定すると、２光子励起蛍光強度ＩFは、蛍光分子の２光子励起スペクトル（または２光子吸収断面積）ｇ(2)(ω)と励起光パルスの第二高調波スペクトル|Ｅ(2)(ω)|2の積で表される。
【００２１】
【数１】
2010014677.tif 000003
【００２２】
ここで、ωは角周波数、Ｅ(ω)及びφ(ω)は、それぞれ励起パルスのスペクトル振幅及びスペクトル位相である。第二高調波のパワースペクトル|Ｅ(2)(ω)|2は、基本波スペクトル|Ｅ(ω)|2やスペクトル位相φ(ω)を変調することにより整形することが可能である。
【００２３】
本発明による蛍光分子の選択的励起、蛍光強度を制御した同時励起のための位相変調関数の作成方法について説明する。本発明では、全ての周波数において第二高調波の光強度が最も高くなる位相関数φc(ω)と、全ての周波数で第二高調波の強度が最も低くなる位相関数φd(ω)を用い、２種類の位相関数φc(ω)とφd(ω)を組み合わせて、スペクトル位相を変調する位相関数を作成する。第二高調波の光強度が最も高くなる位相関数φc(ω)は、フーリエ変換限界パルスのスペクトル位相、すなわち全ての周波数においてφ(ω)＝０である。このとき、パルス幅は最も短く、第二高調波の強度が最も高い。
【００２４】
２種類の蛍光分子１，２のうち蛍光分子１を選択的に励起する場合、次式のように、励起したい蛍光分子１の２光子吸収断面積が大きな周波数帯域Ｒ1には位相関数φc(ω)を用い、蛍光を抑制したい蛍光分子２の２光子吸収断面積が大きな周波数帯域Ｒ2には位相関数φd(ω)を用いる。位相関数φc(ω)が適用された周波数帯域では２光子励起が促進されて強い２光子励起蛍光が発生し、逆に位相関数φd(ω)が適用された周波数帯域では２光子励起が抑制されて２光子励起蛍光が発生しにくくなる。本発明によれば、このようにスペクトル領域を分割する最適周波数を決定するだけで、高コントラストな選択的励起を行うことができる。
【００２５】
【数２】
2010014677.tif 000004
【００２６】
このとき、周波数帯域Ｒ1かＲ2の一方に用いる位相関数に、群遅延時間をもたらす位相φ’(ω−ω0)を与えることにより、周波数帯域Ｒ1とＲ2で構成されるパルスを独立に扱うことを可能とする。すなわち、群遅延を与えることにより、励起パルスを構成している周波数帯域Ｒ1の周波数と周波数帯域Ｒ2の周波数は、蛍光分子と個別に相互作用する。さらに、位相関数φc(ω)に０から１の割合で位相関数φd(ω)を加算した位相関数を与え、その割合αとβ（０≦α，β≦１）を周波数帯域Ｒ1とＲ2で各々調整することにより、２種類の蛍光分子から発生する２光子励起蛍光強度を独立かつ任意に制御することができる。
【００２７】
【数３】
2010014677.tif 000005
【００２８】
このように、群遅延及び２種類の位相関数の加算割合のパラメータにより、２種類の蛍光分子の２光子励起蛍光強度を独立に制御することができる。スペクトル領域の分割数を増し、群遅延を導入することにより、この方法は３種類以上の蛍光分子に対しても容易に拡張することができる。
【００２９】
例えば、３種類の蛍光分子１，２，３が存在するとき、そのうちの蛍光分子１のみを選択的に励起する場合には、励起したい蛍光分子１の２光子吸収断面積が大きな周波数帯域Ｒ1には位相関数φc(ω)を用い、蛍光を抑制したい蛍光分子２と蛍光分子３の２光子吸収断面積が大きな周波数帯域Ｒ2，Ｒ3には位相関数φd(ω)を用いることにより、蛍光分子１だけを選択的に励起することができる。
【００３０】
【数４】
2010014677.tif 000006
【００３１】
また、３種類の蛍光分子の蛍光強度を任意に調節して同時励起するには、次のようにすればよい。すなわち、次式のように、蛍光分子１，２，３の２光子吸収断面積が大きな周波数帯域Ｒ1，Ｒ2，Ｒ3に対して、それぞれ位相関数φc(ω)に０から１の割合で位相関数φd(ω)を加算した位相関数を与え、その割合α，β，γ（０≦α，β，γ≦１）を周波数帯域Ｒ1，Ｒ2，Ｒ3で各々調整すると共に、群遅延時間をもたらす位相φ1’(ω−ω0)，φ2’(ω−ω0)を２つの周波数帯域の位相関数に付加し、遅延時間を異なる値に設定する。
【００３２】
【数５】
2010014677.tif 000007
【００３３】
次に、最も第二高調波光強度が高くなる位相関数φc(ω)と最も低くなる位相関数φd(ω)の取得方法について説明する。位相関数φc(ω)及びφd(ω)は、以下の２つの方法のいずれかによって決定することができる。
【００３４】
（１）試料の代わりに非線形光学結晶を用い、最も第二高調波光強度が高く、または低くなるようにスペクトル位相の最適制御（焼きなまし法、遺伝的アルゴリズム）を行い、それぞれ位相関数φc(ω)とφd(ω)を取得する。
（２）最も第二高調波光強度が高くなるのは、全ての周波数成分で位相がそろったフーリエ変換限界パルスである。パルス測定技術（FROG：Frequency Resolved Optical Gating, R.Trebino et al., J Opt. Soc. Am 10, 1101 (1993) SPIDER：Spectral Interferometry for Direct Electric Field Reconstruction, C. Iaconis et al., Opt. Lett. 23, 79 (1998)など）により、試料位置における光パルスのスペクトル位相φp(ω)を測定し、それを打ち消すような位相（−φp(ω)）を与えてフーリエ変換限界パルスを発生させる（最も第二高調波光強度が高い）。このとき測定光路中に配置した位相変調器により与えた位相関数がφc(ω)=−φp(ω)である。また、第二高調波光強度を低くするには、パルス幅を長くするような位相関数を与えればよく、φd(ω)=Ｄ2(ω-ω0)2/2 がひとつの位相関数である。Ｄ2，ω0は群遅延分散、特定周波数（何でもよい）である。Ｄ2が大きいほどパルス幅が長く、第二高調波強度が小さくなる。位相変調が可能な最大のＤ2は光路中に挿入される位相変調素子のピクセル数で制限され、ピクセル数が多いほどＤ2は大きくできる。ピクセル数が少ない場合には、Ｄ2は大きくできないため、この設計した関数より、（１）の方法によって最適制御により位相関数φd(ω)を取得した法が良い場合がある。
【００３５】
次に、実施例により、本発明を詳細に説明する。本実施例では、試料として、蛍光タンパク質Azuriteで核を、EGFPで細胞質を標識したHeLa細胞を用いた。図４は、AzuriteとEGFPの２光子励起スペクトルを示す図である。
【００３６】
図１は、本発明による２光子励起蛍光観察装置の構成例を示す模式図である。広帯域パルスレーザー光源１１には、繰り返し周波数８０ＭＨｚ、パルス幅８フェムト秒以下でパルス発光するチタンサファイアレーザーを用いた。図２は、チタンサファイアレーザーの発光スペクトルを示し、７００ｎｍから１０４０ｎｍの範囲に広がっている。なお、光源は６６０ｎｍから１１００ｎｍであるが、スペクル位相変調におけるスペクトル分解能を向上させるために、空間位相変調器でスペクトル帯域幅を制限している。空間位相変調器のピクセル数と分光素子と集光素子の焦点距離の関係からスペクル位相変調におけるスペクトル分解能が決まる。レーザー光源１１から射出された光パルスはプリズム対１２、チャープミラー対１３を通ってパルス整形器１４に入射され、パルス整形器１４内の空間光変調器２７によりスペクトルが位相変調される。位相変調された光パルスは、顕微鏡１５へ導入される。プリズム対１２とチャープミラー対１３は、光学部品の分散を補償するためのものである。
【００３７】
パルス整形器１４は、プリズムや回折格子などの分光素子、集光レンズや凹面鏡などの集光素子、及び空間光変調器２７を備える。図１の例では、分光素子としての一対の回折格子２１，２２、集光素子として対向配置された一対の凹面鏡２５，２６を備える。第１の回折格子２１は第１の凹面鏡２５の焦点に設置され、第２の回折格子２２は第２の凹面鏡２６の焦点に設置されている。パルス整形器１４に入射したレーザー光は第１の回折格子２１に入射して分光される。分光されたレーザー光は、ミラー２３によって反射されて第１の凹面鏡２５に入射し、そこで反射されて平行光となり、空間光変調器２７に入射する。空間光変調器２７を通って波長毎に位相変調された平行レーザー光は、第２の凹面鏡２６及びミラー２４で反射されて、第２の回折格子２２上に集光される。第２の回折格子２２で反射されたレーザー光は、パルス整形器１４への入射時と逆の過程を経て１本の細いビームとなり、後段の顕微鏡１５に入射する。
【００３８】
図１の顕微鏡は、蛍光分子が２種類の場合の構成例を示している。顕微鏡１５に入射したレーザー光は、対物レンズ３１により試料容器３２中の試料に集光され、２光子励起過程で試料中の蛍光分子から蛍光を発生させる。試料容器３２を通ったレーザー光はショートパスフィルター３３によって励起光が除去されたのち、ダイクロイックビームスプリッター３４によって２つの光路に分割される。一方の光路のレーザー光は、第１のバンドパスフィルター３５を通った後、光電子増倍管３７で検出される。他方の光路のレーザー光は、第２のバンドパスフィルター３６を通り、光電子増倍管３８で検出される。試料容器３２はステージ上に保持されており、ステージはステージ駆動部３９により入射光に対して垂直方向に２次元移動することができる。第１のバンドパスフィルター３５としては、Azuriteの２光子励起蛍光を透過するように中心波長４４７ｎｍ、半値全幅６０ｎｍのものを用いた。また第２のバンドパスフィルター３８は、EGFPの２光子励起蛍光を透過するように、中心波長５２０ｎｍ、半値全幅３５ｎｍのものを用いた。
【００３９】
図３は、空間光変調器２７の模式図である。空間光変調器２７は、スペクトル分散方向に並べて一次元アレイ状に配列された多数の位相変調素子２８を備え、各位相変調素子２８はそこを通過する光の位相を互いに独立して変調することができる。パルス整形器１４に入射した光パルスは第１の回折格子２１によって分光され、空間光変調器２７の各素子位置はスペクトルを構成する各波長に対応する。すなわち、空間光変調器２７の各位相変調素子２８は、光パルスの対応する波長成分の位相を相互に独立して変調することができる。空間光変調器２７の素子数は多ければ多いほど精密な制御を行うことができる。本実施例では、素子数を１２８とし、１つの位相変調素子２８によって変調される波長範囲を約３ｎｍとした。本実施例では位相変調素子として液晶を利用した素子を用いたが、液晶以外に、変形ミラーを用いて位相を変調する素子、音響光学効果を利用して位相変調を行う素子、電気光学効果を用いて位相を変調する素子など、種々のタイプのものを用いることができる。
【００４０】
顕微鏡１５の２つの光電子増倍管３７，３８の出力信号は制御部１６の画像処理ユニット４３に入力される。また、パルス整形器１４の空間光変調器２７は、制御部１６の駆動制御ユニット４１からの制御信号によって駆動される。さらに、ステージ駆動部３９は、試料中の観察対象物が光ビームによって走査されるように、制御部１６の駆動制御ユニット４１からの制御信号によってステージを２次元駆動する。画像処理ユニット４３によって生成された画像は表示部１８に表示される。また、位相関数φc(ω)とφd(ω)は制御部１６のメモリ４２に記憶されており、駆動制御ユニット４１は入力部１７から入力された条件に従って２つの位相関数φc(ω)，φd(ω)から合成した位相関数を生成し、その合成位相関数によって空間光変調器２７を制御する。
【００４１】
最初に、光源が有する全周波数帯域において励起エネルギーが最も大きくなる位相関数φc(ω)と最も小さくなる位相関数φd(ω)を取得した。ここでは、上記（１）の方法によって２つの位相関数を取得した。具体的には、試料の代わりに非線形光学結晶である厚さ１０μｍのＢＢＯ（β−barium borate）結晶を用い、結晶から発生する第二高調波（ＳＨ）の信号が最も大きくなるように、または最も小さくなるように空間変調器２７の各位相変調素子２８に対して適応制御を行い、位相関数φc(ω)とφd(ω)を取得した。ここで位相関数φc(ω)の取得に当たり、試料内部のレーザー集光点においてφ(ω)＝０を達成するには、光学素子の分散を補償する必要がある。分散補償素子には、プリズム対、回折格子対、位相変調素子を用いたスペクトル位相変調器がある。広帯域パルスを用いた場合には、高次の分散が大きく影響するため、完全に分散を補償するためにはスペクトル位相変調器による調整が必要となる。すなわちスペクトル位相変調器により、φ(ω)＝０となる変調関数を与える必要がある。
【００４２】
図５に、取得した位相関数φc(ω)とφd(ω)を示す。また、このとき発生したＳＨスペクトルを、図６のグラフに示す。図６のスペクトル６１はスペクトル位相を位相関数φc(ω)で変調したときのＳＨスペクトル、スペクトル６２はスペクトル位相を位相関数φd(ω)で変調したときのＳＨスペクトルである。ただし、位相関数φc(ω)はフーリエ限界パルスとなることから、図５では、位相関数φc(ω)のグラフは位相を全て０とし、位相関数φd(ω)のグラフは得られた位相関数φc(ω)と位相関数φd(ω)の差として表している。図６のスペクトル６１に示されるように、空間光変調器２７に位相関数φc(ω)を設定すると、全波長帯域において強め合う干渉によりＳＨ強度が最も高くなる。また、空間光変調器２７に位相関数φd(ω)を設定すると、スペクトル６２に示されるように、全波長帯域において弱め合う干渉によりＳＨ強度が最も低くなっている。 取得した位相関数φc(ω)とφd(ω)は制御部１６のメモリ４２に格納する。
【００４３】
次に、選択的励起について説明する。図５に示した位相関数φc(ω)と位相関数φd(ω)を組み合わせて、短波長側あるいは長波長側のみに励起エネルギーを集中するような位相変調を行った。２光子励起スペクトルでEGFPとAzuriteは、それぞれ基本波の波長９２０ｎｍと７６０ｎｍにピークを有する（図４参照）。これらの波長はレーザースペクトル中で長波長と短波長に対応するため、スペクトル領域を波長λdより長波長と波長λdより短
波長の２つの波長帯域に分割した。波長λdとして、EGFPとAzuriteの２光子励起スペクトルにおいて、両者のピークを与える波長の間にあり、強度が同程度となるような波長を選択した。図４を見るとその波長は約４１０ｎｍであり、基本波の波長で考えるとその２倍の８２０ｎｍとなる。分割波長を波長８２０ｎｍとしたとき、EGFP／Azurite及びAzurite／EGFPに対して２光子励起蛍光強度のコントラスト比は、それぞれ２１：１及び１：５．６となった。このとき、EGFP及びAzuriteの２光子励起蛍光強度は、フーリエ変換限界パルスに対する２光子励起蛍光強度の０．４２倍及び０．２９倍であり、バイオイメージングにとって十分な強度であった。このEGFP／Azurite及びAzurite／EGFPに対するコントラスト比の積は１００を超えている。
【００４４】
図６に示したグラフ６３は、図７（ａ）に示すように、８２０ｎｍ以上の長波長帯域では位相関数φd(ω)を、８２０ｎｍ以下の短波長帯域では位相関数φc(ω)を用いた場合のＳＨスペクトルを示している。ここで、図７の横軸のピクセル番号は空間光変調器２７のピクセル番号を表し、番号が小さいピクセルほど短波長に、番号が大きいピクセルほど長波長に対応する。また、図６に示したグラフ６４は、図７（ｂ）に示すように、８２０ｎｍ以上の長波長帯域では位相関数φc(ω)を、８２０ｎｍ以下の短波長帯域では位相関数φd(ω)を用いた場合のＳＨスペクトルを示している。それぞれ励起エネルギーが特定の波長帯域に集中し、他の波長帯域では励起エネルギーが十分小さいことがわかる。これらのＳＨスペクトルとなる位相変調は、AzuriteとEGFPの選択的励起に応用可能である。
【００４５】
このように、２つの波長帯域で２種類の位相関数を切り替えて適用するに当たっては、入力部１７から２つの波長帯域の境界の波長値及びどの波長帯域にどの位相関数を適用するかの情報を入力する。制御部１６の駆動制御ユニット４１は、その入力を受けてメモリ４２から読み出した２つの位相関数φc(ω)，φd(ω)を境界波長でつなぎ合わせた合成位相関数を発生し、その合成位相関数によってパルス整形器１４の空間光変調器２７を制御し、試料に入射する光パルスのスペクトルの各波長に所望の位相変調を与える。
【００４６】
蛍光タンパク質Azuriteで核を、EGFPで細胞質を標識したHeLa細胞の２光子励起蛍光イメージングを行った。集光点近部のみで発生した２光子励起蛍光を光電子増倍管３７，３８により測定した。このとき、ショートパスフィルター３３を用いて励起光を除去するとともに、中心波長４４７ｎｍ、半値全幅６０ｎｍのバンドパスフィルター３５、及び中心波長５２０ｎｍ、半値全幅３５ｎｍのバンドパスフィルター３６を用いて、AzuriteとEGFPを蛍光波長ごとに測定した。ステージに固定された試料をステージ移動して２次元走査することにより、２チャンネルで同時に２光子励起蛍光像を取得した。
【００４７】
図８（ａ）にAzuriteを選択的に励起した場合、図８（ｂ）にEGFPを選択的に励起した場合の２光子励起蛍光像を示す。Azuriteの選択的励起に当たっては、図７（ａ）に示した位相関数を用いた。また、EGFPの選択的励起に当たっては、図７（ｂ）に示した位相関数を用いた。図８（ａ）にはAzuriteに起因する青色の蛍光像が強く見えているものの、EGFPに起因する緑色の蛍光像がほとんど見えず、また、図８（ｂ）には、Azuriteに起因する青色の蛍光像がほとんどみられず、逆にEGFPに起因する緑色の蛍光像が強く見えており、それぞれ、一方の蛍光像のみが選択的に得られていることがわかる。 
【００４８】
次に、蛍光強度の制御について説明する。AzuriteとEGFPの蛍光強度を独立に制御するような位相変調を行った。８２０ｎｍ以下の波長帯域に群遅延時間１００ｆｓに相当する位相変調を行い、さらに位相関数φc(ω)に位相関数φd(ω)を、８２０ｎｍ以下の波長帯域と８２０ｎｍを超える波長帯域で、それぞれ任意の割合で加算した位相関数を用いた。すなわち、上記式(4)において、Ｒ1を８２０ｎｍ以下の波長帯域、Ｒ2を８２０ｎｍを超える波長帯域とし、加算割合αとβ（０≦α，β≦１）を変化させた位相関数を空間光変調器に設定し、蛍光観察した。空間光変調器に設定した位相関数を図９に示す。
【００４９】
図１０〜１２に、加算割合α，βを変化させることにより得られた２光子励起蛍光像を示す。図中の（ｘ，ｙ）は８２０ｎｍ以下の波長帯域に設定した位相関数のαの値と、８２０ｎｍ以上の波長帯域に設定した位相関数のβの値を表している。図１０はAzuriteの蛍光強度のみを制御した場合の２光子励起蛍光像、図１１はEGFPの蛍光強度のみを制御した場合の２光子励起蛍光像、図１２はAzuriteとEGFPの蛍光強度をそれぞれ同時に制御した場合の２光子励起蛍光像である。各図の上段は図１の光電子増倍管３７で検出したAzuriteに対応する２光子励起蛍光像、中段は光電子増倍管３８で検出したEGFPに対応する２光子励起蛍光像、下段は両者を重ね合わせた像である。
【００５０】
φc(ω)に対するφd(ω)の加算割合を２つの波長帯域において変化させることにより、AzuriteとEGFPの蛍光強度を独立に制御できていることがわかる。標識された分子の濃度差により、一方の蛍光強度が強すぎる場合においても本発明により、同程度の蛍光強度でイメージングを行うことが可能となる。このように、本発明の位相変調法を用いれば、単一の光源のみで複数種類の蛍光分子の２光子励起蛍光強度を独立かつ任意に制御可能となり、試料作成における複数種類の蛍光分子の濃度差を制御する必要が軽減される。また、空間光変調器２７に設定する位相関数を入れ替えることにより、選択的励起と同時励起を瞬時に切り替えながら観察を続けることが可能になる。
【００５１】
上記実施例では、２種類の蛍光分子としてAzuriteとEGFPを用いたため、２つの波長帯域の境界波長を８２０ｎｍとした。他の蛍光分子の場合には、その蛍光分子の２光子励起スペクトルを参照して境界波長を決める必要があるのは勿論である。また、本実施例に示した顕微鏡やその他の光学系の構成や蛍光像の取得の仕方は単なる一例であり、本実施例に示したものに限られない。
【図面の簡単な説明】
【００５２】
【図１】本発明による２光子励起蛍光観察装置の構成例を示す模式図。
【図２】チタンサファイアレーザーの発光スペクトルを示す図。
【図３】空間光変調器の模式図。
【図４】Azurite及びEGFPの２光子励起スペクトルを示す図。
【図５】位相関数φc(ω)，φd(ω)を示す図。
【図６】位相関数を変えて得られるＳＨスペクトルの図。
【図７】空間光変調器に設定した位相関数の例を示す図。
【図８】選択的に励起された２光子励起蛍光像の例を示す図。
【図９】空間光変調器に設定した位相関数の例を示す図。
【図１０】Azuriteの蛍光強度のみを制御した場合の２光子励起蛍光像を示す図。
【図１１】EGFPの蛍光強度のみを制御した場合の２光子励起蛍光像を示す図。
【図１２】AzuriteとEGFPの蛍光強度をそれぞれ同時に制御した場合の２光子励起蛍光像を示す図。
【符号の説明】
【００５３】
１１：広帯域パルスレーザー光源
１２：プリズム対
１３：チャープミラー対
１４：パルス整形器
１５：顕微鏡
１６：制御部
１７：入力部
１８：表示部
２１，２２：回折格子
２３，２４：ミラー
２５，２６：凹面鏡
２７：空間光変調器
２８：位相変調素子
３１：対物レンズ
３２：試料容器
３３：ショートパスフィルター
３４：ビームスプリッター
３５，３６：バンドパスフィルター
３７，３８：光電子増倍管
３９：ステージ駆動部
４１：駆動制御ユニット
４２：メモリ
４３：画像処理ユニット
【図１】
2010014677.tif 000008
【図２】
2010014677.tif 000009
【図３】
2010014677.tif 000010
【図４】
2010014677.tif 000011
【図５】
2010014677.tif 000012
【図６】
2010014677.tif 000013
【図７】
2010014677.tif 000014
【図８】
2010014677.tif 000015
【図９】
2010014677.tif 000016
【図１０】
2010014677.tif 000017
【図１１】
2010014677.tif 000018
【図１２】
2010014677.tif 000019