Properties of isolated intercalators (ethidium (E), daunomycin (D), ellipticine (EL), and 4,6&... more Properties of isolated intercalators (ethidium (E), daunomycin (D), ellipticine (EL), and 4,6'-diaminide-2-phenylindole (DAPI)) and their stacking interactions with adenine...thymine (AT) and guanine...cytosine (GC) nucleic acid base pairs were investigated by means of a nonempirical correlated ab initio method. All intercalators exhibit large charge delocalization, and none of them (including the DAPI dication) exhibits a site with dominant charge. All intercalators have large polarizability and are good electron acceptors, while base pairs are good electron donors. MP2/6-31G*(0.25) stabilization energies of intercalator...base pair complexes are large (E...AT, 22.4 kcal/mol; D...GC, 17.8 kcal/mol; EL...GC, 18.2 kcal/mol; DAPI...GC, 21.1 kcal/mol) and are well reproduced by modified AMBER potential (van der Waals radii of intercalator atoms are enlarged and their energy depths are increased). Standard AMBER potential underestimates binding, especially for DAPI-containing complexes. Because the DAPI dication is the best electron acceptor (among all intercalators studied), this difference is explained by the importance of the charge-transfer term, which is not included in the AMBER potential. For the neutral EL molecule, the standard AMBER force field provides correct results. The Hartree-Fock and DFT/B3LYP methods, not covering the dispersion energy, fail completely to reveal any energy minimum at the potential energy curve of the E...AT complex, and these methods thus cannot be recommended for a study of intercalation process. On the other hand, an approximate version of the DFT method, which was extended to cover London dispersion energy, yields for all complexes very good stabilization energies that are well comparable with referenced ab initio data. Besides the vertical dependence of the interaction, an energy twist dependence of the interaction energy was also investigated by a reference correlated ab initio method and empirical potentials. It is concluded that, despite the cationic (E +1, D +1, DAPI +2) or polar (EL) character of the intercalators investigated, it is the dispersion energy which predominantly contributes to the stability of intercalator...base pair complexes. Any procedure which does not cover dispersion energy is thus not suitable for studying the process of intercalation.
Properties of isolated intercalators (ethidium (E), daunomycin (D), ellipticine (EL), and 4,6&... more Properties of isolated intercalators (ethidium (E), daunomycin (D), ellipticine (EL), and 4,6'-diaminide-2-phenylindole (DAPI)) and their stacking interactions with adenine...thymine (AT) and guanine...cytosine (GC) nucleic acid base pairs were investigated by means of a nonempirical correlated ab initio method. All intercalators exhibit large charge delocalization, and none of them (including the DAPI dication) exhibits a site with dominant charge. All intercalators have large polarizability and are good electron acceptors, while base pairs are good electron donors. MP2/6-31G*(0.25) stabilization energies of intercalator...base pair complexes are large (E...AT, 22.4 kcal/mol; D...GC, 17.8 kcal/mol; EL...GC, 18.2 kcal/mol; DAPI...GC, 21.1 kcal/mol) and are well reproduced by modified AMBER potential (van der Waals radii of intercalator atoms are enlarged and their energy depths are increased). Standard AMBER potential underestimates binding, especially for DAPI-containing complexes. Because the DAPI dication is the best electron acceptor (among all intercalators studied), this difference is explained by the importance of the charge-transfer term, which is not included in the AMBER potential. For the neutral EL molecule, the standard AMBER force field provides correct results. The Hartree-Fock and DFT/B3LYP methods, not covering the dispersion energy, fail completely to reveal any energy minimum at the potential energy curve of the E...AT complex, and these methods thus cannot be recommended for a study of intercalation process. On the other hand, an approximate version of the DFT method, which was extended to cover London dispersion energy, yields for all complexes very good stabilization energies that are well comparable with referenced ab initio data. Besides the vertical dependence of the interaction, an energy twist dependence of the interaction energy was also investigated by a reference correlated ab initio method and empirical potentials. It is concluded that, despite the cationic (E +1, D +1, DAPI +2) or polar (EL) character of the intercalators investigated, it is the dispersion energy which predominantly contributes to the stability of intercalator...base pair complexes. Any procedure which does not cover dispersion energy is thus not suitable for studying the process of intercalation.
Charge transfer in DNA has received much attention in the last few years due to its role in oxida... more Charge transfer in DNA has received much attention in the last few years due to its role in oxidative damage and repair in DNA and also due to possible applications of DNA in nanoelectronics. Despite intense experimental and theoretical efforts, the mechanism underlying long-range hole transport is still unresolved. This is in particular due to the sensitive dependence of charge transfer on the complex structure and dynamics of DNA and the interaction with the solvent, which could not be addressed adequately in the modeling approaches up to now. In this work, we study the factors governing hole transfer in detail, using a DFT-based fragment-orbital method, which allows to compute the charge transfer parameters along multinanosecond molecular dynamics simulations. Environmental effects are captured using a hybrid quantum mechanics-molecular mechanics (QM/MM) coupling scheme. This methodology allows to analyze several factors responsible for charge transfer in DNA in detail. The fluctuation of counterions, strongly counterbalanced by the surrounding water, leads to large oscillations of onsite energies, which govern the energetics of hole propagation along the DNA strand. In contrast, the electronic couplings depend only on DNA conformation and are not affected by the solvent. In particular, the onsite energies are strongly correlated between neighboring nucleobases, indicating that a conformational-gating type of mechanism may be induced by the collective environmental degrees of freedom.
In this work, we derive and test a new automatized strategy to construct repulsive potentials for... more In this work, we derive and test a new automatized strategy to construct repulsive potentials for the self-consistent charge density functional tight-binding (SCC-DFTB) method. This approach allows one to explore the parameter space in a systematic fashion in order to find optimal solutions. We find that due to the limited flexibility of the SCC-DFTB electronic part, not all properties can be optimized simultaneously. For example, the optimization of heats of formation is in conflict with the optimization of vibrational frequencies. Therefore, a special parametrization for vibrational frequencies is derived. It is shown that the performance of SCC-DFTB can be significantly improved using a more elaborate fitting strategy. A new fit for C and H is presented, which results in an average error of 2.6 kcal/mol for heats of formations for a large set of hydrocarbons, indicating that the performance of SCC-DFTB can be systematically improved also for other elements.
Description/Abstract Theory and Applications of Computational Chemistry: The First Forty Years is... more Description/Abstract Theory and Applications of Computational Chemistry: The First Forty Years is a collection of articles on the emergence of computational chemistry. It shows the breadth of theoretical and computational chemistry today and establishes how theory and computation have become increasingly linked as methodologies and technologies have advanced.
We report the parametrization of a density functional tight binding method (DFTB3) for copper in ... more We report the parametrization of a density functional tight binding method (DFTB3) for copper in a spin-polarized formulation. The parametrization is consistent with the framework of 3OB for main group elements (ONCHPS) and can be readily used for biological applications that involve copper proteins/peptides. The key to our parametrization is to introduce orbital angular momentum dependence of the Hubbard parameter and its charge derivative, thus allowing the 3d and 4s orbitals to adopt different sizes and responses to the change of charge state. The parametrization has been tested by applying to a fairly broad set of molecules of biological relevance, and the properties of interest include optimized geometries, ligand binding energies, and ligand proton affinities. Compared to the reference QM level (B3LYP/aug-cc-pVTZ, which is shown here to be similar to the B97-1 and CCSD(T) results, in terms of many properties of interest for a set of small copper containing molecules), our parametrization generally gives reliable structural properties for both Cu(I) and Cu(II) compounds, although several exceptions are also noted. For energetics, the results are more accurate for neutral ligands than for charged ligands, likely reflecting the minimal basis limitation of DFTB3; the results generally outperform NDDO based methods such as PM6 and even PBE with the 6-31+G(d,p) basis. For all ligand types, single-point B3LYP calculations at DFTB3 geometries give results very close (∼1-2 kcal/mol) to the reference B3LYP values, highlighting the consistency between DFTB3 and B3LYP structures. Possible further developments of the DFTB3 model for a better treatment of transition-metal ions are also discussed. In the current form, our first generation of DFTB3 copper model is expected to be particularly valuable as a method that drives sampling in systems that feature a dynamical copper binding site.
Journal of chemical theory and computation, Jan 13, 2015
Despite decades of investigations, the principal mechanisms responsible for the high affinity and... more Despite decades of investigations, the principal mechanisms responsible for the high affinity and specificity of proteins for key physiological cations K(+), Na(+), and Ca(2+) remain a hotly debated topic. At the core of the debate is an apparent need (or lack thereof) for an accurate description of the electrostatic response of the charge distribution in a protein to the binding of an ion. These effects range from partial electronic polarization of the directly ligating atoms to long-range effects related to partial charge transfer and electronic delocalization effects. While accurate modeling of cation recognition by metalloproteins warrants the use of quantum-mechanics (QM) calculations, the most popular approximations used in major biomolecular simulation packages rely on the implicit modeling of electronic polarization effects. That is, high-level QM computations for ion binding to proteins are desirable, but they are often unfeasible, because of the large size of the reactive-...
Semi-empirical quantum mechanical methods traditionally expand the electron density in a minimal,... more Semi-empirical quantum mechanical methods traditionally expand the electron density in a minimal, valence-only electron basis set. The minimal-basis approximation causes molecular polarization to be underestimated, and hence intermolecular interaction energies are also underestimated, especially for intermolecular interactions involving charged species. In this work, the third-order self-consistent charge density functional tight-binding method (DFTB3) is augmented with an auxiliary response density using the chemical-potential equalization (CPE) method and an empirical dispersion correction (D3). The parameters in the CPE and D3 models are fitted to high-level CCSD(T) reference interaction energies for a broad range of chemical species, as well as dipole moments calculated at the DFT level; the impact of including polarizabilities of molecules in the parameterization is also considered. Parameters for the elements H, C, N, O, and S are presented. The Root Mean Square Deviation (RMS...
A linear scaling QM/MM model for studying charge transport in high-mobility molecular semiconduct... more A linear scaling QM/MM model for studying charge transport in high-mobility molecular semiconductors is presented and applied to an anthracene single crystal and a hexabenzocoronene derivative in its liquid crystalline phase. The model includes both intra- and intermolecular electron-phonon couplings, long-range interactions with the environment, and corrections to the self-interaction error of density functional theory. By performing Ehrenfest simulations of the cationic system, hole mobilities are derived and compared to the experiment. A detailed picture of the charge carrier dynamics is given, and the performance of our method is discussed.
Photochromic ligands (PCLs), defined as photoswitchable molecules that are able to endow native r... more Photochromic ligands (PCLs), defined as photoswitchable molecules that are able to endow native receptors with a sensitivity towards light, have become a promising photopharmacological tool for various applications in biology. In general, PCLs consist of a ligand of the target receptor covalently linked to an azobenzene, which can be reversibly switched between two configurations upon light illumination. Gluazo, as a PCL that targets excitatory amino acid receptors, in its dark-adapted trans iso-form was characterized to be a partial agonist of the kainate glutamate receptor GluK2. Application of UV light leads to the formation of the cis form, with remarkedly reduced affinity towards GluK2. The mechanism of the change of ligand affinity induced by the photoisomerization was unresolved. The presented computational study explains how the isomerization of such a PCL affects the structural changes in the target receptor that lead to its activation.
We propose a new adaptive QM/MM method, the size-consistent multipartitioning (SCMP) QM/MM scheme... more We propose a new adaptive QM/MM method, the size-consistent multipartitioning (SCMP) QM/MM scheme, which enables stable and computationally efficient QM/MM simulations. A number of partitionings with identical size of the QM and the MM regions are considered with a new adaptive scheme in order to (1) realize smooth QM/MM switching, (2) introduce a conserved quantity (total energy, Hamiltonian), (3) avoid spurious artificial forces on the QM/MM border, and (4) allow for an efficient parallel implementation. Benchmark simulations performed for "QM water in MM water" show that energy conservation can be significantly improved and the computational efficiency allows treating also larger QM regions, for which previous methods had to face an intractable increase in computer time.
Theory and Applications of Computational Chemistry, 2005
Description/Abstract Theory and Applications of Computational Chemistry: The First Forty Years is... more Description/Abstract Theory and Applications of Computational Chemistry: The First Forty Years is a collection of articles on the emergence of computational chemistry. It shows the breadth of theoretical and computational chemistry today and establishes how theory and computation have become increasingly linked as methodologies and technologies have advanced.
ABSTRACT We present an extension to the recent 3OB parametrization of the Density Functional Tigh... more ABSTRACT We present an extension to the recent 3OB parametrization of the Density Functional Tight Binding Model DFTB3 for biological and organic systems. Parameters for the halogens F, Cl, Br, and I have been developed for use in covalently bound systems and benchmarked on a test set of 106 molecules (the ‘OrgX’ set), using bonding distances, bonding angles, atomization energies, and vibrational frequencies to assess the performance of the parameters. Additional testing has been done with the X40 set of 40 supramolecular systems containing halogens, adding a simple correction for the halogen bonds that are strongly overbound in DFTB3. Furthermore, parameters for Ca, K, and Na as counterions in biological systems have been created. To benchmark geometries as well as ligand binding energies a test set ‘BioMe’ of 210 molecules has been created that cover coordination to various functional groups frequently occurring in biological systems. The new DFTB3/3OB parameter set outperforms DFT calculations with a double-ζ basis set in terms of energies and can reproduce DFT geometries, with some minor deviations in bond distances and angles due to the use of a minimal basis set. http://dx.doi.org/10.1021/ct5009137 (Manuscript not provided due to American Chemical Society publishing policy)
In the last years, we have developed a computationally efficient approximation to density functio... more In the last years, we have developed a computationally efficient approximation to density functional theory, the so called self-consistent charge density functional tight-binding scheme (SCC-DFTB). To extend its applicability to biomolecular structures, this method has been implemented into quantum mechanical/molecular mechanics (QM/MM) and linear scaling schemes and augmented with an empirical treatment of the dispersion forces. We review here applications
Here we shall review our studies of peptides, both in zwitterionic and neutral forms, in the gas ... more Here we shall review our studies of peptides, both in zwitterionic and neutral forms, in the gas phase and aqueous solution. Our benchmarking of the ab initio restricted Hartree–Fock and Møller Plesset theory, density functional theory (DFT) (Becke3LYP), semi-empirical and approximate DFT methods for the description of structure, energetics and vibrational properties are also discussed. In addition, in our modeling
Properties of isolated intercalators (ethidium (E), daunomycin (D), ellipticine (EL), and 4,6&... more Properties of isolated intercalators (ethidium (E), daunomycin (D), ellipticine (EL), and 4,6'-diaminide-2-phenylindole (DAPI)) and their stacking interactions with adenine...thymine (AT) and guanine...cytosine (GC) nucleic acid base pairs were investigated by means of a nonempirical correlated ab initio method. All intercalators exhibit large charge delocalization, and none of them (including the DAPI dication) exhibits a site with dominant charge. All intercalators have large polarizability and are good electron acceptors, while base pairs are good electron donors. MP2/6-31G*(0.25) stabilization energies of intercalator...base pair complexes are large (E...AT, 22.4 kcal/mol; D...GC, 17.8 kcal/mol; EL...GC, 18.2 kcal/mol; DAPI...GC, 21.1 kcal/mol) and are well reproduced by modified AMBER potential (van der Waals radii of intercalator atoms are enlarged and their energy depths are increased). Standard AMBER potential underestimates binding, especially for DAPI-containing complexes. Because the DAPI dication is the best electron acceptor (among all intercalators studied), this difference is explained by the importance of the charge-transfer term, which is not included in the AMBER potential. For the neutral EL molecule, the standard AMBER force field provides correct results. The Hartree-Fock and DFT/B3LYP methods, not covering the dispersion energy, fail completely to reveal any energy minimum at the potential energy curve of the E...AT complex, and these methods thus cannot be recommended for a study of intercalation process. On the other hand, an approximate version of the DFT method, which was extended to cover London dispersion energy, yields for all complexes very good stabilization energies that are well comparable with referenced ab initio data. Besides the vertical dependence of the interaction, an energy twist dependence of the interaction energy was also investigated by a reference correlated ab initio method and empirical potentials. It is concluded that, despite the cationic (E +1, D +1, DAPI +2) or polar (EL) character of the intercalators investigated, it is the dispersion energy which predominantly contributes to the stability of intercalator...base pair complexes. Any procedure which does not cover dispersion energy is thus not suitable for studying the process of intercalation.
Properties of isolated intercalators (ethidium (E), daunomycin (D), ellipticine (EL), and 4,6&... more Properties of isolated intercalators (ethidium (E), daunomycin (D), ellipticine (EL), and 4,6'-diaminide-2-phenylindole (DAPI)) and their stacking interactions with adenine...thymine (AT) and guanine...cytosine (GC) nucleic acid base pairs were investigated by means of a nonempirical correlated ab initio method. All intercalators exhibit large charge delocalization, and none of them (including the DAPI dication) exhibits a site with dominant charge. All intercalators have large polarizability and are good electron acceptors, while base pairs are good electron donors. MP2/6-31G*(0.25) stabilization energies of intercalator...base pair complexes are large (E...AT, 22.4 kcal/mol; D...GC, 17.8 kcal/mol; EL...GC, 18.2 kcal/mol; DAPI...GC, 21.1 kcal/mol) and are well reproduced by modified AMBER potential (van der Waals radii of intercalator atoms are enlarged and their energy depths are increased). Standard AMBER potential underestimates binding, especially for DAPI-containing complexes. Because the DAPI dication is the best electron acceptor (among all intercalators studied), this difference is explained by the importance of the charge-transfer term, which is not included in the AMBER potential. For the neutral EL molecule, the standard AMBER force field provides correct results. The Hartree-Fock and DFT/B3LYP methods, not covering the dispersion energy, fail completely to reveal any energy minimum at the potential energy curve of the E...AT complex, and these methods thus cannot be recommended for a study of intercalation process. On the other hand, an approximate version of the DFT method, which was extended to cover London dispersion energy, yields for all complexes very good stabilization energies that are well comparable with referenced ab initio data. Besides the vertical dependence of the interaction, an energy twist dependence of the interaction energy was also investigated by a reference correlated ab initio method and empirical potentials. It is concluded that, despite the cationic (E +1, D +1, DAPI +2) or polar (EL) character of the intercalators investigated, it is the dispersion energy which predominantly contributes to the stability of intercalator...base pair complexes. Any procedure which does not cover dispersion energy is thus not suitable for studying the process of intercalation.
Charge transfer in DNA has received much attention in the last few years due to its role in oxida... more Charge transfer in DNA has received much attention in the last few years due to its role in oxidative damage and repair in DNA and also due to possible applications of DNA in nanoelectronics. Despite intense experimental and theoretical efforts, the mechanism underlying long-range hole transport is still unresolved. This is in particular due to the sensitive dependence of charge transfer on the complex structure and dynamics of DNA and the interaction with the solvent, which could not be addressed adequately in the modeling approaches up to now. In this work, we study the factors governing hole transfer in detail, using a DFT-based fragment-orbital method, which allows to compute the charge transfer parameters along multinanosecond molecular dynamics simulations. Environmental effects are captured using a hybrid quantum mechanics-molecular mechanics (QM/MM) coupling scheme. This methodology allows to analyze several factors responsible for charge transfer in DNA in detail. The fluctuation of counterions, strongly counterbalanced by the surrounding water, leads to large oscillations of onsite energies, which govern the energetics of hole propagation along the DNA strand. In contrast, the electronic couplings depend only on DNA conformation and are not affected by the solvent. In particular, the onsite energies are strongly correlated between neighboring nucleobases, indicating that a conformational-gating type of mechanism may be induced by the collective environmental degrees of freedom.
In this work, we derive and test a new automatized strategy to construct repulsive potentials for... more In this work, we derive and test a new automatized strategy to construct repulsive potentials for the self-consistent charge density functional tight-binding (SCC-DFTB) method. This approach allows one to explore the parameter space in a systematic fashion in order to find optimal solutions. We find that due to the limited flexibility of the SCC-DFTB electronic part, not all properties can be optimized simultaneously. For example, the optimization of heats of formation is in conflict with the optimization of vibrational frequencies. Therefore, a special parametrization for vibrational frequencies is derived. It is shown that the performance of SCC-DFTB can be significantly improved using a more elaborate fitting strategy. A new fit for C and H is presented, which results in an average error of 2.6 kcal/mol for heats of formations for a large set of hydrocarbons, indicating that the performance of SCC-DFTB can be systematically improved also for other elements.
Description/Abstract Theory and Applications of Computational Chemistry: The First Forty Years is... more Description/Abstract Theory and Applications of Computational Chemistry: The First Forty Years is a collection of articles on the emergence of computational chemistry. It shows the breadth of theoretical and computational chemistry today and establishes how theory and computation have become increasingly linked as methodologies and technologies have advanced.
We report the parametrization of a density functional tight binding method (DFTB3) for copper in ... more We report the parametrization of a density functional tight binding method (DFTB3) for copper in a spin-polarized formulation. The parametrization is consistent with the framework of 3OB for main group elements (ONCHPS) and can be readily used for biological applications that involve copper proteins/peptides. The key to our parametrization is to introduce orbital angular momentum dependence of the Hubbard parameter and its charge derivative, thus allowing the 3d and 4s orbitals to adopt different sizes and responses to the change of charge state. The parametrization has been tested by applying to a fairly broad set of molecules of biological relevance, and the properties of interest include optimized geometries, ligand binding energies, and ligand proton affinities. Compared to the reference QM level (B3LYP/aug-cc-pVTZ, which is shown here to be similar to the B97-1 and CCSD(T) results, in terms of many properties of interest for a set of small copper containing molecules), our parametrization generally gives reliable structural properties for both Cu(I) and Cu(II) compounds, although several exceptions are also noted. For energetics, the results are more accurate for neutral ligands than for charged ligands, likely reflecting the minimal basis limitation of DFTB3; the results generally outperform NDDO based methods such as PM6 and even PBE with the 6-31+G(d,p) basis. For all ligand types, single-point B3LYP calculations at DFTB3 geometries give results very close (∼1-2 kcal/mol) to the reference B3LYP values, highlighting the consistency between DFTB3 and B3LYP structures. Possible further developments of the DFTB3 model for a better treatment of transition-metal ions are also discussed. In the current form, our first generation of DFTB3 copper model is expected to be particularly valuable as a method that drives sampling in systems that feature a dynamical copper binding site.
Journal of chemical theory and computation, Jan 13, 2015
Despite decades of investigations, the principal mechanisms responsible for the high affinity and... more Despite decades of investigations, the principal mechanisms responsible for the high affinity and specificity of proteins for key physiological cations K(+), Na(+), and Ca(2+) remain a hotly debated topic. At the core of the debate is an apparent need (or lack thereof) for an accurate description of the electrostatic response of the charge distribution in a protein to the binding of an ion. These effects range from partial electronic polarization of the directly ligating atoms to long-range effects related to partial charge transfer and electronic delocalization effects. While accurate modeling of cation recognition by metalloproteins warrants the use of quantum-mechanics (QM) calculations, the most popular approximations used in major biomolecular simulation packages rely on the implicit modeling of electronic polarization effects. That is, high-level QM computations for ion binding to proteins are desirable, but they are often unfeasible, because of the large size of the reactive-...
Semi-empirical quantum mechanical methods traditionally expand the electron density in a minimal,... more Semi-empirical quantum mechanical methods traditionally expand the electron density in a minimal, valence-only electron basis set. The minimal-basis approximation causes molecular polarization to be underestimated, and hence intermolecular interaction energies are also underestimated, especially for intermolecular interactions involving charged species. In this work, the third-order self-consistent charge density functional tight-binding method (DFTB3) is augmented with an auxiliary response density using the chemical-potential equalization (CPE) method and an empirical dispersion correction (D3). The parameters in the CPE and D3 models are fitted to high-level CCSD(T) reference interaction energies for a broad range of chemical species, as well as dipole moments calculated at the DFT level; the impact of including polarizabilities of molecules in the parameterization is also considered. Parameters for the elements H, C, N, O, and S are presented. The Root Mean Square Deviation (RMS...
A linear scaling QM/MM model for studying charge transport in high-mobility molecular semiconduct... more A linear scaling QM/MM model for studying charge transport in high-mobility molecular semiconductors is presented and applied to an anthracene single crystal and a hexabenzocoronene derivative in its liquid crystalline phase. The model includes both intra- and intermolecular electron-phonon couplings, long-range interactions with the environment, and corrections to the self-interaction error of density functional theory. By performing Ehrenfest simulations of the cationic system, hole mobilities are derived and compared to the experiment. A detailed picture of the charge carrier dynamics is given, and the performance of our method is discussed.
Photochromic ligands (PCLs), defined as photoswitchable molecules that are able to endow native r... more Photochromic ligands (PCLs), defined as photoswitchable molecules that are able to endow native receptors with a sensitivity towards light, have become a promising photopharmacological tool for various applications in biology. In general, PCLs consist of a ligand of the target receptor covalently linked to an azobenzene, which can be reversibly switched between two configurations upon light illumination. Gluazo, as a PCL that targets excitatory amino acid receptors, in its dark-adapted trans iso-form was characterized to be a partial agonist of the kainate glutamate receptor GluK2. Application of UV light leads to the formation of the cis form, with remarkedly reduced affinity towards GluK2. The mechanism of the change of ligand affinity induced by the photoisomerization was unresolved. The presented computational study explains how the isomerization of such a PCL affects the structural changes in the target receptor that lead to its activation.
We propose a new adaptive QM/MM method, the size-consistent multipartitioning (SCMP) QM/MM scheme... more We propose a new adaptive QM/MM method, the size-consistent multipartitioning (SCMP) QM/MM scheme, which enables stable and computationally efficient QM/MM simulations. A number of partitionings with identical size of the QM and the MM regions are considered with a new adaptive scheme in order to (1) realize smooth QM/MM switching, (2) introduce a conserved quantity (total energy, Hamiltonian), (3) avoid spurious artificial forces on the QM/MM border, and (4) allow for an efficient parallel implementation. Benchmark simulations performed for "QM water in MM water" show that energy conservation can be significantly improved and the computational efficiency allows treating also larger QM regions, for which previous methods had to face an intractable increase in computer time.
Theory and Applications of Computational Chemistry, 2005
Description/Abstract Theory and Applications of Computational Chemistry: The First Forty Years is... more Description/Abstract Theory and Applications of Computational Chemistry: The First Forty Years is a collection of articles on the emergence of computational chemistry. It shows the breadth of theoretical and computational chemistry today and establishes how theory and computation have become increasingly linked as methodologies and technologies have advanced.
ABSTRACT We present an extension to the recent 3OB parametrization of the Density Functional Tigh... more ABSTRACT We present an extension to the recent 3OB parametrization of the Density Functional Tight Binding Model DFTB3 for biological and organic systems. Parameters for the halogens F, Cl, Br, and I have been developed for use in covalently bound systems and benchmarked on a test set of 106 molecules (the ‘OrgX’ set), using bonding distances, bonding angles, atomization energies, and vibrational frequencies to assess the performance of the parameters. Additional testing has been done with the X40 set of 40 supramolecular systems containing halogens, adding a simple correction for the halogen bonds that are strongly overbound in DFTB3. Furthermore, parameters for Ca, K, and Na as counterions in biological systems have been created. To benchmark geometries as well as ligand binding energies a test set ‘BioMe’ of 210 molecules has been created that cover coordination to various functional groups frequently occurring in biological systems. The new DFTB3/3OB parameter set outperforms DFT calculations with a double-ζ basis set in terms of energies and can reproduce DFT geometries, with some minor deviations in bond distances and angles due to the use of a minimal basis set. http://dx.doi.org/10.1021/ct5009137 (Manuscript not provided due to American Chemical Society publishing policy)
In the last years, we have developed a computationally efficient approximation to density functio... more In the last years, we have developed a computationally efficient approximation to density functional theory, the so called self-consistent charge density functional tight-binding scheme (SCC-DFTB). To extend its applicability to biomolecular structures, this method has been implemented into quantum mechanical/molecular mechanics (QM/MM) and linear scaling schemes and augmented with an empirical treatment of the dispersion forces. We review here applications
Here we shall review our studies of peptides, both in zwitterionic and neutral forms, in the gas ... more Here we shall review our studies of peptides, both in zwitterionic and neutral forms, in the gas phase and aqueous solution. Our benchmarking of the ab initio restricted Hartree–Fock and Møller Plesset theory, density functional theory (DFT) (Becke3LYP), semi-empirical and approximate DFT methods for the description of structure, energetics and vibrational properties are also discussed. In addition, in our modeling
Uploads
Papers by Marcus Elstner