Homodimeric complexes of the 90-231 human prion: a multilayered computational study based on FMO/GRID-DRY approach

A computational insight on the aromatic amino acids conjugation with [Cp*Rh(H2O)3]2+ by using the meta-dynamics/FMO3 approach

CONTEXT

Rh(III) complexes demonstrated to exert promising pharmacological effects with potential applications as anti-cancer, anti-bacterial, and antimicrobial agents. One important Rh(III)-ligand is the pentamethylcyclopentadienyl (Cp*) group forming in water the [Cp*Rh(H2O)3]2+ complex. Among of its attractive chemical properties is the ability to react specifically with Tyr amino acid side chain of G-protein-coupled receptor (GPCR) peptides by means of highly chemoselective bioconjugation reaction, at room temperature and at pH 5-6. In this computational work, in order to deepen the mechanism of this chemoselective conjugation, we study the ligand exchange reaction between [Cp*Rh(H2O)3]2+ and three small molecules, namely p-cresol, 3-methylimidazole, and toluene, selected as mimetic of aromatic side chains of tyrosine (Tyr), tryptophan (Trp) and phenylalanine (Phe), respectively. Our outcomes suggest that the high selectivity for Tyr side chain might be related to OH group able to affect both thermodynamic and kinetic of ligand exchange reaction, due to its ability to act as both H bond acceptor and donor. These mechanistic aspects can be used to design new metal drugs containing the [Cp*Rh]2+ scaffold targeting specifically Tyr residues involved in biological/pathological processes such as phosphorylation by means of Tyr-kinase enzyme and protein-protein interactions.

METHODS

The geometry of three encounter complexes and product adducts were optimized at the B3LYP//CPCM/ωB97X-D level of theory, adopting the 6-311+G(d,p) basis set for all non-metal atoms and the LANL2DZ pseudopotential for the Rh atom. Meta-dynamics RMSD (MTD(RMSD)) calculations at GFN2-xTB level of theory were performed in NVT conditions at 298.15 K to investigate the bioconjugation reactions (simulation time: 100 ps; integration step 2.0; implicit solvent model: GBSA). The MTD(RMSD) simulation was performed in two replicates for each encounter complex. Final representative subsets of 100 structures for each run were gained with a sampling rate of 1 ps and analyzed by performing single point calculations using the FMO3 method at RI-MP2/6-311G//PCM[1] level of theory, adopting the MCP-TZP core potential for Rh atom.

Improving the accuracy of the FMO binding affinity prediction of ligand-receptor complexes containing metals

Polarization and charge transfer strongly characterize the ligand-receptor interaction when metal atoms are present, as for the Au(I)-biscarbene/DNA G-quadruplex complexes. In a previous work (J Comput Aided Mol Des2022, 36, 851-866) we used the ab initio FMO2 method at the RI-MP2/6-31G* level of theory with the PCM [1] solvation approach to calculate the binding energy (ΔEFMO) of two Au(I)-biscarbene derivatives, [Au(9-methylcaffein-8-ylidene)2]+ and [Au(1,3-dimethylbenzimidazole-2-ylidene)2]+, able to interact with DNA G-quadruplex motif. We found that ΔEFMO and ligand-receptor pair interaction energies (EINT) show very large negative values making the direct comparison with experimental data difficult and related this issue to the overestimation of the embedded charge transfer energy between fragments containing metal atoms. In this work, to improve the accuracy of the FMO method for predicting the binding affinity of metal-based ligands interacting with DNA G-quadruplex (Gq), we assess the effect of the following computational features: (i) the electron correlation, considering the Hartree-Fock (HF) and a post-HF method, namely RI-MP2; (ii) the two (FMO2) and three-body (FMO3) approaches; (iii) the basis set size (polarization functions and double-ζ vs. triple-ζ) and (iv) the embedding electrostatic potential (ESP). Moreover, the partial screening method was systematically adopted to simulate the solvent screening effect for each calculation. We found that the use of the ESP computed using the screened point charges for all atoms (ESP-SPTC) has a critical impact on the accuracy of both ΔEFMO and EINT, eliminating the overestimation of charge transfer energy and leading to energy values with magnitude comparable with typical experimental binding energies. With this computational approach, EINT values describe the binding efficiency of metal-based binders to DNA Gq more accurately than ΔEFMO. Therefore, to study the binding process of metal containing systems with the FMO method, the adoption of partial screening solvent method combined with ESP-SPCT should be considered. This computational protocol is suggested for FMO calculations on biological systems containing metals, especially when the adoption of the default ESP treatment leads to questionable results.

IRMPD spectroscopy and quantum-chemical simulations of the reaction products of cisplatin with the dipeptide CysGly

The inorganic antineoplastic drug cisplatin was made to react in solution with the dipeptide cysteinylglycine (CysGly), chosen as a functional model of glutathione, and the reaction products were analyzed using electrospray ionization mass spectrometry (ESI-MS). Selected complexes, i.e., the primary substitution product cis-[PtCl(NH3)2(CysGly)]+ and the chelate cis-[PtCl(NH3)(CysGly)]+, were submitted to IR multiple photon dissociation (IRMPD) spectroscopy obtaining their vibrational features. The experimental IR ion spectra were compared with the calculated IR absorptions of different plausible isomeric families, finding CysGly to bind preferentially platinum(II) via its deprotonated thiolic group in the monovalent complex, cis-[PtCl(NH3)2(CysGly)]+, and to evolve in the S,N-bound chelate structure cis-[PtCl(NH3)(CysGly)]+ through the SH and NH2 functionality of the cysteine residue. Moreover, our findings indicate that the platination reaction does not affect the CysGly peptide bond, which remains in its trans configuration. These results provide additional insights into the reactivity of Pt(II)-complexes with glutathione which is involved in cellular cisplatin resistance.

Inorganic Drugs as a Tool for Protein Structure Solving and Studies on Conformational Changes

Inorganic drugs are capable of tight interactions with proteins through coordination towards aminoacidic residues, and this feature is recognized as a key aspect for their pharmacological action. However, the "protein metalation process" is exploitable for solving the phase problem and structural resolution. In fact, the use of inorganic drugs bearing specific metal centers and ligands capable to drive the binding towards the desired portions of the protein target could represent a very intriguing and fruitful strategy. In this context, a theoretical approach may further contribute to solve protein structures and their refinement. Here, we delineate the main features of a reliable experimental-theoretical integrated approach, based on the use of metallodrugs, for protein structure solving.

Structural Reshaping of the Zinc-Finger Domain of the SARS-CoV-2 nsp13 Protein Using Bismuth(III) Ions: A Multilevel Computational Study

The identification of novel therapeutics against the pandemic SARS-CoV-2 infection is an indispensable new address of current scientific research. In the search for anti-SARS-CoV-2 agents as alternatives to the vaccine or immune therapeutics whose efficacy naturally degrades with the occurrence of new variants, the salts of Bi³⁺ have been found to decrease the activity of the Zn²⁺-dependent non-structural protein 13 (nsp13) helicase, a key component of the SARS-CoV-2 molecular tool kit. Here, we present a multilevel computational investigation based on the articulation of DFT calculations, classical MD simulations, and MIF analyses, focused on the examination of the effects of Bi³⁺/Zn²⁺ exchange on the structure and molecular interaction features of the nsp13 protein. Our calculations confirmed that Bi³⁺ ions can replace Zn²⁺ in the zinc-finger metal centers and cause slight but appreciable structural modifications in the zinc-binding domain of nsp13. Nevertheless, by employing an in-house-developed ATOMIF tool, we evidenced that such a Bi³⁺/Zn²⁺ exchange may decrease the extension of a specific hydrophobic portion of nsp13, responsible for the interaction with the nsp12 protein. The present study provides for a detailed, atomistic insight into the potential anti-SARS-CoV-2 activity of Bi³⁺ and, more generally, evidences the hampering of the nsp13–nsp12 interaction as a plausible therapeutic strategy.

The Zn²⁺/Bi³⁺ exchange in the zinc finger domains of SARS-CoV-2 nsp13 hampers the viral machinery