MM2 Force Field Research Articles

Avoiding Electron Spill-Out in QM/MM Calculations on Excited States with Simple Pseudopotentials.

QM/MM calculations of electronic excitations with diffuse basis sets often have large errors due to spill-out of electrons from the quantum subsystem. The Pauli repulsion of the electrons by the environment has to be included to avoid this. We propose transferable atomic all-electron pseudopotentials that can readily be combined with most MM force fields to avoid electron spill-out. QM/MM excitation energies computed with time-dependent Hartree-Fock and the algebraic diagrammatic construction through second-order are benchmarked against supermolecular calculations to validate these new pseudopotentials. The QM/MM calculations with pseudopotentials give accurate results that are stable with augmentation of the basis set with diffuse functions. We show that the largest contribution to residual deviations from full QM calculations is caused by the missing London dispersion interaction.

Topological Analysis of Electron Density in Large Biomolecular Systems.

A great step toward describing the structure of the molecular electron was made in the era of quantum chemical methods. Methods play a very important role in the prediction of molecular properties and in the description of the reactivity of compounds, which cannot be overestimated. There are many works, books, and articles on quantum methods, their applications, and comparisons. At the same time, quantum methods of a high level of theory, which give the most accurate results, are time-consuming, which makes them almost impossible to describe large complex molecular systems, such as macromolecules, enzymes, supramolecular compounds, crystal fragments, and so on. To propose an approach that allows real-time estimation of electron density in large systems, such as macromolecules, nanosystems, proteins. AlteQ approach was applied to the tolopogical analysis of electron density for "substrate - cytochrome" complexes. The approach is based on the use of Slater's type atomic contributions. Parameters of the atomic contributions were found using high resolution X-ray diffraction data for organic and inorganic molecules. Relationships of the parameters with atomic number, ionization potentials and electronegativities were determined. The sufficient quality of the molecular electron structure representation was shown under comparison of AlteQ predicted and observed electron densities. AlteQ algorithm was applied for evaluation of electron structure of "CYP3A4 - substrate" complexes modeled using BiS/MC restricted docking procedure. Topological analysis (similar to Atoms In Molecules (AIM) theory suggested by Richard F.W. Bader) of the AlteQ molecular electron density was carried out for each complex. The determination of (3,-1) bond, (3,+1) ring, (3,+3) cage critical points of electron density in the intermolecular "CYP3A4 - substrate" space was performed. Different characteristics such as electron density, Laplacian eigen values, etc. at the critical points were computed. Relationship of pKM (KM is Michaelis constant) with the maximal value of the second Laplacian eigen value of electron density at the critical points and energy of complex formation computed using MM3 force field was determined. It was shown that significant number of (3,-1) bond critical points are located in the intermolecular space between the enzyme site and groups of substrate atoms eliminating during metabolism processes.

Efficiently Transplanting Potential Energy Interpolation Database between Two Systems: Bacteriochlorophyll Case with FMO and LH2 Complexes.

Generating a reliable potential energy surface (PES) is an important issue for studying the dynamics of any system with computational simulations. Interpolation mechanics/molecular mechanics (IM/MM) based on a PES interpolation scheme is a useful tool in that regard as it provides an accuracy of a quantum chemistry (QC) level while maintaining its computational cost comparable to conventional MM force fields. Despite this benefit, constructing the database for interpolation itself is still challenging and time-consuming. Here, we present a method with which we can construct the IM database of one system based on a preexisting data set for another related system. We adopt the case of constructing bacteriochlorophyll PESs for the light-harvesting 2 (LH2) complex by utilizing already available IM database for the BChls from the Fenna-Matthews-Olson (FMO) complex. In this method, the IM database from FMO is first transplanted to LH2 by considering BChl displacement vectors that take into account the geometry differences induced by the protein scaffolds. From this transplanted primitive database entries, a relatively small number of effective ones are selected by a survival process based on a genetic algorithm such that the IM energies evaluated at geometries in a conveniently collected prediction set can closely match with the reference QC energies. The selection process is expedited by using two different levels of basis sets for the QC calculations. To demonstrate the utility of the PES thus constructed, we carry out 1 ns of IM/MM dynamics simulations with the finally optimized BChl database for LH2. Indeed, the energy profiles of the snapshots are found to be closely matching with the reference QC calculation data, with only ∼0.07 eV of errors in the ground- and excited-state energies and ∼0.008 eV of errors in the transition energies. We also show that properly selecting data points is actually quite important for generating an IM PES toward performing molecular dynamics simulations.

Correlation between single crystal data and molecular mechanics calculation of 3-(2-hydroxy-phenylimino)-1,3-diphenyl-propan-1-one

The Schiff base 3-(2-hydroxy-phenylimino)-1,3-diphenyl-propan-1-one (L), was prepared and characterized by elemental analysis, IR, 1H NMR and mass spectral techniques. The crystal and molecular structures were determined by using the single crystal X-ray diffraction data. The crystal was triclinic with unit-cell dimensions: a = 9.051(3) Å, b = 10.245(3) Å, c = 11.016(3) Å, α = 69.28(2)°, β = 81.66(2)°, γ = 64.25(2)°, space group is Pī and Z = 2. The structure was refined by least squares method. The single crystal was prepared by supersaturating method using tetramethoxysilane as gelling agent. Energy minimization for locating stable conformations and single point energy calculations for comparing the obtained and calculated conformations of the same molecule were carried out using molecular mechanics method (MM2 force field). The bond lengths and the bond angles were calculated and compared with the values obtained from the single crystal analysis. The obtained results were discussed and evaluated. Generally, the calculated parameters are in good agreement with the corresponding X-ray diffraction values. The obtained results reveal that, the amine-imine tautomerism is more preferable than keto-enol form and the intra-hydrogen bond stabilize the amine-imine tautomerism.

QM/MM Calculations for the Cl- + CH3Cl SN2 Reaction in Water Using CM5 Charges and Density Functional Theory.

The prototypical SN2 reaction of chloride ion with methyl chloride has been reinvestigated in aqueous solution using QM/MM methodology featuring MO6-2X/6-31+G(d) calculations with the TIP4P water model, and partial charges were computed with the CM5 method. Though the DFT method yields excellent gas-phase energetics for the reaction, the QM/MM approach is found to yield overestimation of the activation barrier by ca. 12 kcal/mol. The discrepancy is traced to underestimate of the magnitude of the partial charges on the chlorine atoms in the transition structure. When CM1 or CM3 charges based on semiempirical wave functions are used instead, the agreement with experiment is much improved. The findings emphasize the sensitivity of the results of QM/MM calculations to the choice of QM method, the MM force field, and implementation of the QM/MM interface.

Anharmonic Molecular Mechanics: Ab Initio Based Morse Parametrizations for the Popular MM3 Force Field

Methodologies for creating reactive potential energy surfaces from molecular mechanics force-fields are becoming increasingly popular. To date, molecular mechanics force-fields in biochemistry and small molecule organic chemistry tend to use harmonic expressions to treat bonding stretches, which is a poor approximation in reactive and nonequilibirum molecular dynamics simulations since bonds are often displaced significantly from their equilibrium positions. For such applications there is need for a better treatment of anharmonicity. In this contribution, Morse bonding potentials have been extensively parametrized for the atom types in the MM3 force field of Allinger and co-workers using high level CCSD(T)(F12*) energies. To our knowledge this is among the first instances of a comprehensive parametrization of Morse potentials in a popular organic chemistry force field. In the context of molecular dynamics simulations, these data will: (1) facilitate the fitting of reactive potential energy surfaces using empirical valence bond approaches and (2) enable more accurate treatments of energy transfer.

Interfacing the Core-Shell or the Drude Polarizable Force Field With Car-Parrinello Molecular Dynamics for QM/MM Simulations.

We report a quantum mechanics/polarizable–molecular mechanics (QM/p–MM) potential based molecular dynamics (MD) technique where the core–shell (or the Drude) type polarizable MM force field is interfaced with the plane-wave density functional theory based QM force field which allows Car–Parrinello MD for the QM subsystem. In the QM/p-MM Lagrangian proposed here, the shell (or the Drude) MM variables are treated as extended degrees of freedom along with the Kohn–Sham (KS) orbitals describing the QM wavefunction. The shell and the KS orbital degrees of freedom are then adiabatically decoupled from the nuclear degrees of freedom. In this respect, we also present here the Nosé–Hoover Chain thermostat implementation for the dynamical subsystems. Our approach is then used to investigate the effect of MM polarization on the QM/MM results. Especially, the consequence of MM polarization on reaction free energy barriers, defect formation energy, and structural and dynamical properties are investigated. A low point charge polarizable potential (p–MZHB) for pure siliceous systems is also reported here.

Electron momentum spectroscopy of cyclopentane: A molecular dynamical investigation

The experimental results of electron momentum spectroscopy studies on the valence orbitals of cyclopentane have been reinterpreted on the basis of Molecular Dynamical simulations employing the classical MM3 force field. Analysis using the produced molecular structures demonstrates that the pseudorotation and the variations of the ring puckering amplitude due to thermally induced vibrational effects have significant influence on the electron momentum distributions for the six highest-energy orbitals 1e2″, 3e1′ and 3e2′. The conformational average of the structures of C2 symmetry and the structures of Cs symmetry is also an equitable model to interpret the valence electron momentum distributions.

Combined Quantum Mechanics and Molecular Mechanics Studies of Enzymatic Reaction Mechanisms.

The combined quantum mechanics/molecular mechanics (QM/MM) methods have become a valuable tool in computational biochemistry and received versatile applications for studying the reaction mechanisms of enzymes. The approach combines the calculations of the electronic structure of the active site by QM, with modeling of the protein environment using MM force field, which allows the long-range electrostatics and steric effects on the enzyme reactivity to be accounted for. In this review, we review some key theoretical and computational aspects of the method and we also present some applications to particular enzymatic reactions such as tryptophan-7-halogenase, cyclooxygenase-1, and the epidermal growth factor receptor.

Conformational space of 4,4′-methoxypropylstilbene molecule

Quantum-chemical simulation of the ground state [the density function B3LYP/6-31G, B3LYP/6-31G(d), and B3LYP/6-31+G(d,p) and the perturbation theory MP2/6-31G(d) methods] and the transition states [the B3LYP/6-31G(d) method] of 4,4′-methoxypropylstilbene molecule has been performed. Using the Ellinger MM2 force field method, the potentials of internal rotation have been obtained for each rotational degree of freedom of the molecule. The B3LYP simulation has revealed the planarity of the conjugated system and the orthogonal position of the alkyl substituent, whereas the benzene rings have deviated by about 20° with respect to the double bond plane according to the MP2 data. Three transition states of the molecule corresponding to the saddle points of the first and the second orders have been revealed. The stationary points have been identified by means of vibrational analysis.

Producing DFT/MM enzyme reaction trajectories from SCC-DFTB/MM driving forces to probe the underlying electronics of a glycosyltransferase reaction.

The SCC-DFTB/MIO/CHARMM free energy surface for a glycosyltransferase, TcTS, is benchmarked against a DFT/MM reaction trajectory using the same CHARMM MM force field ported to the NWChem package. The popular B3LYP functional, against which the MIO parameter set was parameterized is used to optimize TS structures and run DFT reaction dynamics. A novel approach was used to generate reaction forces from a SCC-DFTB/MIO/CHARMM reaction surface to drive B3LYP/6-31G/MM and B3LYP/6-31G(d)/MM reaction trajectories. Although TS structures compare favorably, differences stemming primarily from a minimal basis set approximation prevented a successful 6-31G(d) FEARCF reaction dynamics trajectory. None the less, the dynamic evolution of the B3LYP/6-31G/MM-computed electron density provided an opportunity to perform NBO analysis along the reaction trajectory. Here, we illustrate that a successful ab initio reaction trajectory is computationally accessible when the underlying potential energy function of the semi-empirical method used to produce driving forces is sufficiently close to the ab initio potential. © 2017 Wiley Periodicals, Inc.

Systematic Quantum Mechanical Region Determination in QM/MM Simulation.

Hybrid quantum mechanical-molecular mechanical (QM/MM) simulations are widely used in enzyme simulation. Over ten convergence studies of QM/MM methods have revealed over the past several years that key energetic and structural properties approach asymptotic limits with only very large (ca. 500-1000 atom) QM regions. This slow convergence has been observed to be due in part to significant charge transfer between the core active site and the surrounding protein environment, which cannot be addressed by improvement of MM force fields or the embedding method employed within QM/MM. Given this slow convergence, it becomes essential to identify strategies for the most atom-economical determination of optimal QM regions and to gain insight into the crucial interactions captured only in large QM regions. Here, we extend and develop two methods for quantitative determination of QM regions. First, in the charge shift analysis (CSA) method, we probe the reorganization of electron density when core active site residues are removed completely, as determined by large-QM region QM/MM calculations. Second, we introduce the highly parallelizable Fukui shift analysis (FSA), which identifies how core/substrate frontier states are altered by the presence of an additional QM residue in smaller initial QM regions. We demonstrate that the FSA and CSA approaches are complementary and consistent on three test case enzymes: catechol O-methyltransferase, cytochrome P450cam, and hen eggwhite lysozyme. We also introduce validation strategies and test the sensitivities of the two methods to geometric structure, basis set size, and electronic structure methodology. Both methods represent promising approaches for the systematic, unbiased determination of quantum mechanical effects in enzymes and large systems that necessitate multiscale modeling.

Metal Complexes of a Novel Schiff Base Based on Penicillin: Characterization, Molecular Modeling, and Antibacterial Activity Study.

A novel Schiff base ligand of type HL was prepared by the condensation of amoxicillin trihydrate and nicotinaldehyde. The metal complexes of Co+2, Ni+2, Cu+2, and Zn+2 were characterized and investigated by physical and spectral techniques, namely, elemental analysis, melting point, conductivity, 1H NMR, IR, UV-Vis spectra, ESR, SEM, and mass spectrometry measurements. They were further analyzed by thermal technique (TGA/DTA) to gain better insight about the thermal stability and kinetic properties of the complexes. Thermal data revealed high thermal stability and nonspontaneous nature of the decomposition steps. The Coats-Redfern method was applied to extract thermodynamic parameters to explain the kinetic behavior. The molar conductance values were relatively low, showing their nonelectrolytic nature. The powder XRD pattern revealed amorphous nature except copper complex (1c) that crystallized in the triclinic crystal system. The EPR study strongly recommends the tetrahedral geometry of 1c. The structure optimization by MM force field calculation through ArgusLab 4.0.1 software program supports the concerned geometry of the complexes. The in vitro antibacterial activity of all the compounds, at their two different concentrations, was screened against four bacterial pathogens, namely, E. coli, P. vulgaris, K. pneumoniae, and S. aureus, and showed better activity compared to parent drug and control drug.

Hydrogen Bonding in Tight Environments: Simulated Force Spectroscopy of Nanoconfined Hydrogen-Bonded Complexes

The single-molecule force spectroscopy of a prototypical class of hydrogen-bonded complexes is computationally investigated. The complexes consist of derivatives of a barbituric acid and a Hamilton receptor that can form up to six simultaneous hydrogen bonds. The force–extension (F–L) isotherms of the host–guest complexes are simulated using classical molecular dynamics and the MM3 force field, for which a refined set of hydrogen bond parameters was developed from MP2 ab initio computations. The F–L curves exhibit peaks that signal conformational changes during elongation, the most prominent of which is in the 60–180 pN range and corresponds to the force required to break the hydrogen bonds. These peaks in the F–L curves are shown to be sensitive to relatively small changes in the chemical structure of the host molecule. Thermodynamic insights into the supramolecular assembly were obtained by reconstructing, from the force measurements, the Helmholtz free energy profile along the extension coordinate and ...

The effects of conformation and zwitterionic tautomerism on the structural and vibrational spectral data of anserine

In this study, the stable conformers of neutral anserine were searched by molecular dynamics simulations and energy minimization calculations using the MM2 force field. Thermochemical calculations at B3LYP/6-31G(d) level of theory followed these preliminary calculations. The results confirmed that neutral anserine has quite a flexible structure and many stable gauche and trans conformers at room temperature. Nevertheless, two are considerably more favourable in energy than the others and expected to dominate the gas-phase and matrix IR spectra of the molecule. The corresponding structural and vibrational spectral data for these two conformers of neutral anserine, whose relative stabilities were also examined by high-accuracy energy calculations carried out using G3MP2B3 method, and for the most stable conformer of anserine in zwitterion form were calculated at B3LYP/6-311++G(d,p) level of theory. The calculated harmonic force constants were refined using the Scaled Quantum Mechanical Force Field (SQM-FF) method and then used to produce the refined wavenumbers, potential energy distributions (PEDs) and IR and Raman intensities. These refined data together with the scaled harmonic wavenumbers obtained using another method, Dual Scale factors (DS), enabled us to correctly analyse the observed IR and Raman spectra of anserine and revealed the effects of conformation and zwitterionic tautomerism on its structural and vibrational spectral data.

Molecular Mechanics (MM4) Studies on Unusually Long Carbon-Carbon Bond Distances in Hydrocarbons.

The carbon-carbon single bond is of central importance in organic chemistry. When the molecular mechanics MM4 force field was developed beginning in the early 1990s, C-C bond lengths were not known very reliably for many important molecules, and bond lengths greater than 1.6 Å were quite poorly known experimentally. Quantum-mechanically computed values could not yet be obtained with useful accuracy in a general way. This paper examines structures now available from experiment and quantum-mechanical computations and extends the fit of the MM4 methodology to include new bond distances as long as 1.71 Å.

Comparing classical approaches with empirical or quantum-mechanically derived force fields for the simulation electronic lineshapes: application to coumarin dyes

The simulation of optical properties in complex and heterogeneous systems, as those involved in photovoltaic devices, requires efficient protocols able to account for the complexity of the environment. An usual procedure is based on molecular dynamics (MD) sampling and subsequent computation of the transition properties, at quantum mechanics/molecular mechanics (QM/MM) level, adopting, e.g., time-dependent density functional theory for the QM part. In this work, we adopt such protocol to simulate the spectra of two structurally related coumarin dyes employed as sensitizers in solar cells, c343 and nkx-2586, and we focus on gas-phase results with the scope to investigate in detail the impact of two key aspects of the protocol: the choice of the force field (FF) and the classical approximation of the spectra. We start from the selection of a FF able to accurately describe the intramolecular interactions in the system, comparing the predictions of a generic empirical force field (namely, GAFF) and a molecule-specific quantum-mechanically derived FF. The MD simulations generated by adopting each FF are analysed in terms of the sampled normal mode distributions, providing a rationale to understand the features (position and width) of the absorption spectra simulated from the distribution of vertical energies. By comparing those simulations with classical, semiclassical and fully quantum formulations of the vibronic lineshapes based on harmonic models, we demonstrate the necessity of high-quality FFs. Notwithstanding this, our results highlight that irrespective of the choice between the two FFs, the largest errors arise from the classical description of the fast nuclear motions (e.g. bond stretchings), thus emphasizing the need to account for quantum effects in order to achieve an accurate simulation of the spectrum. We further analyse the spectra simulated using a number of harmonic models of the initial and final states at either QM or FF levels, pointing out the occurrence of remarkable inaccuracies due to the large differences in the QM and FF normal modes. This problem is inherent to the fixed functional forms associated to any MM FF, and our results indicate a convenient route to avoid it. This knowledge is a prerequisite towards the goal of developing hybrid quantum/classical QM/MM approaches for the calculations of spectra in complex environments.

Simulated Solute Tempering in Fully Polarizable Hybrid QM/MM Molecular Dynamics Simulations.

We successfully apply a solute tempering approach, which substantially reduces the large number of temperature rungs required in conventional tempering methods by solvent charge scaling, to hybrid molecular dynamics simulations combining quantum mechanics with molecular mechanics (QM/MM). Specifically, we integrate a combination of density functional theory (DFT) and polarizable MM (PMM) force fields into the simulated solute tempering (SST) concept. We show that the required DFT/PMM-SST weight parameters can be obtained from inexpensive calculations and that for alanine dipeptide (DFT) in PMM water three rungs suffice to cover the temperature range from 300 to 550 K.

How Photoisomerization Drives Peptide Folding and Unfolding: Insights from QM/MM and MM Dynamics Simulations

AbstractPhotoswitchable azobenzene cross‐linkers can control the folding and unfolding of peptides by photoisomerization and can thus regulate peptide affinities and enzyme activities. Using quantum mechanics/molecular mechanics (QM/MM) methods and classical MM force fields, we report the first molecular dynamics simulations of the photoinduced folding and unfolding processes in the azobenzene cross‐linked FK‐11 peptide. We find that the interactions between the peptide and the azobenzene cross‐linker are crucial for controlling the evolution of the secondary structure of the peptide and responsible for accelerating the folding and unfolding events. They also modify the photoisomerization mechanism of the azobenzene cross‐linker compared with the situation in vacuo or in solution.

How Photoisomerization Drives Peptide Folding and Unfolding: Insights from QM/MM and MM Dynamics Simulations.

Photoswitchable azobenzene cross-linkers can control the folding and unfolding of peptides by photoisomerization and can thus regulate peptide affinities and enzyme activities. Using quantum mechanics/molecular mechanics (QM/MM) methods and classical MM force fields, we report the first molecular dynamics simulations of the photoinduced folding and unfolding processes in the azobenzene cross-linked FK-11 peptide. We find that the interactions between the peptide and the azobenzene cross-linker are crucial for controlling the evolution of the secondary structure of the peptide and responsible for accelerating the folding and unfolding events. They also modify the photoisomerization mechanism of the azobenzene cross-linker compared with the situation in vacuo or in solution.