用于产生化学或生化信号的系统和方法

文档序号:440143阅读:871来源:国知局
专利名称:用于产生化学或生化信号的系统和方法
技术领域
本发明涉及用于产生化学或生物信号的装置和效应,该化学或生物信号能产生一种或多种化学或生物效应剂的作用。
背景技术
在化学和生化领域中可接受的范例之一是化学或生化效应剂,例如分子,通过多种物理化学力,例如离子、电荷或者分散力,或者通过电荷诱导的键的共价的分裂或者形成,与靶系统相互作用。这些作用力可以包括在效应剂或靶系统中的振动或旋转能模式。
因此,例如,当药物分子被施予生物体时,药物的作用包括其与目标成分的相互作用,例如膜,酶或者核酸成分,从而产生或者触发一连串药物的关联反应。同样的,如果将酶加入至生物底物中,该酶能够通过一些类型的空间配位与底物相互作用,系统中存在的能量模式被转换成活化状态或激活状态,这导致共价键的分裂或形成。
这个范例的显而易见的必然结果是在效应物-靶系统中,对目标环境中的效应剂的需求。可是,不知道或者不了解的是该需求是否与效应物的实际存在有关,或者至少对于某些效应物功能,这是否可以归因于表征效应物的能量模式的存在。如果效应物功能可以通过某些特有的能量模式,至少部分地被模拟,它可能通过使系统暴露于表征效应物的某些模式而在靶系统中“模拟”效应剂的作用。如果这样,自然会出现以下问题什么效应物-分子能量模式是有效的,它们怎样可以转换成可测量的信号形式,以及这些信号怎样用来作用于靶系统,更确切地说,至少模仿在靶系统中的分子的某些效应物功能。
通过将靶系统暴露于特征性效应物-分子信号来实现效应物-分子功能,而不需要效应剂的实际存在,这种可能性有许多重要的应用。通过将生物体暴露于药物特异性信号,可以替代通过使用药物治疗生物体,能实现同样的作用。在纳米制造领域,现在可以通过将表征多化合价效应物分子的信号引入该组装系统来催化和促进自组装模式,所述效应物分子能够促进需要的自组装模式。
本发明描述了通过表征效应物分子的低频信号的转换,在已知对效应物应答的系统中获得效应物特异性结果的装置和方法。

发明内容
一方面,本发明包括一种化学或生化制剂在对该制剂产生应答的系统中产生作用的方法。实际的方法中,该系统被放置在电磁转换器的磁场区域中,将低频,时域信号施加于电磁转换器,该信号的特征是在谱分析中呈现大量制剂特异性谱峰。这些谱峰从通过以下方法生成的低频,时域信号的频谱中识别(i)将这种化学或生化制剂放置在兼有磁屏蔽和电磁屏蔽的容器中,并且(ii)在时域信号谱线中有效产生这样可识别的谱峰的噪声水平,将噪声注入样本的同时,记录样本中的低频,时间相关信号。该样本暴露于转换器产生的磁场,以所施加的信号功率,并且持续足以对该系统产生制剂特异性效应的时间周期。
在一个普通的实施例中,施加于转换器的时域信号通过以下步骤产生(i)将这样的化学或生化制剂放置于兼有磁屏蔽和电屏蔽的容器中,并且(ii)在时域信号谱线中有效产生这样可识别的谱峰的噪声水平,将噪声注入样本的同时,记录样本中的低频时间相关信号。该方法中的步骤(ii)可以包括步骤(a)在给定的噪声振幅将噪声注入样本;(b)记录由叠加在注入噪声上的样本源辐射组成的电磁时域信号,以及(c)在选定的噪声水平范围内的多个噪声水平,重复步骤(a)和(b),以及(d)分析通过产生时域信号的谱线而生成的多个时域信号,并且基于谱线中的信息,识别最优的制剂特异性时域信号。
以如下步骤将噪声注入样本的噪声源(iia)包括功率可调的高斯噪声发生器和Helmholz线圈,并且其从噪声发生器接收直到1伏特范围内的所选噪声输出信号。
分析步骤(d)可以用以下步骤之一实现(i)生成直方图,其对于dc至8khz的范围内的所选频率范围上的每一个事件箱(event bin)f,示出每箱(bin)中的事件计数的数量,与超过给定阈值的箱的数量有关的得分;并且根据该得分选择时域信号,其中f是对时域信号采样的、分配给直方图的采样率。
(ii)使时域信号自相关,在dc至8khz范围之内的所选的频率范围内,生成自相关信号的FFT(傅立叶变换),将与超过平均噪声值的峰的数量相关的得分分配给FFT信号,并且根据该得分选择时域信号。并且(iii)计算多个定义的时间段中的每个上的时域信号的一系列傅立叶谱,在所选的在dc和8kHz之间的频率范围内,对傅立叶谱求平均;将与超过平均噪声值的峰的数量相关的得分分配给平均FFT信号,并且根据该得分选择时域信号。
使用与超导量子干涉仪(Squid)耦合的梯度计可以实现记录信号,所述注入包括将噪声注入梯度计。
系统可应答的制剂可以是一种配体,该配体能够通过配体特异性非共价相互作用,在生物学系统中与抗-配体细胞靶产生相互作用。
为使基因或基因簇上调或下调,该制剂可以是能够与基因启动子相互作用的化合物,并且低频时域信号通过记录该化合物的低频时域信号而产生。典型的化合物是l(+)阿拉伯糖,其中基因包括细菌lac操纵子。
为了产生可测量的生物学系统的抑制或生长以及生存,该制剂可以是能够与生物学系统生长和生存所需要的酶相互作用的化合物,以竞争性抑制所述酶。低频时域信号可以通过记录化合物的低频,时域信号产生。一种典型的制剂是甘胺酸磷酸盐(glyphosphate),并且靶酶为植株中的5-烯醇式丙酮莽草酸酯-3-磷酸(EPSP)合酶。另一个典型的制剂是phepropeptin D,靶酶是与真核细胞内的蛋白体结合的蛋白水解酶。
在治疗哺乳动物癌症的使用中,该制剂可以从其中之一选择(a)微管蛋白结合制剂,其中细胞靶是染色体微管纺锤体;(b)蒽环类抗生素(anthracycline),其细胞靶为双链DNA;(c)拓扑异构酶抑制剂,其中细胞靶为拓扑异构酶;(d)抗代谢剂,其细胞靶为细胞代谢所必需的酶;(e)免疫抑制素,其细胞靶为免疫应答细胞;以及(f)抑癌蛋白,其靶为细胞核内的DNA复制结构。典型的微管蛋白结合制剂为紫杉醇或者紫杉醇类似物。典型的抑癌因子为p53。
电磁转换器可以包括定义开口内部的线圈绕组,其中样本被放置在绕组的开口内部中。在另一个具体实施例中,转换器包括将被植入生物学系统的可植入线圈,例如,邻近的脉管区域。
可以选择在有效产生所选择的磁场强度的功率水平施加MIDS信号,该磁场强度在1-200mG(毫高斯)范围之内,优选10-100mG,更优选30-80mG。暴露可以通过间歇的MIDS信号的1小时或更多的周期的暴露来实现,例如在治疗期暴露12小时,不暴露12小时。
另一方面,本发明包括一种用于产生化学或生化制剂对制剂系统的作用的装置,所述系统应答于这种制剂。该装置包括(a)用于存储低频时域信号的存储设备,该信号的特征是在谱分析中呈现大量制剂特异性谱峰,其中制剂特异性谱峰从以下方法产生的低频时域信号的谱线中识别(i)将这样的化学或生化制剂放置于具有磁屏蔽和电磁屏蔽的容器中,并且(ii)当在时域信号谱线中有效产生可识别的谱峰的噪声水平将噪声注入样本的同时,记录来自样本的低频时域信号,并且装置中还包括(b)电磁转换器,其定义了活性磁场(active magnetic field)的区域,该活性磁场在电磁信号被施加于转换器并且样本被放置在其中的时候产生,以及(c)放大器有效地将存储设备连接到转换器上,以将信号施加给转换器,以所施加的信号功率并且持续一段时间周期,该时间周期足以在系统中对该系统产生制剂特异性效应,这样放置在转换器的活性磁场区域内。
电磁转换器可以包括线圈绕组和适于将样本放入其中的开口内部。在一个实施例中,转换器是具有一对对齐的电磁线圈的Helmholz线圈,这一对电磁线圈定义了它们中间的暴露的位置,其中样本,例如对象,被放置在该位置中。在另一个实施例中,转换器包括可植入的线圈。
存储设备可以在远离转换器和放大器的位置,并且信号从远离的位置传送到转换器。
根据另一个方面,提供了化学或生化活性制剂的低频时域信号,其特征是在谱分析中呈现大量制剂特异性谱峰,其中制剂特异性谱峰从以下方法产生的低频时域信号的谱线中识别
(i)将这样的化学或生化制剂放置于具有磁屏蔽和电磁屏蔽的容器中,并且(ii)当在时域信号谱线中有效产生可识别的谱峰的噪声水平将噪声注入样本的同时,记录来自样本的低频时域信号。
该信号本身可以通过以下步骤产生(i)将这样的化学或生化制剂放置于具有磁屏蔽和电磁屏蔽的容器中,并且(ii)当在时域信号谱线中有效产生可识别的谱峰的噪声水平将噪声注入样本的同时,记录来自样本的低频时域信号。以及当将下述本发明的详细描述与附图结合起来阅读时,本发明的这些以及其他的目的和特征将有更完整的体现。


图1是依据本发明的一个实施例形成的分子电磁信令检测装置的一个实施例的立体图;图2是图1所示的法拉第笼及其内容物的放大的详图;以及图3是图1和2中所示的衰减管之一的放大的横断面图;图4是图2所示的法拉第笼及其内容物的横断面图;图5是图1至图4所示的本发明的可选择的实施例的横断面图;图6是这里所述的Helmholtz变压器的支撑线圈的框架的放大详图;图7是可选择的电磁辐射检测系统的图;图8是包括在上述附图的检测系统中的处理单元的图;图9是图8之外另一个可选择的处理单元的图;图10是本系统执行的信号检测和处理的流程图;图11A是第一个样本辐射的谱线;图11B是第二个样本辐射的谱线;图12A和12B是饱和NaCl样本在500-530Hz之间的谱区内的谱图,其通过对非相关的时域样本信号进行傅立叶变换(12A),并且对互相关的样本谱进行傅立叶变换而产生(12B);图13A和13B是烷基醚硫酸盐样本在500-530Hz之间的谱区中的谱图,其通过对非相关的时域样本信号进行傅立叶变换(13A),并且对互相关的样本谱进行傅立叶变换而产生(13B);图14A-14F是样本去离子水(14A),饱和的NaCl溶液(14B),1%NaCl去离子水溶液(14C);饱和的NaBr样本(14D),烷基醚硫酸盐去离子水溶液(14E),以及无样本(14F)在500-530Hz之间的谱区中的谱线;图15A-15F是氨基酸样本溶液1∶100wt每体积(15A)以及增大稀释到1∶10,000(15B),1∶1百万(15C),1∶1亿(15D),1∶1百亿(15E和15F),在500和535Hz之间的谱区中的谱线,其中图15A-15E中的频谱是通过50秒记录和40分钟相关而产生的,并且图15F中的频谱由12小时相关的4:25分记录生成;图16是示意图,显示了分子电磁信令检测装置的可选择的实施例;图17A是图16的可选择实施例的横断面图;图17B是图17A的部分放大图;图18是图17B的横断面等距视图;图19是图9的可选择的处理单元的图;图20显示了用于本发明直方图谱图方法的数据流的高级流程图;图21是依据本发明生成谱图直方图的算法的流程图;以及图22A-22D是取自四个不同的噪声功率水平的样本的频谱的直方图;图23A-23C是用于拍摄显示用于生成和显示谱图直方图的用户界面的电脑屏幕;图24是根据本发明方法的第二实施例的用于识别最优时域信号的步骤的流程图;图25A-25E显示了包含40%除草剂复合物的样本的时域信号部分(25A),来自图25A中的样本的自相关的时域信号的FFT,所述时域信号在噪声水平70.9-dbm(25B),74.8-dbm(25C和25D),以及78.3-dbm(25E)记录;图26是图25中样本的自相关得分和噪声设置的对比图;图27,图24是根据本发明方法的第三实施例识别最优时域信号的步骤的流程图;图28A显示了典型的转换试验中的转换仪器的布置;以及图28B和28C显示了典型的转换试验中使用的转换线圈和容器(28B);以及转换线圈在振动台上的布置(28C);图29A-29B显示了Helmholz线圈转换器中的磁场作为线圈(29A);螺管转换器(29B)和改良的螺管转换器(29C)之间的位置的函数的图;图30A-30E显示了AraL_(+)-可诱导绿色荧光蛋白在E线圈中对不同的刺激的应答,包括MIDS AraL(+)信号;图31显示了应答不同刺激的茎长的箱式图,包括除草剂的MIDS信号(31A),茎细胞长的概率图(31B),以及同样实验中茎长的经验CDF(31C);图32A和32B是甜豌豆芽没有暴露(32A)和暴露(32B)于MIDS除草剂信号的照片;图33显示了不同的刺激对于茎长的影响,包括生长刺激物MIDS信号,表达为茎长应答刺激的箱式图;图34是显示了用生化抑制剂或MIDS信号治疗1.5小时之后20S蛋白体活性的抑制(深色条)和治疗之后20小时化验的抑制(浅色条)的条状图;以及图35A和35B在小鼠中显示了肿瘤重量变化的百分比,作为用肿瘤组织接种后天数的函数,小鼠以40或60mG,单独用赋形剂,紫杉醇,或紫杉醇MIDS信号治疗(Fig.35A);以及同样的4组动物的肿瘤重量变化的百分比(图35B)。
在所有特征部件或技术的表述中,参考标号中,第一次引入部件的附图标号是最显著的。
具体实施例方式
I.定义下列术语具有下述的定义除非另外指出“呈现分子转动(molecular rotation)的样本”是指样本材料,其可以是气体,液体或固体形式(固体金属除外)其中一个或多个分子化合物或原子离子组成或出现在呈现转动的样本中。
“磁屏蔽”是指由于屏蔽材料的导磁性抑制或阻止磁通量的通过的屏蔽。
“电磁屏蔽”是指,例如,标准法拉第电磁屏蔽。
“时域信号”或“时序信号”是指具有随着时间改变的瞬时信号特性的信号。
“样本源辐射”是指来自样本分子运动的磁通量的发射,例如磁场中分子偶极子的旋转。
“高斯噪声”表示具有高斯功率分布的随机噪声。
“稳定的高斯白噪声”表示没有可预报的将来的分量的随机高斯噪声。
“均匀噪声”表示具有恒定振幅的噪声。
“频域谱”是指时域信号的傅立叶频率曲线。
“谱分量”是指在时域信号里的单一或重复的性质,它可以由频率,振幅,和/或相畴测量。典型的谱分量是指存在于频域中的信号。
“相似的样本”,参考第一样本,是指相同的样本或与第一样本具有基本相同样本成分的样本。
“法拉第笼”是指由于不必要的电磁辐射而提供了至地面的电路径的电磁屏蔽结构,借此使电磁环境平静。
“谱特征得分”是指以制剂特异性谱峰的数量和/或振幅为基础的得分,所述数量和/或振幅通过在对已经由适当方法处理过的制剂或样本记录的时域信号中,观察可选择的低频范围,例如,DC到1kHz或DC到8kHz,所述例如在这里描述的三种方法中的一种,从而揭示制剂或样本特有的可辨认的谱特征。
“最优制剂特异性时域信号”是指具有最大值或接近最大值的谱特征得分的时域信号。
“MIDS”或“分子查询和数据系统”是指从样本记录的时域信号,且包含样本相关的在信号的谱分析中显而易见的谱峰。MIDS信号优选最优的或被增强,如下面详细描述。MIDS信号还可以是通过结合和/或过滤个体样本时域信号来构成,如在这里描述的。
II.用于产生和处理低频时域信号的装置在下面详细描述的是一种用于检测,处理,和呈现所关注的样本的低频电磁辐射或信号的系统和方法。在一种实施方式中,一种已知的白噪声或高斯噪声信号被引入至样本中。高斯噪声配置成容许来自样本的电磁辐射被信号检测系统充分检测到。检测的信号组在一起处理以确保重复率和统计关联性。所得到的辐射图案或频谱可以作为一种特殊的物质显示,储存,和/或识别。
下列描述提供了对本发明的完全理解的特定细节和实施本发明实施例的具体细节。然而,本领域技术人员将理解没有这些细节也可以实践本发明。在其它的情况中,公知的结构和功能没有在细节中显示或描述,以避免不必要的模糊对本发明实施方式的描述。
如下列内容中说明的,本发明的实施方式直接提供了一种用于低域值分子电磁信号的可重复检测和记录的装置和方法。磁屏蔽的法拉第笼屏蔽了样本材料和来自外部电磁信号的探测装置。在磁屏蔽的法拉第笼里,线圈注入白噪声或高斯噪声,非铁的托盘固定住样本,以及梯度计检测低域值分子电磁信号。该装置还包括超导量子干涉器件(“SQUID”)和前置放大器。
通过放置样本在磁屏蔽的法拉第笼里接近噪声线圈和梯度计来使用该装置。白噪声被穿过噪声线圈注入并且被调整直至分子电磁信号通过随机共振而增强。增强的分子电磁信号,被来自法拉第笼的外界干扰和噪声线圈产生的磁场所屏蔽,然后被梯度计和SQUID检测和测量。该信号随后被放大且传送至任何适当的记录或测量仪器。
参考图1,显示有屏蔽结构10,其在由外向内的方向包括导电线笼16,它是磁屏蔽,和内部导电线笼18和20,它们提供了电磁屏蔽。在另一个实施方式里,外部磁屏蔽由具有铝-镍合金涂层的固体镀铝材料形成,并且由两个内壁结构提供电磁屏蔽,每一个由固体铝形成。
参考图2,法拉第笼10在顶端是开口的,并且包括侧开口12和14。法拉第笼10还由三个铜网笼16,18和20组成,三者互相紧靠。铜网笼16,18和20中的每一个都通过在每个笼子之间的介电屏障(没有显示)和其它笼电绝缘。
侧开口12和14还包括衰减管22和24,当将笼的内部从外干扰源隔离时提供通向法拉第笼10内部的通道。参考图3,衰减管24是由三个铜网管26,28和30组成,三者互相紧靠。外部的铜网笼16,18和20各自分别电连接至铜网管26,28和30中的一个。衰减管24还由顶盖32覆盖,顶盖有洞34。衰减管22相似地由铜网管26,28和30组成,但不包括顶盖32。
还是参考图2,低密度非铁样本托盘50安装在法拉第笼10内部。安装样本托盘50使得可以穿过衰减管22和侧开口12从法拉第笼10移除样本托盘。三根杆52,其中的每一根在长度上都大于从法拉第笼10的中心垂直轴到衰减管22最外缘的距离,附着在样本托盘50上。三根杆52适于与衰减管22内部曲线相一致,以便于样本托盘50可以通过使杆在衰减管中静止来放置在法拉第笼10的中心。在所示的实施方式中,样本托盘50和杆52可由玻璃纤维环氧树脂制造。对于那些本领域的技术人员而言很容易想到样本托盘50和杆52可由其它非铁材料制造,并且托盘可以通过别的方法安在法拉第笼10里,例如通过一根单独的杆。
还是参考图2,安装在法拉第笼10里且在样本托盘50上面的是低温的杜瓦瓶(Dewar flask)100。在公开的实施方式中,杜瓦瓶100适于安装在法拉第笼10顶部的开口中,并且是由Tristan Technologies公司制造的型号BMD-6的液体氦杜瓦瓶。杜瓦瓶100用玻璃纤维环氧树脂合成物制造。带有非常狭窄视野的梯度计110安装在杜瓦瓶100里适当的位置以便它的视野环绕样本托盘50。在所示的实施方式里,梯度计110是第一级轴检测线圈,标称直径1厘米,有2%的差额,并且由超导体形成。梯度计可以是平面梯度计以外的任意形式的梯度计。梯度计110连接至一种低温直流超导量子干扰设备(“SQUID”)120的输入线圈。在公开的实施方式中,SQUID是由Tristan Technologies公司制造的型号LSQ/20LTS dc的SQUID。那些本领域的技术人员可以意识到在不背离本发明实质和范围时也可以使用高温或交流SQUID。在可选实施方式中,SQUID 120包括噪声抑制线圈124。
当测量磁场时,公开的梯度计110和SQUID 120的组合具有5微特斯拉/√赫兹的灵敏度。
SQUID 120的输出连接到由Tristan Technologies公司制造的型号SP的低温电缆130。低温电缆130能够经受住杜瓦瓶100内部和外部的温度并且将信号从SQUID 120传送至磁通锁定环140,它安装在法拉第笼10和杜瓦瓶100的外部。公开的实施方式中的磁通锁定环140是由Tristan Technologies公司制造的iFL-301-L磁通锁定环。
参考图1,磁通锁定环140还放大且输出从SQUID 120经由高级输出电路142至iMC-303 iMAGSQUID控制器150接收的信号。磁通锁定环140还通过型号CC-60的6米光导纤维合成物连接电缆144连接至SQUID控制器150。光导纤维连接电缆144和SQUID控制器150由TristanTechnologies公司制造。控制器150安装在磁屏蔽笼40的外部。光导纤维连接电缆144承载从SQUID控制器150到磁通锁定环140的控制信号,这进一步降低了对需要测量的信号的电磁干扰的可能性。对于本领域技术人员在不背离本发明实质和范围的情况下可以使用其它的磁通锁定环,连接电缆,以及Squid控制器是显而易见的。
SQUID控制器150还包括高分辨率模拟-数字转换器152,用于输出数字化信号的标准GP-IB总线154,和用于输出模拟信号的BNC连接器156。在所示实施方式里,BNC连接器穿过接插线162连接双线示波器160。
参考图2,当样本托盘完全插进法拉第笼10时,二元Helmholtz变压器60安装在样本托盘50的任一边。在所示实施方式里,Helmholtz变压器60的线圈绕组62和64被设计成在直流电至50千赫范围内操作,具有25千赫的中心频率和8.8兆赫的自振频率。在所示的实施方式里,线圈绕组62和64通常外形上是矩形的且大约8英尺高4英尺宽。可以使用其它的Helmholtz线圈形状,但应该有满足如下要求的形状和大小,即梯度计110和样本托盘50放置在Helmholtz线圈产生的磁场内。线圈绕组62和64中的每一个安装在两个低密度非铁框架66和68中的一个上。框架66和68互相铰链连接并由支架70支撑。框架66和68可滑动地附着在支架70上以允许与杜瓦瓶100的较低部分有关的框架的垂直运动。框架的运动允许对Helmholtz变压器60的线圈绕组62和64进行调节以改变梯度计110所接受的白噪声的振幅。支架70搁置在上面或被粘上,例如,用环氧树脂粘合在法拉第笼10的底部上。在所示的实施方式中,框架66和68以及支架70由玻璃纤维环氧树脂制成。在不背离本发明实质和范围的情况下,可以在样本托盘50周围使用变压器或线圈的其它布置。
参考图4,显示了法拉第笼和其内含物的横断面图,显示了与杜瓦瓶100和法拉第笼10有关的Helmholtz变压器60的绕组62。注意图4中的样本托盘50和样本200的定位。
参考图5,显示了可选择的实施方式,其中Helmholtz线圈绕组62和64按垂直方向固定,并且附加噪声线圈300置于样本托盘50下面。附加噪声线圈300的绕组基本上垂直对着Helmholtz变压器60的垂直绕组62和64,并且附加噪声线圈300的绕组因此基本上与法拉第笼10的底部平行定向。
在这个可选择的实施方式里,噪声将由相同的双绞线(未示出)馈送到噪声线圈300,如供应Helmholtz线圈一样。噪声源由用于提供噪声给Helmholtz线圈的相同的噪声发生器发起。噪声将在经由附加输出连接的噪声发生器处,或经由从输出连接到噪声发生器的平衡分离器被取样。在附加噪声线圈300处的噪声信号的衰减将穿过可调的RF信号衰减电路,其中许多在商业上可用,或者经由合适的定值RF衰减滤波器的串连。
参考图6,可以看见支撑Helmholtz变压器60的线圈的框架的细节;图6的参照点是从图4视图旋转90度,并且省略了法拉第笼10。框架66和68布置为基本在垂直位置显示Helmholtz线圈的线圈绕组并且互相平行。框架66’和68’示出了关于连接框架的铰链连接的轴的框架的旋转,以便布置Helmholtz变压器的线圈绕组成为彼此不平行的关系。
还是参考图1,振幅可调的白噪声发生器80在磁屏蔽笼40的外部,并且由电缆82通过滤波器90电连接到Helmholtz变压器60。参考图3,电缆82穿过侧开口12,衰减管24,并且经洞34穿过顶盖32。电缆82是同轴电缆,其还包括被内外磁屏蔽86和88分别围绕的铜导体84的双绞线。在其它的实施方式中,导体可以是任何非磁性的导电材料,例如银或金。内外磁场屏蔽86和88终止在顶盖32,留下双绞线84横跨在图1中所示的从末端顶盖到Helmholtz变压器60的剩余距离。当如图1所示外部磁屏蔽电连接到磁屏蔽笼40时,内部磁屏蔽86穿过顶盖32电连接到法拉第笼16。
参考图1,白噪声发生器80能产生近似均匀的噪声,其频谱穿过从0到100千赫的频谱。在所示实施方式里,滤波器90滤出在50千赫之上的噪声,但在不背离本发明实质和范围的情况下,可以使用其它频率范围。
白噪声发生器80还可以通过接插线164电连接到双线示波器160的其它输入端。
参考图1,2和3,需要测量的物质样本200被放置在样本托盘50上并且样本托盘放置在法拉第笼10里。在第一个实施方式里,白噪声发生器80用来注入白噪声穿过Helmholtz变压器60。噪声信号在梯度计110中产生感应电压。梯度计110中的感应电压随后被检测到且经SQUID 120放大,来自SQUID的输出通过磁通锁定环140被进一步放大,并且被传送到SQUID控制器150,然后传送到双线示波器160。双线示波器160还可以用于显示由白噪声发生器80生成的信号。
白噪声信号通过改变白噪声发生器80的输出和旋转围绕样本200的Helmholtz变压器60来调整,如图2所示。关于框架66和68的铰链连接的轴的Helmholtz变压器60的旋转改变它关于梯度计110的相位。取决于所期望的相位改变,框架66和68的铰链连接允许当围绕样本托盘50旋转约30到40度时绕组62和64保持彼此平行。铰链连接还允许绕组62和64旋转超出平行多达大约60度,以便改变由Helmholtz变压器60生成的磁场关于梯度计110的信号定相。例如,典型的相位调整包括这种平行以外的方向,即使在特定的环境下可能优选其它方向,以调节不规则形状的样本200。噪声被施加并调整直至被检测的噪声超过分子电磁辐射30到35分贝。在这个噪声水平,噪声呈现出通过公知的随机共振现象的分子电磁信号的特征。当反映经梯度计110检测的信号的示波器迹线不同于反映直接来自白噪声发生器80的信号的迹线时,可观察到所寻找的(sought)随机产物。在可选择的实施方式中,信号可被记录和/或经任何商业上可用的仪器处理。
在一个可选择的实施方式里,检测分子电磁信号的方法还包括穿过SQUID 120的噪声抑制线圈124注入与Helmholtz变压器60处施加的原始噪声信号有180°相位差的噪声。当反映由梯度计110检测的信号的示波器迹线变成非随机的,那么可以观察到所寻找的随机产物。
不管噪声如何注入和调整,当谱峰值增加时随机产物还是可以通过观察确定。谱峰值可由示波器160上的线性曲线或数值,或由其它公知的测量设备观察到。
本发明的实施方式提供了一种用于检测没有外界干扰的极低域值分子电磁信号的方法和装置。它们还提供了那些易被各种信号记录和处理仪器使用的格式的信号的输出。
目前参考图7,显示了上面附图所示的分子电磁辐射检测和处理系统的一种可选择的实施方式。系统700包括检测单元702,其耦合到处理单元704。虽然处理单元704在检测单元702的外部显示,但至少处理单元的一部分可定位在检测单元里面。
检测单元702,如图7横断面图所示,包括多个互相嵌套或同轴的元件。样本室或法拉第笼706嵌套在金属笼708里面。样本室706和金属笼708中的每一个可由铝材料构成。样本室706可以维持为真空并且可以控制温度为预先设定的温度。金属笼708配置成如低通滤波器的功能。
在样本室706和金属笼708之间以及环绕样本室706的是一组平行加热线圈或元件710。一个或多个温度传感器711也位于加热元件710和样本室706附近。例如,四个温度传感器可以置于围绕样本室706外部的不同位置。加热元件710和温度传感器711配置成维持样本室706内部的某一温度。
屏蔽712环绕着金属笼708。屏蔽712配置成为样本室706提供附加磁场屏蔽或隔离。屏蔽712可由导线或其它磁场屏蔽材料组成。当由样本室706和/或金属笼708提供了充分的屏蔽时,屏蔽712是可任意选择的。
围绕屏蔽712的是带G10绝缘体的低温层716。低温剂可以是液体氦。低温层716(也称为低温杜瓦瓶)处于4开尔文温度的操作温度。围绕低温层716的是外部屏蔽718。外部屏蔽718由镍合金组成并且配置成磁屏蔽层。沿着笛卡儿坐标系的三个直角平面,由检测单元702提供的磁屏蔽总数大约-100dB,-100dB和-120dB。
上述各种元件是通过气隙或介电屏障(未示出)彼此电隔离的。还应该明白为了减少叙述,元件没有显示出相对于彼此的大小。
样本架720可手工或用机械放置在样本室706里面。样本架720可以从样本室706的顶端降低,升高,或移除。样本架720由不导入涡电流(Eddycurrent)和内部分子转动显示很少或没有的材料组成。举例来说,样本架720可由高品质玻璃或硬质玻璃(Pyrex)组成。
检测单元702配置成处理固体,液体或气体样本。在检测单元702中可以利用各种样本架。例如,取决于样本的大小,可以利用较大的样本架。举另一个例子来说,当样本对空气起反应时,样本架可配置成密封的或围绕样本形成空气密封。在另一个例子中,当样本处于气态时,样本可以不用样本架720而引入样本室706里面。对这样的样本,样本室706维持真空。在样本室706顶部的真空密封721帮助维持真空和/或调节样本架720。
感测线圈722和感测线圈724,还可以指如检测线圈,分别被设置在样本架720的上面和下面。感测线圈722,724的线圈绕组配置成在直流(DC)至大约50千赫范围内运操作,具有25千赫的中心频率和8.8兆赫的自共振频率。感测线圈722,724是二阶导数形式,并被配置成达到大约100%的耦合。在一种实施方式里,线圈722,724通常外形上是矩形并且由G10接线柱固定在位置上。线圈722,724起二阶导数梯度计的作用。
Helmholtz线圈726和728可垂直放置在屏蔽712和金属笼708之间,如在此处说明的。线圈726和728中的每一个可以彼此独立地上升或下降。线圈726和728,还指白噪声或高斯噪声发生线圈,处于室温和环境温度。由线圈726,728产生的噪声为大约0.10高斯。
来自样本的辐射和线圈722,724之间的耦合程度可以通过改变样本架720相对于线圈722,724的位置,或改变线圈726,728中的一个,或两者相对于样本架720的位置而改变。
处理单元704与线圈722,724,726和728电耦合在一起。处理单元704指定由线圈726,728将白噪声或高斯噪声注入样本。处理单元104也接收线圈722,724的感应电压,它来自和注入的高斯噪声混合的样本的电磁辐射。
参考图8,使用本发明的方面的处理单元包括样本托盘840以允许插入和移除样本842,法拉第笼844和Helmholtz线圈746。SQUID/梯度计检测器组件848被放置在低温杜瓦瓶850里面。磁通锁定环852耦合在SQUID/梯度计检测器组件848和SQUID控制器854之间,SQUID控制器854可以是由Tristan提供的型号iMC-303 iMAG的多路控制器。
模拟噪声发生器856给锁相环858提供噪声信号(如上所述)。锁相环的X-轴输出被供给Helmholz线圈846,并且可能被衰减,例如20dB。锁相环的y-轴输出经信号分离器860分离。y-轴输出的一部分被输入在SQUID的噪声消除线圈,它对于梯度计有单独的输入。y-轴信号的其它部分被输入示波器862中,例如像Tektronix TDS 3000b(例如,3032b型号)的具有傅立叶函数的模拟/数字示波器。就是说,锁相环的x-轴输出驱动Helmholz线圈,并且y-轴输出,以转换的形式,被分离输入SQUID和示波器。因而,锁相环起信号转换器的作用。示波器迹线被用于监控模拟噪声信号,例如,用于测定何时达到足够产生非平稳谱分量的噪声水平。模拟磁带记录器或记录设备864,与控制器854耦合,记录从该设备输出的信号,并且优选是宽频(例如,50kHz)记录器。PC控制器866可以是经由,例如,RS 232端口与控制器854通过接口连接的基于MS Windows的PC。
在图9中,显示了处理单元的另一种实施方式的示意图。双相锁定放大器202配置成提供第一信号(例如,“X”或噪声信号)给线圈726,728和第二信号(例如,“y”或噪声消除信号)给超导量子干扰设备(SQUID)206的噪声消除线圈。放大器202配置成不用外部引用而锁定并且可以是PerkinsElmer型号7265 DSP的锁定放大器。该放大器以“虚拟方式”工作,其中它锁定原始基准频率,然后移除基准频率以允许它自由运转并锁定“噪声”。
模拟噪声发生器200是和放大器202电耦合的。发生器200配置成产生或感应经由放大器202在线圈726,728的模拟高斯白噪声。举例来说,发生器200可以是由General Radio制造的型号1380。
阻抗变换器204在SQUID 206和放大器202之间电耦合。阻抗变换器204配置成提供在SQUID 206和放大器202之间的阻抗匹配。
SQUID 206的噪声消除特征可被开启或关闭。当噪声消除特征开启时,SQUID 206能够消除或使来自被检测辐射的注入噪声分量无效。为了提供噪声消除,到线圈726,728的第一信号是超过所寻找的要被检测的分子电磁辐射20dB或35dB的噪声信号。在这个水平,注入的噪声呈现出通过随机共振的分子电磁信号的特征。到SQUID 206的第二信号是噪声消除信号,并从振幅足以无效SQUID输出端噪声的第一信号转换而来(例如,关于第一信号的相位超出180度)。
SQUID 206是低温直流元件SQUID。举例来说,SQUID 206可以是由Tristan Technologies公司制造的型号LSQ/20LTS直流SQUID。可选择地,可使用高温或交流SQUID。组合的线圈722,724(例如,梯度计)和SQUID206(统称SQUID/梯度计检测组件)具有大约5微特斯拉/√赫兹的磁场测量灵敏度。线圈722,724中的感应电压被SQUID 206检测和放大。SQUID 206的输出电压是大约在0.2-0.8微伏范围里的电压。
SQUID 206的输出是SQUID控制器208的输入。SQUID控制器208配置成控制SQUID 206的操作状态并进一步调节所检测的信号。举例来说,SQUID控制器208可以是由Tristan Technologies公司制造的iMC-303 iMAG多路SQUID控制器。
SQUID控制器208的输出被输入放大器210。放大器210配置成提供在0-100dB范围内的增益。当噪声消除节点在SQUID 206处开启时,提供了大约20dB的增益。当SQUID 206没有提供噪声消除时,提供了大约50dB的增益。
放大的信号被输入记录器或存储设备212。记录器212配置成将模拟放大信号转化成数字信号并存储该数字信号。在一种实施方式里,记录器212存储8600个数据点每赫兹并且能处理2.46Mbits/秒。举例来说,记录器212可以是Sony数字录音磁带(DAT)记录器。使用DAT记录器,原信号或数据集可被传送至第三方以按需求显示或进行特定处理。
低通滤波器214对来自记录器212的数字化数据集进行滤波。低通滤波器214是模拟滤波器并且可以是Butterworth滤波器。截止频率大约在50kHz。
接着带通滤波器216对那些经过滤波的数据集进行滤波。带通滤波器216配置成具有频带宽度在DC到50kHz之间的数字滤波器。带通滤波器216可适用于不同的频带宽度。
带通滤波器216的输出被输入傅立叶变换处理器218。傅立叶变换处理器218配置成把在时域里的数据集转换成在频域里的数据集。傅立叶变换处理器218执行快速傅立叶变换(FFT)类型的转换。
经傅立叶变换的数据集是相关和对比处理器220的输入。记录器212的输出也输入到处理器220。处理器220配置成使数据集和先前记录的数据集相关,确定阈值,以及执行噪声消除(当SQUID 206没有提供噪声消除时)。处理器220的输出是代表样本的分子低频电磁辐射谱的最终数据集。
用户接口(UI)222,例如图形用户接口(GUI),也可以连接到至少滤波器216和处理器220以指定信号处理参数。滤波器216,处理器218,和处理器220可当硬件,软件,或固件应用。例如,滤波器216和处理器218可以应用在一个或多个半导体芯片中。处理器220可以是在计算设备中应用的软件。
放大器以“虚拟方式”工作,其中它锁定至原始基准频率,然后移除基准频率以允许其自由运转和锁定“噪声”。模拟噪声发生器(General Radio生产的,真正的模拟噪声发生器)对Helmholz和噪声消除线圈分别需要20dB和45-dB的衰减。
Helmholz线圈可以具有百分比1/100th差额(balance of 1/100thof a percent)的大约1立方英寸的最佳点(sweet spot)。在可选择的实施方式里,Helmholtz线圈既可以垂直移动,也能旋转移动(围绕垂直通路),以及以圆形形状从平行到展开。在一个实施方式里,SQUID,梯度计,和驱动变换器(控制器)分别具有1.8,1.5和0.3微亨值。Helmholz线圈在最佳点可以具有0.5高斯/安培的灵敏度。
对于随机响应大约需要10到15微伏。通过注入噪声,系统已经提高了SQUID设备的灵敏度。没有噪声的SQUID设备具有大约5毫微微特斯拉的灵敏度。该系统已经能够通过注入噪声和使用随机共振响应改善灵敏度25到35dB,总计近乎增长1,500%。
在接收和记录来自系统的信号之后,计算机,例如主计算机,超级计算机或高性能计算机进行预处理和后处理,这样通过使用Systat Software ofRichmond CA制造的Autosignal软件产品,用于预处理,而用Flexpro软件产品进行后处理。Flexpro是由Dewetron公司提供的数据(统计)分析软件。下列方程式或选项可以在Autosignal和Flexpro产品中使用。
正变换x^k=1NΣn=0N-1xne-2πikn/N]]>逆变换
x^k=1NΣn=0N-1xne-2πikn/N]]>FFT算法最精确的N使用Temperton的质因数FFT(C.Temperton,″Implementation of a Self-Sorting In-Place Prime Factor FFT Algorithm,Journalof Computation Physics,v.58,p.283,1985)。
数据锥度窗口(Data Tapering Windows)[cs4 BHarris min]0.35875-0.48829*cos(2*Pi*i/(n-1))+0.14128*cos(4*Pi*i/(n-1))-0.01168*(6*Pi*i/(n-1)),i=0.n-1[矩形]没有可用的固定形状锥度(示波器)量值sqrt(Re*Re+Im*Im)[Re=实分量,Im=虚分量]振幅2.0*sqrt(Re*Re+Im*Im)/ndb,分贝10.0*log10(Re*Re+Im*Im)平均重复采样(replicates)重复采样是基于符合1e-8分数精度内的X-值的。
基准减法(reference subtraction)基准信号减法(基线噪声)是在Y轴(振幅)沿着X(时间)轴的每一点(通路)上执行的。负数Y值然后归零。
互相关函数使用求和法和积分来计算互相关函数。因为信号是瞬态的,该相关函数使用直接的乘法和积分进行计算。位于源通路(数据序列)之外的计算所需的所有数值被置为0。t<0的点也要被计算。
傅立叶显著性水平蒙特卡洛数据适合于参数模型。其中数据大小N是唯一因数,使用单变量TableCurve 2D参数模型。对于分段FFT,其中段大小和重叠(overlap)是附加的影响,应用三元Chebyshev多项式。这些是在Autosignal下选择的选项。一个可以具有单独分析的数据集,或者可以以重叠方式分析,在重叠方式中数据集一将被分析,然后数据集一的后半部分和数据集二的前半部分,然后数据集二,然后后半部分。
图10显示了系统100执行的信号检测和处理的流程图。当样本是关注样本时,执行至少四个信号检测或数据运行(run)第一个数据运行在没有样本的时间t1,第二个数据运行在有样本的时间t2,第三个数据运行在有样本的时间t3,以及第四个数据运行在没有样本的时间t4。执行和采集来自多于一个数据运行的数据集增加了最后的(例如,相关的)数据集的准确性。在四个数据运行中,系统100的参数和条件是保持不变的(例如,温度,放大的倍数,线圈的位置,噪声信号,等等)。
在方框300,适当的样本(或如果它是第一或第四数据运行,没有样本),放置在系统100里。所给样本,在没有注入噪声的情况下,在振幅等于或小于大约0.001微特斯拉时发出DC-50kHz范围内的电磁辐射。为了捕获这些低辐射,在方框301注入白高斯噪声。
在方框302,线圈722,724检测表示样本辐射和注入噪声的感应电压。感应电压包括作为数据运行持续时间的时间的函数的连串的电压值(振幅和相位)。数据运行可持续2-20分钟并且因此,与数据运行对应的数据集包括2-20分钟的作为时间的函数的电压值。
在方框304,当感应电压正在被检测时注入的噪声被消除了。当SQUID206的噪声消除特征关闭时此方框被省略了。
在方框306,数据集的电压值放大达20-50dB,这取决于噪声消除是否在方框304发生。在方框308,放大的数据集经历模拟-数字(A/D)转换并且存储在记录器212中。数字化的数据集可包括数百万的数据行。
在将所获取的数据集存储之后,在方框310执行检查以判断是否已经出现至少四个对于样本的数据运行(例如,已经获取至少四个数据集)。如果已经获得对于所给样本的四个数据集,那么在方框312发生低通滤波。否则,启动下一个数据运行(返回方框300)。
在对数字化数据集进行低通滤波(方框312)和带通滤波(方框314)之后,数据集被转换成位于傅立叶变换方框316处的频域。
接下来,在方框318,相似数据集在每一个数据点彼此相关。例如,与第一数据运行相应的第一数据集(例如,基线或环境噪声数据运行)和与第四数据运行相应的第四数据集(例如,另一个噪声数据运行)彼此互相关。如果在给定频率处的第一数据集的振幅值和在给定频率处的第四数据集的振幅值相同,那么,所给定频率的相关值或数值将是1.0。可选择地,相关值的范围可以设定在0-100之间。这样的相关或比较也会发生在第二和第三数据运行(例如,样本数据运行)。因为所获取的数据集被存储,当剩余的数据运行完成时它们可在稍后被访问。
当SQUID 206没有提供噪声消除时,那么预定的阈值水平被应用于每个相关的数据集以消除统计上无关的相关值。可以使用许多种阈值,取决于数据运行的长度(数据运行越长,获取的数据的准确性越高)和样本实际辐射谱与其它类型样本的可能相似性。除阈值水平之外,相关性被平均。阈值和平均相关性的使用导致注入的噪声分量在所得到的相关数据集中变得很小。
如果在SQUID 206提供了噪声消除,那么阈值和平均相关性的使用不是必需的。
一旦两个样本数据集已经被定义成相关的样本数据集并且两个噪声数据集已经被定义成相关的噪声数据集,那么相关的噪声数据集被从相关的样本数据集减去。得到的数据集是最终的数据集(例如,代表样本辐射谱的数据集)(方框320)。
因为每赫兹具有8600数据点以及最终的数据集能够具有DC-50kHz频率范围内的数据点,最终的数据集可包括数百百万行的数据。每一行的数据可包括频率、振幅、相位和相关值。
在图11A和11B中,显示了样本辐射谱的例子。图11A显示了与饱和氯化钠溶液样本的频谱相应的傅立叶曲线400。图11B显示了与酶样本的频谱相应的傅立叶曲线500。
参考图16,上述系统的另一个可选择实施方式现在描述为系统1600。总的来说,在这里描述的代替方式和可选择的实施方式本质上和前面描述的实施方式相似,并且相同的参考号通常指共同的元件和功能。只有在结构或操作中的明显差别会详细描述。
二阶导数梯度计显示为1602,靶样本放置在上下配对线圈之间。样本两侧的两个内部线圈互相补偿,而两个外部线圈(顶部和底部线圈)各自互相补偿,并且与两个内部线圈相对。这样的排列允许从样本提取更多的信号和改进噪声抑制。
虽然下面在附图中显示和描述更多细节,系统1600采用了同心元件序列和沿着中心轴延伸进杜瓦瓶的排列。步进电动机1604允许将样本轴向地放置在这种同心元件的排列里面。具体地,样本可以放置在梯度计1602中间所期望的位置。
同样地,测微调整机构1606,例如机械测微计或步进电动机,允许Helmholtz线圈可以相对系统中的元件对准(例如样本和梯度计)。象Helmholtz线圈这样的调整辅助系统1600的制造和校准,也允许在系统里面磁场的精确对准,例如提供关于梯度计1602的均匀磁场。它也可能对提供磁场梯度的场偏移(field off set)或改变有用,从而产生更好的随机结果,以抵消系统中的噪声,或提供其它益处。
图17A,17B,和18更清楚地显示出系统1600里面元件的同心排列,其中样本管轴向地延伸穿过低通滤波金属屏蔽1802(例如不锈钢合金)的中心以使2kHz以下的信号通过。外部磁场(MU)屏蔽围绕着梯度计,Helmholtz线圈和样本。系统1600的排列关于图通常是一目了然的。
随机白噪声发生器,型号1381,由General Radio制造,如上所述,可以被由Noise/Com制造的可编程高斯白噪声发生器取代。这样的发生器采用两个输出端,一个从另一个转化而来。一个输出端可被连接至Helmholtz线圈,同时另一个(转化的)输出端连接至上述的SQUID噪声消除线圈。
同样地,如图19所示,上述Tektronix数字示波器可以被由StanfordResearch Systems制造的双通道动态信号分析器1902,型号SR785取代。这个信号分析器可以通过采样多重时域信号并对经过多频域FFT’s的它们求平均来处理输入信号。这可能导致所有非随机信号分量的全频谱频域记录。可能发生的其它变化包括用数字通用光盘(DVD)记录器1904替代数字录音磁带存储系统。此外,可以使用由Keithley制造的数据采集板1906,型号3801,它运行软件以生成直方图,如下所述。
在图19中所示的可选择实施方式里,噪声消除线圈1908连接在梯度计和SQUID之间。(虽然显示了一阶导数梯度计,也可以使用二阶导数梯度计,例如在图16中所示的。)虽然没有在图19中显示,但反向(inverted)的噪声通路(关于施加给Helmholtz线圈的噪声反向)可适用于噪声消除线圈1908(且可首次穿过阻抗变换器,其衰减噪声,例如,45dB)。在可选择的实施方式里,未示出,噪声消除线圈可放置在SQUID 120里面,SQUID输入线圈和输出线圈之间。
III.产生最优时域信号的方法根据本发明的一方面,已经发现在对于给定样本获得的低频时域信号中,样本相关的频谱特征可通过记录样本在噪声水平范围上的时域信号而实现最优化,噪声水平是在信号记录期间注入样本的噪声的功率增益。记录的信号随后被处理以显示频谱信号特征,并且选择具有最优频谱特征得分的时域信号,如下详述的。最优或接近最优的时域信号的选择是有用的,因为已经发现,并且根据本发明,用最优时域信号转换化学或生物学系统可产生比用非最优时域信号给出更强和更加可预测的响应。从另一种方式来看,当靶系统由样本信号转换时,选择最优(或接近最优)时域信号有益于获得可靠,可检测的样本效应。
一般而言,在其上典型地记录时域信号的注入噪声水平的范围是在约0到1伏特之间,典型地,或可选择地,注入的噪声优选超过要检测的所寻找的分子电磁辐射约30到35分贝之间,例如,在70-80-dbm范围内。被记录的样本数量,即,在其上时域信号被记录的噪声水平的区间的数量可以从10-100或更多变化,典型地,并且无论如何,在足够小的区间上,以便于识别最优信号。例如,噪声发生器水平的功率增益可在50 20mV区间上变化。如将要在下面看见的,当信号的频谱特征得分相对于注入噪声水平绘图时,曲线显示在几个不同噪声水平上延伸的波峰,当噪声水平增幅适当小时。
本发明思考了用于计算所记录的时域信号的频谱特征得分的三种不同方法。这些是(1)直方图箱方法(histogram bin method),(2)产生自相关信号的FFT,和(3)将FFT平均,以及这些方法中的每一个在下面详细叙述。
即使没有特别地描述,也会明白每一种方法可以以手工方式执行,其中用户估计频谱特征得分所基于的频谱,所述每种方法使得为下一次记录进行噪声水平调整,并测定当达到峰值得分时,或它可以以自动化的或半自动化的方式执行,其中噪声水平的连续增长和/或频谱特征得分的评估,由驱动计算机的程序执行。
A.产生频谱信息的直方图方法图20是在产生频谱信息的直方图方法中的高级数据流程图。从SQUID获取的数据(框2002)或存储的数据(框2004)被当作16bit WAV数据(框2006)保存,并被转换成双精度浮点数据(框2008)。转换后的数据可被保存(框2010)或显示为原波形(框2012)。转换后的数据随后传送到下述参考图21的并由框2014标记的傅立叶分析表明的算法。直方图可显示在2016。可选择地,以及将在下面描述,转换后的数据可传送到两个附加算法中的一个,用于识别在时域信号中的频谱特征。
参考图21,直方图算法的总流程图采用离散采样的时域信号并使用傅立叶分析将其转换成频域频谱,以便进一步分析。时域信号从ADC(模拟/数字转换器)获取,并存储在2102指示的缓存器中。样本长是SampleDuration秒,并以每秒SampleRate个样本进行采样,因而提供SampleCount(SampleDuration*SampleRate)个样本。如Nyquist定义的,可从信号恢复的FrquencyRange被定义为SampleRate的一半。因而,如果时序信号以每秒10000样本被采样,FrquencyRange为0Hz到5kHz。可使用的一种傅立叶算法是基数2的实数快速傅立叶变换(RFFT),它具有可选择的2直到216的功率频域分辨率(FFTSize)。只要FrquencyRange维持在8kHz或以下,选择8192个FFTSize,从而提供足够的分辨率以使每赫兹至少具有一频谱箱(frequency bin)。SampleDuration应该足够长以便SampleCount>(2*)FFTSize*10以确保可靠的结果。
因为FFT只能在某个时间作用于FFTSize个样本上,程序必须在样本上连续地执行FFT,并将结果一起平均以得到最终频谱。如果选择对于每个FFT略过(skip)FFTSize个样本,那么会引起1/FFTSize^0.5的统计误差。如果,然而,选择由一半FFTSize输入的FFT覆盖,那么这种误差将被减少到1/(0.81*2*FFTSize)^0.5。这使误差从0.0110485435减少至0.0086805556。关于误差和相关分析的附加信息一般而言,参见Bendat &Piersol,″Engineering Applications of Correlation and Spectral Analysis″,1993。
在所给窗口上执行FFT之前,数据锥度滤波器可应用于避免由于采样混叠导致的频谱遗漏。该滤波器可选自作为举例的Rectangular(没有滤波器),Hamming,Hanning,Bartlett,Blackman和Blackman/Harris之间。
在典型方法中,如在框2104中所示的,我们已经选择8192用于变量FFTSize,它是我们一次操作的时域样本数量,也是FFT输出的离散频率数量。注意FFTSize=8192是分辨率,或者是在由采样率指定的范围内的箱数。变量n,指示执行了多少离散RFFT’s(实数FFTs),其通过将SampleCount除以FFTSize*2,FFT箱数,来设定。为了进行产生切合实际的结果的算法,数n应该至少为10到20(尽管其它值是可能的),其中更多可能优选拾取较弱信号。这意味着对于所给SampleRate和FFTSize,SampleDuration必须足够长。计数器m,它从0数到n,被初始化为0,并且如框2104所示的。
程序首先确立了三个缓存器用于FFTSize直方图箱的缓存器2108,其将在每个箱频率累积计数;用于在每个箱频率的平均功率的缓存器2110,以及含有每个m的FFTSize个拷贝样本的缓存器2112。
程序初始化直方图和阵列(框2113),且在2114将波形数据的FFTSize个样本复制入缓存器2112,并对该波形数据(框2115)执行RFFT。FFT被归一化,使得最高振幅是1(框2116),且从归一化信号确定所有FFTSize个箱的平均功率(框2117)。对于每个箱频率,在这个频率的来自FFT的归一化值被添加到缓存器2108的每个箱(框2118)。
在框2119中程序随后察看在每个箱频率的功率,相对于从上面计算的平均功率。如果功率在平均功率的特定因数小正数(在0和1之间)内,那么它被计数并且相应的箱在16在直方图缓存器增加。否则它被丢弃。
注意进行比较的平均功率仅用于此FFT实例。增强的,即使较慢的算法可能采用两条途径通过数据并在设定直方图水平之前计算所有时间的平均值。与小正数的比较帮助表示对频率箱足够显著的功率值。或者在扩展的方面,采取小正数的方程式帮助回答这个问题,“在这个频率在这个时间有信号吗?”如果答案是有,它可能导致两种情况中的一种(1)仅这一次陷入这个箱的平稳噪声,或者(2)实数低水平周期信号,它近乎每次发生。因而,直方图计数将清除噪声冲击(hit),并且提高低水平信号冲击。因此,平均和小正数因数允许选择被认为显著的最小功率水平。
计数器m在框2120增值,并且上述过程被重复用于每个WAV数据的n集合直到m等于n(框2121)。在每个循环,在2118,每个箱的平均功率被添加到相关联的箱,并且当满足在2114的功率振幅条件时,每个直方图箱增加一。
当已经考虑过数据的所有n个周期时,每个箱里的平均功率通过将在每个箱中的总累积平均功率除以n来确定,显示周期的总数(框2122)和结果(框2123)。除了结构噪声存在的地方,例如,DC=0或在60Hz的倍数,每个箱的平均功率将是一些相对较低的数值。这被表示在图22A-D(在400,600,700,和900mV产生的直方图)显示的曲线中。图22A-22D的曲线只显示了一部分的直方图箱,即从7953Hz到8533Hz的频谱。如图22A和22B所示,没有分别在注入噪声的400mV或600mV可见的随机事件。然而,如图22C所示,在700mV,可见的随机事件是明显的。此后,如图22D所示,在900mV,随机事件消失了。
在每个箱中,由上述步骤产生的直方图包含,在那个频率的功率超过(小正数*整个FFT输出量的平均功率)的次数在0和n之间的计数。如果箱计数由于非结构化噪声增大,那么噪声将会随着时间的过去被分配在所有频率箱,因此在给定箱中总数不是很多。如果在所给频率有一致的信号,将存在于每个n时间段,并因此具有接近n的箱计数。大振幅的噪声,例如60赫兹和它的谐波都具有高箱计数以及高平均功率。我们能求这些频率,和我们所感兴趣的具有较低平均功率,但高箱计数的那些之间的微分。
图22A-22D显示了由在四种不同噪声功率输入的方法产生的直方图。如图所示,该程序可将在每个频率的平均功率显示为垂直线条。直方图箱计数可表示为连接的上部线。如果功率被认为是“低”(例如少于平均值/3),并且直方图具有特定的计数,那么连线可在功率线条峰值和直方图线条峰值之间变成可见的。由连线突出显示的箱可能作为低能量分子频谱的候选。
从图22A-22D和上述内容可以理解,在产生有意义的直方图中使用了两个值得注意的设定,即,直方图,其显示了与被询问的样本相关的随机共振效应。首先是供给样本的噪声功率水平(该情况下,高斯白噪声)。如果功率水平太低,那么噪声水平不足以产生随机共振并且箱直方图仅反映噪声。如果功率输入太高,对于每个箱计算的平均功率水平较高并且随机事件无法分辨。从这次研究中,最优噪声水平是在700mV左右,尽管真正的最优条件还可通过将该方法应用于信号来确定,该信号以在许多较小的在大约650和750mV之间的增量的信号增益来记录。
在该方法中的频谱特征得分通过计数随机事件的数量来确定,所述随机事件高于与这样的值相应的箱计数,该值在没有随机事件发生时是统计学上高于该箱的平均噪声的值。在图22A-22C所示的曲线中,该平均箱计数等于或略高于沿频谱轴线分布的显而易见的随机峰值,具体地,如图22A-22C中所示的。在最优噪声增益处(图22A),可观测到许多明显高于该水平的箱峰值,这些峰值可在选定的频率区间上计数,例如DC-1kHz或DC-8kHz,由此确定相应时域信号的频谱特征得分。
在该方法中的关键设置是噪声增益和小正数值。该值决定了将被用于区分平均值上的事件的功率值。在值为1处,由于功率不可能高于平均功率,因此没有事件被检测到。当小正数趋近0时,事实上每个值被放在箱中。在0和1之间,典型地在给出占结构化噪声总箱计数约20-50%的箱计数的值处,小正数具有最大“频谱特征(spectral character)”,这意味着随机共振事件将能够从纯噪声中最大受益。
因此,在实施本发明时,可以在噪声输入中系统地增加功率增益,例如在0和1V之间50mV的增量,和在每个功率设置,调整小正数直到观测到具有明确定义的峰值的直方图。在例如被处理样本表现为20秒的时间区间处,用于每一不同功率和小正数的全部处理时间将约为25秒。当观测到明确定义的信号时,可重新定义功率设置或小正数或其两者,直到产生最优直方图,这意味着产生了具有最大可分辨峰值数量的直方图。
在该算法中,在低频由于噪声的全面出现(如环境噪声),许多箱将被充填,与之相关的直方图表现为低频。从而,系统可以简单地忽略低于给定频率(例如低于1kHz)的箱,但是在更高的频率仍然表现出足够的箱值,以确定样本间独有的信号签名(signal signature)。
任选地,由于小正数变量的目的是适应在每个周期确定的不同平均功率水平,因此程序自身能够采用预先定义的使平均功率水平与小正数的最优值相关的函数来自动调整小正数。
类似地,在每一功率设置中,程序能够比较峰高,并自动调整噪声功率设置,直到在直方图中可观测到最优峰高或特征。
虽然小正数值对于所有频率可以是固定值,但其也可通过采用频率相关值(frequency-dependent value)来考虑,以调整用于可在低频,例如DC-1,000,观测到的更高值的平均能量。频率相关小正数因数可通过例如对大量低频FFT区域求平均,并确定小正数值来确定,该小正数值可调整平均值至与在更高频率观测到的那些可比较的值。
参考图23A-23C,显示了用于产生直方图的用户界面的例子。滑动条2302决定样本波形片段的长度,例如长达300-600秒,并可使用户在波形内有效拖动。框2304允许用户设定Nyquist频率,例如5、10或20kHz,并还提供临近的重置按钮。滑动条2306使用户可以移动直方图基线,而60Hz检验框2308使用户可以用垂线(如图23C中所示)识别60Hz箱和所有相关的60Hz的谐波。当选定获取按钮2312时,软件从样本中生成或获取波形,例如图23B中所示。当选定fft按钮2310时,软件生成直方图(histogramplot),如图23C中所示。
B.自相关信号的FFT在第二个用于确定频谱特征得分的常规方法中,在选定噪声记录的时域信号是自相关的,采用自相关信号的快速傅立叶变换(FFT)来生成频谱特征图,即,在该频率域中的信号图。然后使用FFT对在选定频率范围内,例如DC-1kHz或DC-8kHz,高于平均噪声水平的频谱信号数量评分。
图24是根据该第二个实施方式,在对记录的时域信号评分时,执行步骤的流程图。时域信号如上所述被采样、数字化和滤波(框402),其在噪声水平上的增益设置为初始水平,如在404。图25A显示了典型的样本化合物时域信号,在该情况下除草剂甘胺酸磷酸盐(RoundupR),在这里显示的片段在时间区间14.08-14.16秒上被采用。然后采用标准自相关算法对时域信号进行自相关,如在408,并采用标准FFT算法生成该自相关函数的FFT,如在410。
采用FFT图,例如图25B-25D中所示的,通过计数统计学上大于在自相关FFT中观测到的平均噪声的谱峰数量对该图评分,如在414。通过416和406逻辑关系,重复该过程,直至记录了峰值得分,即直到随着噪声增益增加,给定信号得分开始下降。在418记录下峰值得分,且程序或用户从422处的时域信号文档中选择对应于峰值得分的信号(框420)。
在图25B-25D中的一系列自相关FFT曲线阐明了该方法包括的信号分析。在70.9-dbm(图22B)噪声水平,观测到极少高于背景噪声的峰值(最高表示60周期噪声)。在74.8-dbm最优噪声水平(图25C和25D),其在相同噪声水平表示不同的记录,通过DC-8kHz频率范围,观测到大量统计学上高于平均噪声的峰值。这些峰值中的一些在78.3-dbm的更高噪声增益处不显著或者已经消失。
当这些信号的频谱特征得分作为噪声设置的函数绘制出来时,如图26所示,在约75-dbm的噪声设置处观测到峰值得分。从该图中,选择与峰值得分对应的时域信号。
如上,该实施方式可以以手工方式实施,其中用户手工调整噪声设定增量,手动分析(计数峰值)FFT频谱图,并采用峰值得分识别一个或更多最优时域信号。任选地,该步骤的一个或更多方面可以自动实施。
C.平均FFT在另一个确定谱峰值得分的实施方式中,在每个噪声增益处平均多个如10-20个时域信号的FFT,以生成谱峰图,并如上计算得分。
图27是根据该第三实施方式在对记录的时域信号评分中,执行的步骤的流程图。时域信号如上被采样、数字化和滤波(框424),其在噪声水平上的增益设置为初始水平,如在426。然后在每一噪声增益处,在428,程序生成一系列时域信号的FFT,并且这些图在430被平均化。采用平均后的FFT图,通过计数统计学上大于在平均后的FFT中观测到的平均噪声的谱峰的数量进行评分,如在432、434。通过436和437的逻辑,重复该操作,直到记录下峰值得分,即直到随着噪声增益增加,给定信号得分开始下降。在438记录峰值得分,并且程序或用户从442处的时域信号文档中选择对应于峰值得分的信号(框440)。
如上,该方法可以以手工、半自动或全自动模式实施。
D.其他谱分析方法出于识别制剂特异性谱峰的目的,这部分简要的考虑了用于时域信号谱分析的其它方法。如上所述,系统利用在随机共振试验中得到的声音文件作为输入,输出频率、振幅,和成分的正弦曲线(content sinusoid)的相位。系统可采用软件程序,称为“寻峰器(peakfinder)”,其依次采用其它软件包,例如Octave和Pd,前述两种都是开放源代码的(open-source)和目前得到支持的软件平台。
另外,可使用两个环境变量PF_TMP,其指定临时目录;和PF_BASE,其指定寻峰器(peakfinder)文件夹的位置。如果没有PF_BASE,peakfinder.sh脚本会试图根据自身调用推定它(假定其作为绝对路径名被调用)。输入文件为立体声声音文件,采用44100的标准采样率。文件格式可以是“wav”、“au”或“aiff”,为16、24或32比特样本帧。输出文件为指定一个正弦曲线的ASCII文件。例如

这里,第一字段是以基本分析频率为单位的频率,解释见下,第二字段是以赫兹为单位的频率,第三是输入声音文件本身单位中正弦曲线的峰值,第四和第五是正弦曲线的余弦和正弦分量的振幅,复数振幅的实数和虚数部分。当然,该值可由实分量和虚分量推出。第一字段没有物理意义,其用于调试目的。
在白噪声中用于确定单一正弦曲线振幅和频率的技术是极大似然(ML)方法,其已被扩展为多重正弦曲线。该方法假定正弦曲线数量预先已知。找到非预定的正弦曲线数量的问题在数学上很难处理,但可以通过假定所讨论的正弦曲线在频率上充分独立的方式来处理。此外,需要一种方法来区分正弦曲线的存在和不存在。
通过研究白噪声中的单一正弦曲线开始以下分析,并发展到分析多重正弦曲线和非白(如粉色)噪声问题。给定以下所测量的信号x[n],n=0,……,N,该(离散时间)非归一化傅立叶变换定义为FT{x[n]}(k)=Σn=0N-1e-2πink/Nx[n],]]>其中k是以该分析的基础频率为单位的频率;每样本2π/N弧度。k不必是整数;在实际操作中k的额外数值可根据需要用零-填充(zero-padding)信号填充。在假定存在单一正弦曲线时,其最可能频率由下式给出k=arg max|FT{x[n}(k)|.
换句话说,最佳估计值就是简单的k值,其能够使傅立叶变换量最大。然后,系统确定该k的估计值是否对应一个真实的正弦曲线或仅为随机波动。为此,需要分析无效假设以确定x[n]是否仅含白噪声,例如具有平均值0和RMS振幅σ。在每一个点k的傅立叶变换都是N个独立随机变量的和,每一个等于样本x[n]的单位量复数倍,因此傅立叶变换每一点的平均数仍然为0,标准偏差为 如果独立噪声样本的修减(tail behavior)能够被良好反映(例如对于高斯噪声或均衡噪声),得到的随机变量FT{x[n]}(k)对采用的N值(在106量级上),将非常接近高斯噪声。因此超出约 的可能性极小。
另一方面,具有峰值振幅a和频率k(通常以2π/N为单位)的实数值正弦曲线具有aN/2的傅立叶变换量。为得到 的量,我们只需要a至少为 该方法用零填充记录的声音文件(介于2和4的因数之间,取决于两个指数中的下一个(next power of two)),然后报告超过振幅阈值的峰值。峰值定义为与其周围相比,相对于给定k值具有更大量,并具有至少20个临近k值的量值的一半(约为20π/NHz,或1/3Hz频带,对一分钟的样本)。
如果存在几个正弦曲线,如果所提供的它们的频率彼此相差大于20π/N,上述方法应该可以分别分辨出它们;在计算出的傅立叶变换上,每个正弦曲线的影响在远离谱峰以k频率为单位的振幅上下降了2/3πk。
为补偿非白噪声信号,需要估算测量信号的频谱包络。噪声可假定为是在每一狭窄频率范围内(如上的20π/N)的局部白噪声,并根据该选定的频率范围逐渐改变σ值。另外一个问题是要确定注入的噪声样本能否从试验的测量的输出中减去。在这种情况下,通过简单与上述两者相关的可测量传递函数,即使是非线性的,也可使用该传递函数的估计从所测的信号中移除大量噪声。其也可提高方法的灵敏度。
E.复合信号如前文记载中所示的,该系统使用户可以创建可用于治疗疾病的波形,或者在生物学系统中诱发反应。可以得到由两种或更多化合物生成的波形或谱系。然后这两个信号可以被合并以形成单一信号,该合并的信号具有上述两个单一信号的特性。例如如果该两个源信号涉及两种具有不同治疗特性的不同化合物,那么其结果是,得到的合并信号应具有该两种化合物合并的治疗特性。然后可对该合并信号操作,从而移除已发现与生物学系统中副作用或负反应有关的频率分量。在已知感兴趣的频谱分量的振幅和频率条件下,该合并信号可很容易从公知的合成器技术生成。
任选地,如果该两种化合物在生物学系统中产生类似应答,那么可以比较从上述化合物生成的两种信号,以识别与产生生物学作用相关的共有频率分量。随后可以生成只包含与生物学作用有关的那些频率分量的第三信号。从而,可以比较某些缓解疼痛药物的信号,来识别其共有频率组分,随后生成用于传输、存储或应用在生物学系统中的最终信号。事实上,该系统允许构建出新的信号,其不直接基于从一种或多种化合物生成的信号。相反,该系统允许生成只在需要的频率具有峰值的信号,在该频率处这种峰值具有在生物学系统中需要的结果。因此,这种合成信号与已经存在的化合物无关。
IV.转换靶样本该部分记载了用优化低频时域信号转换样本的设备和方法(A部分),和在一些生物学样本上实施的转换试验(B部分)。用于转换的样本表现出被很好地表征,并对化学或生化制剂具有容易探测的应答,且该转换信号是该化学或生化制剂的优化时域信号。
A.转换装置和方案图28A显示了根据本发明的使用制剂特异性信号转换样本的仪器的设计图。该特殊设计可以容纳五种不同样本,包括保存在转换线圈内和暴露于电磁信号中的的三种样本444、446和448,用作控制的样本450,和用作化学诱导(chemical-induction)对照的样本452。如将要在图28C中看到的,在诱导期,典型地将样本在相同摇动、温度和湿度条件下,置于混合台上并保持在那里。
制剂特异性信号的转换通过“播放”样本的优化制剂特异性信号来实现,该信号记录于CD上,并通过前置放大器456和音频放大器458播放在CD上记录的454。如图所示,该信号通过不同通道提供给电磁线圈444和446。在一个实施方式中,使用Sony Model CDP CE375 CD播放器。播放器的通道1与Adcom Pre Amplifier Model GFP 750的CD输入1连接。通道2与Pre Amplifier Model GFP 750的CD输入2连接。CD可被记录为播放来自每个通道的同样的信号。任选地,CD也可记录为播放来自每个通道的不同的信号。样本448的线圈作为在每次试验中的控制端,主要用于产生白噪声场。例如GR模拟噪声发生器为该线圈提供白高斯噪声源。任选地,该线圈可用于经第二Crown放大器播放任何预记录的转换信号。
CD播放器和前置放大器间的电缆是标准RCA音频连接线缆(6英尺)。Adcom Pre Amplifier Model GFP 750接收从Sony CD播放器输出的信号,并充分放大以驱动Crown放大器。任选地,还可以配置前置放大器接收来自其他源的转换信号,例如来自Sound Blaster PC board。前置放大器和Crown放大器之间的电缆是标准RCA音频连接电缆,其在每条电缆(3英尺)的一端装配有到1/4电话插头的RCA用于连接Crown放大器。该Crown放大器型号Micro-Tech 2400,立体声系统(每条通道1000瓦)收到来自前置放大器的信号,并将信号水平充分提升至能够驱动转换线圈。连接Crown放大器和转换线圈的电缆为在每端带有香蕉(banana)插头的14型标准铜音频电缆。香蕉插头使用固定螺钉机械连接于电缆上。
化学诱导控制样本置于带有混合台的VWR紧凑培养器内,并维持其温度与转换线圈一致。
图28B显示了样本转换仪器466,例如由任意图28A中样本444、446和448表示的那些。该仪器包括装有电磁铁470的室468,和用于监测室内条件,例如温度的各种探测器。电磁铁位于基底474上,包括常规环形铁磁体芯和线绕组。
在一个实施方式中,线圈由American Magnetics设计和生产,以提供线圈间均衡的性能。每个线圈由416转#8规格的正方形铜磁体线组成,用搪瓷覆盖,并有约2″的空气芯。每个线圈在11赫兹在10安培RMS在10瓦RMS且温度上升不超过15摄氏度下,可在中心产生约1500高斯。无诱导对照组。将无诱导对照组样本置于带有混合台的VWR紧凑培养器内,并维持其温度与转换线圈一致。
线圈由两英寸高四英寸宽的PVC支撑体474支撑。将一系列1 7/8英寸的OD PVC管切成不同的垂直长度,以使得能够将样本容易的定位在每个线圈中央。样本由线圈顶部插入,并停留在PVC定位管(未示出)中。线圈和基底用快速干燥的环氧树脂固定在外壳上。线圈柱和RF输入连接器间的连接通过用铝钉、螺母和垫片将4英寸长14规格铜标准线机械连接于线圈而形成。其他端用60%焊料焊接于RF输入连接器。RF探测器由焊接在插入式(male)BNC连接器绝缘中柱的单一6英寸长12规格铜线构成。RF探测器在每一端均使用由BNC连接器安装的RG 6同轴电缆,连接在Stanford Research Systems的型号SR 785 2的通道动态信号分析器上。
使用Sensatronics型号E4的温度监视器监视全部三个线圈和培养器的温度。将传感器分接于每个屏蔽外壳的外壁和培养器内壁。(探测器为Sensatronics Stanford的型号的温度探测器)探测器电缆也由Sensatronics提供,以匹配探测器和监视器。
图28C显示了混合台的布局,其中三个样本室,例如室466,由单独混合台支持,例如台486,上述全部设备均由支持台488支撑,该支持台位于屏蔽室或外壳490内,在转换试验期间维持屏蔽室或外壳为恒定的温度和湿度。
该屏蔽外壳的例子是大小为10×10×10英寸,由125英寸6061 T6铝板制成,例如由Rowe Air Manufacturing,Marysville,Washington生产的。
在图29A阐述的第二个常规实施方式中,转换线圈为Helmholz线圈。理想的Helmholz线圈由具有相同半径的两个共轴循环电流回路构成,彼此分离一个半径。也就是说,该回路为I分离的,且I=r。由线圈产生并以特斯拉为单位的在Helmholtz线圈轴线上任意点产生的磁场Bx由下面的方程给出,其中磁场方向与回路平面垂直。
Bx=μ0i2r[1(r2+r+54)3/2+1(r2-r+54)3/2]]]>在该方程中,m0是磁导率常数(1.26×10-6H/m),i为导线中的电流,单位为安培,r为电流回路的半径,单位为米,g是x/r的比值,其中x为从Helmholtz线圈中心的轴线距离,r为线圈半径。
如图29A所示,由每个线圈产生的磁场叠加,从而在线圈中央得到了相对大的均匀磁场。线圈也可以为正方形或矩形,其中根据公知原理,线圈分离以在线圈间形成均匀磁场。
图29B和29C显示了适用于本发明的可选择转换线圈。图29B中的转换器494是长螺线管,例如长达数英尺。螺线管内的磁场平行于螺线管轴线,并在螺线管内恒定,在螺线管外趋于零(在接近无限长的螺线管内)。该有限长的线圈只有在其中心附近才具有充分均匀的磁场。因此,通过将样本放置于线圈的中央,当用MIDS信号为线圈供能时,在样本上可以形成基本均匀的磁场。
通过向螺线管上增加额外匝数,例如在图29C的螺线管496上增加额外的匝数500,可增加线圈末端的额外磁场强度,以补偿在其末端的线圈磁场下降。
在另一个其他实施方式中,转换线圈可以是很小的可植入铁磁体线圈,在这种情况下,可以是一种脉管内限制线圈(vascular stent coil),其可通过附着于线圈相对两端的电极或通过远程感应系统接收转换信号,在所述远程感应系统中,电磁体被放置在与患者胸腔相对的体表附近,且该信号被诱导传输到植入线圈上。
B.转换研究以下B1-B6部分是对用于阐述本发明的方法的样本/制剂系统的详细说明,该系统包括(1)具有可由L-(+)阿拉伯糖(+)诱导的lac操纵子的阿拉伯糖诱导型细菌系统;(2)茎长增长可由于除草剂甘胺酸磷酸盐的存在而受到抑制的甜豌豆植株;(3)茎长增长可通过使用植株激素赤霉酸受到刺激的甜豌豆植株;(4)涉及蛋白体活性的phepropetin D抑制,和在没有信号的情况下活性恢复的体外系统;(5)在小鼠动物模型系统中,通过暴露于紫杉醇治疗癌症;和(6)在小鼠动物模型系统中,通过暴露于肿瘤抑制剂肽p53的癌症治疗。
获得各样本材料的MIDS信号的详细方法以及以MIDS信号转换生物学系统的详细方法分别详细的记载于下面的例1和2中。
B1.通过L(+)阿拉伯糖MIDS诱导细菌lac操纵子细菌中的阿拉伯糖操纵子是严谨调节的系统,由基因产物控制下的启动子、araC及其同类诱导物、L-阿拉伯糖、糖组成。异源蛋白质表达可受具有阿拉伯糖和araC基因的质粒的严密控制,其既是启动子的正调节剂也是负调节剂。诱导机制涉及在连接阿拉伯糖之后,araC基因与DNA元素的结合特性的变构变化。
用于本研究的化学传感器系统已知对L(+)阿拉伯糖形式具有选择性,而D(-)阿拉伯糖不能诱导araC-PBAD启动子依赖型诱导和报道分子(reporter)GFP uv蛋白质的表达,尽管它们具有相同经验式和类似的结构。因此,该系统和相关的化合物提供了一组理想的工具,用以评估经再现恢复的这些类似的分子特征的分子发射信号(MIDS)而转换的任何基因表达的特异性。
如以下例1中所详述的那样产生L(+)阿拉伯糖和D(-)阿拉伯糖的最优信号。简要地,通过改变注入的噪声水平和用于信号恢复的材料的物理浓度,优化用于随机共振产物的信号。衍生自等摩尔的阿拉伯糖异构体溶液的测量结果的优化SR信号产物的比较存在显著区别。恢复的MIDS信号被认为存在显著区别,并对于两种不同的阿拉伯糖异构体而言是截然不同的。这些结果表明,恢复的分子发射信号对这两种具有相同化学组成和经验式的分子的化学结构很敏感。这两种异构体的区别仅在于与己糖环连接的四个羟基的顺-反排列不同。
利用例2中详述的转换装置和方案,将细菌样本暴露于L(+)阿拉伯糖和D(+)阿拉伯糖的MIDS信号。简要地,在一个典型的试验中,用1-3滴pBAD-GFPuv转化的JM109过夜培养物接种50ml新鲜的LB/amp,所述JM109被培养于含100-500μg/ml的氨卞青霉素钠的Luria肉汤中。将其充分混合,将其中10ml转移至25ml的三角瓶中,并添加0.1ml参比化学诱导剂(D或L阿拉伯糖)的100×原液至所需终浓度。将这些接种的培养物放置在转换线圈内,或置于振动台上用作对照。多数情况下,转换需要通过过夜培养进行。
将五个独立实验的OD测量数据绘制于图30A-30E中。图30A中的图表示GFP应答于MIDS信号和化学制剂本身的诱导。如所见,MIDS L(+)阿拉伯糖信号的诱导效果,和0.2mM L(+)阿拉伯糖相当,并远远大于对D(-)阿拉伯糖MIDS信号或D(-)阿拉伯糖化合物的响应。在不存在MIDS信号或阿拉伯糖化合物时,只可观察到很小的诱导。图30B的数据表明对这两种MIDS信号有类似的效果。在这头两个试验中,只使用了两个线圈。非MIDS控制器独立地置于一个培养器内而不摇动,将该培养器的温度与线圈温度进行手工匹配。将线圈在使用的电流和电压下电阻加热。从而在L(+)和D(-)MIDS样本间采用相关的差值测量。该对照仅表明了该系统正确的化学应答,不能用于定量。
绘制于图30C-30F中的三个试验均显示使用L(+)MIDS信号的诱导远远大于使用D(-)MIDS信号或白噪声控制。这些数据证明了从同源化学诱导剂L(+)-阿拉伯糖而非无活性的化学异构体D(-)阿拉伯糖恢复的MIDS信号对araC-PBAD操纵子的特异性转换。由于仅在立体化学上不同,参比化学制剂在组成上是相同的,因此该系统表现出是转换特异性的良好模型。显然,该生物学系统对结构具有选择性。在进行的试验中,连续过夜再现得到了非常一致的15-20%的L相对于D MIDS信号的相对诱导。在一个试验中,在Surestha等的化学诱导方案中所使用的3.5小时的较短的诱导期,得到了7.8%的明显的相对诱导,接近于背景水平。
在实验条件下,生长和诱导均需相当长的时间。这些数据与MIDS再现可影响操纵子开启通路,以实现基因的特异性诱导的观念相一致。可能的作用位点可能位于与araC蛋白质相互作用的L(+)阿拉伯糖MIDS中,以诱导其模拟L(+)阿拉伯糖结合状态,然后可导致基因诱导。这一模型不需要形成任何化学键,但可能变构地影响该系统。
B2.通过甘胺酸磷酸盐MIDS抑制茎细胞生长单聚体植株酶5-烯醇式丙酮莽草酸酯-3-磷酸(EPSP)合酶(EPSPS)是烯醇式丙酮转移酶类中两种酶中的一种。配体结合使酶依据诱导契合(induced-fit)机制,由开放状态转化为紧缩的闭合状态。EPSP合酶参与到莽草酸途径中,通过利用磷酸烯醇式丙酮酸(PEP)将莽草酸-3-磷酸(S3P)转化为5-乙醇式丙酮-3-莽草酸磷酸,其是在细胞中产生的大量芳香族化合物包括芳香族氨基酸的前体。已经报道在该途径中产生的化合物构成了多达35%或更多的植株干重。在哺乳动物、鱼、爬行动物、鸟和昆虫中不存在这种酶的事实,使之成为抗生素和除草剂的良好靶。
合成化合物草甘膦(glyphosate、N-phosphonomethylglycine),除草剂Roudup中的活性成分作为该酶的竞争性抑制剂,有效地关闭了芳香族氨基酸生物合成以及衍生于这些氨基酸的其他芳香族化合物的合成。甘胺酸磷酸盐是一种结合于PEP-结合位点的过渡态类似物(磷酸盐和甲酸盐离子模拟PEP含氧碳正离子的活性基团)。草甘膦呈现出对EPSP合酶的高特异性,不完全结合于UDP-N-乙酰氨基葡萄糖烯醇式丙酮转移酶(MurA),烯醇丙酮转移酶类中仅有的另一种酶。
采用以下方法部分中详述的系统仪器和方案,如上述部分II和III中记载的那样获得甘胺酸磷酸盐的优化的固有发射信号(MIDS)。在一个实验中,将优化的甘胺酸磷酸盐MIDS信号施加于潮湿介质支持的豌豆芽,在37℃,100%湿度下施加超过72小时的信号。各种控制包括(1)无信号,(2)化学Roundup,和(3)白噪声信号。暴露给刺激物1小时后测量茎长,每组中含有5-50个萌芽。
图31显示了四个实验组的每一组的茎长的条形图,其中每个条的上边界线和下边界线分别代表茎长的最大值和最小值,中间线代表该组茎长的平均值。如图所示,无论化学甘胺酸磷酸盐还是甘胺酸磷酸盐MIDS信号均有效地抑制了植株中茎的生长,然而白噪声控制对茎生长的影响很小或没有影响。
图32A和32B提供了更多甘胺酸磷酸盐MIDS信号对豌豆芽样本的影响的图形证据。图32A是在没有任何化学或MIDS信号处理的条件下,生长72小时的豌豆芽的图。相同生长期的豌豆芽,但是暴露给刺激物1小时后,显示出明显的受干扰和异常生长的迹象。
B3.通过赤霉素MIDS刺激茎细胞生长时序电磁信号模拟特异性生物活性分子的生物活性能力也已经利用赤霉酸-3(GA-3)在甜豌豆芽上得到证明。赤霉酸(赤霉酸钾,也称为MegaGrow)是一种天然存在的植株生长调节剂,其可导致多种效应包括在某些情况下刺激种子萌芽。GA-3天然存在于许多物种种子中,可通过在桶中培养赤酶菌(Gibberella Fujikuroi)真菌培养物,然后提取和提纯GA-3进行商业化生产。在GA-3溶液中预浸渍的种子可以在很多情况下,使得多种类型的高度休眠的种子快速萌芽,否则这些种子将需要冷处理,后熟或老化,或其他延长的预处理。
将大量生长调节剂赤霉酸(.001%)放置在1.5cc Pyrex样本瓶中,并临近第一衍生物超导磁场梯度计(derivative super conducting gradiometer)放置。频率范围为DC-2000Hz的白高斯噪声电流经16匝Helmholtz线圈传播,产生与梯度计轴线正交的B-场(B-Field)定向的磁场。假定该噪声信号促进了用梯度计测量的赤霉酸的分子电磁(磁)辐射的随机信号特征的形成。如上述,以WAV格式记录和存储该随机信号。给每个信号提供独有的控制数字。
对上述记录的生长调节剂的信号进行后处理,以识别只属于激动剂的频谱,并通过使用IFFT转化为WAV文件。创建了五个研究组,每组均由从Sun Grown Inc.,San Diego获得的活体甜豌豆芽组成。以促进甜豌豆萌芽的方式将豆芽放在湿纸巾上准备各研究组。
这五组是(1)无影响对照组(在适当的光、水和通风条件下自然发芽);(2)生长调节剂GA-3组(用生长调节剂处理萌芽);(3)经噪声线圈施加白高斯白噪声;(4)生长调节剂辐射组(对萌芽应用从步骤1记录的生长调节剂);和(5)生长调节剂辐射组(噪声移除-IFFT,即将记录的生长调节剂信号进行后处理以仅包括激动剂独有的频谱,然后再显萌芽)。
信号源为步骤1中产生的WAV文件,将其转化为模拟信号,并使用配置的低频RF放大器放大,以向手性缠绕(hand wound)的4英寸(ID)特斯拉线圈提供具有足够强度和质量的信号;将组3、4和5安置在转换器上方。采用适当的电磁屏蔽使各组彼此隔离。提供适当的光、热量和通风,以促进所有组健康生长。除非更早地得到动态结果,否则每组研究的持续时间为72小时。转换器是18英寸直径的改进型Helmholz,使萌芽位于线圈核心的中央。
研究结果绘制于图33中,表明了上述1-4组的平均茎长。如图所示,相对于自由生长对照组(组1)、白噪声(组3),甚至赤霉酸(组2),赤霉酸的MIDS信号(组4)的平均茎长显著增加。
B4.通过phepropeptin D MIDS抑制蛋白体的活性使用温和去垢剂(SDS)激活20S蛋白体酶。通过被激活的酶对底物Suc-Leu-Leu-Val-Tyr-AMC的剪切释放出荧光AMC分子,其可通过在380nm激发AMC,并在460nm检测荧光对该分子进行监测。蛋白体活性被phepropeptin D抑制。
在第一个研究中,按照上面以及实施例中记载的完全相同的方法,测定了phepropeptin D及其抑制蛋白体活性的MIDS信号的能力,其通过测定底物Suc-Leu-Leu-Val-Tyr-AMC的转化出现的荧光信号测定。
20S蛋白体检测试剂盒(EMD Biosciences Cat#539158)适于使用终容积200μl的试管类型。采用例3中详细记载的方法建立9个样本。阳性对照样本#2设定为700,其他所有测量均相对于该样本进行。化学抑制剂将20S蛋白体的活性抑制到小于未受抑制的活性的20%。当处理1.5小时后进行检测时,衍生自phepropeptin D样本的MIDS场使20S蛋白体酶活性下降至小于36%。
在后续研究中,将样本在室温下过夜培养(20小时),以确定是否底物会继续转化为荧光产物。如果MIDS信号抑制20S蛋白体且所述酶从场中移除后仍然保留活性,则相对于其他含有与化学抑制剂的样本及无抑制剂的对照,暴露给MIDS的样本的活性增加了。
下表中给出的结果表明,如通过对荧光计的灵敏水平下降所指示的那样,阳性对照持续上升。相对于持续增加的阳性对照,化学抑制剂的活性变得更加明显。MIDS样本是唯一一个在20小时后活性增加的样本。后处理信号表明1小时30分钟培养后最初探测到的抑制,在将样本从场中移除后不再存在。下表中给出的结果绘制于图34中。
%归一化的荧光1小时30分钟 20小时后阴性对照 0 0阳性对照 100 100ALLN 4.97479 2.645832PheD 18.1260516.97173PheD+MIDS 34.6386627.7141MIDS 38.4537866.3626
B5.利用紫杉醇MIDS治疗癌症该研究在人乳腺肿瘤的标准小鼠异种移植模型中,评估了紫杉醇转换分子发射信号的生长抑制潜力。紫杉醇是一种已在临床上获证的细胞毒性制剂,其通过与由聚合微管蛋白组成的细胞骨架元素的微管蛋白亚单位非共价结合以避免其分解,从而阻止细胞分裂并引发编程性细胞死亡而特异性地发挥作用。因此,紫杉醇对最初结合之后具有增加效应的蛋白质单体具有别构效应,从而影响较大的大分子结构。细胞毒素效应在快速分裂的癌细胞中更加明显。紫杉醇的作用机制在与靶微管蛋白分子结合后,不需要形成或破坏共价键。试验设计和条件详述于例4中。
参考在图35A中绘制的数据,单个动物肿瘤大小的对比显示,在研究结束后(第36天)肿瘤生长在统计学上被显著抑制。在统计分析前,将从所有组中移除了超出标准偏差的两倍的范围的值定义为异常值(outlier)。在第36天,紫杉醇、40mG和60mG治疗分别抑制肿瘤生长43%、36%和38%。将40mG和60mG治疗与标准的化疗剂紫杉醇对比,所揭示的紫杉醇在这些治疗间没有差别,这表明在移除异常值后,40mG和60mG治疗表现出与紫杉醇治疗类似的功效。
在第36天的肿瘤重量的对比揭示,在紫杉醇对照组和40mG治疗组中,肿瘤重量也同样在统计上被明显抑制。紫杉醇、40mG和60mG分别抑制生长43%、36%和26%(图35B)。总之,40mG和60mG治疗均在MDA-MB-435人乳腺肿瘤模型中显示出类似的功效,而没有有害临床体征。36天后,将肿瘤大小和肿瘤重量与载体对照比较,两种治疗均获得了统计学显著性。然而,在统计分析中包括来自60mG治疗的两个异常值和来自40mG治疗的一个异常值,将会消除其统计学显著性。40mG和60mG治疗与紫杉醇治疗的对比,均显示没有显著区别,这说明治疗36天后具有类似的功效。
B6.通过p53 MIDS抑制肿瘤生长采用MDA-MB-468人乳腺癌细胞系(line)测试Ant-p53pep的作用,所述细胞系含有内源p53-273His突变基因。该突变已作为p53基因中两个最常见突变点之一被报道,占高达所有p53突变的18%。Ant-p53pep对细胞成活力的影响通过台盼蓝排除进行初步检测。细胞成活力计数表明,在MDA-MB-468细胞中,具有30μm Ant-p53pep的处理5小时后活细胞减少了10倍。相反,用30μM肽处理5小时后,在良性wt p53乳腺上皮细胞系,MCF10-2A的成活力没有任何影响。
编程性(apoptotic)细胞死亡通过Annexin V与磷脂酰丝氨酸的结合被进一步证实。在许多细胞类型中,由于编程性细胞死亡的早期的膜不对称性缺失,磷脂酰丝氨酸被暴露给细胞外环境。用30μM Ant-p53pep经各种时间处理MDA-MB-468细胞,然后用Annexin V-FITC培养,并用流式细胞仪进行分析。结果表明与未处理对照相比,约35%的细胞在暴露给Ant-p53pep 30分钟后与Annexin V结合。
也在MDA-MB-231和MCF7细胞中观察到Annexin V结合,但在27sk细胞中没有观察到。MDA-MB-468细胞表现出特有的与编程性细胞死亡有关的形态学变化,包括膜出泡、细胞皱缩和核碎裂。肽处理的MCF 10-2A细胞的形态与对照组细胞类似。此外,用Ant-p53pep处理p53无效/无效乳腺癌细胞系、MDA-MB-453,没有表现出sub-G1DNA含量或Annexin V结合的积聚。其他两个具有无效p53情形的细胞系,H1299(肺腺癌系)和Saos2(骨肉瘤系)也没有被Ant-p53pep诱导至编程性细胞死亡。
此外,用30μM Ant-p53pep处理MDA-MB-468细胞2小时,通过蛋白质印迹(western blof)产生了FLICE(p26)特有的截短的活性形式。这些数据说明,在Ant-p53pep处理的MDA-MB-468细胞中,编程性细胞死亡经Fas/APO-1信令途径调解,并导致FLICE活化。
在MIDS研究中,将培养的MCF10-2A细胞暴露于记录的用于Ant-p53肽溶液的MIDS信号。于60mG,将细胞暴露给1小时开1小时关的MIDS信号,共处理18小时。处理结束时,检查培养盘中的细胞数。28小时后观察到的细胞数相对于对照组(没有MIDS信号)观察到的细胞数显著下降(约为对照组细胞的50%)。
C.MIDS制剂和化学或生化靶的的选择上面记载的方法证明,当已知在细胞环境下,制剂本身与抗配体靶的非共价作用情况下,本发明对于很宽范围的配体生化制剂MIDS信号是可操作的。本发明的研究中的配体/抗配体靶的对包括与细菌系统中的lac操纵子相互作用的L(+)阿拉伯糖;均结合并抑制它们各自靶酶的两种酶抑制剂(甘胺酸磷酸盐和phepropeptin D);与染色体微管蛋白结构结合以抑制自组装微管蛋白结构的解离的紫杉醇;和肿瘤抑制剂肽p53,其与DNA以及可能的蛋白质靶的结合位点相互作用以诱发参与阻止细胞分裂的因子(例如p21)的表达,或涉及编程性细胞死亡的多种蛋白质(包括bax、Fas、FasL和DR5)的表达。
在这些情况中的每一情况中,实际配体通过非共价作用与其抗配体靶的相互作用,非共价作用可包括离子、疏水、分散,和/或氢键。基于这些观察结果,可以合理地推测到,其它系统也对配体MIDS信号显示出类似的应答,在所述其它系统中MIDS信号代表其生物学活动涉及与生物靶的(例如蛋白质和/或DNA靶的)非共价作用的配体。
到目前为止,已经研究的一种配体/生物系统没有表现出统计学上有意义的MIDS信号应答。该情况涉及癌细胞应答系统,在该系统中MIDS信号是对于肿瘤坏死因子相关性编程性细胞死亡诱导的配体(TRAIL)蛋白质记录的低频时域信号。已知该蛋白质通过在细胞中诱导编程性细胞死亡途径抑制某些类型的癌细胞。该蛋白质的活性形式是一种具有配位锌(Zn)原子的三聚体,所述锌原子表现为参与了配体-受体相互作用,例如TRAIL蛋白与DR4和DR5细胞受体结合。众所周知,无论TRAIL化合物(具有配位锌原子的三聚物,且总分子量约为60K)多么复杂,MIDS信号不能反映所述三聚物或结合锌特征或者溶液中的TRAIL成分不能形成必需的配体复合物,均是导致不能观察到MIDS效应的原因。虽然如此,其结果对于选择用于实施所述方法的适当的效应物化合物具有指导意义,这将在下面研究。然而,TRAIL MIDS信号的结果确实证明了,将生物样本暴露于无活性MIDS信号的过程本身影响了系统的表现,因为在细胞暴露于该MIDS信号的过程中没有观察到癌细胞生长的有意义的变化。
已经研究了涉及共价键形成的一个配体/应答系统,但是没有明确的结果。该系统涉及将MIDS信号应用于抗细菌靶系统的β-内酰胺抗生素氨卞青霉素。已知β-内酰胺抗生素通过对细菌转肽酶中的丝氨酸残基共价修饰发挥至少部分作用。由于期望MIDS信号包括远低于共价键形成或破坏所需能量的相互作用的能量,因此如果其破坏或形成共价键所需机制不能在靶中引起制剂特异性反应的MIDS信号将是不被期望的。
基于以上讨论的结果,下列指导原则将会有助于选择适于用作影响生物学系统的的MIDS信号的配体效应物化合物。
1.配体或效应物化合物应该通过非共价相互作用(例如,离子,疏水,分散,和/或氢键相互作用)与生物学靶相互作用;2.配体或效应物化合物可以是一个相对较大的生物分子(如p53抑癌蛋白),但是产生MIDS信号的分子或复合物的溶液形式应为其生物学活性形式。因此,如果该生物学活性形式包括二聚体或三聚体结构,或配位金属键,则当记录MIDS信号时这些活性形式必须占主要地位。代表性的化合物是相对较小的有机分子或肽,例如,分子量小于1kD;3.特别优选的化合物是那些克服已知的生物系统中化合物自身限制的化合物。特别地,如果化合物具有较差的药物动力学参数(例如,由于在靶肿瘤细胞中的多重药物抗性,或诱导P450药物代谢途径,或不能穿过血-脑屏障),则该化合物的MIDS信号输送可有效地克服这些限制;4.相似地,对于由于与化合物在体内蓄积或与药物分解成毒性代谢产物相关的非特异性毒性效应而具有不需要的或禁止的副作用的化合物,MIDS信号可望成为有益的替代。
一旦选择了制剂或样本,则会优选在与靶的相容的介质(例如水介质)中,以及在与靶系统中有效的浓度相似或更高的浓度下产生最优MIDS信号。MIDS信号随后被用于转换靶系统,在一个选定水平(例如,在40-60mG之间)播放信号一段足以产生制剂特异性效应的时间,例如24小时。观察到制剂特异性效应时,MIDS信号被识别为适用于本发明的转换方法。如果没有观察到制剂特异性效应,则可进一步优化该MIDS信号,例如,通过产生仅含有所述样本最显著的频谱分量的合成信号。
在识别适合的化合物或样本之后,可通过将靶系统(例如,体内或体外生物系统)以多个时间间隔中的每一个(例如,持续暴露或在一定时期内(典型地在12小时至一周或更多的时间内)交替暴露)暴露给处于多个信号水平(例如,10至100mG之间)中的每一个的MIDS信号,使所述化合物的(或样本的)MIDS信号效应最优化,,以确立最可行的处理剂量和方案。
以下实例阐明了本发明的各种方法,但决非对本发明的范围的限制。
例1 MIDS信号获取方法测量时,样本材料处于70至74华氏度之间的室温。样本温度的大范围变化可以使辐射的频率根据温度的升高和降低变得更高或更低。如果样本在高于或低于室温时被记录,应当记录测量前和测量后材料的温度,并注意加热或冷却的方法。
如果待测量的材料暴露于室内空气,应当记录测量时的大气湿度。如果待测量的材料暴露于室内空气,应当记录测量时大气压。如果样本被加压,压力以磅/平方英寸或英寸汞柱记录。
进行分子测量前后应当记录周围环境的电磁干扰。在后处理中,该数据将从样本数据中减去。
A.样本制备样本大小0.8 to 1.5cc样本容器1cc或2cc平底三角瓶。
样本存储样本被安全保存以用于不同条件下的多次测量。将样本储存在原先的三角样本瓶中,用无活性的螺纹瓶盖盖紧。
B.SQUID(超导量子干涉设备)参数设置冷却使所有低温成分都达4开尔文的工作温度。
启动SQUID和SQUID控制器被启动并完成内部检测以确保正确操作。
调谐将SQUID调至最优设置(0.4-0.8微伏输出电压,无白噪声注入)。
增益100XDC偏移量(offset)ZERO带宽标准(50kHz)滤波器3Hz(高通)C.噪声发生器信号类型模拟高斯白噪声,或者均匀噪声(等幅)输出电压从最小值到最大值开路输出电压3伏rms最小。施加的输出电压最初设置为零输出电压。
X通道连接至Helmholtz线圈用于噪声注入。
Y通道(反向的)连接至与梯度计串联的噪声消除线圈,在SQUID输入线圈之前(可选的)。
输出阻抗50ohmsD.谱分析通过PCI数据采集板获取从MIDS输出的模拟信号(分子查询和数据系统)并且存储。WavGrab也被设计成与Noise Com型号UFX 9837的白噪声发生器串行连接。根据每一个后来的记录,当施加于Helm Holtz线圈的噪声的量增加时,样本被多次记录。一系列的记录从而形成不同噪声水平的单一靶样本。
可替换地,可以使用用于信号获取的人工方法。
E.随机生成将样本台从SQUID检测器中移走(没有样本)。将小量增益施加给噪声发生器通道X。将噪声发生器通道Y的输出调整,以在SQUID输出端产生可能的最深的0值(deepest null possible)(任选白噪声消除)。噪声发生器的主增益控制减小至零。通过将样本槽滑至梯度计下的位置将样本嵌入检测器。监测平均傅立叶显示并注意基线的频谱(在感兴趣的带宽内)。当监测基线频谱中的改变的平均傅立叶显示时,主增益控制(噪声发生器)被进一步改进。
当出现的明显的谱峰不是典型的基线频谱时,中断噪声增益的前进(advancement)并记录一段时序的辐射记录。
F.后处理采集单分子靶的一些原始时序记录并以WAV格式存储。每一个记录代表在一定白噪声振幅范围上的一系列记录测量值之一。标记从单分子靶获得的所有原始时序记录用于批处理。该批中的每一个原始时序记录是自相关的。
自相关函数被用于以下两个目的A.检测数据中的非随机性。
B.如果数据不是随机的,确定一个合适的时序模式。
将自相关时序与批量原始时序记录存储在一起用于进一步处理。该批中的每个自相关时序通过使用傅立叶变换被转换至频域。对于每个频域变换,计算Y轴穿过X轴(频率)中的所有数据点的RMS平均值(自相关得分)。对于每个频域变换,所有超出RMS平均值的Y轴异常值被列成表显示。列成表的批量中每一个时序记录的异常值总量被写成显示记录名称以及异常值计数的展开图(sheet)。最高的异常值得分的时序记录代表最佳的记录信号。
例2转换方法和条件转换研究通过两组线圈完成。在第一个系统中,一对线圈包括由30规格的缠绕在25/8″外径的PVC芯的3″高的切面周围的铝线组成,缠绕成35X启动线圈的绕组的数量。这些线圈被放置在铝制法拉第箱(faraday box)5 1/2″×4″×8 3/4″(W×H×D)内以隔离用于交叉RF信号的样本。将两个箱放置在定轨振动台上通风转换培养。样本和线圈的温度用Craftsman型号82327非接触型红外摄氏温度计监测。(实验时,配对的RF白噪声控制是不适用的)。
在三个以1″内径围起的28G绝缘铜的4″高的stock 4000-ft线轴的第二个设置中,塑料芯被直接用作转换线圈(″Mega线圈″)并且当第三个控制线圈被随机白噪声发生器驱动的时候,两个线圈被MIDS信号驱动。在这种情况中,全部三个线圈均在法拉第箱中隔离,并且通过在适于定轨振动器的6-ft板上隔离2ft独立安装。这确保全部三个样本全都用一个振动器以同样的速度摇动。将热敏电阻温度探测器也放置在每个最靠近样本的线圈的中心以实时读取温度。两组MIDS信号线圈设置均在最高输入电压设置时被驱动,直至它们的温度自缠绕的金属线电阻发热平衡于94-97附近。将样本放置于线圈的中心以经历最大以及最均匀的RF通量(电流设置的限制与采用控制RF通量的驱动电压加热相关)。
实施例3通过phepropeptin D MIDS抑制蛋白体活性20S蛋白体酶检测试剂盒(EMD Biosciences Cat# 539158)适合于在200μl终体积规格的试管中使用。使用下列步骤准备九个样本1)稀释反应缓冲液(20×),将100μl用880μl的HPLC级水稀释,并加入20μl的100×激活缓冲液,稀释至终体积1ml。
2)将100μl经稀释的2×反应缓冲液加入每个无酶的对照试管。
3)将20S蛋白体酶加入剩余的溶液中,将3.2μl加入至800μl的2×反应缓冲液中,分别取100μl加入剩下的7个试管。
4)向每个试管加入100μl水或水与抑制剂。
5)在一个独立的试管中,用96μl HPLC级水稀释4μl的底物,并将10μl/试管的稀释底物加入到每一个反应试管中至终体积200μl。
6)通过反复移液混合样本,将反应物在37℃培养2小时。暴露于MIDS场的样本在helmholtz线圈中培养,部分淹没在的传播60mG ACPhepropeptin D场的37℃水浴中。
7)通过荧光计(Turner Designs,Sunnyvale,CA-Model TD-700)量化并记录荧光。。
8)在图上以星号表示具有两个重复采样的样本。
所述样本为124.91(Disc) 无酶对照组(阴性对照)2697.9(Disc) 无抑制剂(阳性对照)325.83(DiSc) 化学抑制剂ALLN4105.6(Disc) 化学抑制剂Phepropeptin D579.35(Disc) 化学抑制剂Phepropeptin D+MIDS PhepropeptinD(60mG)6251.1(Disc) MIDS Phepropeptin D(60mG)7519.3(Disc) 无抑制剂(阳性对照)821.96(Disc) 无酶对照组(阴性对照)9209.7(Disc) MIDS Phepropeptin D(60mG)该研究的结果在前面已经述及。
实施例4通过紫杉醇MIDS降低肿瘤生长在玻璃样本架(Kimble自动注射瓶)上测量MIDS,通过(1)预制的临床用紫杉醇药水瓶中的纯净溶液,使用一个1ml的结核菌素注射器,以及22G注射针头,从多计量瓶中抽出0.7ml(无菌保藏并完全溶解,不往回注射),然后(2)将该0.7ml的紫杉醇制剂用于溶解额外的1mg药水瓶的紫杉醇,以增加紫杉醇与Cremophor的相对量。该研究被设计用于评价转换技术抑制体内人肿瘤小鼠模型异种移植模型中的肿瘤生长的功效,并将其功效与标准化学疗法治疗的进行比较,标准化学疗法治疗为将紫杉醇,5百万MDA-MB-435细胞,以悬浮液侧面皮下注射给无胸腺的裸鼠。从细胞移植的同一天(第一天)开始治疗。动物被指定为4组中的一组。在研究持续期间,每周2次监测各动物的肿瘤生长。一组采用紫杉醇治疗,两组采用转换治疗,还有一组采用紫杉醇载体治疗作为对照。
小鼠是无胸腺的BALB/c裸小鼠nu/nu。植入5百万乳腺癌肿瘤系MDA-MB-435(源自NCI的NCL标准,获自美国典型培养物保藏中心ATCC(Rockville,MD))。乳腺导管癌的生长不依赖雌激素。在研究持续期间,常规治疗组(各10只小鼠)只用载体、紫杉醇15mg/kg、腹腔注射每周两次共2周,MIDS转换组(各11只小鼠)在两个不同的能量水平。

将小鼠限制在直径2英尺的直角圆笼中,通过Tristan技术制成的螺管线圈,调节标准小鼠笼使得小鼠始终暴露于紫杉醇的MIDS再现。将一个治疗组中的所有小鼠放在一个笼子中,并在重复再现期间将其保持在大转换线圈的中心圆腔区域之内,同时喂食和给水。期间,测量肿瘤的体积,并在清理笼的时候,将笼从线圈中滑出。这样导致约占研究持续时间的80-90%的工作时间连续暴露。
虽然本发明描述了某些具体的实施方式和应用,但各种变形、修改和其他应用如所附权利要求所包含的一样包含在本发明的范围内。
权利要求
1.一种化学或生化制剂在对该制剂应答的系统中产生作用的方法,包含(a)将所述系统置于电磁转换器的磁场区域中,(b)将低频时域信号施加于所述转换器,该信号的特征是在谱分析中呈现多个制剂特异性谱峰,其中,所述制剂特异性谱峰从通过以下步骤产生的低频时域信号的谱图中识别(i)将这种化学或生化制剂置于具有磁屏蔽和电磁屏蔽的容器中,并且(ii)以在所述时域信号的谱图中有效产生这种可识别的谱峰的噪声水平,将噪声注入样本时,记录来自所述样本的低频时间相关信号,并且(c)当施加所述信号时,以所施加的信号功率将所述样本暴露于所述转换器产生的所述磁场,并且暴露足以在所述系统中产生对于所述系统的制剂特异性作用的周期。
2.如权利要求1所述的方法,其中施加于所述转换器的所述时域信号由以下步骤产生(i)将这种化学或生化制剂置于具有磁屏蔽和电磁屏蔽的容器中,并且(ii)以在所述时域信号的谱图中有效产生这种可识别的谱峰的噪声水平,将噪声注入样本时,记录来自所述样本的低频时间相关信号,并且。
3.如权利要求2所述的方法,其中,步骤(ii)包括如下步骤(a)将噪声以给定的噪声振幅注入所述样本;(b)记录由重叠在所注入噪声上的样本源辐射组成的电磁时域信号,以及(c)在所选的噪声水平范围内,以多个噪声水平中的每个,重复步骤(a)和(b),以及(d)分析所述多个时域信号,所述时域信号是通过产生所述时域信号的谱图,并且基于所述谱线中的信息,识别最优的制剂特异性时域信号而生成的。
4.如权利要求3所述的方法,其中在步骤(iia)中,注入所述样本的噪声的噪声源包括功率可调的高斯噪声发生器和Helmholz线圈,并且其从所述噪声发生器接收所选的达1伏特的范围内的噪声输出信号。
5.如权利要求3所述的方法,其中,所述分析步骤(d)可以由以下步骤之一实现(i)生成直方图,对于DC至8kHz范围内的所选频率范围上的每个事件箱f,所述直方图示出了每箱(bin)中的事件计数的数量,其中f是用于对所述时域信号采样的样本率,将与高过给定阈值的箱的数量相关的得分分配给所述直方图,并且根据所述得分选择时域信号;(ii)使所述时域信号自相关,在DC至8kHz范围内的所选频率范围上,生成所述自相关信号的FFT,将与高过平均噪声值的峰值的数量相关的得分分配给所述FFT信号,并且根据所述得分选择时域信号;并且(iii)计算在DC和8kHz之间的所选频率范围内的多个所定义的时间周期中的每个上的所述时域信号的一系列傅立叶谱,对所述傅立叶谱求平均;将与高过平均噪声值的峰值的数量相关的得分分配给所述平均后的FFT信号,并且根据该记分选择时域信号。
6.如权利要求3所述的方法,其中所述分析步骤(d)可以由以下步骤之一实现(i)存储所述样本在样本持续时间上的时域信号;(ii)选择用于对所述时域信号进行采样的采样率F,其中F*T是样本总计数S,F是以采样率F采样的所述时域信号的实数快速傅立叶变换的频域分辨率的大约两倍,且S>(2)f*n,其中n至少为10,(iii)从所存储的时域信号中选择S/n个样本,并对所述样本执行实数快速傅立叶变换(RFFT),(iv)归一化所述RFFT信号并计算所述信号的平均功率,(v)将事件计数置于f所选频率事件箱中的每个中,其中,获得相应的所选频率处的所测量的功率大于等于平均功率*ε,其中0<ε<1,并且选择相应的所选频率处的所测量的功率大于等于平均功率*ε,使得置于事件箱中的计数的总数约为该箱中的最大可能箱计数的20-50%之间。(vi)重复步骤(iii-v)N>2次,并且(vii)生成直方图,对于所选频率范围上的每个事件箱f,所述直方图示出了每个箱中的事件计数的数量。
7.如权利要求6所述的方法,还包括,在用于归一化所述RFFT信号的步骤(iv)中,包括将来自所述RFFT的所述归一化的功率值置于f相应频率功率箱中,以及在步骤(viii)中(a)用n除置于所述f功率箱中的每个中的累积值,从而得到每个箱中的平均功率,并且(b)在所述直方图上显示每个箱中的所述平均功率。
8.如权利要求2所述的方法,其中,所述记录使用与SQUID耦合的梯度计执行,所述注入包括将噪声注入所述梯度计。
9.如权利要求1所述的方法,其中,所述制剂是能够在生物系统中通过配体特异性、非共价键相互作用与抗配体细胞靶相互作用的配体。
10.如权利要求9所述的方法,所述制剂是能与一种或多种基因的启动子相互作用的化合物,由于所述化合物的存在上调或下调所述启动子,所述低频时域信号通过记录所述化合物的低频时域信号来产生,并且以信号振幅执行所述施加,且施加的时间足以产生所述一个或多个基因的可测量的上调或下调。
11.如权利要求10所述的方法,其中,所述制剂是L(+)阿拉伯糖,并且所述基因包括细菌lac操纵子基因。
12.如权利要求9所述的方法,所述制剂是能与酶相互作用的化合物,从而竞争性地抑制所述酶,该酶是生物系统的生长和生存所需要的,通过记录所述化合物的低频时域信号来产生所述低频时域信号,以信号振幅执行所述施加,且施加的时间足以产生对所述生物系统的可测量的抑制或生长和生存。
13.如权利要求12所述的方法,其中,所述制剂是甘胺酸磷酸盐,以及所述靶酶是植株中的5-烯醇式丙酮莽草酸酯-3-磷酸(EPSP)合酶。
14.如权利要求12所述的方法,其中,所述制剂是phepropeptin D,以及所述靶酶是在真核细胞内与蛋白体相关的蛋白水解酶。
15.如权利要求9所述的方法,用于治疗哺乳动物对象中的癌,其中,所述制剂选自包括由以下组成的组(a)微管蛋白结合制剂,其中所述细胞靶是染色体微管纺锤体;(b)蒽环类抗生素,其中所述细胞靶为双链DNA;(c)拓扑异构酶抑制剂,其中所述细胞靶为拓扑异构酶;(d)抗代谢剂,其中所述细胞靶为细胞代谢所需要的酶;(e)免疫抑制素,其中所述细胞靶为免疫应答细胞;以及(f)抑癌蛋白,其中所述靶为细胞核内的DNA复制结构。
16.如权利要求15所述的方法,其中,所述制剂是紫杉醇或紫杉醇类似物。
17.如权利要求15所述的方法,其中,所述制剂是p53肿瘤抑制因子。
18.如权利要求1所述的方法,其中,所述电磁转换器包括定义了内部开口的线圈绕组,并且所述暴露包括将所述样本置于所述绕组的所述内部开口处。
19.如权利要求1所述的方法,其中,所述电磁转换器包括可植入的线圈,并且所述转换器在所述暴露之前植入生物系统。
20.如权利要求1所述的方法,其中,以有效产生在10-100mG范围之内的所选磁场强度的所选功率水平,施加所述信号。
21.如权利要求1所述的方法,其中所述暴露通过1小时或更长时间的交替的,间歇的暴露来实现。
22.使化学或生化制剂在对该制剂应答的系统中产生作用的装置,包括(a)用于存储低频时域信号的存储设备,该信号的特征是在谱分析中呈现多个制剂特异性谱峰,其中,所述制剂特异性谱峰从通过以下步骤产生的低频时域信号的谱图中识别(i)将这种化学或生化制剂置于有磁屏蔽和电磁屏蔽的容器中,并且(ii)以在所述时域信号的谱图中有效地产生这种可识别谱峰的噪声水平,将噪声注入所述样本时,记录来自所述样本的低频时间相关信号,并且(b)定义了活性磁场区域的电磁转换器,所述活性磁场在电磁信号被施加给所述转换器时产生,并且所述样本置于所述活性磁场中,以及(c)用于将所述存储设备连接到所述转换器上的放大器,通过所述转换器在所述活性磁场区域中的这种放置方式,用于以所施加的信号功率将所述信号施加给所述转换器,并且施加的周期足以在所述系统中产生对所述系统的制剂特异性作用。
23.如权利要求22所述的装置,其中,所述电磁转换器包括线圈绕组和适于放置所述样本的内部开口。
24.如权利要求22所述的装置,其中,所述电磁转换器是有一对对齐的电磁线圈的Helmholz线圈,该线圈中间定义了暴露的操作台,并且所述暴露包括将所述样本置于所述操作台内。
25 如权利要求22所述的装置,其中,所述电磁转换器包括可植入的线圈。
26.如权利要求22所述的装置,其中,所述存储设备处于远离所述转换器和放大器的位置,并且所述信号从这种远端位置被传送到所述转换器。
27.一种化学或生物活性制剂的低频时域信号,其特征是在谱分析中呈现多个制剂特异性谱峰,其中所述制剂特异性谱峰是从通过以下步骤产生的低频时域信号的谱图中识别的(i)将这种化学或生化制剂置于有磁屏蔽和电磁屏蔽的容器中,并且(ii)以在所述时域信号的谱图中有效地产生这种可识别的谱峰的噪声水平,将噪声注入所述样本时,记录来自所述样本的低频时间相关信号。
28.如权利要求27所述的信号,由以下步骤产生(i)将这种化学或生化制剂置于有磁屏蔽和电磁屏蔽的容器中,并且(ii)以在所述时域信号的谱图中有效地产生这种可识别的谱峰的噪声水平,将噪声注入所述样本时,记录来自所述样本的低频时间相关信号,并且。
29.如权利要求28所述的信号,其在步骤(ii)中通过以下步骤产生(a)以给定的噪声振幅,将噪声注入所述样本;(b)记录由重叠在所注入噪声上的样本源辐射组成的电磁时域信号,并且(c)在所选的噪声水平范围内以多个噪声水平中的每个重复步骤(a)和(b),并且(d)分析所述多个时域信号,所述时域信号是通过产生所述时域信号的谱图,并且基于所述谱图中的信息,识别最优的制剂特异性时域信号而生成的。
30.如权利要求27所述的信号,其包括在谱分析中两种或多种不同的化合物的制剂特异性谱峰。
全文摘要
公开了一种使化学或生化制剂在对该制剂应答的系统上产生作用的方法和装置。在实际的本发明的方法中,该系统放置在电磁转换器的磁场区域内,将低频率、时域信号施加于电磁转换器,该信号的特征是在谱分析中呈现大量制剂特异性谱峰。这些谱峰从如下产生的低频,时域信号的谱图中识别(i)放置这样一种化学或生化制剂于兼有磁屏蔽和电磁屏蔽的容器内,并且(ii)在时域信号谱图中有效产生可识别谱峰的噪声水平注射噪声进入样本,记录来自样本的低频,时域信号。将样品暴露于转换器产生的磁场中,在适用的信号功率和充足的时间内,在系统中产生对系统的制剂特异性作用。
文档编号C12M1/42GK101031796SQ200580025199
公开日2007年9月5日 申请日期2005年7月27日 优先权日2004年7月27日
发明者约翰·布特斯, 贝内特·M·巴特斯, 托马斯·钟, 罗宾·瓦尔盖斯, 帕特里克·诺顿 申请人:纳提维斯公司
网友询问留言 已有0条留言
  • 还没有人留言评论。精彩留言会获得点赞!
1