Molecular Modeling and Nuclear Overhauser Enhancement Spectroscopy (NOESY): Tools for Studying the Regioselective Bromination of 3-Bromoanisole

The use of molecular modeling for predicting chemical reactivity has been highly successful in the industrial and academic research communities. For this reason, increased emphasis has been placed on molecular modeling in the undergraduate curriculum. In the described experiment, the bromination of 3-bromoanisole, students are encouraged to use molecular modeling software as a tool for predicting chemical reactivity. Besides introducing students to molecular modeling, this experiment incorporates the use of nontraditional, less hazardous reagents and solvents for electrophilic aromatic bromination reactions. Lastly, nuclear Overhauser enhancement spectroscopy (NOESY) is introduced as a tool for structural elucidation. Although there are a number of aspects to this experiment, two 3-hour laboratory periods are sufficient because the results from semiempirical (AMI) geometry optimizations, which are complete in seconds, were almost identical to the higher order, more time-intensive ab initio (3-21G*) calculations. In addition, the experimental time was greatly shortened by the discovery that catalytic HCl(aq) reduces the reaction time from 5 hours to 18 minutes.     

Using Molecular Modeling in the Organic Chemistry Course for Majors

Molecular modeling provides a way to correlate theoretical concepts with experimental data; therefore, we have introduced organic chemistry students to molecular modeling early in the first semester. This approach provides students with additional skills for clarifying chemical and theoretical concepts by means of demonstrations in the classroom and hands-on tutorial modules. In this manner the impact of the active-learning process is increased. In addition, this tool allows us to further enhance laboratory experiments already developed using a guided-inquiry approach and to design new experiments. Chemical concepts such as conformational analysis, stereochemistry, IR spectra, molecular and electronic properties, molecular orbitals, and chemical reactivity are emphasized through this approach.     

VRML在分子对称性教学中的应用

以椅式环己烷为例,介绍了如何应用分子模拟软件制作含有分子对称元素的VRML模型的过程。利用VRML的交互性可以帮助学生直观地掌握分子的三维结构以及对称元素的分布,切实提高了分子对称性的教学效果。     

结构化学课程中波函数3D图的VRML实现

Using three-dimensional (3D) graphing and molecular modeling softwares, the procedures for presenting the VRML 3D wave functions have been discussed. The interactivities of VRML can help students to observe the characters of wave functions intuitively, which also improve the teaching results of quantum chemistry and electronic structures effectively.     

Integrating Molecular Modeling with Semiempirical Quantum Mechanical Calculations into the Upper-Level Inorganic Chemistry Course

The analysis of chemical bonding and reactivity from the perspective of molecular orbital theory is challenging for students at the undergraduate level. In an attempt to improve the instruction of this material in my upper-level inorganic chemistry course I developed a series of computational experiments using a molecular modeling program that can perform semiempirical quantum mechanical calculations. These exercises explore the chemistry of molecular systems through an analysis of the variation in the attractive and repulsive forces in the system as a function of structure or composition. The exercises challenge the analysis skills of the students by requiring them to consider how two or more factors contribute to the properties of the system. Examples of exercises that demonstrate different types of computational experiments are given. These sample exercises examine the structure of simple molecules, the reactivity of Lewis acids, and the bonding in transition metal complexes.     

Complex molecular assemblies at hand via interactive simulations

Studying complex molecular assemblies interactively is becoming an increasingly appealing approach to molecular modeling. Here we focus on interactive molecular dynamics (IMD) as a textbook example for interactive simulation methods. Such simulations can be useful in exploring and generating hypotheses about the structural and mechanical aspects of biomolecular interactions. For the first time, we carry out low‐resolution coarse‐grain IMD simulations. Such simplified modeling methods currently appear to be more suitable for interactive experiments and represent a well‐balanced compromise between an important gain in computational speed versus a moderate loss in modeling accuracy compared to higher resolution all‐atom simulations. This is particularly useful for initial exploration and hypothesis development for rare molecular interaction events. We evaluate which applications are currently feasible using molecular assemblies from 1900 to over 300,000 particles. Three biochemical systems are discussed: the guanylate kinase (GK) enzyme, the outer membrane protease T and the soluble N‐ethylmaleimide‐sensitive factor attachment protein receptors complex involved in membrane fusion. We induce large conformational changes, carry out interactive docking experiments, probe lipid–protein interactions and are able to sense the mechanical properties of a molecular model. Furthermore, such interactive simulations facilitate exploration of modeling parameters for method improvement. For the purpose of these simulations, we have developed a freely available software library called MDDriver. It uses the IMD protocol from NAMD and facilitates the implementation and application of interactive simulations. With MDDriver it becomes very easy to render any particle‐based molecular simulation engine interactive. Here we use its implementation in the Gromacs software as an example. © 2009 Wiley Periodicals, Inc. J Comput Chem, 2009     

Probing molecular shape. 1. Conformational studies of 5-hydroxyhexahydropyrimidine and related compounds

Understanding the factors that determine molecular shape enables scientists to begin to understand and tailor molecular properties and reactivity. Many biomolecules and bioactive compounds contain aliphatic heterocyclic rings whose conformations play a major role in their biological activity. The interplay of a number of factors, both steric and electronic, is examined for 5-hydroxyhexahydropyrimidine (1) and related compounds with use of spectroscopy and molecular modeling.     

The Dienone–Phenol Rearrangement of Santonin: A Comprehensive Laboratory Experiment for Advanced Undergraduate Organic Chemistry Students

A comprehensive laboratory experiment suitable for advanced undergraduate organic chemistry students has been designed. The experiment is based on the dienone-phenol rearrangement reaction of the sesquiterpene santonin to give -desmotroposantonin acetate. It challenges students to solve an earlier controversial stereochemical problem. The students carry out the reaction and analyze spectroscopic data to determine the stereochemistry of the starting material and the product. In addition, they perform simple molecular modeling calculations, which enable them to rationalize the stereochemical outcome of the transformation and discuss the mechanism of the dienone-phenol rearrangement and related rearrangements reported in the literature.     

Al3+ induced planarization, conformational arrest and metallochromic shift in a pyrimidine dione dye: insight from integrated hybrid quantum-classical calculations

In order to explore the possibilities of simulating metallochromism by modern molecular modeling, we apply a sequential hybrid quantum-classical approach to a prototype metallochromic system-the Al(3+) ion and pyrimidinedione (PY) dye complex. The complex shows several structural features with relevance for the metallochromism: the PY dye exhibits conformers with dynamical transitions between twisted structures, which are inhibited by the addition of the metal ion leading to planarization and a conformational arrest: the Al(3+) ion behaves like a structure-modifier for both intra and intermolecular degrees of freedom and with respect to the intermolecular solvation shell structure. The sequential approach that we have employed uses DFT/MM molecular dynamics for structure modeling and TDDFT/PCM for property modeling. The computed metallochromic shift between PY and the Al(PY)(3+) complex in DMSO solvent is obtained in excellent agreement with experiment. The results infer optimism for future use of such modeling techniques to design metallochromic indicators.     

Macromodel—an integrated software system for modeling organic and bioorganic molecules using molecular mechanics

An integrated molecular modeling system for designing and studying organic and bioorganic molecules and their molecular complexes using molecular mechanics is described. The graphically controlled, atom-based system allows the construction, display and manipulation of molecules and complexes having as many as 10,000 atoms and provides interactive, state-of-the-art molecular mechanics on any subset of up to 1,000 atoms. The system semiautomates the graphical construction and analysis of complex structures ranging from polycyclic organic molecules to biopolymers to mixed molecular complexes. We have placed emphasis on providing effective searches of conformational space by a number of different methods and on highly optimized molecular mechanics energy calculations using widely used force fields which are supplied as external files. Little experience is required to operate the system effectively and even novices can use it to carry out sophisticated modeling operations. The software has been designed to run on Digital Equipment Corporation VAX computers interfaced to a variety of graphics devices ranging from inexpensive monochrome terminals to the sophisticated graphics displays of the Evans & Sutherland PS300 series.     

Ligand-guided optimization of CXCR4 homology models for virtual screening using a multiple chemotype approach

CXCR4 is a G-protein coupled receptor for CXCL12 that plays an important role in human immunodeficiency virus infection, cancer growth and metastasization, immune cell trafficking and WHIM syndrome. In the absence of an X-ray crystal structure, theoretical modeling of the CXCR4 receptor remains an important tool for structure–function analysis and to guide the discovery of new antagonists with potential clinical use. In this study, the combination of experimental data and molecular modeling approaches allowed the development of optimized ligand-receptor models useful for elucidation of the molecular determinants of small molecule binding and functional antagonism. The ligand-guided homology modeling approach used in this study explicitly re-shaped the CXCR4 binding pocket in order to improve discrimination between known CXCR4 antagonists and random decoys. Refinement based on multiple test-sets with small compounds from single chemotypes provided the best early enrichment performance. These results provide an important tool for structure-based drug design and virtual ligand screening of new CXCR4 antagonists.     

Visualizing The Future of Molecular Graphics

Abstract

The field of computer graphics has played an important role in the advancement of structural molecular biology and in the development of structure-based drug design. This article will provide a brief background on the development of this technology, and then focus on the current trends and future directions in molecular graphics and how they will impact the practice of molecular modeling and design. Specific areas that will be covered include: 1) the development of surface and volume based representations of molecular properties and interactions; 2) new approaches to modeling flexible and multi-component structures, and 3) the impact of object-oriented graphics-based programming and the rapidly growing use of network based computing.     

MBO(N)D: A multibody method for long‐time molecular dynamics simulations

A modeling approach that can significantly speed up the dynamics simulation of large molecular systems is presented herein. A multigranular modeling approach, whereby different parts of the molecule are modeled at different levels of detail, is enabled by substructuring. Substructuring the molecular system is accomplished by collecting groups of atoms into rigid or flexible bodies. Body flexibility is modeled by a truncated set of body‐based modes. This approach allows for the elimination of the high‐frequency harmonic motion while capturing the low‐frequency anharmonic motion of interest. This results in the use of larger integration step sizes, substantially reducing the computational time required for a given dynamic simulation. The method also includes the use of a multiple time scale (MTS) integration scheme. Speed increases of 5‐ to 30‐fold over atomistic simulations have been realized in various applications of the method. © 2000 John Wiley & Sons, Inc. J Comput Chem 21: 159–184, 2000     

Bayesian inference of conformational state populations from computational models and sparse experimental observables

We present a Bayesian inference approach to estimating conformational state populations from a combination of molecular modeling and sparse experimental data. Unlike alternative approaches, our method is designed for use with small molecules and emphasizes high‐resolution structural models, using inferential structure determination with reference potentials, and Markov Chain Monte Carlo to sample the posterior distribution of conformational states. As an application of the method, we determine solution‐state conformational populations of the 14‐membered macrocycle cineromycin B, using a combination of previously published sparse Nuclear Magnetic Resonance (NMR) observables and replica‐exchange molecular dynamic/Quantum Mechanical (QM)‐refined conformational ensembles. Our results agree better with experimental data compared to previous modeling efforts. Bayes factors are calculated to quantify the consistency of computational modeling with experiment, and the relative importance of reference potentials and other model parameters. © 2014 Wiley Periodicals, Inc.     

Accurate and interpretable computational modeling of chemical mutagenicity

We describe a method for modeling chemical mutagenicity in terms of simple rules based on molecular features. A classification model was built using a rule-based ensemble method called RuleFit, developed by Friedman and Popescu. We show how performance compares favorably against literature methods. Performance was measured through the use of cross-validation and testing on external test sets. All data sets used are publicly available. The method automatically generated transparent rules in terms of molecular structure that agree well with known toxicology. While we have focused on chemical mutagenicity in demonstrating this method, we anticipate that it may be more generally useful in modeling other molecular properties such as other types of chemical toxicity.     

Homology model directed alignment selection for comparative molecular field analysis: Application to photosystem II inhibitors

The use of a computational docking protocol in conjunction with a protein homology model to derive molecular alignments for Comparative Molecular Field Analysis (CoMFA) was examined. In particular, the DOCK program and a model of the herbicidal target site, photosystem II (PSII), was used to derive alignments for two PSII inhibitor training sets, a set of benzo- and napthoquinones and a set of butenanilides. The protein design software in the QUANTA molecular modeling package was used to develop a homology model of spinach PSII based on the reported amino acid sequence and the X-ray crystal structure of the purple bacterium reaction center. The model is very similar to other reported PSII protein homology models. DOCK was then used to derive alignments for CoMFA modeling by docking the inhibitors in the PSII binding pocket. The molecular alignments produced from docking yielded highly predictive CoMFA models. As a comparison, the more traditional atom-atom alignments of the same two training sets failed to produce predictive CoMFA models. The general utilities of this application for homology model refinement and as an alternative scoring method are discussed.     

Development of force field parameters for molecular simulation of polylactide

Polylactide is a biodegradable polymer that is widely used for biomedical applications, and it is a replacement for some petroleum based polymers in applications that range from packaging to carpeting. Efforts to characterize and further enhance polylactide based systems using molecular simulations have to this point been hindered by the lack of accurate atomistic models for the polymer. Thus, we present force field parameters specifically suited for molecular modeling of PLA. The model, which we refer to as PLAFF3, is based on a combination of the OPLS and CHARMM force fields, with modifications to bonded and nonbonded parameters. Dihedral angle parameters were adjusted to reproduce DFT data using newly developed CMAP dihedral cross terms, and the model was further adjusted to reproduce experimentally resolved crystal structure conformations, melt density, volume expansivity, and the glass transition temperature of PLA. We recommend the use of PLAFF3 in modeling PLA in its crystalline or amorphous states and have provided the necessary input files required for the publicly available molecular dynamics code GROMACS.     

Synthesis,properties and molecular modeling of functional hyperbranched polymers and dendrimers

The design, synthesis, properties and molecular modeling of fully conjugated dendritic molecules and conjugated hyperbranched polymers are described. It has been shown that conjugated hyperbranched molecules are much more soluble than their linear analogues while maintaining all the properties characteristic of conjugated polymers. It was found that the use of polymeric solid support in hyperbranched polymerization allows to control molecular weight and degree of branching (DB). The molecular modeling of hyperbranched conjugated molecules reveals that hyperbranched structure of conjugated molecules affects significantly neither their stability nor the conjugation. On the other hand the terminal groups affect appreciably the electronic structure of conjugated hyperbranched molecules.     

Determining absolute configuration in flexible molecules: a case study

Assigning absolute configuration of molecules continues to be a major problem. Determining absolute configuration in conformationally flexible systems is challenging, even for experts. Here, we present a case study in which we use a combination of molecular modeling, solution NMR, and X-ray crystallography to illustrate why it is difficult to use solution methods alone for configuration assignment. For the case examined, a comparison of calculated and experimental optical rotatory dispersion (ORD) data provides the most straightforward way to assign the absolute configuration.