Cisplatin Binding to DNA Oligomers from Hybrid Car-Parrinello/Molecular Dynamics Simulations

What is the Reason That the Pharmacological Future of Chemotherapeutics in the Treatment of Lung Cancer Could Be Most Closely Related to Nanostructures? Platinum Drugs in Therapy of Non-Small and Small Cell Lung Cancer and Their Unexpected, Possible Interactions. The Review

Over the course of several decades, anticancer treatment with chemotherapy drugs for lung cancer has not changed significantly. Unfortunately, this treatment prolongs the patient's life only by a few months, causing many side effects in the human body. It has also been proven that drugs such as Cisplatin, Carboplatin, Oxaliplatin and others can react with other substances containing an aromatic ring in which the nitrogen atom has a free electron group in its structure. Thus, such structures may have a competitive effect on the nucleobases of DNA. Therefore, scientists are looking not only for new drugs, but also for new alternative ways of delivering the drug to the cancer site. Nanotechnology seems to be a great hope in this matter. Creating a new nanomedicine would reduce the dose of the drug to an absolute minimum, and thus limit the toxic effect of the drug; it would allow for the exclusion of interactions with competitive compounds with a structure similar to nucleobases; it would also permit using the so-called targeted treatment and bypassing healthy cells; it would allow for the introduction of other treatment options, such as radiotherapy directly to the cancer site; and it would provide diagnostic possibilities. This article is a review that aims to systematize the knowledge regarding the anticancer treatment of lung cancer, but not only. It shows the clear possibility of interactions of chemotherapeutics with compounds competitive to the nitrogenous bases of DNA. It also shows the possibilities of using nanostructures as potential Platinum drug carriers, and proves that nanomedicine can easily become a new medicinal product in personalized medicine.

Predicting 195 Pt NMR Chemical Shifts in Water-Soluble Inorganic/Organometallic Complexes with a Fast and Simple Protocol Combining Semiempirical Modeling and Machine Learning

Water-soluble Pt complexes are the key components in medicinal chemistry and catalysis. The well-known cisplatin family of anticancer drugs and industrial hydrosylilation catalysts are two leading examples. On the molecular level, the activity mechanisms of such complexes mostly involve changes in the Pt coordination sphere. Using ¹⁹⁵ Pt NMR spectroscopy for operando monitoring would be a valuable tool for uncovering the activity mechanisms; however, reliable approaches for the rapid correlation of Pt complex structure with ¹⁹⁵ Pt chemical shifts are very challenging and not available for everyday research practice. While NMR shielding is a response property, molecular 3D structure determines NMR spectra, as widely known, which allows us to build up 3D structure to ¹⁹⁵ Pt chemical shift correlations. Accordingly, we present a new workflow for the determination of lowest-energy configurational/conformational isomers based on the GFN2-xTB semiempirical method and prediction of corresponding chemical shifts with a Machine Learning (ML) model tuned for Pt complexes. The workflow was designed for the prediction of ¹⁹⁵ Pt chemical shifts of water-soluble Pt(II) and Pt(IV) anionic, neutral, and cationic complexes with halide, NO₂ ^- , (di)amino, and (di)carboxylate ligands with chemical shift values ranging from -6293 to 7090 ppm. The model offered an accuracy (normalized root-mean-square deviation/RMSD) of 1.08 %/145.02 ppm on the held-out test set.

Will the Interactions of Some Platinum (II)-Based Drugs with B-Vitamins Reduce Their Therapeutic Effect in Cancer Patients? Comparison of Chemotherapeutic Agents such as Cisplatin, Carboplatin and Oxaliplatin-A Review

Pt (II) derivatives show anti-cancer activity by interacting with nucleobases of DNA, thus causing some spontaneous and non-spontaneous reactions. As a result, mono- and diaqua products are formed which further undergo complexation with guanine or adenine. Consequently, many processes are triggered, which lead to the death of the cancer cell. The theoretical and experimental studies confirm that such types of interactions can also occur with other chemical compounds. The vitamins from B group have a similar structure to the nucleobases of DNA and have aromatic rings with single-pair orbitals. Theoretical and experimental studies were performed to describe the interactions of B vitamins with Pt (II) derivatives such as cisplatin, oxaliplatin and carboplatin. The obtained results were compared with the values for guanine. Two levels of simulations were implemented at the theoretical level, namely, B3LYP/6-31G(d,p) with LANL2DZ bases set for platinum atoms and MN15/def2-TZVP. The polarizable continuum model (IEF-PCM preparation) and water as a solvent were used. UV-Vis spectroscopy was used to describe the drug-nucleobase and drug-B vitamin interactions. Values of the free energy (ΔG_r) show spontaneous reactions with mono- and diaqua derivatives of cisplatin and oxaliplatin; however, interactions with diaqua derivatives are more preferable. The strength of these interactions was also compared. Carboplatin products have the weakest interaction with the studied structures. The presence of non-covalent interactions was demonstrated in the tested complexes. A good agreement between theory and experiment was also demonstrated.

How can the cisplatin analogs with different amine act on DNA during cancer treatment theoretically?

Cisplatin is a widely used anti-cancer drug which inhibits the replication and polymerization of DNA molecule while showing some side effects and drug resistance. For this reason, to enhance its therapeutic index, researchers have synthesized several thousand analogs and tested their properties. In this project, several cisplatin analogs were designed to theoretically study the biological activity and lipophilicity effects on amine changes. The amines of the cisplatin molecule were substituted with aliphatic amines in different analogs. Computational methods such as molecular dynamics simulation, molecular docking, and molecular mechanics Poisson-Boltzmann surface area analysis were performed to investigate the binding of six cisplatin derivatives with DNA. The binding affinity and potential interactions of these drugs with double-strand DNA were analyzed. The stability effect of these drugs was investigated via root-mean-square deviation and root-mean-square fluctuation analysis, which showed that some analogs can break base-pair interaction at the end of DNA and reduced the stability of DNA. Also, the results revealed that the hydrogen bond is one of the most important factors in the binding of cisplatin's adduct to DNA. Molecular mechanics Poisson-Boltzmann surface area analysis indicated that electrostatic and van der Waals interactions are the most important deriving forces to the binding of cisplatin's drug to DNA. Finally, data revealed that cisplatin and the cis-dichloro-dimethylamine-platin tendency for binding to DNA are greater than that of other analogs.

Frontiers of metal-coordinating drug design

Introduction: The occurrence of metal ions in biomolecules is required to exert vital cellular functions. Metal-containing biomolecules can be modulated by small-molecule inhibitors targeting their metal-moiety. As well, the discovery of cisplatin ushered the rational discovery of metal-containing-drugs. The use of both drug types exploiting metal-ligand interactions is well established to treat distinct pathologies. Therefore, characterizing and leveraging metal-coordinating drugs is a pivotal, yet challenging, part of medicinal chemistry.Area covered: Atomic-level simulations are increasingly employed to overcome the challenges met by traditional drug-discovery approaches and to complement wet-lab experiments in elucidating the mechanisms of drugs' action. Multiscale simulations, allow deciphering the mechanism of metal-binding inhibitors and metallo-containing-drugs, enabling a reliable description of metal-complexes in their biological environment. In this compendium, the authors review selected applications exploiting the metal-ligand interactions by focusing on understanding the mechanism and design of (i) inhibitors targeting iron and zinc-enzymes, and (ii) ruthenium and gold-based anticancer agents targeting the nucleosome and aquaporin protein, respectively.Expert opinion: The showcased applications exemplify the current role and the potential of atomic-level simulations and reveal how their synergic use with experiments can contribute to uncover fundamental mechanistic facets and exploit metal-ligand interactions in medicinal chemistry.

Does the Affinity of Cisplatin to B-Vitamins Impair the Therapeutic Effect in the Case of Patients with Lung Cancer-consuming Carrot or Beet Juice?

BACKGROUND

Cisplatin (CisPt) has a well-recognized anticancer activity by interacting with DNA and inducing programmed cell death. However, theoretical studies performed on the molecular level suggest that such nonspecific interactions can also take place with many competitive compounds, such as vitamins containing aromatic rings with lone-pair orbitals.

OBJECTIVE

This work is a theoretical study on the initial Pt-N₇(N₁) bond formation with vitamins from B group and their comparison with values characterizing native purines.

METHODS

Geometries of studied structures were optimized with an aid of Gaussian 09 using the B3LYP functional with the 6-31G** basis set. Atomic orbitals of platinum were represented by the lanl2dz basis. Solvation free energies were evaluated by a self-consistent reaction field (SCRF) approach. A dielectric constant of 78 for water was used in the PCM continuum model computations along with radii Bondii.

RESULTS

The affinities of mono-aqua cis~[Pt(NH₃)₂Cl(H₂O)]+ and di-aqua cis~ [Pt(NH₃)₂(H₂O)₂]²⁺ derivatives of Cisplatinum toward compounds belonging to the group of eight B vitamins were studied and compared to interactions with canonical purines. All the values of ΔGr unambiguously indicate that reactions with cisPt-diaqua are more preferable, but the comparison of ΔGr values obtained for compounds from vitamin B group and the ones characterizing complexes created by Guanine molecules indicates higher affinity of cisPt monomers toward purines.

CONCLUSION

Based on the observations, the regular intake of vitamin-rich beetroot or carrot juices is strongly discouraged during anticancer therapy using CisPt drug. To confirm the results of the performed computational study, detailed clinical trials should be performed.

All-atom simulations to studying metallodrugs/target interactions

Metallodrugs are extensively used to treat and diagnose distinct disease types. The unique physical-chemical properties of metal ions offer tantalizing opportunities to tailor effective scaffolds for selectively targeting specific biomolecules. Modern experimental techniques have collected a large body of structural data concerning the interactions of metallodrugs with their biomolecular targets, although being unable to exhaustively assess the molecular basis of their mechanism of action. In this scenario, the complementary use of accurate computational methods allows uncovering the minutiae of metallodrugs/targets interactions and their underlying mechanism of action at an atomic-level of detail. This knowledge is increasingly perceived as an invaluable requirement to rationally devise novel and selective metallodrugs. Building on literature studies, selected largely from the last 2 years, this compendium encompasses a cross-section of the current role, advances, and challenges met by computer simulations to decipher the mechanistic intricacies of prototypical metallodrugs.

Theoretical study of copper binding to GHK peptide

We report ligand field molecular mechanics, density functional theory and semi-empirical studies on the binding of Cu(II) to GlyHisLys (GHK) peptide. Following exhaustive conformational searching using molecular mechanics, we show that relative energy and geometry of conformations are in good agreement between GFN2-xTB semi-empirical and B3LYP-D DFT levels. Conventional molecular dynamics simulation of Cu-GHK shows the stability of the copper-peptide binding over 100 ps trajectory. Four equatorial bonds in 3N1O coordination remain stable throughout simulation, while a fifth in apical position from C-terminal carboxylate is more fluxional. We also show that the automated conformer and rotamer search algorithm CREST is able to correctly predict the metal binding position from a starting point consisting of separated peptide, copper and water.

Mechanistic Insight into the Binding Profile of DCVJ and α-Synuclein Fibril Revealed by Multiscale Simulations

Parkinson's disease (PD) is a serious neurodegenerative disease and is characterized by abnormal α-synuclein (α-syn) accumulation in Lewy bodies (LB) and Lewy neurites (LN), which makes α-syn an important imaging target for PD. An imaging probe that quantifies fibrillar α-syn can enhance the clinical diagnosis of PD and can also be used to evaluate the efficacy of therapeutics aimed at reducing the abnormal aggregation of the α-syn fibril in the brain. In this paper, we study the binding profile of fibrillar α-syn with a fluorescent probe 4-(dicyanovinyl)julolidine (DCVJ), which is being explored for identifying α-syn imaging agents. A multiscale simulation workflow including molecular docking, molecular dynamics, metadynamics, and QM/MM calculations was implemented. We find that DCVJ can bind to multiple sites of α-syn which are located either at the surface or in the core. Free energy calculations using implicit solvent models reveal that the most favorable binding mode for DCVJ is associated with the core binding site and is further confirmed by metadyamics simulation. Besides, a dynamic binding pathway is discovered, which reveals that DCVJ binds gradually into the core of the fibril passing through several intermediate states. The conformational arrest of the dicyano vinyl group in the fibrillar environment could explain the reason behind the fibril-specific fluorescence of DCVJ. Furthermore, based on hybrid QM/MM calculations, the molecular geometry of the dicyano vinyl group is found to be environment specific which explains why DCVJ serves as a staining agent for such fibrillar-like environments. Our results could be helpful for elucidating the binding mechanism of imaging tracers with the fibrillar form of α-syn and explain their fibrillar-specific optical properties, a knowledge that in turn can be used to guide the design and development of compounds with higher affinity and selectivity for α-syn using structure-based strategies.

Square planar or octahedral after all? The indistinct solvation of platinum(ii)

The solvation structures of Pd(ii) and Pt(ii) are typically reduced to the well-known square-planar structural motif, although it has been shown, in both experimental and theoretical investigations, that these solutes demonstrate the affinity to bind ligand molecules at elongated distance in axial coordination sites.

Conformational polymorphism or structural invariance in DNA photoinduced lesions: implications for repair rates

DNA photolesions constitute a particularly deleterious class of molecular defects responsible for the insurgence of a vast majority of skin malignant tumors. Dimerization of two adjacent thymines or cytosines mostly gives rise to cyclobutane pyrimidine dimers (CPD) and pyrimidine(6-4)pyrimidone 64-PP as the most common defects. We perform all-atom classical simulations, up to 2 μs, of CPD and 64-PP embedded in a 16-bp duplex, which reveal the constrasted behavior of the two lesions. In particular we evidence a very limited structural deformation induced by CPD while 64-PP is characterized by a complex structural polymorphism. Our simulations also allow to unify the contrasting experimental structural results obtained by nuclear magnetic resonance or Förster Resonant Energy Transfer method, showing that both low and high bent structures are indeed accessible. These contrasting behaviors can also explain repair resistance or the different replication obstruction, and hence the genotoxicity of these two photolesions.

The Role of the DNA Backbone in Minor-Groove Ligand Binding

Molecular recognition between ligands and nucleic acids plays a key role in therapeutic activity. Some molecules interact with DNA in a nonbonded manner through intercalation or through the DNA grooves. The recognition of minor-groove binders is attributed to a set of hydrogen-bonding interactions between the binders and the hydrogen-bond-acceptor groups on the groove floor and walls. It is commonly considered that interactions with the sugar groups of the DNA backbone are insignificant and do not contribute to the binding affinity or the specificity. However, our group has found that the deoxyribose rings have a central function in the recognition and the intercalation of metal complexes into DNA. Herein, we determined the specific interactions between the binder CGP 40215A and the minor-groove atoms, based on the local properties of electron density. We found that specific interactions between the deoxyribose moiety within the backbone of DNA and the binder are essential for molecular recognition, and they are responsible for one third of the interaction energy.

Molecular dynamics simulation and free energy analysis of the interaction of platinum-based anti-cancer drugs with DNA

Cisplatin and oxaliplatin are two widely-used anti-cancer drugs which covalently bind to a same location in DNA strands. Platinum agents make intrastrand and interstrand cross-links with the N7 atoms of guanine nucleotides which prevent DNA from polymerization by causing a distortion in the double helix. Molecular dynamics simulations and free energy calculations were carried out to investigate the binding of two platinum-based anti-cancer drugs with DNA. We compared the binding of these drugs which differ in their carrier ligands, and hence their potential interactions with DNA. When a platinum agent binds to nucleotides, it causes a high amount of deformation in DNA structure. To find the extent of deformation, torsion angles and base pair and groove parameters of DNA were considered. These parameters were compared with normal B-DNA which was considered as the undamaged DNA. The formation of hydrogen bonds between drugs and DNA nucleotides was examined in solution. It was shown that oxaliplatin forms more hydrogen bonds than cisplatin. Our results confirm that the structure of the platinated DNA rearranges significantly and cisplatin tries to deform DNA more than oxaliplatin. The binding free energies were also investigated to understand the affinities, types and the contributions of interactions between drugs and DNA. It was concluded that oxaliplatin tendency for binding to DNA is more than cisplatin in solvent environment. The binding free energy was calculated based on the MM/PBSA and MM/GBSA methods and the results of QM/MM calculations verified them.

Interactions of the "piano-stool" [ruthenium(II)(η(6) -arene)(quinolone)Cl](+) complexes with water; DFT computational study

Full optimizations of stationary points along the reaction coordinate for the hydration of several quinolone Ru(II) half-sandwich complexes were performed in water environment using the B3PW91/6-31+G(d)/PCM/UAKS method. The role of diffuse functions (especially on oxygen) was found crucial for correct geometries along the reaction coordinate. Single-point (SP) calculations were performed at the B3LYP/6-311++G(2df,2pd)/DPCM/saled-UAKS level. In the first part, two possible reaction mechanisms-associative and dissociative were compared. It was found that the dissociative mechanism of the hydration process is kinetically slightly preferred. Another important conclusion concerns the reaction channels. It was found that substitution of chloride ligand (abbreviated in the text as dechlorination reaction) represents energetically and kinetically the most feasible pathway. In the second part the same hydration reaction was explored for reactivity comparison of the Ru(II)-complexes with several derivatives of nalidixic acid: cinoxacin, ofloxacin, and (thio)nalidixic acid. The hydration process is about four orders of magnitude faster in a basic solution compared to neutral/acidic environment with cinoxacin and nalidixic acid as the most reactive complexes in the former and latter environments, respectively. The explored hydration reaction is in all cases endergonic; nevertheless the endergonicity is substantially lower (by ∼6 kcal/mol) in basic environment. © 2016 Wiley Periodicals, Inc.

Empirical force field for cisplatin based on quantum dynamics data: case study of new parameterization scheme for coordination compounds

Parameterization of molecular complexes containing a metallic compound, such as cisplatin, is challenging due to the unconventional coordination nature of the bonds which involve platinum atoms. In this work, we develop a new methodology of parameterization for such compounds based on quantum dynamics (QD) calculations. We show that the coordination bonds and angles are more flexible than in normal covalent compounds. The influence of explicit solvent is also shown to be crucial to determine the flexibility of cisplatin in quantum dynamics simulations. Two empirical topologies of cisplatin were produced by fitting its atomic fluctuations against QD in vacuum and QD with explicit first solvation shell of water molecules respectively. A third topology built in a standard way from the static optimized structure was used for comparison. The later one leads to an excessively rigid molecule and exhibits much smaller fluctuations of the bonds and angles than QD reveals. It is shown that accounting for the high flexibility of cisplatin molecule is needed for adequate description of its first hydration shell. MD simulations with flexible QD-based topology also reveal a significant decrease of the barrier of passive diffusion of cisplatin accross the model lipid bilayer. These results confirm that flexibility of organometallic compounds is an important feature to be considered in classical molecular dynamics topologies. Proposed methodology based on QD simulations provides a systematic way of building such topologies.

Asymmetric environmental effects on the structure and vibrations of cis-[Pt(NH3)2Cl2] in condensed phases

We report the structural and vibrational properties of anticancer drug cisplatin (cis-[Pt(NH3)2Cl2]) in gas phase, in solid phase, and in aqueous solution using density functional theory (DFT) calculations, quantum mechanical/molecular mechanical (QM/MM) molecular dynamics, and effective normal modes analysis. In contrast with the gas-phase case, asymmetric hydrogen bonding environments are found in both solid phase and aqueous solution. It is shown that the discrepancy of the molecular geometry between previous gas phase calculations and the X-ray crystal structure can be resolved by considering intermolecular hydrogen bonds in the calculations of solid phase. In addition, our simulations in solid phase and aqueous solution reveal that asymmetric environmental effects lead to several spectral features observed in experiments, such as the blue-shift in the N-H stretching region and the frequency splitting of NH3 symmetric deformation modes. Furthermore, a similar decoupling and localization of several vibrational modes of cisplatin is found in solid phase and aqueous solution, in comparison to those of O-H stretching modes of water molecules in liquid water [ J. Phys. Chem. Lett. 2013 , 4 ( 19 ), 3245 - 3250 ].

QM/MM molecular dynamics studies of metal binding proteins

Mixed quantum-classical (quantum mechanical/molecular mechanical (QM/MM)) simulations have strongly contributed to providing insights into the understanding of several structural and mechanistic aspects of biological molecules. They played a particularly important role in metal binding proteins, where the electronic effects of transition metals have to be explicitly taken into account for the correct representation of the underlying biochemical process. In this review, after a brief description of the basic concepts of the QM/MM method, we provide an overview of its capabilities using selected examples taken from our work. Specifically, we will focus on heme peroxidases, metallo-β-lactamases, α-synuclein and ligase ribozymes to show how this approach is capable of describing the catalytic and/or structural role played by transition (Fe, Zn or Cu) and main group (Mg) metals. Applications will reveal how metal ions influence the formation and reduction of high redox intermediates in catalytic cycles and enhance drug metabolism, amyloidogenic aggregate formation and nucleic acid synthesis. In turn, it will become manifest that the protein frame directs and modulates the properties and reactivity of the metal ions.

NMR solution structure of rat aβ(1-16): toward understanding the mechanism of rats' resistance to Alzheimer's disease

In an attempt to reveal the mechanism of rats' resistance to Alzheimer's disease, we determined the structure of the metal-binding domain 1-16 of rat β-amyloid (rat Aβ(1-16)) in solution in the absence and presence of zinc ions. A zinc-induced dimerization of the domain was detected. The zinc coordination site was found to involve residues His-6 and His-14 of both peptide chains. We used experimental restraints obtained from analyses of NMR and isothermal titration calorimetry data to perform structure calculations. The calculations employed an explicit water environment and a simulated annealing molecular-dynamics protocol followed by quantum-mechanical/molecular-mechanical optimization. We found that the C-tails of the two polypeptide chains of the rat Aβ(1-16) dimer are oriented in opposite directions to each other, which hinders the assembly of rat Aβ dimers into oligomeric aggregates. Thus, the differences in the structure of zinc-binding sites of human and rat Aβ(1-16), their ability to form regular cross-monomer bonds, and the orientation of their hydrophobic C-tails could be responsible for the resistance of rats to Alzheimer's disease.

Theoretical analysis of trans-[PtCl2(NH3)(thiazole)] and trans-[PtCl2(thiazole)2] binding to biological targets — Factors influence binding kinetics and adduct stability

Full reaction energy profiles for trans-[PtCl2(NH3)(thiazole)] and trans-[PtCl2(thiazole)2] binding to sulfur- and nitrogen-containing biorelevant ligands were constructed by the density functional theory (DFT) method. Calculated results demonstrate that trans-platinum complexes can interact with biological targets, affording cis and trans products via very similar transition states. For different substituents, sulfur-containing ligands constitute kinetically preferred targets for platination, whereas the platination of nitrogen-containing ligands is more favorable thermodynamically. This is consistent with previous experimental studies. Calculated results also suggest that the trans effect, the influence of the ligand, the size of the ligand, and hydrogen bonding play important roles in binding kinetics and stabilizing adducts.

Cisplatin interaction with amino acids cysteine and methionine from gas phase to solutions with constant pH

This work is focused on the computational studies of reactions of hydrated forms of cisplatin with sulphur-containing amino acids cysteine and methionine. First, the appropriate model for solvation of the examined complexes was searched for. The suggested procedure employs the B3LYP density functional, 6-311++G(2df,2pd) basis set with Stuttgart-Dresden pseudopotentials on heavy atoms, the D-PCM solvation model and the UAKS cavity which uses more realistic NPA partial charges instead of formal partial charges for platinum ligands. In the second part this model is applied to the evaluation of the Legendre transformed reaction Gibbs free energy of cisplatin with cysteine and methionine in solution at constant pH.