Computation of Electromagnetic Properties of Molecular Ensembles

We outline a methodology for efficiently computing the electromagnetic response of molecular ensembles. The methodology is based on the link that we establish between quantum-chemical simulations and the transfer matrix (T-matrix) approach, a common tool in physics and engineering. We exemplify and analyze the accuracy of the methodology by using the time-dependent Hartree-Fock theory simulation data of a single chiral molecule to compute the T-matrix of a cross-like arrangement of four copies of the molecule, and then computing the circular dichroism of the cross. The results are in very good agreement with full quantum-mechanical calculations on the cross. Importantly, the choice of computing circular dichroism is arbitrary: Any kind of electromagnetic response of an object can be computed from its T-matrix. We also show, by means of another example, how the methodology can be used to predict experimental measurements on a molecular material of macroscopic dimensions. This is possible because, once the T-matrices of the individual components of an ensemble are known, the electromagnetic response of the ensemble can be efficiently computed. This holds for arbitrary arrangements of a large number of molecules, as well as for periodic or aperiodic molecular arrays. We identify areas of research for further improving the accuracy of the method, as well as new fundamental and technological research avenues based on the use of the T-matrices of molecules and molecular ensembles for quantifying their degrees of symmetry breaking. We provide T-matrix-based formulas for computing traditional chiro-optical properties like (oriented) circular dichroism, and also for quantifying electromagnetic duality and electromagnetic chirality. The formulas are valid for light-matter interactions of arbitrarily-high multipolar orders.     

The distributed diagonal force decomposition method for parallelizing molecular dynamics simulations

Parallelization is an effective way to reduce the computational time needed for molecular dynamics simulations. We describe a new parallelization method, the distributed-diagonal force decomposition method, with which we extend and improve the existing force decomposition methods. Our new method requires less data communication during molecular dynamics simulations than replicated data and current force decomposition methods, increasing the parallel efficiency. It also dynamically load-balances the processors' computational load throughout the simulation. The method is readily implemented in existing molecular dynamics codes and it has been incorporated into the CHARMM program, allowing its immediate use in conjunction with the many molecular dynamics simulation techniques that are already present in the program. We also present the design of the Force Decomposition Machine, a cluster of personal computers and networks that is tailored to running molecular dynamics simulations using the distributed diagonal force decomposition method. The design is expandable and provides various degrees of fault resilience. This approach is easily adaptable to computers with Graphics Processing Units because it is independent of the processor type being used.     

计算机辅助蛋白质分子理性设计:从肌红蛋白到一氧化氮还原酶

计算机辅助蛋白质分子理性设计在解决化学及生物学重要问题中被证实十分有效。在NOR本身三维结构未知的情况下通过计算机分子模拟,使用肌红蛋白(Mb)作为蛋白质分子模型,设计了结构功能型一氧化氮还原酶(NOR),所设计的NOR蛋白质模型——Fe_BMb一年后被天然NOR的晶体结构所证实。本文综述了设计Fe_BMb, I107E Fe_BMb以及Fe_BMb(-His)的研究过程及其设计合理性,评述了通过使用计算机分子模拟,获得Mb处于双组氨酸配位的非天然状态的原子层次结构信息,而这些信息很难通过实验方法来获得。计算机辅助蛋白质分子理性设计的广泛应用将会为生物体系提供更深刻的内涵。     

计算机化学模拟—分子构象识别的新方法

简介了几种利用计算机图形技术研究化合物分子构象的新方法。重点介绍了分子力学计算方法中的系统搜索、随机搜索方法和分子动力学计算方法中的模拟淬火、模拟退火等新技术 ,为药物分子设计中受体与配体分子构象的识别提供了合理可行的方法。     

Aerolysin单分子界面的构建及高选择性单分子检测

提出了一种基于Aerolysin膜蛋白质分子构建单分子界面的方法, 运用蛋白质工程技术对单分子界面进行定点修饰, 所建立方法灵活、 可控且重复性好. 采用Poly(dA)₄为探针分子对修饰后的单分子界面进行了表征, 结果表明, 在孔口处的Arginine修饰影响了寡聚核苷酸的选择性. 为进一步理解Aerolysin单分子界面及合理设计功能性单分子界面提供了参考.     

Protein‐specific force field derived from the fragment molecular orbital method can improve protein–ligand binding interactions

Accurate computational estimate of the protein–ligand binding affinity is of central importance in rational drug design. To improve accuracy of the molecular mechanics (MM) force field (FF) for protein–ligand simulations, we use a protein‐specific FF derived by the fragment molecular orbital (FMO) method and by the restrained electrostatic potential (RESP) method. Applying this FMO‐RESP method to two proteins, dodecin, and lysozyme, we found that protein‐specific partial charges tend to differ more significantly from the standard AMBER charges for isolated charged atoms. We did not see the dependence of partial charges on the secondary structure. Computing the binding affinities of dodecin with five ligands by MM PBSA protocol with the FMO‐RESP charge set as well as with the standard AMBER charges, we found that the former gives better correlation with experimental affinities than the latter. While, for lysozyme with five ligands, both charge sets gave similar and relatively accurate estimates of binding affinities. © 2013 Wiley Periodicals, Inc.     

Steered molecular dynamics simulations of protein-ligand interactions

Molecular recognition and specific protein-ligandinteractions are central to many biochemical processes,such as enzyme catalysis, assembly of organelles, en-ergy transduction, signaling, diverse control functions,and replication, expression and storage of the geneticmaterial[1]. Moreover, protein-ligand interactions pro-vide the mechanism of many drug therapies and un-derstanding of such interactions is thus significant forrational drug design[1,2]. For the experimental studiesof protein-ligan…     

Efficient quantification of the importance of contacts for the dynamical stability of proteins

Understanding the stability of the native state and the dynamics of a protein is of great importance for all areas of biomolecular design. The efficient estimation of the influence of individual contacts between amino acids in a protein structure is a first step in the reengineering of a particular protein for technological or pharmacological purposes. At the same time, the functional annotation of molecular evolution can be facilitated by such insight. Here, we use a recently suggested, information theoretical measure in biomolecular design - the Kullback-Leibler-divergence - to quantify and therefore rank residue-residue contacts within proteins according to their overall contribution to the molecular mechanics. We implement this protocol on the basis of a reduced molecular model, which allows us to use a well-known lemma of linear algebra to speed up the computation. The increase in computational performance is around 10(1)- to 10(4)-fold. We applied the method to two proteins to illustrate the protocol and its results. We found that our method can reliably identify key residues in the molecular mechanics and the protein fold in comparison to well-known properties in the serine protease inhibitor. We found significant correlations to experimental results, e.g., dissociation constants and Φ values.     

Precision Molecular Threading/Dethreading

The general principles guiding the design of molecular machines based on interlocked structures are well known. Nonetheless, the identification of suitable molecular components for a precise tuning of the energetic parameters that determine the mechanical link is still challenging. Indeed, what are the reasons of the “all‐or‐nothing” effect, which turns a molecular “speed‐bump” into a stopper in pseudorotaxane‐based architectures? Here we investigate the threading and dethreading processes for a representative class of molecular components, based on symmetric dibenzylammonium axles and dibenzo[24]crown‐8 ether, with a joint experimental–computational strategy. From the analysis of quantitative data and an atomistic insight, we derive simple rules correlating the kinetic behaviour with the substitution pattern, and provide rational guidelines for the design of modules to be integrated in molecular switches and motors with sophisticated dynamic features.     

PRO_LIGAND: An approach to de novo molecular design. 6. Flexible fitting in the design of peptides

Summary This paper describes the further development of the functionality of our in-house de novo design program, PRO_LIGAND. In particular, attention is focussed on the implementation and validation of the directed tweak method for the construction of conformationally flexible molecules, such as peptides, from molecular fragments. This flexible fitting method is compared to the original method based on libraries of prestored conformations for each fragment. It is shown that the directed tweak method produces results of comparable quality, with significant time savings. By removing the need to generate a set of representative conformers for any new library fragment, the flexible fitting method increases the speed and simplicity with which new fragments can be included in a fragment library and also reduces the disk space required for library storage. A further improvement to the molecular construction process within PRO_LIGAND is the inclusion of a constrained minimisation procedure which relaxes fragments onto the design model and can be used to reject highly strained structures during the structure generation phase. This relaxation is shown to be very useful in simple test cases, but restricts diversity for more realistic examples. The advantages and disadvantages of these additions to the PRO_LIGAND methodology are illustrated by three examples: similar design to an alpha helix region of dihydrofolate reductase, complementary design to the active site of HIV-1 protease and similar design to an epitope region of lysozyme.     

Identification of Isomeric N-Glycans by Conformer Distribution Fingerprinting using Ion Mobility Mass Spectrometry

Glycans possess unparalleled structural complexity arising from chemically similar monosaccharide building blocks, configurations of anomeric linkages and different branching patterns, potentially giving rise to many isomers. This level of complexity is one of the main reasons that identification of exact glycan structures in biological samples still lags behind that of other biomolecules. Here, we introduce a methodology to identify isomeric N-glycans by determining gas phase conformer distributions (CDs) by measuring arrival time distributions (ATDs) using drift-tube ion mobility spectrometry-mass spectrometry. Key to the approach is the use of a range of well-defined synthetic glycans that made it possible to investigate conformer distributions in the gas phase of isomeric glycans in a systematic manner. In addition, we have computed CD fingerprints by molecular dynamics (MD) simulation, which compared well with experimentally determined CDs. It supports that ATDs resemble conformational populations in the gas phase and offer the prospect that such an approach can contribute to generating a library of CCS distributions (CCSDs) for structure identification.     

Molecular processes in an aqueous cluster

A cluster of 216 water molecules in a vacuum is simulated (molecular-dynamic trajectory is 40 ps); the average kinetic energy of water molecule movements corresponds to 300 K. The HLC model of data representation in time intervals of 0.082 ps is used to investigate the molecular scale processes. Volumetric and surface movements of water molecules are distinguished by the rank distributions based on the molecular movement character. The method for calculating volumetric and surface self-diffusion coefficients of molecules in a cluster is proposed. The cooperative character of movements of water molecules in a 5.4 Å-thick surface layer of the cluster is revealed.     

Development of a unique 3D interaction model of endogenous and synthetic peripheral benzodiazepine receptor ligands

Different classes of Peripheral-type Benzodiazepine Receptor (PBR) ligands were examined and common structural elements were detected and used to develop a rational binding model based on energetically allowed ligand conformations. Two lipophilic regions and one electrostatic interaction site are essential features for high affinity ligand binding, while a further lipophilic region plays an important modulator role. A comparative molecular field analysis, performed over 130 PBR ligands by means of the GRID/GOLPE methodology, led to a PLS model with both high fitting and predictive values (r² = 0.898, Q² = 0.761). The outcome from the 3D QSAR model and the GRID interaction fields computed on the putative endogenous PBR ligands DBI (Diazepam Binding Inhibitor) and TTN (Tetracontatetraneuropeptide) was used to identify the amino acids most probably involved in PBR binding. Three amino acids, bearing lipophilic side chains, were detected in DBI (Phe49, Leu47 and Met46) and in TTN (Phe33, Leu31 and Met30) as likely residues underlying receptor binding. Moreover, a qualitative comparison of the molecular electrostatic potentials of DBI, TTN and selected synthetic ligands indicated also similar electronic properties. Convergent results from the modeling studies of synthetic and endogenous ligands suggest a common binding mode to PBRs. This may help the rational design of new high affinity PBR ligands.     

MCDOCK: A Monte Carlo simulation approach to the molecular docking problem

Prediction of the binding mode of a ligand (a drug molecule) to its macromolecular receptor, or molecular docking, is an important problem in rational drug design. We have developed a new docking method in which a non-conventional Monte Carlo (MC) simulation technique is employed. A computer program, MCDOCK, was developed to carry out the molecular docking operation automatically. The current version of the MCDOCK program (version 1.0) allows for the full flexibility of ligands in the docking calculations. The scoring function used in MCDOCK is the sum of the interaction energy between the ligand and its receptor, and the conformational energy of the ligand. To validate the MCDOCK method, 19 small ligands, the binding modes of which had been determined experimentally using X-ray diffraction, were docked into their receptor binding sites. To produce statistically significant results, 20 MCDOCK runs were performed for each protein–ligand complex. It was found that a significant percentage of these MCDOCK runs converge to the experimentally observed binding mode. The root-mean-square (rms) of all non-hydrogen atoms of the ligand between the predicted and experimental binding modes ranges from 0.25 to 1.84 Å for these 19 cases. The computational time for each run on an SGI Indigo2/R10000 varies from less than 1 min to 15 min, depending upon the size and the flexibility of the ligands. Thus MCDOCK may be used to predict the precise binding mode of ligands in lead optimization and to discover novel lead compounds through structure-based database searching.     

分子印迹技术的回顾、现状与展望

分子印迹技术是近年来集高分子合成、分子设计、分子识别、仿生生物工程等众多学科优势发展起来的一门边缘学科分支。分子印迹聚合物由于具有与天然抗体同样的识别性能和与高分子同样的抗腐蚀性能的双重优点,因而广泛应用于生物工程、临床医学、环境监测、食品工业等众多领域。本文回顾了分子印迹技术近十多年来的发展过程,总结了目前的研究现状,并展望了分子印迹技术未来的发展趋势。     

Influence of Intermolecular Vibrations on the Electronic Coupling in Organic Semiconductors: The Case of Anthracene and Perfluoropentacene

We have performed classical molecular dynamics simulations and quantum‐chemical calculations on molecular crystals of anthracene and perfluoropentacene. Our goal is to characterize the amplitudes of the room‐temperature molecular displacements and the corresponding thermal fluctuations in electronic transfer integrals, which constitute a key parameter for charge transport in organic semiconductors. Our calculations show that the thermal fluctuations lead to Gaussian‐like distributions of the transfer integrals centered around the values obtained for the equilibrium crystal geometry. The calculated distributions have been plugged into Monte‐Carlo simulations of hopping transport, which show that lattice vibrations impact charge transport properties to various degrees depending on the actual crystal structure.     

Low‐molecular‐weight organogelators as functional materials for oil spill remediation

Oil spills from tankers are one of the major types of man‐made disasters that impact the marine environment, and they have been shown to have long‐lasting effects. On prevention of the spread of oil through rapid cleanup of spills, low‐molecular‐weight organogelators have received much attention because of their ability to tune their properties through rational design. In this mini‐review, I present a brief summary of studies focused on the remediation of oil spills via a chemical method, which involves the use of low‐molecular‐weight organogelators that form organogels with fuel oils or organic solvents. Moreover, recent attempts to create new improved molecular organogels composed of commercially available simple organogelators via a mixing induced enhancement method for solidifying oil are also discussed. In addition, polymer organogelators for oil spills are discussed in relation to low‐molecular weight gelators. Copyright © 2015 John Wiley & Sons, Ltd.     

Performance of multi-configurational calculations for a 1,4-bis(phenylethynyl)benzene derivative conjugated molecule.

The theoretical challenge of finding a single method that quantitatively reproduces both the experimental low-lying excitation energies and the torsional barrier of a prototypical conjugated molecule, which could act as a molecular wire, has been addressed here. The results indicate that this goal can be reasonably achieved when multi-reference perturbation theory up to second order (MRMP2) based on a complete active space self-consistent field (CASSCF) wave function using large active spaces is used. The results obtained were also used to compare with less expensive Kohn-Sham (KS) density functional theory (DFT) calculations when applied to these properties. The results obtained with BLYP and B3LYP exchange-correlation functionals indicate that quantitative agreement with all the experimental data cannot be obtained with this methodology, with a clear dependence on the exchange-correlation form selected. We thus encourage a careful testing of pure and hybrid density functionals whenever KS DFT is used for the rational design of conjugated materials for charge conduits.     

Multicomponent reaction design in the quest for molecular complexity and diversity

Multicomponent reactions have become increasingly popular as tools for the rapid generation of small-molecule libraries. However, to ensure sufficient molecular diversity and complexity, there is a continuous need for novel reactions. Although serendipity has always played an important role in the discovery of novel (multicomponent) reactions, rational design strategies have become much more important over the past decade. In this Review, we present an overview of general strategies that allow the design of novel multicomponent reactions. The challenges and opportunities for the future will be discussed.