Computational insight into the interaction of oxaliplatin with insulin

TALAIA: a 3D visual dictionary for protein structures

MOTIVATION

Graphical analysis of the molecular structure of proteins can be very complex. Full-atom representations retain most geometric information but are generally crowded, and key structural patterns can be challenging to identify. Non-full-atom representations could be more instructive on physicochemical aspects but be insufficiently detailed regarding shapes (e.g. entity beans-like models in coarse grain approaches) or simple properties of amino acids (e.g. representation of superficial electrostatic properties). In this work, we present TALAIA a visual dictionary that aims to provide another layer of structural representations.TALAIA offers a visual grammar that combines simple representations of amino acids while retaining their general geometry and physicochemical properties. It uses unique objects, with differentiated shapes and colors to represent amino acids. It makes easier to spot crucial molecular information, including patches of amino acids or key interactions between side chains. Most conventions used in TALAIA are standard in chemistry and biochemistry, so experimentalists and modelers can rapidly grasp the meaning of any TALAIA depiction.

RESULTS

We propose TALAIA as a tool that renders protein structures and encodes structure and physicochemical aspects as a simple visual grammar. The approach is fast, highly informative, and intuitive, allowing the identification of possible interactions, hydrophobic patches, and other characteristic structural features at first glance. The first implementation of TALAIA can be found at https://github.com/insilichem/talaia.

Successes and challenges in multiscale modelling of artificial metalloenzymes: the case study of POP-Rh2 cyclopropanase

Molecular modelling applications in metalloenzyme design are still scarce due to a series of challenges. On top of that, the simulations of metal-mediated binding and the identification of catalytic competent geometries require both large conformational exploration and simulation of fine electronic properties. Here, we demonstrate how the incorporation of new tools in multiscale strategies, namely substrate diffusion exploration, allows taking a step further. As a showcase, the enantioselective profiles of the most outstanding variants of an artificial Rh₂-based cyclopropanase (GSH, HFF and RFY) developed by Lewis and co-workers (Nat. Commun., 2015, 6, 7789 and Nat. Chem., 2018, 10, 318-324) have been rationalized. DFT calculations on the free-cofactor-mediated process identify the carbene insertion and the cyclopropanoid formation as crucial events, the latter being the enantiodetermining step, which displays up to 8 competitive orientations easily altered by the protein environment. The key intermediates of the reaction were docked into the protein scaffold showing that some mutated residues have direct interaction with the cofactor and/or the co-substrate. These interactions take the form of a direct coordination of Rh in GSH and HFF and a strong hydrophobic patch with the carbene moiety in RFY. Posterior molecular dynamics sustain that the cofactor induces global re-arrangements of the protein. Finally, massive exploration of substrate diffusion, based on the GPathFinder approach, defines this event as the origin of the enantioselectivity in GSH and RFY. For HFF, fine molecular dockings suggest that it is likely related to local interactions upon diffusion. This work shows how modelling of long-range mutations on the catalytic profiles of metalloenzymes may be unavoidable and software simulating substrate diffusion should be applied.

Introducing the penicillamine moiety into a metallopeptide mimicking the NiSOD enzyme: electronic and kinetic effects

The novel NiSOD related metallopeptide incorporates penicillamine moiety in the active center which alters both the electronic and kinetic features.

Hybrid Cyclobutane/Proline-Containing Peptidomimetics: The Conformational Constraint Influences Their Cell-Penetration Ability

A new family of hybrid β,γ-peptidomimetics consisting of a repetitive unit formed by a chiral cyclobutane-containing trans-β-amino acid plus a N^α-functionalized trans-γ-amino-l-proline joined in alternation were synthesized and evaluated as cell penetrating peptides (CPP). They lack toxicity on the human tumoral cell line HeLa, with an almost negligible cell uptake. The dodecapeptide showed a substantial microbicidal activity on Leishmania parasites at 50 µM but with a modest intracellular accumulation. Their previously published γ,γ-homologues, with a cyclobutane γ-amino acid, showed a well-defined secondary structure with an average inter-guanidinium distance of 8–10 Å, a higher leishmanicidal activity as well as a significant intracellular accumulation. The presence of a very rigid cyclobutane β-amino acid in the peptide backbone precludes the acquisition of a defined conformation suitable for their cell uptake ability. Our results unveiled the preorganized charge-display as a relevant parameter, additional to the separation among the charged groups as previously described. The data herein reinforce the relevance of these descriptors in the design of CPPs with improved properties.

Rationalizing the Decavanadate(V) and Oxidovanadium(IV) Binding to G-Actin and the Competition with Decaniobate(V) and ATP

The experimental data collected over the past 15 years on the interaction of decavanadate(V) (V₁₀O₂₈^6-; V₁₀), a polyoxometalate (POM) with promising anticancer and antibacterial action, with G-actin, were rationalized by using several computational approaches (docking, density functional theory (DFT), and molecular dynamics (MD)). Moreover, a comparison with the isostructural and more stable decaniobate(V) (Nb₁₀O₂₈^6-; Nb₁₀) was carried out. Four binding sites were identified, named α, β, γ, and δ, the site α being the catalytic nucleotide site located in the cleft of the enzyme at the interface of the subdomains II and IV. It was observed that the site α is preferred by V₁₀, whereas Nb₁₀ is more stable at the site β; this indicates that, differently from other proteins, G-actin could contemporaneously bind the two POMs, whose action would be synergistic. Both decavanadate and decaniobate induce conformational rearrangements in G-actin, larger for V₁₀ than Nb₁₀. Moreover, the binding mode of oxidovanadium(IV) ion, V^IVO²⁺, formed upon the reduction of decavanadate(V) by the -SH groups of accessible cysteine residues, is also found in the catalytic site α with (His161, Asp154) coordination; this adduct overlaps significantly with the region where ATP is bound, accounting for the competition between V₁₀ and its reduction product V^IVO²⁺ with ATP, as previously observed by EPR spectroscopy. Finally, the competition with ATP was rationalized: since decavanadate prefers the nucleotide site α, Ca²⁺-ATP displaces V₁₀ from this site, while the competition is less important for Nb₁₀ because this POM shows a higher affinity for β than for site α. A relevant consequence of this paper is that other metallodrug-protein systems, in the absence or presence of eventual inhibitors and/or competition with molecules of the organism, could be studied with the same approach, suggesting important elements for an explanation of the biological data and a rational drug design.

Integrated experimental/computational approaches to characterize the systems formed by vanadium with proteins and enzymes

An integrated instrumental/computational approach to characterize metallodrug–protein adducts at the molecular level is reviewed. A series of applications are described, focusing on potential vanadium drugs with a generalization to other metals.

BioMetAll: Identifying Metal-Binding Sites in Proteins from Backbone Preorganization

With a large amount of research dedicated to decoding how metallic species bind to proteins, in silico methods are interesting allies for experimental procedures. To date, computational predictors mostly work by identifying the best possible sequence or structural match of the target protein with metal-binding templates. These approaches are fundamentally focused on the first coordination sphere of the metal. Here, we present the BioMetAll predictor that is based on a different postulate: the formation of a potential metal-binding site is related to the geometric organization of the protein backbone. We first report the set of convenient geometric descriptors of the backbone needed for the algorithm and their parameterization from a statistical analysis. Then, the successful benchmark of BioMetAll on a set of more than 90 metal-binding X-ray structures is presented. Because BioMetAll allows structural predictions regardless of the exact geometry of the side chains, it appears extremely valuable for systems whose structures (either experimental or theoretical) are not optimal for metal-binding sites. We report here its application on three different challenging cases: (i) the modulation of metal-binding sites during conformational transition in human serum albumin, (ii) the identification of possible routes of metal migration in hemocyanins, and (iii) the prediction of mutations to generate convenient metal-binding sites for de novo biocatalysts. This study shows that BioMetAll offers a versatile platform for numerous fields of research at the interface between inorganic chemistry and biology and allows to highlight the role of the preorganization of the protein backbone as a marker for metal binding. BioMetAll is an open-source application available at https://github.com/insilichem/biometall.

Host-Guest Complexation of Oxaliplatin and Para-Sulfonatocalix[n]Arenes for Potential Use in Cancer Therapy

P-sulfonatocalix[n]arenes have demonstrated a great potential for encapsulation of therapeutic drugs via host-guest complexation to improve solubility, stability, and bioavailability of encapsulated drugs. In this work, guest-host complexes of a third-generation anticancer drug (oxaliplatin) and p-4-sulfocalix[n]arenes (n = 4 and 6; p-SC4 and p-SC6, respectively) were prepared and investigated, using 1H NMR, UV, Job's plot analysis, and DFT calculations, for use as cancer therapeutics. The peak amplitude of the prepared host-guest complexes was linearly proportional to the concentration of oxaliplatin in the range of 1.0 × 10-5 M-1 to 2.1 × 10-4 M-1. The reaction stoichiometry between either p-SC4 or p-SC6 and oxaliplatin in the formed complexes was 1:1. The stability constants for the complexes were 5.07 × 104 M-1 and 6.3 × 104 M-1. These correspond to complexation free energy of -6.39 and -6.52 kcal/mol for p-SC4 and p-SC6, respectively. Complexation between oxaliplatin and p-SC4 or p-SC6 was found to involve hydrogen bonds. Both complexes exhibited enhanced biological and high cytotoxic activities against HT-29 colorectal cells and MCF-7 breast adenocarcinoma compared to free oxaliplatin, which warrants further investigation for cancer therapy.

Chiral Cyclobutane-Containing Cell-Penetrating Peptides as Selective Vectors for Anti-Leishmania Drug Delivery Systems

Two series of new hybrid γ/γ-peptides, γ-CC and γ-CT, formed by (1S,2R)-3-amino-2,2,dimethylcyclobutane-1-carboxylic acid joined in alternation to a Nα-functionalized cis- or trans-γ-amino-l-proline derivative, respectively, have been synthesized and evaluated as cell penetrating peptides (CPP) and as selective vectors for anti-Leishmania drug delivery systems (DDS). They lacked cytotoxicity on the tumoral human cell line HeLa with a moderate cell-uptake on these cells. In contrast, both γ-CC and γ-CT tetradecamers were microbicidal on the protozoan parasite Leishmania beyond 25 μM, with significant intracellular accumulation. They were conjugated to fluorescent doxorubicin (Dox) as a standard drug showing toxicity beyond 1 μM, while free Dox was not toxic. Intracellular accumulation was 2.5 higher than with Dox-TAT conjugate (TAT = transactivator of transcription, taken as a standard CPP). The conformational structure of the conjugates was approached both by circular dichroism spectroscopy and molecular dynamics simulations. Altogether, computational calculations predict that the drug-γ-peptide conjugates adopt conformations that bury the Dox moiety into a cavity of the folded peptide, while the positively charged guanidinium groups face the solvent. The favorable charge/hydrophobicity balance in these CPP improves the solubility of Dox in aqueous media, as well as translocation across cell membranes, making them promising candidates for DDS.

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.