利用晶体结构预测与Form 4结构解析重新审视利托那韦多晶型景观：对历史危机的再评估与预警

《Communications Chemistry》：Predicting the ritonavir crisis by revisiting the polymorph landscape with crystal structure prediction and form 4 structure solution

本研究针对药物多晶型风险预警难题，通过晶体结构预测(CSP)与三维电子衍射(3D-ED)技术，首次解析了利托那韦Form 4的晶体结构并重构其能量景观。研究表明现代CSP技术可提前10年预警Form 2的稳定性风险，揭示构象无序对晶型稳定的关键作用，为药物固态研发提供了重要风险评估范式。

　　

在药物研发史上，利托那韦的"多晶型灾难"堪称里程碑事件——这种抗HIV药物在上市后，其有效晶型Form 1意外转变为更稳定但溶解度较低的Form 2，导致药物生物利用度急剧下降，不仅造成2.5亿美元的经济损失，更危及患者的治疗。这一事件深刻暴露了药物多晶型风险对药品安全性的重大威胁，也促使制药行业加强固态形式筛选。然而，一个关键问题始终悬而未决：如果当时具备现代计算技术，能否预测这场危机？此外，利托那韦的实验多晶型景观至今仍不完整，特别是Form 4的结构一直未被解析。

发表在《Communications Chemistry》的这项研究，通过结合前沿计算方法与实验技术，对利托那韦多晶型风险进行了全面回顾性评估。研究团队运用晶体结构预测(CSP)技术重构了利托那韦的晶体能量景观，并首次通过三维电子衍射(3D-ED)技术解析了Form 4的晶体结构，填补了该领域的重要知识空白。

关键技术方法包括：采用GRACE 2.8软件进行柔性分子CSP计算，覆盖21个手性空间群；利用TRHu(ST)方法计算室温自由能；结合毛细管结晶、X射线粉末衍射(XRPD)和3D-ED解析Form 4结构；应用交叉相关函数法匹配计算与实验数据。

结果与讨论

Form 4的结构解析

研究团队通过优化结晶条件，利用毛细管结晶技术成功获得了Form 1和Form 4的混合物。由于Form 4呈极细的针状形态且对电子束敏感，传统单晶X射线衍射难以解析其结构。研究创新性地结合3D-ED、XRPD和CSP生成的结构，最终成功解析了Form 4的完整晶体结构。晶体学数据显示Form 4属于P2₁空间群，晶胞参数为a=14.13 ?，b=5.16 ?，c=26.47 ?，β=96.67°，体积为1915.0 ?3。

结构分析表明，Form 4与不稳定的Form 1和Form 3类似，均含有反式氨基甲酸酯构型，而与稳定的Form 2的顺式构型不同。其氢键网络分析显示，Form 4的酰胺基团形成连续氢键链，但醇基与酰胺氧原子形成氢键，这与Form 1的醇基-噻唑氮原子氢键有所不同。特别重要的是，3D-ED结构解析揭示了Form 4边缘异丙基存在构型无序，两种构型的占据率分别为58%和42%，这一特征与Form 1相似，但与有序的Form 2形成鲜明对比。

晶体能量景观

CSP研究构建的利托那韦晶体能量景观显示，与Form 2匹配的结构在300K时是最稳定的预测多晶型，比Form 1低1.62 kJ·mol^-1，与实验测得的2-3 kJ·mol^-1稳定性差异一致。考虑计算误差范围，排名第三的预测多晶型比Form 2不稳定5.51 kJ·mol^-1，仅有0.2%的概率可能比Form 2更稳定，表明利托那韦不太可能存在比Form 2更稳定的未知多晶型。

从风险评估角度，若在利托那韦开发初期进行CSP分析，即可预测到比Form 1更稳定的全球能量最小值（即后来的Form 2）的存在风险。结构比较显示这一预测结构含有罕见构象但更完整的氢键网络，是动力学受阻多晶型的典型特征。溶解度计算表明Form 2的溶解度比Form 1低约两倍，与实际观察到的溶解度下降幅度相符。

研究还量化了无序对晶型稳定性的影响：无序使Form 1相对于有序Form 2稳定了0.19 kJ·mol^-1，而使Form 4稳定了0.58 kJ·mol^-1。虽然这一稳定效应不足以改变全局稳定性排序，但在能量接近的多晶型系统中可能起决定性作用。

对于已知的Form 3（YIGPIO06），计算显示其在300K时排名第22位，比全局最小值不稳定15.42 kJ·mol^-1，远大于实验溶解度比值推算的约8 kJ·mol^-1差异。这种差异可能由于Form 3存在Z'=4（四个不对称单元分子）且所有分子均存在无序，而计算仅优化了主要无序组分，未充分考虑构型熵的稳定化效应。

对制药行业的意义

这项研究证实，药物开发史上最著名的多晶型灾难是可以通过现代CSP技术预测的。对于利托那韦这样的案例，CSP计算成本与潜在损失（2.5亿美元）相比微不足道。CSP作为实验研究的补充，可帮助判断已知晶型是热力学稳定还是动力学产物，但当前技术仍存在局限：对成核生长动力学的预测能力不足；Z'=1以上的晶体结构预测仍不现实（Form 3的Z'=4结构无法预测）；溶剂化物的预测更为复杂。

研究同时证明了3D-ED在微晶结构解析中的价值。当无法获得适合单晶X射线衍射的晶体时，3D-ED结合计算方法可提供准确的结构解决方案，这对风险评估至关重要，因为许多通过高通量筛选发现的晶型仅能获得XRPD图谱，而无法解析结构确认相对稳定性。

结论

本研究通过状态参数计算工作流证明，利托那韦灾难本可被预测：基于计算结果的风险评估很可能避免其发生。同时，通过3D-ED、XRPD和计算建模的结合，首次解析了利托那韦Form 4的结构，证实其热力学不稳定而不构成开发风险。这种结构解析和无序确定方法对于无法通过单晶X射线衍射确定多晶型结构的风险评估案例具有重要意义。研究还强调了确定晶体结构中无序现象的重要性，因为无序对多晶型稳定性的影响在风险评估中不容忽视。

这项工作不仅为利托那韦多晶型事件提供了完整的科学解释，也为制药行业提供了实用的风险评估框架和方法论支持，彰显了计算预测与实验表征相结合在药物固态研究中的强大潜力。