Characterization of cisplatin/membrane interactions by QM/MM energy decomposition analysis

QM/MM-Based Energy Decomposition Analysis Method for Large Systems

In this work, a QM/MM-based EDA method, called GKS-EDA(QM/MM), is proposed. As an extension of GKS-EDA, this scheme divides the total interaction energy into electrostatic, exchange-repulsion, polarization, and correlation/dispersion terms. GKS-EDA(QM/MM) can be applied to describe the interactions of large-scale systems combined with various QM/MM platforms. By using the examples of a hydrated hydronium ion complex in water solution, the barnase-barstar complex, and MMP-13-pyrimidinetrione in a metalloprotein, the capability of GKS-EDA(QM/MM) for various interactions in large systems is validated.

Translocation Processes of Pt(II)-Based Drugs through Human Breast Cancer Cell Membrane: In Silico Experiments

Breast cancer is one of the most frequent modalities of cancer worldwide, with notable mortality. The medication based on platinum drugs (cisplatin (cddp), carboplatin (cpx), and oxaliplatin (oxa)) is a conventional chemotherapy despite severe side effects and the development of drug resistance. In order to provide a deeper molecular description of the influx and efflux processes of platinum drugs through breast cancer tissues, this study focuses on molecular dynamics (MD) simulations of the passive translocation process through a realistic plasma membrane prototype of human breast cancer cell (c_memb). The results showed that the permeation events were mainly mediated by neutral lipids (DOPC, DOPE, and cholesterol), producing a low and temporary membrane deformation. The drug insertion in the region of polar heads was the most favorable stage of the translocation mechanism, especially for cddp and oxa with potential wells of -8.6 and -9.8 kcal mol-1, respectively. However, the potentials of mean force (PMF) revealed unfavorable kinetics for the permeation of these drugs through lipid tails, with energy barriers of 28.3 (cddp), 32.2 (cpx), and 30.4 kcal mol-1 (oxa). The low permeability coefficients (P) of cpx and oxa, which were 3 and 1 orders of magnitude inferior than for cddp, resulted from the high energy barriers for their traslocation processes through the membrane. The obtained results provide a more accurate picture of the permeation of Pt(II)-based drugs through breast cancer cells, which may be relevant for the design and evaluation of new platinum complexes.

MoBioTools: A toolkit to setup quantum mechanics/molecular mechanics calculations

We present a toolkit that allows for the preparation of QM/MM input files from a conformational ensemble of molecular geometries. The package is currently compatible with trajectory and topology files in Amber, CHARMM, GROMACS and NAMD formats, and has the possibility to generate QM/MM input files for Gaussian (09 and 16), Orca (≥4.0), NWChem and (Open)Molcas. The toolkit can be used in command line, so that no programming experience is required, although it presents some features that can also be employed as a python application programming interface. We apply the toolkit in four situations in which different electronic-structure properties of organic molecules in the presence of a solvent or a complex biological environment are computed: the reduction potential of the nucleobases in acetonitrile, an energy decomposition analysis of tyrosine interacting with water, the absorption spectrum of an azobenzene derivative integrated into a voltage-gated ion channel, and the absorption and emission spectra of the luciferine/luciferase complex. These examples show that the toolkit can be employed in a manifold of situations for both the electronic ground state and electronically excited states. It also allows for the automatic correction of the active space in the case of CASSCF calculations on an ensemble of geometries, as it is shown for the azobenzene derivative photoswitch case.

Effect of the QM Size, Basis Set, and Polarization on QM/MM Interaction Energy Decomposition Analysis

Herein, an Energy Decomposition Analysis (EDA) scheme extended to the framework of QM/MM calculations in the context of electrostatic embeddings (QM/MM-EDA) including atomic charges and dipoles is applied to assess the effect of the QM region size on the convergence of the different interaction energy components, namely, electrostatic, Pauli, and polarization, for cationic, anionic, and neutral systems interacting with a strong polar environment (water). Significant improvements are found when the bulk solvent environment is described by a MM potential in the EDA scheme as compared to pure QM calculations that neglect bulk solvation. The predominant electrostatic interaction requires sizable QM regions. The results reported here show that it is necessary to include a surprisingly large number of water molecules in the QM region to obtain converged values for this energy term, contrary to most cluster models often employed in the literature. Both the improvement of the QM wave function by means of a larger basis set and the introduction of polarization into the MM region through a polarizable force field do not translate to a faster convergence with the QM region size, but they lead to better results for the different interaction energy components. The results obtained in this work provide insight into the effect of each energy component on the convergence of the solute-solvent interaction energy with the QM region size. This information can be used to improve the MM FFs and embedding schemes employed in QM/MM calculations of solvated systems.

Benefits of hybrid QM/MM over traditional classical mechanics in pharmaceutical systems

Hybrid quantum mechanics/molecular mechanics (QM/MM) is one of the most reliable approaches for accurately modeling and studying the complex pharmaceutical discovery system. Classical mechanics has significantly accelerated the drug discovery process in the past decade. However, the current challenge is the large pool of false positives, which require extensive validation. Hybrid QM/MM is an effective solution for accurately studying ligand binding, structural mechanisms, free energy evaluation, and spectroscopic characterization. This article highlights the methodological details relevant to cost-effective hybrid QM/MM methods. This approach, combined with traditional pharmacoinformatics methods, could be a reliable strategy to balance the cost and accuracy of the calculations.

On the Permeation of Polychlorinated Dibenzodioxins and Dibenzofurans through Lipid Membranes: Classical MD and Hybrid QM/MM-EDA Analysis

The permeation of dioxin-like pollutants, namely, chlorinated dibenzodioxins and dibenzofurans, through lipid membranes has been simulated using classic molecular dynamics (CMD) combined with the umbrella sampling approach. The most toxic forms of chlorinated dibenzodioxin and dibenzofuran, 2,3,7,8-tetrachloro-p-dibenzodioxin (TCDD) and 2,3,7,8-tetrachlorodibenzofuran (TCDF), and a dioleyl-phosphatidylcholine (DOPC) lipid membrane of 50 Å wide have been chosen for our study. The free energy profile shows the penetration process is largely favoured thermodynamically (ΔG ≈ -12 kcal/mol), with a progressively decrease of the free energy until reaching the energy minima at distances of 8 Å and 9.5 Å from the centre of the membrane for, respectively, TCDD and TCDF. At the centre of the membrane, both molecules display subtle local maxima with free energy differences of 0.5 and 1 kcal/mol with respect to the energy minima for TCDD and TCDF, respectively. Furthermore, the intermolecular interactions between the molecules and the lipid membrane have been characterized at the minima and the local maxima using hybrid quantum mechanics/molecular mechanics energy decomposition analysis (QM/MM-EDA). Total interaction energies of -17.5 and -16.5 kcal/mol have been found at the energy minima for TCDD and TCDF, respectively. In both cases, the dispersion forces govern the molecule-membrane interactions, no significant changes have been found at the local maxima, in agreement with the classical free energy profile. The small differences found in the results obtained for TCDD and TCDF point out that the adsorption and diffusion processes through the cell membrane are not related to the different toxicity shown by these pollutants.