课前准备--小分子药物的结构分析

作者，Evil Genius

接下来为我们的生化小课，CADD计算机辅助药物设计，分子对接/分子动力学做好准备。

两种主要类型的药物是小分子药物和生物制剂。

小分子药物通常通过化学合成，并且比生物制剂简单。它们被设计为针对与疾病过程相关的特定分子，并且通常具有定义的结构，允许它们与这些靶点高特异性地结合。

生物制剂通常是由活细胞产生的大型复杂分子，如蛋白质,它们被设计为与人体内的特定受体或通路相互作用。

官能团是小分子药物设计的关键结构特征之一。如羟基(-OH)、羧基(-COOH)和氨基(-NH2)等官能团的存在可以使小分子与特定的酶(受体)相互作用并在体内赋予特定的药理学特性。

例如，乙醇分子中羟基的存在可以使它们与大脑中的受体相互作用，导致众所周知的饮酒效应。

立体化学是小分子药物设计的一个关键方面，指一个分子中原子的三维布局以及这种布局如何影响分子与生物靶点的相互作用。例如，抗炎药异布比酮的活性形式具有特定的立体化学，这使其能够与环加氧酶结合并抑制其活性。

分子量在小分子药物设计中也是一个重要的考虑因素，它可以影响药物的药代动力学和药效动力学。例如，低分子量的药物可以被体内快速吸收和排泄,而高分子量的药物可能有更长的半衰期和需要更低的剂量。

分子的大小和形状可以影响其物理和化学性质。例如，较大的分子往往比小分子具有更高的沸点和熔点。

分子中特定化学键的存在可以影响其结构和性质。例如，共价键是通过两个原子之间共享价壳电子对形成的，IV到VI族的原子结合以完成价壳层八电子。

分子的电荷特性，如偶极矩和电荷分布，可以影响其与其他分子的相互作用及其反应活性。通过了解这些结构特征的重要性，研究人员可以设计小分子药物候选化合物，优化其药理学特性并减少副作用。

计算方法在设计新型候选药物和优化其性质方面起着至关重要的作用。近年来，计算方法在药物设计中的应用迅速增加，有许多工具可支持这一过程

分子对接是小分子药物设计中使用的主要计算方法之一。

分子对接涉及预测小分子(称为配体)与较大的靶点(如蛋白质或核酸)之间的结合亲和力和取向。

对接算法使用能量最小化和评分函数来评估配体-蛋白质复合物的稳定性并对潜在的药物候选化合物进行排名。

对接过程有助于研究人员深入了解各种生物过程(如药物-蛋白质相互作用、蛋白质-蛋白质相互作用和酶促反应)的分子机制。

随着计算机技术和生命科学领域的融合,分子对接也取得了巨大进步，显著增强了这项技术的精确性和有效性。

(i)自由能微扰(FEP)和热力学积分(Ti)等ALCHEMICAL自由能方法已被开发来计算配体与蛋白质的结合自由能。这些方法可以比传统对接方法提供更准确的结合亲和力预测。

(ii)开发了新的对接工具,如HADDOCK和RosettaDock,可以将蛋白质柔性引入对接模拟。这些工具可以预测配体结合导致的蛋白质构象变化,从而提高对接预测的准确性。

(iii)机器学习已被纳入对接工具以提高其预测准确性。

虚拟筛选

虚拟筛选是使用分子对接或其他技术对大量化合物库进行筛选,以识别潜在的药物候选化合物。

·节省成本和时间，帮助研究人员缩小进一步实验分析的化合物数量。

·用于识别可能与靶向蛋白质结合的小分子。

· 实现工具:GOLD、GLIDE和AutoDockVina等。AutoDock Vina是免费的且能提供良好结果的用于虚拟筛选的工具。

虚拟筛选结果的准确性取决于蛋白质结构和化合物库的质量。

虚拟筛选过程不考虑化合物的药代动力学和毒理学特性。

通过虚拟筛选识别的化合物需要实验验证,以确定其疗效和安全性。

分子动力学模拟

MD模拟是一种强大的计算工具，用于研究体系中原子和分子的物理运动和相互作用。

MD模拟按时间积分牛顿运动方程,以获得体系中原子/分子的运动,从而在原子水平上提供该体系宏观行为的定量和定性信息。

·MD模拟提供了研究系统“演化”动态的“电影”。

MD模拟由经典运动方程的数值步进解决组成对于一个简单的原子系统:

1.每个粒子上的力由原子间相互作用确定,这通常由描述系统势能作为原子位置函数的数学函数来表示(即力场)。

2.原子间势能函数可以使用量子力学、经验势函数或两者的组合来获得。

(i)深度学习算法越来越多地应用于MD模拟。例如，原子卷积网络已用于直接从分子的原子坐标预测分子的势能,这可以显著提高MD模拟的准确性。

(i)开发了增强取样技术,如metadynamics,以克服传统MD模拟的局限性。这些技术可以探索更大的构象空间，并提供系统能量面更准确的代表。

(iii)粗粒度模型通过将原子集合在一起来简化系统的表示，从而降低了MD模拟的计算成本。这些模型可以模拟传统全原子模型无法模拟的更大系统和更长时间尺度。

(iv)基于GPU的工具,如GROMACS和AMBER,已被开发以加速MD模拟。这些工具比基于CPU的工具执行模拟更快,使得可以研究更大和更复杂的系统。

量子力学与分子力学的计算

量子力学/分子力学(QM/MM)计算是一种强大的计算方法,用于在原子和分子水平研究化学和生化系统。QM/MM方法将量子力学计算的精确性与分子力学模拟的速度相结合，来研究化学和生化系统,这使其成为研究溶液或酶中的化学和生化系统的有价值工具。QM/MM计算可以提供准确的反应能和反应路径计算。

在传统的分子力学和量子力学计算中,仅考虑给定冻结配置系统内分子和分子间相互作用，通常不带显式溶剂化。

QM/MM方法可以通过将QM区域视为溶质(或关键区域)和MM区域视为溶剂(或非关键区域)来考虑溶剂化效应。

QM/MM方法由Warshel和Levitt在1976年首次提出。从那时起，它已被广泛用于研究各种化学和生化系统,包括酶、反应机制和蛋白质-配体结合相互作用。

QM/MM方法能够模拟仅用量子力学方法很难有效处理的较大系统。

计算方法的局限性

蛋白质灵活性:传统的对接算法通常将蛋白质视为刚性体,这是一个非常大的简化。蛋白质高度动态并可以采用一系列构象。考虑这种灵活性是对接算法的一个重大挑战。

评分函数:对接算法中用于预测配体与蛋白质的结合亲和力的评分函数通常不准确。它们通常基于分子间相互作用的简化模型，不能完全捕获这些相互作用的复杂性。

溶剂效应:许多对接算法没有充分考虑溶剂的影响。溶剂可以在相互作用中发挥关键作用，忽略这些效应会导致不准确的预测。

时间尺度:分子动力学模拟可以模拟分子的动态行为，但它们受可以模拟的时间范围的限制。大多数模拟在纳秒到微秒尺度上,而许多生物过程发生在毫秒到秒尺度上。

计算成本:对接和分子动力学模拟计算量都很大。这限制了可以研究的系统的大小和可以执行的模拟的长度。

实验数据的缺乏:对接和分子动力学模拟的预测通常缺乏实验数据进行验证。这使得很难评估这些方法的准确性。

预测诱导配合的困难:许多对接算法在预测配体与蛋白质结合时发生的构象变化(称为诱导配合)方面遇到困难。

预测非靶向效应的局限:虽然计算方法在预测药物与其预期靶点的相互作用方面正在改进,但预测潜在的非靶向效应仍然是一个重大挑战。

生活很好，有你更好。