Thermodynamically–Weighted Conformational Ensemble of Cyclic RGD Peptidomimetics from NOE Data

Molecular dynamics simulations of cRGD-conjugated PEGylated TiO2 nanoparticles for targeted photodynamic therapy

The conjugation of high-affinity cRGD-containing peptides is a promising approach in nanomedicine to efficiently reduce off-targeting effects and enhance the cellular uptake by integrin-overexpressing tumor cells. Herein we utilize atomistic molecular dynamics simulations to evaluate key structural-functional parameters of these targeting ligands for an effective binding activity towards α_Vβ₃ integrins. An increasing number of cRGD ligands is conjugated to PEG chains grafted to highly curved TiO₂ nanoparticles to unveil the impact of cRGD density on the ligand's presentation, stability, and conformation in an explicit aqueous environment. We find that a low density leads to an optimal spatial presentation of cRGD ligands out of the "stealth" PEGylated layer around the nanosystem, favoring a straight upward orientation and spaced distribution of the targeting ligands in the bulk-water phase. On the contrary, high densities favor over-clustering of cRGD ligands, driven by a concerted mechanism of enhanced ligand-ligand interactions and reduced water accessibility over the ligand's molecular surface. These findings strongly suggest that the ligand density modulation is a key factor in the design of cRGD-targeting nanodevices to maximize their binding efficiency into over-expressed α_Vβ₃ integrin receptors.

Acid-base and lipophilic properties of peptide nucleic acid derivatives

The first combined experimental and theoretical study on the ionization and lipophilic properties of peptide nucleic acid (PNA) derivatives, including eleven PNA monomers and two PNA decamers, is described. The acidity constants (pK_a) of individual acidic and basic centers of PNA monomers were measured by automated potentiometric pH titrations in water/methanol solution, and these values were found to be in agreement with those obtained by MoKa software. These results indicate that single nucleobases do not change their pK_a values when included in PNA monomers and oligomers. In addition, immobilized artificial membrane chromatography was employed to evaluate the lipophilic properties of PNA monomers and oligomers, which showed the PNA derivatives had poor affinity towards membrane phospholipids, and confirmed their scarce cell penetrating ability. Overall, our study not only is of potential relevance to evaluate the pharmacokinetic properties of PNA, but also constitutes a reliable basis to properly modify PNA to obtain mimics with enhanced cell penetration properties.

•

The first study on acid-base and lipophilic properties of peptide nucleic acids (PNA).

•

pK_a of acid-base centers of PNA evaluated by potentiometric method and MoKa prediction.

•

NMR experiments provide additional information on the protonation of PNA monomers.

•

Lipophilicity of PNA monomers and oligomers is investigated by IAM chromatography.

•

This study can lay the basis of evaluating the pharmacokinetic properties of PNA.

Managing Experimental 3D Structures in the Beyond-Rule-of-5 Chemical Space: The Case of Rifampicin

The beyond-Rule-of-5 (bRo5) chemical space is a source of new oral drugs and includes large and flexible compounds. Because of their size and conformational variability, bRo5 molecules assume different privileged conformations in the compartments of human body, i. e., they can exhibit chameleonic properties. The elucidation of the ensemble of 3D structures explored by such molecules under different conditions is therefore critical to check the role played by chameleonicity to modulate cell permeability. Here we characterized the conformational ensembles of rifampicin, a bRo5 drug, in polar and nonpolar solvents and in the solid state. We performed NMR experiments, analyzed their results with a novel algorithm and set-up a pool of ad hoc in silico strategies to investigate crystallographic structures retrieved from the CSD. Moreover, a polarity descriptor often related to permeability (SA-3D-PSA) was calculated for all the conformers and its variation with the environment analyzed. Results showed that the conformational behavior of rifampicin in solution and in the solid state is not superposable. The identification of dynamic intramolecular hydrogen bonds can be assessed by NMR spectroscopy but not by X-ray structures. Moreover, SA-3D-PSA revealed that dynamic IMHBs do not provide rifampicin with chameleonic properties. Overall, this study highlights that the peculiarity of rifampicin, which is cell permeable probably because of the presence of static IMHBs but is devoid of any chameleonic behavior, can be assessed by a proper analysis of experimental 3D structures.

Identification of Compounds Targeting HuD. Another Brick in the Wall of Neurodegenerative Disease Treatment

ELAV-like (ELAVL) RNA-binding proteins play a pivotal role in post-transcriptional processes, and their dysregulation is involved in several pathologies. This work was focused on HuD (ELAVL4), which is specifically expressed in nervous tissues, and involved in differentiation and synaptic plasticity mechanisms. HuD represents a new, albeit unexplored, candidate target for the treatment of several relevant neurodegenerative diseases. The aim of this pioneering work was the identification of new molecules able to recognize and bind HuD, thus interfering with its activity. We combined virtual screening, molecular dynamics (MD), and STD-NMR techniques. Starting from around 51 000 compounds, four promising hits eventually provided experimental evidence of their ability to bind HuD. Among the selected best hits, folic acid was found to be the most interesting one, being able to well recognize the HuD binding site. Our results provide a basis for the identification of new HuD interfering compounds which may be useful against neurodegenerative syndromes.

Cyclic RGD and isoDGR Integrin Ligands Containing cis-2-amino-1-cyclopentanecarboxylic (cis-β-ACPC) Scaffolds

Integrin ligands containing the tripeptide sequences Arg-Gly-Asp (RGD) and iso-Asp-Gly- Arg (isoDGR) were actively investigated as inhibitors of tumor angiogenesis and directing unit in tumor-targeting drug conjugates. Reported herein is the synthesis, of two RGD and one isoDGR cyclic peptidomimetics containing (1S,2R) and (1R,2S) cis-2-amino-1-cyclopentanecarboxylic acid (cis-β-ACPC), using a mixed solid phase/solution phase synthetic protocol. The three ligands were examined in vitro in competitive binding assays to the purified α_vβ₃ and α₅β₁ receptors using biotinylated vitronectin (α_vβ₃) and fibronectin (α₅β₁) as natural displaced ligands. The IC50 values of the ligands ranged from nanomolar (the two RGD ligands) to micromolar (the isoDGR ligand) with a pronounced selectivity for α_vβ₃ over α₅β₁. In vitro cell adhesion assays were also performed using the human skin melanoma cell line WM115 (rich in integrin α_vβ₃). The two RGD ligands showed IC₅₀ values in the same micromolar range as the reference compound (cyclo[RGDfV]), while for the isoDGR derivative an IC₅₀ value could not be measured for the cell adhesion assay. A conformational analysis of the free RGD and isoDGR ligands by NMR (VT-NMR and NOESY experiments) and computational studies (MC/EM and MD), followed by docking simulations performed in the α_Vβ₃ integrin active site, provided a rationale for the behavior of these ligands toward the receptor.

Hydrogen Bonding in a l-Glutamine-Based Polyamidoamino Acid and its pH-Dependent Self-Ordered Coil Conformation

This paper reports on synthesis, acid-base properties, and self-structuring in water of a chiral polyamidoamino acid, M-l-Gln, obtained from the polyaddition of N,N'-methylenebisacrylamide with l-glutamine, with the potential of establishing hydrogen bonds through its prim-amide pendants. The M-l-Gln showed pH-responsive circular dichroism spectra, revealing ordered conformations. Structuring was nearly insensitive to ionic strength but sensitive to denaturing agents. The NMR diffusion studies were consistent with a population of unimolecular nanoparticles thus excluding aggregation. The M-l-Gln had the highest molecular weight and hydrodynamic radius among all polyamidoamino acids described. Possibly, transient hydrogen bonds between l-glutamine molecules and M-l-Gln growing chains facilitated the polyaddition reaction. Theoretical modeling showed that M-l-Gln assumed pH-dependent self-ordered coil conformations with main chain transoid arrangements reminiscent of the protein hairpin motif owing to intramolecular dipole moments and hydrogen bonds. The latter were most numerous at the isoelectric point (pH 4.5), where they mainly involved even topologically distant main chain amide N-H and side chain amide C=O brought to proximity by structuring. Hydrogen bonds at pH 4.5 were also suggested by variable temperature NMR. The 2D NOESY experiments at pH 4.5 confirmed the formation of compact structures through the analysis of the main chain/side chain hydrogen contacts, in line with MD simulations.

Insight into GEBR-32a: Chiral Resolution, Absolute Configuration and Enantiopreference in PDE4D Inhibition

Alzheimer's disease is the most common type of dementia, affecting millions of people worldwide. One of its main consequences is memory loss, which is related to downstream effectors of cyclic adenosine monophosphate (cAMP). A well-established strategy to avoid cAMP degradation is the inhibition of phosphodiesterase (PDE). In recent years, GEBR-32a has been shown to possess selective inhibitory properties against PDE type 4 family members, resulting in an improvement in spatial memory processes without the typical side effects that are usually correlated with this mechanism of action. In this work, we performed the HPLC chiral resolution and absolute configuration assignment of GEBR-32a. We developed an efficient analytical and semipreparative chromatographic method exploiting an amylose-based stationary phase, we studied the chiroptical properties of both enantiomers and we assigned their absolute configuration by 1H-NMR (nuclear magnetic resonance). Lastly, we measured the IC50 values of both enantiomers against both the PDE4D catalytic domain and the long PDE4D3 isoform. Results strongly support the notion that GEBR-32a inhibits the PDE4D enzyme by interacting with both the catalytic pocket and the regulatory domains.

Determination of Structural Ensembles of Flexible Molecules in Solution from NMR Data Undergoing Spin Diffusion

Spin diffusion is a formidable problem when interpreting NMR data of chemical compounds. We developed a method to reconstruct the conformational ensemble of flexible molecules displaying spin diffusion, which minimizes the subjective bias in the interpretation of experimental data and which can be used routinely to obtain sets of structures with the correct thermodynamic weights. We showed in the case of a flexible molecule that the correct conformational ensemble is quite different from that obtained with standard methods.

Exploring E-cadherin-peptidomimetics interaction using NMR and computational studies

Cadherins are homophilic cell-cell adhesion molecules whose aberrant expression has often been shown to correlate with different stages of tumor progression. In this work, we investigate the interaction of two peptidomimetic ligands with the extracellular portion of human E-cadherin using a combination of NMR and computational techniques. Both ligands have been previously developed as mimics of the tetrapeptide sequence Asp1-Trp2-Val3-Ile4 of the cadherin adhesion arm, and have been shown to inhibit E-cadherin-mediated adhesion in epithelial ovarian cancer cells with millimolar potency. To sample a set of possible interactions of these ligands with the E-cadherin extracellular portion, STD-NMR experiments in the presence of two slightly different constructs, the wild type E-cadherin-EC1-EC2 fragment and the truncated E-cadherin-(Val3)-EC1-EC2 fragment, were carried out at three temperatures. Depending on the protein construct, a different binding epitope of the ligand and also a different temperature effect on STD signals were observed, both suggesting an involvement of the Asp1-Trp2 protein sequence among all the possible binding events. To interpret the experimental results at the atomic level and to probe the role of the cadherin adhesion arm in the dynamic interaction with the peptidomimetic ligand, a computational protocol based on docking calculations and molecular dynamics simulations was applied. In agreement with NMR data, the simulations at different temperatures unveil high variability/dynamism in ligand-cadherin binding, thus explaining the differences in ligand binding epitopes. In particular, the modulation of the signals seems to be dependent on the protein flexibility, especially at the level of the adhesive arm, which appears to participate in the interaction with the ligand. Overall, these results will help the design of novel cadherin inhibitors that might prevent the swap dimer formation by targeting both the Trp2 binding pocket and the adhesive arm residues.

Classical cadherins are the main adhesive proteins at the intercellular junctions and play an essential role in tissue morphogenesis and homeostasis. A large number of studies have shown that cadherin aberrant expression and/or dysregulation often correlate with pathological processes, such as tumor development and progression. Notwithstanding the emerging role played by cadherins in a number of solid tumors, the rational design of small inhibitors targeting these proteins is still in its infancy, likely due to the challenges posed by the development of small drug-like molecules that modulate protein-protein interactions and to the structural complexity of the various cadherin dimerization interfaces that constantly form and disappear as the protein moves along its highly dynamic and reversible homo-dimerization trajectory. In this work, we study the interaction of two small molecules with the extracellular portion of human E-cadherin using a combination of spectroscopic and computational techniques. The availability of molecules interfering in the cadherin homophilic interactions could provide a useful tool for the investigation of cadherin function in tumors, and potentially pave the way to the development of novel alternative diagnostic and therapeutic interventions in cadherin-expressing solid tumors.

A Combined NMR-Computational Study of the Interaction between Influenza Virus Hemagglutinin and Sialic Derivatives from Human and Avian Receptors on the Surface of Transfected Cells

The development of small-molecule inhibitors of influenza virus Hemagglutinin could be relevant to the opposition of the diffusion of new pandemic viruses. In this work, we made use of Nuclear Magnetic Resonance (NMR) spectroscopy to study the interaction between two derivatives of sialic acid, Neu5Ac-α-(2,6)-Gal-β-(1⁻4)-GlcNAc and Neu5Ac-α-(2,3)-Gal-β-(1⁻4)-GlcNAc, and hemagglutinin directly expressed on the surface of recombinant human cells. We analyzed the interaction of these trisaccharides with 293T cells transfected with the H5 and H1 variants of hemagglutinin, which thus retain their native trimeric conformation in such a realistic environment. By exploiting the magnetization transfer between the protein and the ligand, we obtained evidence of the binding event, and identified the epitope. We analyzed the conformational features of the glycans with an approach combining NMR spectroscopy and data-driven molecular dynamics simulations, thus obtaining useful information for an efficient drug design.

Investigating the Interaction of Cyclic RGD Peptidomimetics with αVβ₆ Integrin by Biochemical and Molecular Docking Studies

The interaction of a small library of cyclic RGD (Arg-Gly-Asp) peptidomimetics with α_Vβ₆ integrin has been investigated by means of competitive solid phase binding assays to the isolated receptor and docking calculations in the crystal structure of the α_Vβ₆ binding site. To this aim, a rigid receptor-flexible ligand docking protocol has been set up and then applied to predict the binding mode of the cyclic RGD peptidomimetics to α_Vβ₆ integrin. Although the RGD interaction with α_Vβ₆ recapitulates the RGD binding mode observed in α_Vβ₃, differences between the integrin binding pockets can strongly affect the ligand binding ability. In general, the peptidomimetics exhibited IC₅₀ values for integrin α_Vβ₆ (i.e., the concentration of compound required for 50% inhibition of biotinylated fibronectin binding to isolated α_Vβ₆ integrin) in the nanomolar range (77–345 nM), about 10–100 times higher than those for the related α_Vβ₃ receptor, with a single notable ligand displaying a low nanomolar IC₅₀ value (2.3 nM). Insights from the properties of the binding pocket combined with the analysis of the docking poses provided a rationale for ligand recognition and selectivity.

Principles of protein structural ensemble determination

The biological functions of protein molecules are intimately dependent on their conformational dynamics. This aspect is particularly evident for disordered proteins, which constitute perhaps one-third of the human proteome. Therefore, structural ensembles often offer more useful representations of proteins than individual conformations. Here, we describe how the well-established principles of protein structure determination should be extended to the case of protein structural ensembles determination. These principles concern primarily how to deal with conformationally heterogeneous states, and with experimental measurements that are averaged over such states and affected by a variety of errors. We first review the growing literature of recent methods that combine experimental and computational information to model structural ensembles, highlighting their similarities and differences. We then address some conceptual problems in the determination of structural ensembles and define future goals towards the establishment of objective criteria for the comparison, validation, visualization and dissemination of such ensembles.