Development of rare-event sampling techniques for predicting structures and free energies of crystal polymorphs and oligopeptides

纽约大学的Mark塔克曼获得了化学系化学理论、模型和计算方法项目的奖励，以开发预测分子晶体结构的方法和软件。 该奖项由CISE/ACI软件重用风险基金共同资助。 在材料科学中，有序的分子阵列形成被称为分子晶体的结构，在制药，电子和国防工业中发挥着重要作用。通常，关键问题是应该为特定应用制造哪些晶体。 值得注意的是，最广泛使用的药物分子晶体之一阿司匹林基本上是偶然发现的。通常，在晶体工程中，有必要筛选潜在候选化合物的大型数据库。不幸的是，在实验室中制造和表征分子晶体通常耗时且昂贵，使得通过这样的数据库进行试错的方法不切实际。如果能够应用系统的、有针对性的方法，还能发现多少更重要的分子晶体系统？理论和计算在原则上可以快速预测分子晶体结构及其性质，在创建这样一种有针对性的方法方面具有独特的优势。然而，需要的是用于进行这些预测的鲁棒算法。 塔克曼小组开发了计算技术和软件，用于预测给定化合物可以形成的晶体结构，并根据称为自由能的热力学性质对其进行排名，自由能在科学界被认为是这种排名的适当品质因数，但仍然是一个难以捉摸的属性来确定。塔克曼小组还采用这些算法来研究被称为寡肽的氨基酸短链的构象偏好，以探索这些重要的生物分子在免疫原性和新型药物设计中的作用。 塔克曼和他的同事们从事许多软件活动，包括开发用于晶体结构预测的计算机包，提高分子动力学软件PINY-MD的效率，并继续为许多社区软件代码贡献软件。 在这个项目中开发的所有软件都提供给更广泛的研究社区。固态分子材料的基本性质通常受到其晶体结构细节和多晶型物存在的强烈影响。 这些结构的实验测定是昂贵和耗时的，这使得理论和计算的作用越来越重要。 类似地，小寡肽的生物化学功能，从免疫原性到抑制，受到它们在不同环境中的平衡构象的影响。 在复杂系统中，结构的计算预测是具有挑战性的，因为在粗糙的势能景观上存在所谓的稀有事件采样问题，这是在试图研究许多复杂系统的平衡热力学和动力学时出现的。能量表面的粗糙度是指构象和结构变化的高障碍的存在。 塔克曼小组已经提出开发鲁棒的基于自由能的增强的采样算法和软件，用于克服在寡肽的晶体结构预测和构象采样中出现的罕见事件问题，从而允许以有效的方式识别有利的结构并进行分级。在所提出的方法中，自由能景观表示在选择一组集体变量（CV），旨在区分这些系统中的不同结构图案。CV首先受到新的表面导航技术，以确定最小值和鞍点，统称为“地标”的景观，然后有针对性的增强采样，以产生的地标的自由能排名。这些新技术被应用于预测刚性和柔性有机小分子的晶体结构和多晶型物，研究免疫原性肽与主要组织相容性复合物结合的构象自由能景观，并了解机械力对<$-发夹肽解折叠机制的影响。通过塔克曼小组组织的黑客马拉松，软件开发将得到加速。 通过在纽约大学全球校园举办的讲习班，对学生进行罕见事件方法的教育。 最后，塔克曼集团通过在纽约市有业务的全国性组织接触代表性不足的群体，以帮助设计和参与与STEM相关的教育活动。

Mark Tuckerman of New York University is supported by an award from the Chemical Theory, Models and Computational Methods program in the Chemistry Division to develop methods and software for the prediction of molecular crystal structure. This award is cofunded by the CISE/ACI Software Reuse Venture Fund. In the science of materials, ordered arrays of molecules forming structures known as molecular crystals play an essential role in the pharmaceutical, electronics, and defense industries. Often, the crucial question is which crystals should be made for a particular application. It is worth noting that one of the most widely used pharmaceutical molecular crystals, aspirin, was discovered essentially by accident. Typically, in crystal engineering, it is necessary to screen large databases of potential candidate compounds. Unfortunately, making and characterizing molecular crystals in the laboratory is generally time consuming and costly, rendering a trial-and-error approach through such a database impractical. How many more important molecular crystal systems might be discovered if a systematic, targeted approach could be applied? Theory and computation, which can, in principle, rapidly predict molecular crystal structures and their properties, are uniquely poised to play a key role in creating such a targeted approach. What is needed, however, are robust algorithms for making these predictions. The Tuckerman group develops computational techniques and software for predicting the crystal structures a given compound can form and ranking them according to a thermodynamic property known as free energy, which has been recognized in the scientific community as the proper figure of merit for such a ranking but has remained an elusive property to determine. The Tuckerman group also adapts these algorithms for studying the conformational preferences of short chains of amino acids known as oligopeptides in order to explore the role these important biological molecules play in immunogenicity and the design of new classes of pharmaceuticals. Tuckerman and his coworkers are engaged in many software activities including developing a computer package for crystal structure prediction, improving the efficiency of their molecular dynamics software, PINY-MD and continuing to contribute software to many community software codes. All of the software developed in this project is made available to the broader research community. The basic properties of molecular materials in the solid state are often strongly influenced by the details of their crystal structures and the existence of polymorphs. Experimental determination of these structures is costly and time-consuming, which places increased importance on the role of theory and computation. Similarly, the biochemical function of small oligopeptides, from immunogenicity to inhibition, is affected by their equilibrium conformations in different environments. Computational prediction of structure in complex systems such as these is challenging due to the so-called rare-event sampling problem on a rough potential energy landscape, which arises when attempting to study the equilibrium thermodynamics and kinetics of many complex systems. Roughness on an energy surface refers to the existence of high barriers to conformational and structural changes. The Tuckerman group has proposed to develop robust free-energy based enhanced sampling algorithms and software for overcoming the rare-event problem that arises in the crystal structure prediction and conformational sampling of oligopeptides, thereby allowing favored structures to be identified and thermodynamically ranked in an efficient manner. In the proposed methods, the free energy landscape is expressed in terms of select set of collective variables (CVs) designed to distinguish the different structural motifs in these systems. The CVs are first be subject to new surface navigation techniques in order to identify the minima and saddles points, collectively referred to as "landmarks" on the landscape, and then targeted for enhanced sampling in order to produce the free energy ranking of the landmarks. The new techniques are applied to predict the crystal structures and polymorphs of both rigid and flexible small organic molecules, to study the conformational free energy landscape of an immunogenic peptide binding to the major histocompatibility complex, and to understand the influence of mechanical force on the unfolding mechanism of â-hairpin peptide. Software creation will be accelerated via hackathons organized by the Tuckerman group. Education of students in rare-event methods is aided through workshops organized at New York University's global campus sites. Finally, the Tuckerman group reaches out to underrepresented groups via national organizations having a presence in New York City in order to help devise and participate in STEM-related educational activities.