An insight of early PrP-E200K aggregation by combined molecular dynamics/fragment molecular orbital approaches

Combining the Fragment Molecular Orbital and GRID Approaches for the Prediction of Ligand-Metalloenzyme Binding Affinity: The Case Study of hCA II Inhibitors

Polarization and charge-transfer interactions play an important role in ligand-receptor complexes containing metals, and only quantum mechanics methods can adequately describe their contribution to the binding energy. In this work, we selected a set of benzenesulfonamide ligands of human Carbonic Anhydrase II (hCA II)-an important druggable target containing a Zn2+ ion in the active site-as a case study to predict the binding free energy in metalloprotein-ligand complexes and designed specialized computational methods that combine the ab initio fragment molecular orbital (FMO) method and GRID approach. To reproduce the experimental binding free energy in these systems, we adopted a machine-learning approach, here named formula generator (FG), considering different FMO energy terms, the hydrophobic interaction energy (computed by GRID) and logP. The main advantage of the FG approach is that it can find nonlinear relations between the energy terms used to predict the binding free energy, explicitly showing their mathematical relation. This work showed the effectiveness of the FG approach, and therefore, it might represent an important tool for the development of new scoring functions. Indeed, our scoring function showed a high correlation with the experimental binding free energy (R2 = 0.76-0.95, RMSE = 0.34-0.18), revealing a nonlinear relation between energy terms and highlighting the relevant role played by hydrophobic contacts. These results, along with the FMO characterization of ligand-receptor interactions, represent important information to support the design of new and potent hCA II inhibitors.

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.

Preventive or promotive effects of PRNP polymorphic heterozygosity on the onset of prion disease

The polymorphic heterozygosity of PRNP at codon 129 or 219 prevents the onset of sporadic Creutzfeldt-Jakob disease (sCJD). We investigated the association between polymorphic genotypes at codon 129 or 219 and comprehensive prion disease onset using non-CJD as a reference. EK heterozygotes at codon 219, versus EE homozygotes, showed a preventive effect on the extensive prion diseases-sCJD, genetic CJD (gCJD) with V180I or M232R mutation, and Gerstmann-Straussler-Scheinker disease with P102L mutation. No preventive effect was observed for E200K-gCJD and dura-grafted CJD (dCJD) in 129 MV and 219 EK heterozygotes. It was suggested that unlike other prion diseases, E200K-gCJD may not benefit from the preventive effect of 219 EK heterozygosity because complementary electrostatic interactions between PrP molecules at K200 and E219 might make homodimer formation easier. Comparison of sCJD and dCJD indicates that 219 EK heterozygosity strongly inhibits de novo synthesis of PrP^Sc (initial PrP^Sc formation), but does not inhibit accelerated propagation of existing PrP^Sc.

Effects of the pathological E200K mutation on human prion protein: A computational screening and molecular dynamics approach

The human prion protein gene (PRNP) is mapped to the short arm of chromosome 20 (20pter-12). Prion disease is associated with mutations in the prion protein-encoding gene sequence. Earlier studies found that the mutation G127V in the PRNP increases protein stability. In contrast, the mutation E200K, which has the highest mutation rate in the prion protein, causes Creutzfeldt-Jakob disease (CJD) in humans and induces protein aggregation. We aimed to identify the structural mechanisms of E200k and G127V mutations causing CJD. We used a variety of bioinformatic algorithms, including SIFT, PolyPhen, I-Mutant, PhD-SNP, and SNP& GO, to predict the association of the E200K mutation with prion disease. MD simulation is performed, and graphs for root mean square deviation, root mean square fluctuation, radius of gyration, DSSP, principal component analysis, porcupine, and free energy landscape are generated to confirm and prove the stability of the wild-type and mutant protein structures. The protein is analyzed for aggregation, and the results indicate more fluctuations in the protein structure during the simulation owing to the E200K mutation; however, the G127V mutation makes the protein structure stable against aggregation during the simulation.

Human prion protein: exploring the thermodynamic stability and structural dynamics of its pathogenic mutants

Human familial prion diseases are known to be associated with different single-point mutants of the gene coding for prion protein with a primary focus at several locations of the globular domain. We have identified 12 different single-point pathogenic mutants of human prion protein (HuPrP) with the help of extensive perturbations/mutation technique at multiple locations of HuPrP sequence related to potentiality towards conformational disorders. Among these, some of the mutants include pathogenic variants that corroborate well with the literature reported proteins while majority include some unique single-point mutants that are either not explicitly studied early or studied for variants with different residues at the specific position. Primarily, our study sheds light on the unfolding mechanism of the above mentioned mutants in depth. Besides, we could identify some mutants under investigation that demonstrates not only unfolding of the helical structures but also extension and generation of the β-sheet structures and or simultaneously have highly exposed hydrophobic surface which is assumed to be linked with the production of aggregate/fibril structures of the prion protein. Among the identified mutants, Q212E needs special attention due to its maximum exposure of hydrophobic core towards solvent and E200Q is found to be important due to its maximum extent of β-content. We are also able to identify different respective structural conformations of the proteins according to their degree of structural unfolding and those conformations can be extracted and further studied in detail. Communicated by Ramaswamy H. Sarma.

The FMO2 analysis of the ligand-receptor binding energy: the Biscarbene-Gold(I)/DNA G-Quadruplex case study

In this work, the ab initio fragment molecular orbital (FMO) method was applied to calculate and analyze the binding energy of two biscarbene-Au(I) derivatives, [Au(9-methylcaffein-8-ylidene)₂]⁺ and [Au(1,3-dimethylbenzimidazol-2-ylidene)₂]⁺, to the DNA G-Quadruplex structure. The FMO2 binding energy considers the ligand-receptor complex as well as the isolated forms of energy-minimum state of ligand and receptor, providing a better description of ligand-receptor affinity compared with simple pair interaction energies (PIE). Our results highlight important features of the binding process of biscarbene-Au(I) derivatives to DNA G-Quadruplex, indicating that the total deformation-polarization energy and desolvation penalty of the ligands are the main terms destabilizing the binding. The pair interaction energy decomposition analysis (PIEDA) between ligand and nucleobases suggest that the main interaction terms are electrostatic and charge-transfer energies supporting the hypothesis that Au(I) ion can be involved in π-cation interactions further stabilizing the ligand-receptor complex. Moreover, the presence of polar groups on the carbene ring, as C = O, can improve the charge-transfer interaction with K⁺ ion. These findings can be employed to design new powerful biscarbene-Au(I) DNA-G quadruplex binders as promising anticancer drugs. The procedure described in this work can be applied to investigate any ligand-receptor system and is particularly useful when the binding process is strongly characterized by polarization, charge-transfer and dispersion interactions, properly evaluated by ab initio methods.

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 Bi3+ have been found to decrease the activity of the Zn2+-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 Bi3+/Zn2+ exchange on the structure and molecular interaction features of the nsp13 protein. Our calculations confirmed that Bi3+ ions can replace Zn2+ 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 Bi3+/Zn2+ 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 Bi3+ and, more generally, evidences the hampering of the nsp13-nsp12 interaction as a plausible therapeutic strategy.

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

The molecular interaction properties and aggregation capabilities disclosed by PrP-E200K, a pathogenic mutant of the human prion protein, were investigated in detail using multilayered computational approaches. In a previous work, we reported that the electrostatic complementarity between region1 (negative) and region3 (positive) has been assumed to lead to a head-to tail interaction between 120 and 231 PrP-E200K units and to initiation of the aggregation process. In this work, we extended the PrP-E200K structure by including the unstructured 90-120 segment which was found to assume different conformations. Plausible models of 90-231 PrP-E200K dimers were calculated and analyzed in depth to identify the nature of the involved protein-protein interactions. The unstructured 90-120 segment was found to extend the positively charged region3 involved in the association of PrP-E200K units which resulted to be driven by hydrophobic interactions. The combination of molecular dynamics, protein-protein docking, grid-based mapping, and fragment molecular orbital approaches allowed us to provide a plausible mechanism of the early state of 90-231 PrP-E200K aggregation, considered a preliminary step of amyloid conversion.

Hampering the early aggregation of PrP-E200K protein by charge-based inhibitors: a computational study

A multilayered computational workflow was designed to identify a druggable binding site on the surface of the E200K pathogenic mutant of the human prion protein, and to investigate the effect of the binding of small molecules in the inhibition of the early aggregation of this protein. At this purpose, we developed an efficient computational tool to scan the molecular interaction properties of a whole MD trajectory, thus leading to the characterization of plausible binding regions on the surface of PrP-E200K. These structural data were then employed to drive structure-based virtual screening and fragment-based approaches to the seeking of small molecular binders of the PrP-E200K. Six promising compounds were identified, and their binding stabilities were assessed by MD simulations. Therefore, analyses of the molecular electrostatic potential similarity between the bound complexes and unbound protein evidenced their potential activity as charged-based inhibitors of the PrP-E200K early aggregation.

Characterization of PD-L1 binding sites by a combined FMO/GRID-DRY approach

The programmed cell death protein 1 (PD-1) and its ligand, PD-L1, constitute an important co-inhibitory immune checkpoint leading to downregulation of immune system. Tumor cells developed a strategy to trigger PD-1/PD-L1 pathway reducing the T cell anticancer activity. Anti-PD-L1 small drugs, generally with improved pharmacokinetic and technological profiles than monoclonal antibodies, became an attractive research topic. Nevertheless, still few works have been published on the chemical features of possible binding sites. In this work, we applied a novel computational protocol based on the combination of the ab initio Fragment Molecular Orbital (FMO) method and a newly developed GRID-DRY approach in order to characterize the PD-L1 binding sites, starting from PD-1/PD-L1 and PD-L1/BMS-ligands (Bristol-Mayers Squibb ligands) complexes. The FMO method allows the calculation of the pair-residues as well as the ligand-residues interactions with ab initio accuracy, whereas the GRID-DRY approach is an effective tool to investigate hydrophobic interactions, not easily detectable by ab initio methods. The present GRID-DRY protocol is able to determine the energy contributions of each ligand atoms to each hydrophobic interaction, both qualitatively and quantitatively. We were also able to identify the three specific hot regions involved in PD-1/PD-L1 protein-protein interaction and in PD-L1/BMS-ligand interactions, in agreement with preceding theoretical/experimental results, and to suggest a specific pharmacophore for PD-L1 inhibitors.