COBRAMM 2.0 — A software interface for tailoring molecular electronic structure calculations and running nanoscale (QM/MM) simulations

We present a new version of the simulation software COBRAMM, a program package interfacing widely known commercial and academic software for molecular modeling. It allows a problem-driven tailoring of computational chemistry simulations with effortless ground and excited-state electronic structure computations. Calculations can be executed within a pure QM or combined quantum mechanical/molecular mechanical (QM/MM) framework, bridging from the atomistic to the nanoscale. The user can perform all necessary steps to simulate ground state and photoreactions in vacuum, complex biopolymer, or solvent environments. Starting from ground-state optimization, reaction path computations, initial conditions sampling, spectroscopy simulation, and photodynamics with deactivation events, COBRAMM is designed to assist in characterization and analysis of complex molecular materials and their properties. Interpretation of recorded spectra range from steady-state to time-resolved measurements. Various tools help the user to set up the system of interest and analyze the results.

[1]  W. Thiel,et al.  Adaptive time steps in trajectory surface hopping simulations. , 2016, The Journal of chemical physics.

[2]  S. Grimme,et al.  A COMBINATION OF KOHN-SHAM DENSITY FUNCTIONAL THEORY AND MULTI-REFERENCE CONFIGURATION INTERACTION METHODS , 1999 .

[3]  Francesco Aquilante,et al.  Multiconfigurational Second-Order Perturbation Theory with Frozen Natural Orbitals Extended to the Treatment of Photochemical Problems. , 2015, Journal of chemical theory and computation.

[4]  Michael J. Bearpark,et al.  A direct method for the location of the lowest energy point on a potential surface crossing , 1994 .

[5]  H. Wilks,et al.  Enzymology below 200 K: The kinetics and thermodynamics of the photochemistry catalyzed by protochlorophyllide oxidoreductase , 2002, Proceedings of the National Academy of Sciences of the United States of America.

[6]  L. Serrano-Andrés,et al.  Computation of conical intersections by using perturbation techniques. , 2005, The Journal of chemical physics.

[7]  A. Krylov,et al.  Calculations of non-adiabatic couplings within equation-of-motion coupled-cluster framework: Theory, implementation, and validation against multi-reference methods. , 2018, The Journal of chemical physics.

[8]  M. Olivucci,et al.  Comparison of the isomerization mechanisms of human melanopsin and invertebrate and vertebrate rhodopsins , 2014, Proceedings of the National Academy of Sciences.

[9]  C. Hunter,et al.  Identification and characterization of the product release steps within the catalytic cycle of protochlorophyllide oxidoreductase. , 2004, Biochemistry.

[10]  U. Rothlisberger,et al.  Mixed Quantum Mechanical/Molecular Mechanical Molecular Dynamics Simulations of Biological Systems in Ground and Electronically Excited States. , 2015, Chemical reviews.

[11]  C. Hunter,et al.  Protochlorophyllide oxidoreductase: "dark" reactions of a light-driven enzyme. , 2003, Biochemistry.

[12]  R. Improta,et al.  UV-induced long-lived decays in solvated pyrimidine nucleosides resolved at the MS-CASPT2/MM level. , 2018, Physical chemistry chemical physics : PCCP.

[13]  O. Weingart,et al.  Spectral properties and isomerisation path of retinal in C1C2 channelrhodopsin. , 2015, Physical chemistry chemical physics : PCCP.

[14]  M. W. van der Kamp,et al.  Combined quantum mechanics/molecular mechanics (QM/MM) methods in computational enzymology. , 2013, Biochemistry.

[15]  H. Senn,et al.  QM/MM Methods for Biological Systems , 2006 .

[16]  M. Olivucci,et al.  Relationship between Excited State Lifetime and Isomerization Quantum Yield in Animal Rhodopsins: Beyond the One-Dimensional Landau-Zener Model. , 2018, The journal of physical chemistry letters.

[17]  Giulio Cerullo,et al.  Ab Initio Simulations of Two-Dimensional Electronic Spectra: The SOS//QM/MM Approach , 2014 .

[18]  M. Garavelli,et al.  Modelling vibrational coherence in the primary rhodopsin photoproduct. , 2012, The Journal of chemical physics.

[19]  Bernard R. Brooks,et al.  ENZYME MECHANISMS WITH HYBRID QUANTUM AND MOLECULAR MECHANICAL POTENTIALS.I. THEORETICAL CONSIDERATIONS , 1996 .

[20]  Chengyin Wu,et al.  Decay Pathways of Pyrimidine Bases: From Gas Phase to Solution , 2008 .

[21]  J. Watts.  An Introduction to Equation-of-Motion and Linear-ResponseCoupled-Cluster Methods for Electronically Excited Statesof Molecules , 2008 .

[22]  Arieh Warshel,et al.  Bicycle-pedal model for the first step in the vision process , 1976, Nature.

[23]  Jerzy Leszczynski,et al.  Radiation Induced Molecular Phenomena in Nucleic Acids , 2008 .

[24]  W. Thiel,et al.  Implementation of surface hopping molecular dynamics using semiempirical methods , 2008 .

[25]  Julian D. Gale,et al.  The General Utility Lattice Program (GULP) , 2003 .

[26]  S. Mukamel,et al.  UV-Light-Induced Vibrational Coherences: The Key to Understand Kasha Rule Violation in trans-Azobenzene. , 2018, The journal of physical chemistry letters.

[27]  D. Truhlar,et al.  Non-Born-Oppenheimer Liouville-von Neumann Dynamics. Evolution of a Subsystem Controlled by Linear and Population-Driven Decay of Mixing with Decoherent and Coherent Switching. , 2005, Journal of chemical theory and computation.

[28]  T. Martínez,et al.  Conical intersection dynamics in solution: the chromophore of Green Fluorescent Protein. , 2004, Faraday discussions.

[29]  Hans Lischka,et al.  Newton‐X: a surface‐hopping program for nonadiabatic molecular dynamics , 2014 .

[30]  Qi Ou,et al.  First-order derivative couplings between excited states from adiabatic TDDFT response theory. , 2015, The Journal of chemical physics.

[31]  D. Truhlar,et al.  QM/MM: what have we learned, where are we, and where do we go from here? , 2007 .

[32]  L. Verlet.  Computer "Experiments" on Classical Fluids. I. Thermodynamical Properties of Lennard-Jones Molecules , 1967 .

[33]  D. Arendt.  Evolution of eyes and photoreceptor cell types. , 2003, The International journal of developmental biology.

[34]  Walter Thiel,et al.  Photoinduced nonadiabatic dynamics of pyrimidine nucleobases: on-the-fly surface-hopping study with semiempirical methods. , 2009, The journal of physical chemistry. B.

[35]  J. Cadet,et al.  Benzophenone photosensitized DNA damage. , 2012, Accounts of chemical research.

[36]  N. Scrutton,et al.  Stepwise Hydride Transfer in a Biological System: Insights into the Reaction Mechanism of the Light‐Dependent Protochlorophyllide Oxidoreductase , 2018, Angewandte Chemie.

[37]  Walter Thiel,et al.  OMx-D: semiempirical methods with orthogonalization and dispersion corrections. Implementation and biochemical application. , 2008, Physical chemistry chemical physics : PCCP.

[38]  P. M. Rodger,et al.  DL_POLY: Application to molecular simulation , 2002 .

[39]  R. Mathies,et al.  Wavepacket splitting and two-pathway deactivation in the photoexcited visual pigment isorhodopsin. , 2014, Angewandte Chemie.

[40]  S. Mukamel,et al.  Two-Dimensional Electronic Spectroscopy of Benzene, Phenol, and Their Dimer: An Efficient First-Principles Simulation Protocol. , 2015, Journal of chemical theory and computation.

[41]  Alexey A. Sokol,et al.  ChemShell—a modular software package for QM/MM simulations , 2014 .

[42]  Weitao Yang,et al.  Electron transport through single conjugated organic molecules: basis set effects in ab initio calculations. , 2007, The Journal of chemical physics.

[43]  Roland Lindh,et al.  The ultrafast photoisomerizations of rhodopsin and bathorhodopsin are modulated by bond length alternation and HOOP driven electronic effects. , 2011, Journal of the American Chemical Society.

[44]  S. Hammes-Schiffer,et al.  Proton transfer in solution: Molecular dynamics with quantum transitions , 1994 .

[45]  Felipe Zapata,et al.  Molcas 8: New capabilities for multiconfigurational quantum chemical calculations across the periodic table , 2016, J. Comput. Chem..

[46]  W. Hase,et al.  On the dynamics of state selected unimolecular reactions: Chloroacetylene dissociation and predissociation , 1977 .

[47]  M. Barbatti,et al.  Nonadiabatic Photodynamics of a Retinal Model in Polar and Nonpolar Environment , 2013, The journal of physical chemistry. A.

[48]  Paul K. Weiner,et al.  Ground states of molecules , 1972 .

[49]  Z. Lan,et al.  An On-the-Fly Surface-Hopping Program JADE for Nonadiabatic Molecular Dynamics of Polyatomic Systems: Implementation and Applications. , 2015, Journal of chemical theory and computation.

[50]  H. C. Andersen.  Rattle: A “velocity” version of the shake algorithm for molecular dynamics calculations , 1983 .

[51]  J. Tully,et al.  Trajectory Surface Hopping Approach to Nonadiabatic Molecular Collisions: The Reaction of H+ with D2 , 1971 .

[52]  I. Rivalta,et al.  The effect of solvent relaxation in the ultrafast time-resolved spectroscopy of solvated benzophenone , 2018, Photochemical & photobiological sciences : Official journal of the European Photochemistry Association and the European Society for Photobiology.

[53]  S. Mukamel.  Principles of Nonlinear Optical Spectroscopy , 1995 .

[54]  R. Lindh,et al.  Low-cost evaluation of the exchange Fock matrix from Cholesky and density fitting representations of the electron repulsion integrals. , 2007, The Journal of chemical physics.

[55]  F. Maseras,et al.  QM/MM methods in inorganic chemistry. , 2008, Dalton transactions.

[56]  A. Clarke,et al.  Protochlorophyllide oxidoreductase: A homology model examined by site‐directed mutagenesis , 2001, Proteins.

[57]  I. V. van Stokkum,et al.  Protochlorophyllide excited-state dynamics in organic solvents studied by time-resolved visible and mid-infrared spectroscopy. , 2010, The journal of physical chemistry. B.

[58]  J. Tully.  Molecular dynamics with electronic transitions , 1990 .

[59]  Walter Thiel,et al.  Ground States of Molecules. 38. The MNDO Method. Approximations and Parameters , 1977 .

[60]  Tsutomu Kouyama,et al.  Crystal structure of squid rhodopsin , 2008, Nature.

[61]  J. Herbert,et al.  Analytic derivative couplings for spin-flip configuration interaction singles and spin-flip time-dependent density functional theory. , 2014, The Journal of chemical physics.

[62]  Francesco Aquilante,et al.  Cholesky Decomposition-Based Multiconfiguration Second-Order Perturbation Theory (CD-CASPT2): Application to the Spin-State Energetics of Co(III)(diiminato)(NPh). , 2008, Journal of chemical theory and computation.

[63]  Oliver Weingart,et al.  The role of HOOP-modes in the ultrafast photo-isomerization of retinal models , 2008 .

[64]  I. Rivalta,et al.  Resolving the Singlet Excited State Manifold of Benzophenone by First-Principles Simulations and Ultrafast Spectroscopy. , 2018, Journal of chemical theory and computation.

[65]  Hinayat Sawhney.  IMPLEMENTATION AND APPLICATIONS OF DATA MIN , 2017 .

[66]  Mickaël G. Delcey,et al.  Benzophenone Ultrafast Triplet Population: Revisiting the Kinetic Model by Surface-Hopping Dynamics , 2016, The journal of physical chemistry letters.

[67]  M. Alexandre,et al.  Enzyme activation and catalysis: characterisation of the vibrational modes of substrate and product in protochlorophyllide oxidoreductase. , 2011, Physical chemistry chemical physics : PCCP.

[68]  T. Begley,et al.  Protochlorophyllide reductase. 1. Determination of the regiochemistry and the stereochemistry of the reduction of protochlorophyllide to chlorophyllide , 1989 .

[69]  Frank Neese,et al.  The ORCA program system , 2012 .

[70]  R. Improta,et al.  Resolving Ultrafast Photoinduced Deactivations in Water-Solvated Pyrimidine Nucleosides. , 2017, The journal of physical chemistry letters.

[71]  R. van Grondelle,et al.  The First Catalytic Step of the Light-driven Enzyme Protochlorophyllide Oxidoreductase Proceeds via a Charge Transfer Complex* , 2006, Journal of Biological Chemistry.

[72]  T. Martínez,et al.  Ab Initio Multiple Spawning: Photochemistry from First Principles Quantum Molecular Dynamics , 2000 .

[73]  Mixed quantumñclassical dynamics , 1998 .

[74]  Bern Kohler,et al.  Internal conversion to the electronic ground state occurs via two distinct pathways for pyrimidine bases in aqueous solution , 2007, Proceedings of the National Academy of Sciences.

[75]  Q. Cui,et al.  A critical evaluation of different QM/MM frontier treatments with SCC-DFTB as the QM method. , 2005, The journal of physical chemistry. B.

[76]  A. Klamt.  Conductor-like Screening Model for Real Solvents: A New Approach to the Quantitative Calculation of Solvation Phenomena , 1995 .

[77]  J. Park,et al.  Analytical Derivative Coupling for Multistate CASPT2 Theory. , 2017, Journal of chemical theory and computation.

[78]  John M Herbert,et al.  Simultaneous benchmarking of ground- and excited-state properties with long-range-corrected density functional theory. , 2008, The Journal of chemical physics.

[79]  T. Martínez,et al.  Dynamical Correlation Effects on Photoisomerization: Ab Initio Multiple Spawning Dynamics with MS-CASPT2 for a Model trans-Protonated Schiff Base. , 2016, The journal of physical chemistry. B.

[80]  Roland Lindh,et al.  Analytical State-Average Complete-Active-Space Self-Consistent Field Nonadiabatic Coupling Vectors: Implementation with Density-Fitted Two-Electron Integrals and Application to Conical Intersections. , 2016, Journal of chemical theory and computation.

[81]  G. Granucci,et al.  Critical appraisal of the fewest switches algorithm for surface hopping. , 2007, The Journal of chemical physics.

[82]  Jianpeng Ma,et al.  CHARMM: The biomolecular simulation program , 2009, J. Comput. Chem..

[83]  Bipin K. Shah,et al.  The S2 → S1 Internal Conversion of Benzophenone and p-Iodobenzophenone1 , 2004 .

[84]  W. Thiel,et al.  Nonradiative deexcitation dynamics of 9H-adenine: an OM2 surface hopping study. , 2008, The journal of physical chemistry. A.

[85]  Walter Thiel,et al.  Implementation of a general multireference configuration interaction procedure with analytic gradients in a semiempirical context using the graphical unitary group approach , 2003, J. Comput. Chem..

[86]  Hans-Joachim Werner,et al.  Explicitly Correlated Local Electron Correlation Methods , 2017 .

[87]  S. Mukamel,et al.  Bidimensional electronic spectroscopy on indole in gas phase and in water from first principles , 2014 .

[88]  J. Herbert,et al.  Time-resolved infrared spectroscopy of the lowest triplet state of thymine and thymidine. , 2008, Chemical physics.

[89]  Edward F. Valeev,et al.  Explicitly correlated R12/F12 methods for electronic structure. , 2012, Chemical reviews.

[90]  Marco Garavelli,et al.  A tunable QM/MM approach to chemical reactivity, structure and physico-chemical properties prediction , 2007 .

[91]  L. Chambers.  Practical methods of optimization (2nd edn) , by R. Fletcher. Pp. 436. £34.95. 2000. ISBN 0 471 49463 1 (Wiley). , 2001, The Mathematical Gazette.

[92]  Giulio Cerullo,et al.  Modeling the high-energy electronic state manifold of adenine: Calibration for nonlinear electronic spectroscopy. , 2015, The Journal of chemical physics.

[93]  Chris Oostenbrink,et al.  Testing of the GROMOS Force-Field Parameter Set 54A8: Structural Properties of Electrolyte Solutions, Lipid Bilayers, and Proteins , 2013, Journal of chemical theory and computation.

[94]  R. Improta,et al.  Quantum Mechanical Studies on the Photophysics and the Photochemistry of Nucleic Acids and Nucleobases. , 2016, Chemical reviews.

[95]  R. Friesner,et al.  Ab initio quantum chemical and mixed quantum mechanics/molecular mechanics (QM/MM) methods for studying enzymatic catalysis. , 2005, Annual review of physical chemistry.

[96]  Walter Thiel,et al.  Semiempirical Quantum-Chemical Orthogonalization-Corrected Methods: Benchmarks of Electronically Excited States. , 2016, Journal of chemical theory and computation.

[97]  M. Levitt,et al.  Theoretical studies of enzymic reactions: dielectric, electrostatic and steric stabilization of the carbonium ion in the reaction of lysozyme. , 1976, Journal of molecular biology.

[98]  K. Morokuma,et al.  A NEW ONIOM IMPLEMENTATION IN GAUSSIAN98. PART I. THE CALCULATION OF ENERGIES, GRADIENTS, VIBRATIONAL FREQUENCIES AND ELECTRIC FIELD DERIVATIVES , 1999 .

[99]  J. Ponder,et al.  An efficient newton‐like method for molecular mechanics energy minimization of large molecules , 1987 .

[100]  R. Mathies,et al.  Conical intersection dynamics of the primary photoisomerization event in vision , 2010, Nature.

[101]  A. Bottoni,et al.  Product formation in rhodopsin by fast hydrogen motions. , 2011, Physical chemistry chemical physics : PCCP.

[102]  Giulio Cerullo,et al.  Disentangling Peptide Configurations via Two-Dimensional Electronic Spectroscopy: Ab Initio Simulations Beyond the Frenkel Exciton Hamiltonian , 2014, The journal of physical chemistry letters.

[103]  O. Weingart.  Combined Quantum and Molecular Mechanics (QM/MM) Approaches to Simulate Ultrafast Photodynamics in Biological Systems , 2017 .

[104]  Walter Thiel,et al.  QM/MM study of the absorption spectra of DsRed.M1 chromophores , 2010, J. Comput. Chem..

[105]  S. Mukamel,et al.  Towards Accurate Simulation of Two-Dimensional Electronic Spectroscopy , 2018, Topics in Current Chemistry.

[106]  Nonadiabatic Dynamics: Mean-Field and Surface Hopping , 2002 .

[107]  Giulio Cerullo,et al.  Tracking conformational dynamics of polypeptides by nonlinear electronic spectroscopy of aromatic residues: a first-principles simulation study. , 2014, Chemphyschem : a European journal of chemical physics and physical chemistry.

[108]  P. Szalay,et al.  Analytic evaluation of the nonadiabatic coupling vector between excited states using equation-of-motion coupled-cluster theory. , 2009, The Journal of chemical physics.

[109]  H. Schlegel,et al.  Finding minima, transition states, and following reaction pathways on ab initio potential energy surfaces , 2005 .

[110]  J. Tully.  Mixed quantum–classical dynamics , 1998 .

[111]  Walter Thiel,et al.  Beyond the MNDO model: Methodical considerations and numerical results , 1993, J. Comput. Chem..

[112]  R. Fletcher.  Practical Methods of Optimization , 1988 .

[113]  Klaus Schulten,et al.  Photochemical reaction dynamics of the primary event of vision studied by means of a hybrid molecular simulation. , 2009, Biophysical journal.

[114]  Berk Hess,et al.  GROMACS: High performance molecular simulations through multi-level parallelism from laptops to supercomputers , 2015 .

[115]  Hans Lischka,et al.  Optimization of mixed quantum-classical dynamics: Time-derivative coupling terms and selected couplings , 2009 .

[116]  Walter Thiel,et al.  Benchmark of Electronically Excited States for Semiempirical Methods: MNDO, AM1, PM3, OM1, OM2, OM3, INDO/S, and INDO/S2. , 2010, Journal of chemical theory and computation.