Steric Effects Govern the Photoactivation of Phytochromes

Phytochromes constitute a superfamily of photoreceptor proteins existing in two forms that absorb red (Pr) and far‐red (Pfr) light. Although it is well‐known that the conversion of Pr into Pfr (the biologically active form) is triggered by a Z→E photoisomerization of the linear tetrapyrrole chromophore, direct evidence is scarce as to why this reaction always occurs at the methine bridge between pyrrole rings C and D. Here, we present hybrid quantum mechanics/molecular mechanics calculations based on a high‐resolution Pr crystal structure of Deinococcus radiodurans bacteriophytochrome to investigate the competition between all possible photoisomerizations at the three different (AB, BC and CD) methine bridges. The results demonstrate that steric interactions with the protein are a key discriminator between the different reaction channels. In particular, it is found that such interactions render photoisomerizations at the AB and BC bridges much less probable than photoisomerization at the CD bridge.     

应用QM与ABEEM/MM结合的方法研究NAMI-A的结构性质

应用量子力学(QM)与ABEEM浮动电荷力场(ABEEM/MM)相结合的方法研究了抗癌药物NAMI-A在水溶液中的结构性质. 所有的结构优化都是在DFT的B3LYP方法下采用6-31G(d,p)和LanL2DZ基组完成的, 没有加入任何限制性条件. 结果表明, 优化得到的NAMI-A构型受不同环境及方法的影响均有变化. 与气相中得到的构型相比, QM/MM迭代优化得到构型要比PCM的构型变化更明显. QM/MM (ABEEM/MM)迭代优化得到的NAMI-A构型比QM/MM (OPLS-AA)的变化要小. 总之, 溶剂通过极化效应对NAMI-A结构、电荷分布及径向分布函数等性质均有影响, 客观地处理极化效应才能正确地反映QM区的性质.     

Hybrid QM/MM studies on energetics of malonaldehyde in condensed phase

Energetics of eight enol isomers of a malonaldehyde (MA) at the ground state in CCl₄ solvent environment have been investigated using a hybrid quantum mechanical (QM)/molecular mechanics (MM) method. It is found that relative energies of the isomers slightly change due to interactions between MA and surrounding atoms. In an isolated environment all eight isomers have stable planar structure. On the other hand, most of the isomers have nonplanar structure in CCl₄ whose interactions with the solute molecules are, however, weak. Mainly, structural changes are found in the H? O? C?C dihedral angle ?_{H? O? C?C}, i.e., a hydrogen atom that is connected directly with an oxygen atom, is located at a nonplanar position, and other atoms remain almost planar. Vertical excitation energies of low‐lying excited states at the resultant optimized structure of each isomer are evaluated. The vertical excitation energies in CCl₄ are almost the same as those in the isolated environment, but some changes were found in the triplet excitation states. © 2004 Wiley Periodicals, Inc. Int J Quantum Chem, 2004     

Is inhibition process better described with MD(QM/MM) simulations? The case of urokinase type plasminogen activator inhibitors

Urokinase plasminogen activator (uPA) is an enzyme involved in cancer growth and metastasis. Therefore, the design of inhibitors of uPA is of high therapeutic value, and several chemical families have been explored, even if none has still emerged, emphasizing the need of a rationalized approach. This work represents a complete computational study of uPA complexed with five inhibitors, which present weak similarities. Molecular dynamics simulations in explicit solvent were conducted, and structural analyses, along with molecular mechanics (MM)/Poisson-Boltzmann surface area free energies estimations, yield precious structure-activity relationships of these inhibitors. Besides, we realized supplemental QM/MM computations that improved drastically the quality of our models providing original information on the hydrogen bonds and charge transfer effects, which are, most often, neglected in other studies. We suggest that these simulations and analyses could be reproduced for other systems involving protein/ligand molecular recognitions.     

Geometry optimization with QM/MM,ONIOM, and other combined methods. I. Microiterations and constraints

Hybrid energy methods such as QM/MM and ONIOM, that combine different levels of theory into one calculation, have been very successful in describing large systems. Geometry optimization methods can take advantage of the partitioning of these calculations into a region treated at a quantum mechanical (QM) level of theory and the larger, remaining region treated by an inexpensive method such as molecular mechanics (MM). A series of microiterations can be employed to fully optimize the MM region for each optimization step in the QM region. Cartesian coordinates are used for the MM region and are chosen so that the internal coordinates of the QM region remain constant during the microiterations. The coordinates of the MM region are augmented to permit rigid body translation and rotation of the QM region. This is essential if any atoms in the MM region are constrained, but it also improves the efficiency of unconstrained optimizations. Because of the microiterations, special care is needed for the optimization step in the QM region so that the system remains in the same local valley during the course of the optimization. The optimization methodology with microiterations, constraints, and step-size control are illustrated by calculations on bacteriorhodopsin and other systems.