Dynamic self-assembly of supramolecular catalysts from precision macromolecules.

We show the emergence of strong catalytic activity at low concentrations in dynamic libraries of complementary sequence-defined oligomeric chains comprising pendant functional catalytic groups and terminal recognition units. In solution, the dynamic constitutional library created from pairs of such complementary oligomers comprises free oligomers, self-assembled di(oligomeric) macrocycles, and a virtually infinite collection of linear poly(oligomeric) chains. We demonstrate, on an exemplary catalytic system requiring the cooperation of no less than five chemical groups, that supramolecular di(oligomeric) macrocycles exhibit a catalytic turnover frequency ca. 20 times larger than the whole collection of linear poly(oligomers) and free chains. Molecular dynamics simulations and network analysis indicate that self-assembled supramolecular di(oligomeric) macrocycles are stabilized by different interactions, among which chain end pairing. We mathematically model the catalytic properties of such complex dynamic libraries with a small set of physically relevant parameters, which provides guidelines for the synthesis of oligomers capable to self-assemble into functionally-active supramolecular macrocycles over a larger range of concentrations.

Sequence-defined antibody-recruiting macromolecules.

Antibody-recruiting molecules represent a novel class of therapeutic agents that mediate the recruitment of endogenous antibodies to target cells, leading to their elimination by the immune system. Compared to single-ligand copies, macromolecular scaffolds presenting multiple copies of an antibody-binding ligand offer advantages in terms of increased complex avidity. In this study, we describe the synthesis of sequence-defined macromolecules designed for antibody recruitment, utilising dinitrophenol (DNP) as a model antibody-recruiting motif. The use of discrete macromolecules gives access to varying the spacing between DNP motifs while maintaining the same chain length. This characteristic enables the investigation of structure-dependent binding interactions with anti-DNP antibodies. Through solid-phase thiolactone chemistry, we synthesised a series of oligomers with precisely localised DNP motifs along the backbone and a terminal biotin motif for surface immobilisation. Utilising biolayer interferometry analysis, we observed that oligomers with adjacent DNP motifs exhibited enhanced avidity for anti-DNP antibodies. Molecular modelling provided insights into the structures and dynamics of the various macromolecules, shedding light on the accessibility of the ligands to the antibodies. Overall, our findings highlight that the use of sequence-defined macromolecules can contribute to our understanding of structure-activity relationships and provide insights for the design of novel antibody-recruiting therapeutic agents.

Towards a Rational Design of Antibody-Recruiting Molecules through a Computational Microscopy View of their Interactions with the Target Antibody.

Antibody recruiting molecules (ARMs) are an innovative class of chimeric molecules, consisting of an antibody-binding ligand (ABL) and a target-binding ligand (TBL). ARMs mediate ternary complex formation between a target cell of interest for elimination and endogenous antibodies that are present in human serum. Clustering of fragment crystallizable (Fc) domains on the surface of antibody-bound cells mediate destruction of the target cell by innate immune effector mechanisms. ARMs are typically designed by conjugating small molecule haptens to a (macro)molecular scaffold, without considering the structure of the respective anti-hapten antibody. Here we report on a computational molecular modeling methodology that allows for studying the close contacts between ARMs and the anti-hapten antibody, considering (1) the spacer length between ABL and TBL; (2) the number of ABL and TBL, and (3) the molecular scaffold onto which these are positioned. Our model predicts the difference in binding modes of the ternary complex and predicts which ARMs are optimal recruiters. Avidity measurements of the ARM-antibody complex and ARM-mediated antibody recruitment to cell surfaces in vitro confirmed these computational modeling predictions. This kind of multiscale molecular modelling holds potential for design of drug molecules that rely on antibody binding for their mechanism of action.

Uncovering the antimalarial potential of toad venoms through a bioassay-guided fractionation process.

Malaria remains to date one of the most devastating parasitic diseases worldwide. The fight against this disease is rendered more difficult by the emergence and spread of drug-resistant strains. The need for new therapeutic candidates is now greater than ever. In this study, we investigated the antiplasmodial potential of toad venoms. The wide array of bioactive compounds present in Bufonidae venoms has allowed researchers to consider many potential therapeutic applications, especially for cancers and infectious diseases. We focused on small molecules, namely bufadienolides, found in the venom of Rhinella marina (L.). The developed bio-guided fractionation process includes a four solvent-system extraction followed by fractionation using flash chromatography. Sub-fractions were obtained through preparative TLC. All samples were characterized using chromatographic and spectrometric techniques and then underwent testing on in vitro Plasmodium falciparum cultures. Two strains were considered: 3D7 (chloroquine-sensitive) and W2 (chloroquine-resistant). This strategy highlighted a promising activity for one compound named resibufogenin. With IC50 values of (29 ± 8) µg/mL and (23 ± 1) µg/mL for 3D7 and W2 respectively, this makes it an interesting candidate for further investigation. A molecular modelling approach proposed a potential binding mode of resibufogenin to Plasmodium falciparum adenine-triphosphate 4 pump as antimalarial drug target.

Gaining Insights on the Interactions of a Class of Decorated (2-([2,2'-Bipyridin]-6-yl)phenyl)platinum Compounds with c-Myc Oncogene Promoter G-Quadruplex and Other DNA Structures.

Organometallic molecules offer some of the most promising scaffolds for interaction with G-quadruplex nucleic acids. We report the efficient synthesis of a family of organoplatinum(II) complexes, featuring a 2-([2,2'-bipyridin]-6-yl)phenyl tridentate (N∧ N∧ C) ligand, that incorporates peripheral side-chains aiming at enhancing and diversifying its interaction capabilities. These include a di-isopropyl carbamoyl amide, a morpholine ethylenamide, two enantiomeric proline imides and an oxazole. The binding affinities of the Pt-complexes were evaluated via UV-vis and fluorescence titrations, against 5 topologically-distinct DNA structures, including c-myc G-quadruplex, two telomeric (22AG) G-quadruplexes, a duplex (ds26) and a single-stranded (polyT) DNA. All compounds exhibited binding selectivity in favour of c-myc, with association constants (Ka ) in the range of 2-5×105  M-1 , lower affinity for both folds of 22AG and for ds26 and negligible affinity for polyT. Remarkable emission enhancements (up to 200-fold) upon addition of excess DNA were demonstrated by a subset of the compounds with c-myc, providing a basis for optical selectivity, since optical response to all other tested DNAs was low. A c-myc DNA-melting experiment showed significant stabilizing abilities for all compounds, with the most potent binder, the morpholine-Pt-complex, exhibiting a ΔTm >30 °C, at 1 : 5 DNA-to-ligand molar ratio. The same study implied contributions of the diverse side-chains to helix stabilization. To gain direct evidence of the nature of the interactions, mixtures of c-myc with the four most promising compounds were studied via UV Resonance Raman (UVRR) spectroscopy, which revealed end-stacking binding mode, combined with interactions of side-chains with loop nucleobase residues. Docking simulations were conducted to provide insights into the binding modes for the same four Pt-compounds, suggesting that the binding preference for two alternative orientations of the c-myc G-quadruplex thymine 'cap' ('open' vs. 'closed'), as well as the relative contributions to affinity from end-stacking and H-bonding, are highly dependent on the nature of the interacting Pt-complex side-chain.

Revealing the Organization of Catalytic Sequence-Defined Oligomers via Combined Molecular Dynamics Simulations and Network Analysis.

Similar to biological macromolecules such as DNA and proteins, the precise control over the monomer position in sequence-defined polymers is of paramount importance for tuning their structures and properties toward achieving specific functions. Here, we apply molecular network analysis on three-dimensional structures issued from molecular dynamics simulations to decipher how the chain organization of trifunctional catalytic oligomers is influenced by the oligomer sequence and the length of oligo(ethylene oxide) spacers. Our findings demonstrate that the tuning of their primary structures is crucial for favoring cooperative interactions between the catalytic units and thus higher catalytic activities. This combined approach can assist in establishing structure-property relationships, leading to a more rational design of sequence-defined catalytic oligomers via computational chemistry.

Hierarchical Self-Assembly and Multidynamic Responsiveness of Fluorescent Dynamic Covalent Networks Forming Organogels.

Smart stimuli-responsive fluorescent materials are of interest in the context of sensing and imaging applications. In this project, we elaborated multidynamic fluorescent materials made of a tetraphenylethene fluorophore displaying aggregation-induced emission and short cysteine-rich C-hydrazide peptides. Specifically, we show that a hierarchical dynamic covalent self-assembly process, combining disulfide and acyl-hydrazone bond formation operating simultaneously in a one-pot reaction, yields cage compounds at low concentration (2 mM), while soluble fluorescent dynamic covalent networks and even chemically cross-linked fluorescent organogels are formed at higher concentrations. The number of cysteine residues in the peptide sequence impacts directly the mechanical properties of the resulting organogels, Young's moduli varying 2500-fold across the series. These materials underpinned by a nanofibrillar network display multidynamic responsiveness following concentration changes, chemical triggers, as well as light irradiation, all of which enable their controlled degradation with concomitant changes in spectroscopic outputs─self-assembly enhances fluorescence emission by ca. 100-fold and disassembly quenches fluorescence emission.

Teaching photosensitizers a new trick: red light-triggered G-quadruplex alkylation by ligand co-localization.

We propose a bimolecular approach for G-quadruplex alkylation, using a pro-reactive furan-containing ligand, activated by red-light irradiation of a proximate G4-binding photosensitizer. G4- over dsDNA alkylation can be achieved selectively and proves high-yielding at low ligand excess. HPLC and modelling studies allowed identifying potential residues involved in the alkylation.

Binding Mode Multiplicity and Multiscale Chirality in the Supramolecular Assembly of DNA and a π-Conjugated Polymer.

Water-soluble π-conjugated polymers are increasingly considered for DNA biosensing. However, the conformational rearrangement, supramolecular organization and dynamics upon interaction with DNA have been overlooked, which prevents the rational design of such detection tools. To elucidate the binding of a cationic polythiophene (CPT) to DNA with atomistic resolution, we performed molecular simulations of their supramolecular assembly. Comparison of replicated simulations show a multiplicity of CPT binding geometries that contribute to the wrapping of CPT around DNA. The different binding geometries are stabilized by both electrostatic interactions between CPT lateral cations and DNA phosphodiesters and van der Waals interactions between the CPT backbone and the DNA grooves. Simulated circular dichroism (CD) spectra show that the induced CD signal stems from a conserved geometrical feature across the replicated simulations, i. e. the presence of segments of syn configurations between thiophene units along the CPT chain. At the macromolecular scale, we inspected the different shapes related to the CPT binding modes around the DNA through symmetry metrics. Altogether, molecular dynamics (MD) simulations, model Hamiltonian calculations of the CD spectra, and symmetry indices provide insights into the origin of induced chirality from the atomic to the macromolecular scale. Our multidisciplinary approach points out the hierarchical aspect of CPT chiral organization induced by DNA.

A Cationic Tetraphenylethene as a Light-Up Supramolecular Probe for DNA G-Quadruplexes.

Guanine-quadruplexes (G4s) are targets for anticancer therapeutics. In this context, human telomeric DNA (HT-DNA) that can fold into G4s sequences are of particular interest, and their stabilization with small molecules through a visualizable process has become a challenge. As a new type of ligand for HT-G4, we designed a tetraimidazolium tetraphenylethene (TPE-Im) as a water-soluble light-up G4 probe. We study its G4-binding properties with HT-DNA by UV-Visible absorption, circular dichroism and fluorescence spectroscopies, which provide insights into the interactions between TPE-Im and G4-DNA. Remarkably, TPE-Im shows a strong fluorescence enhancement and large shifts upon binding to G4, which is valuable for detecting G4s. The association constants for the TPE-Im/G4 complex were evaluated in different solution conditions via isothermal titration calorimetry (ITC), and its binding modes were explored by molecular modeling showing a groove-binding mechanism. The stabilization of G4 by TPE-Im has been assessed by Fluorescence Resonance Energy Transfer (FRET) melting assays, which show a strong stabilization (ΔT 1/2 around +20°C), together with a specificity toward G4 with respect to double-stranded DNA.

Interaction of POPC, DPPC, and POPE with the µ opioid receptor: A coarse-grained molecular dynamics study.

The µ opioid receptor (µOR), which is part of the G protein-coupled receptors family, is a membrane protein that is modulated by its lipid environment. In the present work, we model µOR in three different membrane systems: POPC (1-palmitoyl-2-oleoyl-sn-glycero-3-phosphocholine), POPE (1-palmitoyl-2-oleoyl-sn-glycero-3-phosphoethanolamine), and DPPC (1, 2-dipalmitoyl-sn-glycero-3-phosphocholine) through 45 µs molecular dynamics (MD) simulations at the coarse-grained level. Our theoretical studies provide new insights to the lipid-induced modulation of the receptor. Particularly, to characterize how µOR interacts with each lipid, we analyze the tilt of the protein, the number of contacts occurring between the lipids and each amino acid of the receptor, and the µOR-lipid interface described as a network graph. We also analyze the variations in the number and the nature of the protein contacts that are induced by the lipid structure. We show that POPC interacts preferentially with helix 1 (H1) and helices H5-H6, POPE, with H5-H6 and H6-H7, and DPPC, with H4 and H6. We demonstrate how each of the three lipids shape the structure of the µOR.

Detection of the Enzymatic Cleavage of DNA through Supramolecular Chiral Induction to a Cationic Polythiophene.

Water-soluble π-conjugated polymers are increasingly envisioned in biosensors, in which their unique optical and electronic properties permit a highly sensitive detection of biomolecular targets. In particular, cationic π-conjugated polymers are attractive for DNA sensing technologies, through the use of the fluorescence signals either in physiological solutions or in thin films. However, in the context of enzymatic activity assays, fluorescence-based methods require covalently labeling DNA with a dye or an antibody and are limited to short time scale due to dye photobleaching. In this frame, we report here a novel possible approach to probe the cleavage of DNA by a restriction enzyme, in continuous and without covalently labeled DNA substrate. This is achieved by exploiting unique chiroptical signals arising from the chiral induction of DNA to a poly[3-(6'-(trimethylphosphonium)hexyl)thiophene-2,5-diyl] upon interaction. The cleavage of DNA by HpaI, an endonuclease enzyme, is monitored through circular dichroism (CD) signals in the spectral range where the polymer absorbs light, i.e., far away from the spectral ranges of both DNA and the enzyme. We compare the results to a conventional noncontinuous assay by polyacrylamide gel electrophoresis, and we demonstrate that induced CD signals are effective in probing the enzymatic activity. By means of molecular dynamics simulations and calculations of CD spectra, we bring molecular insights into the structure of DNA/polymer supramolecular complexes before and after the cleavage of DNA. We show that the cleavage of DNA modifies the dynamics and the organization of the polymer backbone induced by the DNA helix. Altogether, our results provide detailed spectroscopic and structural insights into the enzymatic cleavage of DNA in interaction with a π-conjugated polymer, which could be helpful for developing chiroptical detection tools to monitor the catalytic activity in real time.

Binding Modes and Selectivity of Ruthenium Complexes to Human Telomeric DNA G-Quadruplexes.

Metal complexes constitute an important class of DNA binders. In particular, a few ruthenium polyazaaromatic complexes are attractive as "light switches" because of their strong luminescence enhancement upon DNA binding. In this paper, a comprehensive study on the binding modes of several mononuclear and binuclear ruthenium complexes to human telomeric sequences, made of repeats of the d(TTAGGG) fragment is reported. These DNA sequences form G-quadruplexes (G4s) at the ends of chromosomes and constitute a relevant biomolecular target in cancer research. By combining spectroscopy experiments and molecular modelling simulations, several key properties are deciphered: the binding modes, the stabilization of G4 upon binding, and the selectivity of these complexes towards G4 versus double-stranded DNA. These results are rationalized by assessing the possible deformation of G4 and the binding free energies of several binding modes via modelling approaches. Altogether, this comparative study provides fundamental insights into the molecular recognition properties and selectivity of Ru complexes towards this important class of DNA G4s.

Self-assembly and chiroptical properties in supramolecular complexes of adenosine phosphates and guanidinium-bispyrene.

Supramolecular systems that respond to the hydrolysis of adenosine phosphates (APs) are attractive for biosensing and to fabricate bioinspired self-assembled materials. Here, we report on the formation of supramolecular complexes between an achiral guanidinium derivative bearing two pyrene moieties, with each of the three adenosine phosphates: AMP, ADP, and ATP. By combining results from circular dichroism spectroscopy and molecular modeling simulations, we explore the induced chirality, the dynamics of the complexes, and the interactions at play, which altogether provide insights into the supramolecular self-assembly between APs and the guanidinium-bispyrene. Finally, we identify the chiroptical signals of interest in mixtures of the guanidinium derivative with the three APs in different proportions. This study constitutes a basis to evolve toward a chiroptical detection of the hydrolysis of APs based on organic supramolecular probes.

Photomodulation of DNA-Templated Supramolecular Assemblies.

A new type of DNA ligand that contains a phosphate-binding group and a photoresponsive azobenzene moiety is reported. When the azobenzene is in trans configuration, the ligand binds to the minor groove of a double-stranded DNA, whereas it partially desorbs upon trans-cis isomerisation with light. The ability to photoswitch the ligand upon interaction with DNA is evidenced by (chir)optical signatures, and deciphered by the differences of binding geometry, stability, and dynamics of the DNA/ligand complexes for the two isomers. We exploit these properties to photomodulate DNA-templated self-assembly, through the incorporation of another π-stacking DNA ligand, which together with the photoresponsive ligand form mixed supramolecular complexes along DNA. Our study demonstrates that well-designed photoresponsive DNA binders can be used to modulate multicomponent supramolecular DNA assemblies.

Toward a new and noninvasive diagnostic method of papillary thyroid cancer by using peptide vectorized contrast agents targeted to galectin-1.

The incidence of papillary thyroid cancer has increased these last decades due to a better detection. High prevalence of nodules combined with the low incidence of thyroid cancers constitutes an important diagnostic challenge. We propose to develop an alternative diagnostic method to reduce the number of useless and painful thyroidectomies using a vectorized contrast agent for magnetic resonance imaging. Galectin-1 (gal-1), a protein overexpressed in well-differentiated thyroid cancer, has been targeted with a randomized linear 12-mer peptide library using the phage display technique. Selected peptides have been conjugated to ultrasmall superparamagnetic particles of iron oxide (USPIO). Peptides and their corresponding contrast agents have been tested in vitro for their specific binding and toxicity. Two peptides (P1 and P7) were selected according to their affinity toward gal-1. Their binding has been revealed by immunohistochemistry on human thyroid cancer biopsies, and they were co-localized with gal-1 by immunofluorescence on TPC-1 cell line. Both peptides induce a decrease in TPC-1 cells' adhesion to gal-1 immobilized on culture plates. After coupling to USPIO, the peptides preserved their affinity toward gal-1. Their specific binding has been corroborated by co-localization with gal-1 expressed by TPC-1 cells and by their ability to compete with anti-gal-1 antibody. The peptides and their USPIO derivatives produce no toxicity in HepaRG cells as determined by MTT assay. The vectorized contrast agents are potential imaging probes for thyroid cancer diagnosis. Moreover, the two gal-1-targeted peptides prevent cancer cell adhesion by interacting with the carbohydrate-recognition domain of gal-1.

One-Pot Self-Assembly of Peptide-Based Cage-Type Nanostructures Using Orthogonal Ligations.

The designed arrangement of biomolecular entities within monodisperse nanostructures is an important challenge toward functional biomaterials. We report herein a method for the formation of water-soluble peptide-based cages using orthogonal ligation reactions-acylhydrazone condensation and thiol-maleimide addition. The results show that using preorganized cyclic peptides and heterobifunctional spacers as building blocks and a set of orthogonal and chemoselective ligation reactions enable cage formation in one pot from six components and through eight reactions. Molecular modelling simulations reveal the structural dynamics of these structures. Finally, we exploited the reactional dynamics of the acylhydrazone by demonstrating the controlled dissociation of the cage through directed component exchange.

Multiscale design of coarse-grained elastic network-based potentials for the µ opioid receptor.

Despite progress in computer modeling, most biological processes are still out of reach when using all-atom (AA) models. Coarse-grained (CG) models allow classical molecular dynamics (MD) simulations to be accelerated. Although simplification of spatial resolution at different levels is often investigated, simplification of the CG potential in itself has been less common. CG potentials are often similar to AA potentials. In this work, we consider the design and reliability of purely mechanical CG models of the µ opioid receptor (µOR), a G protein-coupled receptor (GPCR). In this sense, CG force fields (FF) consist of a set of holonomic constraints guided by an elastic network model (ENM). Even though ENMs are used widely to perform normal mode analysis (NMA), they are not often implemented as a single FF in the context of MD simulations. In this work, various ENM-like potentials were investigated by varying their force constant schemes and connectivity patterns. A method was established to systematically parameterize ENM-like potentials at different spatial resolutions by using AA data. To do so, new descriptors were introduced. The choice of conformation descriptors that also include flexibility information is important for a reliable parameterization of ENMs with different degrees of sensitivity. Hence, ENM-like potentials, with specific parameters, can be sufficient to accurately reproduce AA MD simulations of µOR at highly coarse-grained resolutions. Therefore, the essence of the flexibility properties of µOR can be captured with simple models at different CG spatial resolutions, opening the way to mechanical approaches to understanding GPCR functions. Graphical Abstract All atom structure, residue interaction network and coarse-grained elastic network models of the µ opioid receptor (µOR).

On the modularity of the intrinsic flexibility of the µ opioid receptor: a computational study.

The µ opioid receptor (µOR), the principal target to control pain, belongs to the G protein-coupled receptors (GPCRs) family, one of the most highlighted protein families due to their importance as therapeutic targets. The conformational flexibility of GPCRs is one of their essential characteristics as they take part in ligand recognition and subsequent activation or inactivation mechanisms. It is assessed that the intrinsic mechanical properties of the µOR, more specifically its particular flexibility behavior, would facilitate the accomplishment of specific biological functions, at least in their first steps, even in the absence of a ligand or any chemical species usually present in its biological environment. The study of the mechanical properties of the µOR would thus bring some indications regarding the highly efficient ability of the µOR to transduce cellular message. We therefore investigate the intrinsic flexibility of the µOR in its apo-form using all-atom Molecular Dynamics simulations at the sub-microsecond time scale. We particularly consider the µOR embedded in a simplified membrane model without specific ions, particular lipids, such as cholesterol moieties, or any other chemical species that could affect the flexibility of the µOR. Our analyses highlighted an important local effect due to the various bendability of the helices resulting in a diversity of shape and volume sizes adopted by the µOR binding site. Such property explains why the µOR can interact with ligands presenting highly diverse structural geometry. By investigating the topology of the µOR binding site, a conformational global effect is depicted: the correlation between the motional modes of the extra- and intracellular parts of µOR on one hand, along with a clear rigidity of the central µOR domain on the other hand. Our results show how the modularity of the µOR flexibility is related to its pre-ability to activate and to present a basal activity.