Maytansinol Functionalization: Towards Useful Probes for Studying Microtubule Dynamics.

Maytansinoids are a successful class of natural and semisynthetic tubulin binders, known for their potent cytotoxic activity. Their wider application as cytotoxins and chemical probes to study tubulin dynamics has been held back by the complexity of natural product chemistry. Here we report the synthesis of long-chain derivatives and maytansinoid conjugates. We confirmed that bulky substituents do not impact their high activity or the scaffold's binding mode. These encouraging results open new avenues for the design of new maytansine-based probes.

Maytansinol Functionalization: Towards Useful Probes for Studying Microtubule Dynamics.

Invited for the cover of this issue are the groups of Professors Passarella and Pieraccini at the University of Milan, in collaboration with some of the members of TubInTrain consortium. The image depicts work with the elements of nature, in particular the destabilising effect of maytansinol (the constellation) on microtubules (the trees). Read the full text of the article at 10.1002/chem.202203431.

Resonance Raman and Visible Micro-Spectroscopy for the In-Vivo and In-Vitro Characterization of Anthocyanin-Based Pigments in Blue and Violet Flowers: A Comparison with HPLC-ESI- MS Analysis of the Extracts.

Microanalysis techniques based on resonance Raman and reflection visible spectroscopy have been applied to the characterization of pigments responsible for the blue or violet coloration in flowers; in particular of Lobelia erinus, Campanula portenschlagiana, Cineraria, Viola tricolor, Anemone coronaria, Agapanthus, Platycodon, Salvia farinacea, Plumbago capensis, Ceratostigma plumbaginoides, Commelina communis and Salvia patens. The spectroscopic methods were applied both in vivo on the flower petals and in vitro on extracts obtained through a procedure based on SPE (solid-phase extraction) optimized for minimal quantities of vegetable raw material. Different patterns obtained for the Raman spectra have been correlated, also on the basis of density functional theory (DFT) calculations, with different schemes of substitution of the benzopyrilium nucleus of the anthocyanins and with various possible forms of copigmentation responsible for the stabilization of the blue color. The results obtained were verified by comparison with the analysis of the extracts by HPLC-ESI-MS (liquid chromatography-electrospray ionization-mass spectrometry).

Maytansinol Derivatives: Side Reactions as a Chance for New Tubulin Binders.

Maytansinol is a valuable precursor for the preparation of maytansine derivatives (known as maytansinoids). Inspired by the intriguing structure of the macrocycle and the success in targeted cancer therapy of the derivatives, we explored the maytansinol acylation reaction. As a result, we were able to obtain a series of derivatives with novel modifications of the maytansine scaffold. We characterized these molecules by docking studies, by a comprehensive biochemical evaluation, and by determination of their crystal structures in complex with tubulin. The results shed further light on the intriguing chemical behavior of maytansinoids and confirm the relevance of this peculiar scaffold in the scenario of tubulin binders.

Molecular Similarity Perception Based on Machine-Learning Models.

Molecular similarity is an impressively broad topic with many implications in several areas of chemistry. Its roots lie in the paradigm that 'similar molecules have similar properties'. For this reason, methods for determining molecular similarity find wide application in pharmaceutical companies, e.g., in the context of structure-activity relationships. The similarity evaluation is also used in the field of chemical legislation, specifically in the procedure to judge if a new molecule can obtain the status of orphan drug with the consequent financial benefits. For this procedure, the European Medicines Agency uses experts' judgments. It is clear that the perception of the similarity depends on the observer, so the development of models to reproduce the human perception is useful. In this paper, we built models using both 2D fingerprints and 3D descriptors, i.e., molecular shape and pharmacophore descriptors. The proposed models were also evaluated by constructing a dataset of pairs of molecules which was submitted to a group of experts for the similarity judgment. The proposed machine-learning models can be useful to reduce or assist human efforts in future evaluations. For this reason, the new molecules dataset and an online tool for molecular similarity estimation have been made freely available.

Exploring Orthogonality between Halogen and Hydrogen Bonding Involving Benzene.

The concept of orthogonality between halogen and hydrogen bonding, brought out by Ho and coworkers some years ago, has become a widely accepted idea within the chemists' community. While the original work was based on a common carbonyl oxygen as acceptor for both interactions, we explore here, by means of M06-2X, M11, ωB97X, and ωB97XD/aug-cc-PVTZ DFT calculations, the interdependence of halogen and hydrogen bonding with a shared π-electron system of benzene. The donor groups (specifically NCBr and H2O) were placed on either or the same side of the ring, according to a double T-shaped or a perpendicular geometry, respectively. The results demonstrate that the two interactions with benzene are not strictly independent on each other, therefore outlining that the orthogonality between halogen and hydrogen bonding, intended as energetical independence between the two interactions, should be carefully evaluated according to the specific acceptor group.

On-resin multicomponent 1,3-dipolar cycloaddition of cyclopentanone-proline enamines and sulfonylazides as an efficient tool for the synthesis of amidino depsipeptide mimics.

Depsipeptides are biologically active peptide derivatives that possess a high therapeutic interest. The development of depsipeptide mimics characterized by a chemical diversity could lead to compounds with enhanced features and activity. In this work, an on-resin multicomponent procedure for the synthesis of amidino depsipeptide mimics is described. This approach exploits a metal-free 1,3-dipolar cycloaddition of cyclopentanone-proline enamines and sulfonylazides. In this reaction, the obtained primary cycloadduct undergoes a ring opening and molecular rearrangement giving access to a linear sulfonyl amidine functionalized with both a peptide chain and a diazoalkane. The so-obtained diazo function "one pot" reacts with the carboxylic group of N-Fmoc-protected amino acids leading to amidino depsipeptide mimics possessing a C4 aliphatic chain. An important advantage of this procedure is the possibility to easily obtain amidino-functionalized derivatives that are proteolytically stable peptide bond bioisosteres. Moreover, the conformational freedom given by the alkyl chain could promote the obtainment of cyclic depsipeptide with a stabilized secondary structure as demonstrated with both in silico calculations and experimental conformational studies. Finally, labeled depsipeptide mimics can be also synthesized using a fluorescent sulfonylazide in the multicomponent reaction.

X-ray Constrained Spin-Coupled Wavefunction: a New Tool to Extract Chemical Information from X-ray Diffraction Data.

The X-ray constrained wavefunction (XCW) approach is a reliable and widely used method of quantum crystallography that allows the determination of wavefunctions compatible with X-ray diffraction data. So far, all the existing XCW techniques have been developed in the framework of molecular orbital theory and, consequently, provide only pictures of the "experimental" electronic structures that are far from the traditional chemical perception. Here a new strategy is proposed that, by combining the XCW philosophy with the spin-coupled method of valence bond theory, enables direct extraction of traditional chemical information (e.g., weights of resonance structures) from X-ray diffraction measurements. Preliminary results have shown that the new technique is really able to efficiently capture the effects of the crystal environment on the electronic structure, and can be considered as a new useful tool to perform chemically sound analyses of the X-ray diffraction data.

On the molecular optical nonlinearity of halogen-bond-forming azobenzenes.

We study hyper-Rayleigh scattering and computed molecular hyperpolarizability in a series of azobenzene chromophores in chloroform and dimethylformamide as solvents. The chromophores form halogen or hydrogen bonds of varying strength with dimethylformamide molecules, differently from what is expected for chloroform. We show that hyperpolarizability is unaffected or sligthly lower with the azobenzene forming the strongest halogen bond. Solid supramolecular polymers with the same chromophores have previously demonstrated clearly higher second-order nonlinear responses when a halogen-bond-accepting polymer is used, the larger increase being associated with the stronger halogen bond. The present study proves that the higher optical nonlinearity in polymers lies in the better ordering of the chromophores instead of changes in molecular hyperpolarizability, highlighting the unique properties of halogen bonding in supramolecular chemistry.

Assessment of DFT Functionals for QTAIM Topological Analysis of Halogen Bonds with Benzene.

Halogen bonding, a noncovalent interaction between a halogen atom and a nucleophilic site, is receiving a growing attention in the chemical community stimulating a large number of theoretical investigations. The density functional theory (DFT) approach revealed to be one of the most suitable methods owing to its accuracy and low computational cost. We report here a detailed analysis of the performance of an extensive set of DFT functionals in reproducing accurate binding energies and topological properties for the halogen-bonding interaction of either NCX or PhX molecules (X = F, Cl, Br, I) with the aromatic system of benzene in the T-shaped configuration. It was found that the better performance for both sets of properties is provided by a small subset of functionals able to take into account, implicitly or explicitly (by inclusion of an additive pairwise potential), the dispersion contribution, that is, ωB97X, M06-2X, M11, mPW2PLYP-D, and B2PLYP-D3.

Halogen bonds with benzene: an assessment of DFT functionals.

The performance of an extensive set of density functional theory functionals has been tested against CCSD(T) and MP2 results, extrapolated to the complete basis set (CBS) limit, for the interaction of either DCl or DBr (D = H, HCC, F, and NC) with the aromatic system of benzene. It was found that double hybrid functionals explicitly including dispersion, that is, B2PLYPD and mPW2PLYPD, provide the better agreement with the CCSD(T)/CBS results on both energies and equilibrium geometry, indicating the importance of dispersive contributions in determining this interaction. Among the less expensive functionals, the better performance is provided by the ωB97X and M062X functionals, while the ωB97XD and B97D functionals are shown to work very well for bromine complexes but not so well for chlorine complexes.

Quinazolinecarboline alkaloid evodiamine as scaffold for targeting topoisomerase I and sirtuins.

This paper reports the synthesis of a series of evodiamine derivatives. We assayed the ability to inhibit cell growth on three human tumour cell lines (H460, MCF-7 and HepG2) and we evaluated the capacity to interfere with the catalytic activity of topoisomerase I both by the relaxation assay and the occurrence of the cleavable complex. Moreover, whose effect on sirtuins 1, 2 and 3 was investigated. Finally, molecular docking analyses were performed in an attempt to rationalize the biological results.

Halogen-bonding interactions with π systems: CCSD(T), MP2, and DFT calculations.

Halogen bonding is a noncovalent interaction between a halogen atom and a nucleophilic site. Interactions involving the π electrons of aromatic rings have received, up to now, little attention, despite the large number of systems in which they are present. We report binding energies of the interaction between either NCX or PhX (X = F, Cl, Br, I) and the aromatic benzene system as determined with the coupled cluster with perturbative triple excitations method [CCSD(T)] extrapolated at the complete basis set limit. Results are compared with those obtained by Møller-Plesset perturbation theory to second order (MP2) and density functional theory (DFT) calculations by using some of the most common functionals. Results show the important role of DFT in studying this interaction.

Metadynamics study of a ß-hairpin stability in mixed solvents.

Understanding the molecular mechanisms that allow some organisms to survive in extremely harsh conditions is an important achievement that might disclose a wide range of applications and that is constantly drawing the attention of many research fields. The high adaptability of these living creatures is related to the presence in their tissues of a high concentration of osmoprotectants, small organic, highly soluble molecules. Despite osmoprotectants having been known for a long time, a full disclosure of the machinery behind their activity is still lacking. Here we describe a computational approach that, taking advantage of the recently developed metadynamics technique, allows one to fully describe the free energy surface of a small ß-hairpin peptide and how it is affected by an osmoprotectant, glycine betaine (GB) and for comparison by urea, a common denaturant. Simulations led to relevant thermodynamic information, including how the free energy difference of denaturation is affected by the two cosolvents; unlike urea, GB caused a considerable increase of the folded basin stability, which transposes into a higher melting temperature. NMR experiments confirmed the picture derived from the theoretical study. Further molecular dynamics simulations of selected conformations allowed investigation into deeper detail the role of GB in folded state protection. Simulations of the protein in GB solutions clearly showed an excess of osmoprotectant in the solvent bulk, rather than in the protein domain, confirming the exclusion from the protein surface, but also highlighted interesting features on its interactions, opening to new scenarios besides the classic "indirect mechanism" hypothesis.

Halogen bonding in ligand-receptor systems in the framework of classical force fields.

Halogen bond is an important non-covalent interaction which is receiving a growing attention in the study of protein-ligand complexes. Many drugs are halogenated molecules and it has been recently shown that many halogenated ligands establish halogen bonds with biomolecules. As the halogen bond nature is due to an anisotropy of the electrostatic potential around halogen atoms, it is not possible to use traditional force fields based on a set of atom-centred charges to study halogen bonds in biomolecules. We show that the introduction of pseudo-atoms on halogens permits us to correctly describe the anisotropy of the electrostatic potential and to perform molecular dynamics simulations on complexes of proteins with halogenated ligands that reproduce experimental values. The results are compared with crystallographic data and with hybrid quantum mechanics/molecular mechanics calculations.

Vinblastine perturbation of tubulin protofilament structure: a computational insight.

Tubulin is a heterodimeric protein whose self assembly leads to the formation of protofilaments and of more complex structures called microtubules, key components of the cytoskeleton which have a fundamental role in the cell division process. Due to its biological function, tubulin is the target of many antitumoral molecules that exert their action on proliferating tumoral cells. Among these drugs, vinblastine has been widely used in therapy for a long time, albeit its mechanism of interaction with tubulin has remained elusive until recently. Vinblastine acts as a microtubule destabilizing agent and induces the formation of curved or ring-shaped tubulin polymers instead of linear protofilaments in vitro. In this paper we compare, using molecular dynamics simulations and free energy calculations, the network of interactions that allow the assembly of model linear protofilaments with those present in curved tubulin polymers complexed with vinblastine. It is shown that vinblastine, wedging between tubulin heterodimers, actually mediates part of the interactions between them and acts by crosslinking the two proteins, leading to the observed curved polymers rather than to their disassembly.

Modelling of short synthetic antifreeze peptides: Insights into ice-pinning mechanism.

Organisms living in icy environments produce antifreeze proteins to control ice growth and recrystallization. It has been proposed that these molecules pin the surface of ice crystals, thus inducing the formation of a curved surface that arrests crystal growth. Such proteins are very appealing for many potential applications in food industry, material science and cryoconservation of organs and tissues. Unfortunately, their structural complexity has seriously hampered their practical use, while efficient and accessible synthetic analogues are highly desirable. In this paper, we used molecular dynamics based techniques to model the interaction of three short antifreeze synthetic peptides with an ice surface. The employed protocols succeeded in reproducing the ice pinning action of antifreeze peptides and the consequent ice growth arrest, as well as in distinguishing between antifreeze and control peptides, for which no such effect was observed. Principal components analysis of peptides trajectories in different simulation settings permitted to highlight the main structural features associated to antifreeze activity. Modeling results are highly correlated with experimentally measured properties, and insights on ice-peptide interactions and on conformational patterns favoring antifreeze activity will prompt the design of new and improved antifreeze peptides.

The Origin of the σ-Hole in Halogen Atoms: a Valence Bond Perspective.

A detailed Valence Bond-Spin Coupled analysis of a series of halogenated molecules is here reported, allowing to get a rigorous ab initio demonstration of the qualitative models previously proposed to explain the origin of halogen bonding. The concepts of σ-hole and negative belt observed around the halogen atoms in the electrostatic potential maps are here interpreted by analysis of the relevant Spin Coupled orbitals.

Simulating Multiple Substrate-Binding Events by γ-Glutamyltransferase Using Accelerated Molecular Dynamics.

γ-Glutamyltransferase (GGT) is an enzyme that uses γ-glutamyl compounds as substrates and catalyzes their transfer to a water molecule or an acceptor substrate with varied physiological function in bacteria, plants, and animals. Crystal structures of GGT are known for different species and in different states of the chemical reaction; however, the structural dynamics of the substrate binding to the catalytic site of GGT are unknown. Here, we modeled Escherichia coli GGT's glutamine binding by using a swarm of accelerated molecular dynamics (aMD) simulations. Characterization of multiple binding events identified three structural binding motifs composed of polar residues in the binding pocket that govern glutamine binding into the active site. Simulated open and closed conformations of a lid-loop protecting the binding cavity suggest its role as a gating element by allowing or blocking substrates entry into the binding pocket. Partially open states of the lid-loop are accessible within thermal fluctuations, while the estimated free energy cost of a complete open state is 2.4 kcal/mol. Our results suggest that both specific electrostatic interactions and GGT conformational dynamics dictate the molecular recognition of substrate-GGT complexes.

X-ray constrained spin-coupled technique: theoretical details and further assessment of the method.

One of the well-established methods of modern quantum crystallography is undoubtedly the X-ray constrained wavefunction (XCW) approach, a technique that enables the determination of wavefunctions which not only minimize the energy of the system under examination, but also reproduce experimental X-ray diffraction data within the limit of the experimental errors. Initially proposed in the framework of the Hartree-Fock method, the strategy has been gradually extended to other techniques of quantum chemistry, but always remaining limited to a single-determinant ansatz for the wavefunction to extract. This limitation has been recently overcome through the development of the novel X-ray constrained spin-coupled (XCSC) approach [Genoni et al. (2018). Chem. Eur. J. 24, 15507-15511] which merges the XCW philosophy with the traditional spin-coupled strategy of valence bond theory. The main advantage of this new technique is the possibility of extracting traditional chemical descriptors (e.g. resonance structure weights) compatible with the experimental diffraction measurements, without the need to introduce information a priori or perform analyses a posteriori. This paper provides a detailed theoretical derivation of the fundamental equations at the basis of the XCSC method and also introduces a further advancement of its original version, mainly consisting in the use of molecular orbitals resulting from XCW calculations at the Hartree-Fock level to describe the inactive electrons in the XCSC computations. Furthermore, extensive test calculations, which have been performed by exploiting high-resolution X-ray diffraction data for salicylic acid and by adopting different basis sets, are presented and discussed. The computational tests have shown that the new technique does not suffer from particular convergence problems. Moreover, all the XCSC calculations provided resonance structure weights, spin-coupled orbitals and global electron densities slightly different from those resulting from the corresponding unconstrained computations. These discrepancies can be ascribed to the capability of the novel strategy to capture the information intrinsically contained in the experimental data used as external constraints.