Influence of chemical modifications of the crystallophore on protein nucleating properties and supramolecular interactions network.

Crystallophores are lanthanide complexes that have demonstrated outstanding induction of crystallization for various proteins. This article explores the effect of tailored modifications of the crystallophore first generation and their impact on the nucleating properties, and protein crystal structures. Through high-throughput crystallization experiments and dataset analysis, we evaluated the effectiveness of these variants, in comparison to the first crystallophore generation G1. In particular, the V1variant, featuring a propanol pendant arm, demonstrated the ability to produce new crystallization conditions for the proteins tested (hen-egg white lysozyme, proteinase K and thaumatin). Structural analysis performed in the case of hen egg-white lysozyme along with Molecular Dynamics simulations, highlights V1's unique behavior, taking advantage of the flexibility of its propanol arm to explore different protein surfaces and form versatile supramolecular interactions.

One touch is all it takes: the supramolecular interaction between ubiquitin and lanthanide complexes revisited by paramagnetic NMR and molecular dynamics.

The supramolecular interaction between lanthanide complexes and proteins is at the heart of numerous chemical and biological studies. Some of these complexes have demonstrated remarkable interaction properties with proteins or peptides in solution and in the crystalline state. Here we have used the paramagnetism of lanthanide ions to characterize the affinity of two lanthanide complexes for ubiquitin. As the interaction process is dynamic, the acquired NMR data only reflect the time average of the different steps. We have used molecular dynamics (MD) simulations to get a deeper insight into the detailed interaction scenario at the microsecond scale. This NMR/MD approach enabled us to establish that the tris-dipicolinate complex interacts specifically with arginines and lysines, while the crystallophore explores the protein surface through weak interactions with carboxylates. These observations shed new light on the dynamic interaction properties of these complexes, which will ultimately enable us to propose a crystallization mechanism.

Situating the phosphonated calixarene-cytochrome C association by molecular dynamics simulations.

Protein-calixarenes binding plays an increasingly central role in many applications, spanning from molecular recognition to drug delivery strategies and protein inhibition. These ligands obey a specific bio-supramolecular chemistry, which can be revealed by computational approaches, such as molecular dynamics simulations. In this paper, we rely on all-atom, explicit-solvent molecular dynamics simulations to capture the electrostatically driven association of a phosphonated calix-[4]-arene with cytochome-C, which critically relies on surface-exposed paired lysines. Beyond two binding sites identified in direct agreement with the x-ray structure, the association has a larger structural impact on the protein dynamics. Then, our simulations allow a direct comparison to analogous calixarenes, namely, sulfonato, similarly reported as "molecular glue." Our work can contribute to a robust in silico predictive tool to assess binding sites for any given protein of interest for crystallization, with the specificity of a macromolecular cage whose endo/exo orientation plays a role in the binding.

DNA-Histone Cross-Link Formation via Hole Trapping in Nucleosome Core Particles.

Radical cations (holes) produced in DNA by ionizing radiation and other oxidants yield DNA-protein cross-links (DPCs). Detailed studies of DPC formation in chromatin via this process are lacking. We describe here a comprehensive examination of DPC formation within nucleosome core particles (NCPs), which are the monomeric component of chromatin. DNA holes are introduced at defined sites within NCPs that are constructed from the bottom-up. DPCs form at DNA holes in yields comparable to those of alkali-labile DNA lesions that result from water trapping. DPC-forming efficiency and site preference within the NCP are dependent on translational and rotational positioning. Mass spectrometry and the use of mutant histones reveal that lysine residues in histone N-terminal tails and amino termini are responsible for the DPC formation. These studies are corroborated by computational simulation at the microsecond time scale, showing a wide range of interactions that can precede DPC formation. Three consecutive dGs, which are pervasive in the human genome, including G-quadruplex-forming sequences, are sufficient to produce DPCs that could impact gene expression.

Free energy landscape and thermodynamics properties of novel mutations in PncA of pyrazinamide resistance isolates of Mycobacterium tuberculosis.

Pyrazinamide (PZA) is one of the first-line antituberculosis therapy, active against non-replicating Mycobacterium tuberculosis (Mtb). The conversion of PZA into pyrazinoic acid (POA), the active form, required the activity of pncA gene product pyrazinamidase (PZase) activity. Mutations occurred in pncA are the primary cause behind the PZA resistance. However, the resistance mechanism is important to explore using high throughput computational approaches. Here we aimed to explore the mechanism of PZA resistance behind novel P62T, L120R, and V130M mutations in PZase using 200 ns molecular dynamics (MD) simulations. MD simulations were performed to observe the structural changes for these three mutants (MTs) compared to the wild types (WT). Root means square fluctuation, the radius of gyration, free energy landscape, root means square deviation, dynamic cross-correlation motion, and pocket volume were found in variation between WT and MTs, revealing the effects of P62T, L120R, and V130M. The free energy conformational landscape of MTs differs significantly from the WT system, lowering the binding of PZA. The geometric shape complementarity of the drug (PZA) and target protein (PZase) further confirmed that P62T, L120R, and V130M affect the protein structure. These effects on PZase may cause vulnerability to convert PZA into POA.Communicated by Ramaswamy H. Sarma.

Self-assembly of achiral building blocks into chiral cyclophanes using non-directional interactions.

The diastereoselective assembly of achiral constituents through a single spontaneous process into complex covalent architectures bearing multiple stereogenic elements still remains a challenge for synthetic chemists. Here, we show that such an extreme level of control can be achieved by implementing stereo-electronic information on synthetic organic building blocks and templates and that non-directional interactions (i.e., electrostatic and steric interactions) can transfer this information to deliver, after self-assembly, high-molecular weight macrocyclic species carrying up to 16 stereogenic elements. Beyond the field of supramolecular chemistry, this proof of concept should stimulate the on-demand production of highly structured polyfunctional architectures.

Probing the dynamical interaction of the para-sulfonato-calix[4]arene with an antifungal protein.

Calixarenes are hallmark molecules in supramolecular chemistry as hosts for small ligands. They have also conversely proved their interest as ligands toward assisted co-crystallization of proteins. These functionalized macrocycles target positively-charged residues, and notably surface-exposed lysines, with a site-selectivity finely characterized experimentally, but that remains to be assessed. Relying on a tailored molecular dynamics simulations protocol, we explore the association of para-sulfonato-calix[4]arenes with an antifungal protein, as a small yet most competitive system with 13 surface-exposed lysines. Our computational approach probes de novo the electrostatically-driven interaction, ruled out by a competition with salt bridges, corroborating the presence of two main binding sites probed by X-ray. The attach-pull-release (APR) method provides a very good assessment of the overall binding free energy measured experimentally (-6.42 ± 0.5 vs. -5.45 kcal mol-1 by isothermal titration calorimetry). This work also probes dynamic modifications upon ligand binding, and our computational protocol could be generalized to situate the supramolecular forces ruling out the calixarene-assisted co-crystallization of proteins.

Rationalizing the environment-dependent photophysical behavior of a DNA luminescent probe by classical and non-adiabatic molecular dynamics simulations.

Environment-sensitive fluorescent nucleoside analogs are of utmost importance to investigate the structure of nucleic acids, their intrinsic flexibility, and sequence-specific DNA- and RNA-binding proteins. The latter play indeed a key role in transcription, translation as well as in the regulation of RNA stability, localization and turnover, and many other cellular processes. The sensitivity of the embedded fluorophore to polarity, hydration, and base stacking is clearly dependent on the specific excited-state relaxation mechanism and can be rationalized combining experimental and computational techniques. In this work, we elucidate the mechanisms leading to the population of the triplet state manifold for a versatile nucleobase surrogate, namely the 2-thienyl-3-hydroxychromone in gas phase, owing to non-adiabatic molecular dynamics simulations. Furthermore, we analyze its behavior in the B-DNA environment via classical molecular dynamics simulations, which evidence a rapid extrusion of the adenine facing the 2-thienyl-3-hydroxychromone nucleobase surrogate. Our simulations provide new insights into the dynamics of this family of chromophores, which could give rise to an integrated view and a fine tuning of their photochemistry, and namely the role of excited-state intramolecular proton transfer for the rational design of the next generation of fluorescent nucleoside analogs.

A new G-quadruplex-specific photosensitizer inducing genome instability in cancer cells by triggering oxidative DNA damage and impeding replication fork progression.

Photodynamic therapy (PDT) ideally relies on the administration, selective accumulation and photoactivation of a photosensitizer (PS) into diseased tissues. In this context, we report a new heavy-atom-free fluorescent G-quadruplex (G4) DNA-binding PS, named DBI. We reveal by fluorescence microscopy that DBI preferentially localizes in intraluminal vesicles (ILVs), precursors of exosomes, which are key components of cancer cell proliferation. Moreover, purified exosomal DNA was recognized by a G4-specific antibody, thus highlighting the presence of such G4-forming sequences in the vesicles. Despite the absence of fluorescence signal from DBI in nuclei, light-irradiated DBI-treated cells generated reactive oxygen species (ROS), triggering a 3-fold increase of nuclear G4 foci, slowing fork progression and elevated levels of both DNA base damage, 8-oxoguanine, and double-stranded DNA breaks. Consequently, DBI was found to exert significant phototoxic effects (at nanomolar scale) toward cancer cell lines and tumor organoids. Furthermore, in vivo testing reveals that photoactivation of DBI induces not only G4 formation and DNA damage but also apoptosis in zebrafish, specifically in the area where DBI had accumulated. Collectively, this approach shows significant promise for image-guided PDT.

How SARS-CoV-2 Alters the Regulation of Gene Expression in Infected Cells.

Nonstructural accessory proteins in viruses play a key role in hijacking the basic cellular mechanisms, which is essential to promote the virus survival and evasion of the immune system. The immonuglobulin-like open reading frame 8 (ORF8) protein expressed by SARS-CoV-2 accumulates in the nucleus and may influence the regulation of the gene expression in infected cells. In this contribution, by using microsecond time-scale all-atom molecular dynamics simulations, we unravel the structural bases behind the epigenetic action of ORF8. In particular, we highlight how the protein is able to form stable aggregates with DNA through a histone tail-like motif, and how this interaction is influenced by post-translational modifications, such as acetylation and methylation, which are known epigenetic markers in histones. Our work not only clarifies the molecular mechanisms behind the perturbation of the epigenetic regulation caused by the viral infection but also offers an unusual perspective which may foster the development of original antivirals.

Molecular Basis of the pH-Controlled Maturation of the Tick-Borne Encephalitis Flavivirus.

Flaviviruses are enveloped viruses causing high public concerns. Their maturation spans several cellular compartments having different pH. Thus, complex control mechanisms are in place to avoid premature maturation. Here we report the dynamical behavior at neutral and acidic pH of the precursor of the membrane fusion protein E of tick-borne encephalitis, showing the different stabilizations of the E dimer and the role played by the small fusion-assisting protomer (pr). The comprehension, at atomic resolution, of the fine regulation of viral maturation will be fundamental to the development of efficient strategies against emerging viral threats.

Correction: Light-induced in situ chemical activation of a fluorescent probe for monitoring intracellular G-quadruplex structures.

Correction for 'Light-induced in situ chemical activation of a fluorescent probe for monitoring intracellular G-quadruplex structures' by Marco Deiana et al., Nanoscale, 2021, 13, 13795-13808, https://doi.org/10.1039/D1NR02855C.

Dynamics and energetics of PCBP1 binding to severely oxidized RNA.

Oxidatively generated lesions such as 8-oxo-7, 8-dihydroguanine (8-oxoG) on RNA strands constitute a hallmark marker of the oxidative stress in the cell. Poly-C binding protein 1 (PCBP1) is able to specifically recognize severely damaged RNA strands containing two 8-oxoG lesions separated by five nucleobases, which trigger a signaling pathway leading to cell apoptosis. We apply an in silico protocol based on microsecond timescale all-atom classical molecular dynamics simulations associated with conformational and energy analyses to unveil the specific recognition mechanism at a molecular level. By comparing the RNA and protein behavior for sequences with six different damage profiles, our results highlight an allosteric mechanism, allowing a stronger binding of the oxidized guanine at position 9 only if another 8-oxoG lesion is present at position 15, in full agreement with experiments. We assess the role of lysine K23 and the additional ketone group of the oxidized guanine, thanks to computational site-directed mutagenesis.

Ni-Centered Coordination-Induced Spin-State Switching Triggered by Electrical Stimulation.

We herein report the synthesis and magnetic properties of a Ni(II)-porphyrin tethered to an imidazole ligand through a flexible electron-responsive mechanical hinge. The latter is capable of undergoing a large amplitude and fully reversible folding motion under the effect of electrical stimulation. This redox-triggered movement is exploited to force the axial coordination of the appended imidazole ligand onto the square-planar Ni(II) center, resulting in a change in its spin state from low spin (S = 0) to high spin (S = 1) proceeding with an 80% switching efficiency. The driving force of this reversible folding motion is the π-dimerization between two electrogenerated viologen cation radicals. The folding motion and the associated spin state switching are demonstrated on the grounds of NMR, (spectro)electrochemical, and magnetic data supported by quantum calculations.

Capturing the Recognition Dynamics of para-Sulfonato-calix[4]arenes by Cytochrome c: Toward a Quantitative Free Energy Assessment.

Calix[n]arenes' selective recognition of protein surfaces covers a broad range of timely applications, from controlling protein assembly and crystallization to trapping partially disordered proteins. Here, the interaction of para-sulfonated calix-[4]-arenes with cytochrome c is investigated through all-atom, explicit water molecular dynamics simulations which allow characterization of two binding sites in quantitative agreement with experimental evidence. Free energy calculations based on the MM-PBSA and the attach-pull-release (APR) methods highlight key residues implicated in the recognition process and provide binding free energy results in quantitative agreement with isothermal titration calorimetry. Our study emphasizes the role of MD simulations to capture and describe the "walk" of sulfonated calix-[4]-arenes on the cytochrome c surface, with the arginine R13 as a pivotal interacting residue. Our MD investigation allows, through the quasi-harmonic multibasin (QHMB) method, probing an allosteric reinforcement of several per-residue interactions upon calixarene binding, which suggests a more complex mode of action of these supramolecular auxiliaries.

How Fragile We Are: Influence of Stimulator of Interferon Genes (STING) Variants on Pathogen Recognition and Immune Response Efficiency.

The stimulator of interferon genes (STING) protein is a cornerstone of the human immune response. Its activation by cGAMP in the presence of cytosolic DNA stimulates the production of type I interferons and inflammatory cytokines. In the human population, several STING variants exist and exhibit dramatic differences in their activity, impacting the efficiency of the host defense against infections. Understanding the molecular mechanisms of these variants opens perspectives for personalized medicine treatments against diseases such as viral infections, cancers, or autoinflammatory diseases. Through microsecond-scale molecular modeling simulations, contact analyses, and machine learning techniques, we reveal the dynamic behavior of four STING variants (wild type, G230A, R293Q, and G230A/R293Q) and rationalize the variability of efficiency observed experimentally. Our results show that the decrease in STING activity is linked to a stiffening of key structural elements of the binding cavity together with changes in the interaction patterns within the protein.

Reactivity of Singlet Oxygen with DNA, an Update.

The reactivity of singlet oxygen with DNA constituents and in particular with the guanine base has been studied during more than four decades but the exact mechanisms by which such a reactive oxygen species reacts with DNA is still a matter of debate. In this review article, a summary of the data that were obtained from several laboratories and using complementary experimental and theoretical approaches are presented. Reaction mechanisms of 1 O2 with guanine and its oxidation product 8-oxo7,8-dihydroguanine are presented both at the nucleoside level and when the base is inserted into DNA since significant differences have been observed. Efforts have been made to propose tentative mechanisms to explain the conflicting results that were sometimes reported and hypotheses have been put forward to tentatively explain still contradictory observations.

The Behavior of Triplet Thymine in a Model B-DNA Strand. Energetics and Spin Density Localization Revealed by ab initio Molecular Dynamics Simulations.

Among the naturally occurring nucleobases, thymine presents the lowest triplet state, hence it represents a hotspot for energy transfer and photosensitization. In turn, the population of the triplet state may lead to thymine dimerization and hence to the production of toxic DNA lesions and has been the subject of intensive theoretical and experimental investigations. Relying on QM/MM molecular dynamics simulations, we have sought to situate the energy of the lowest triplet state of thymine embedded in a B-DNA environment. The energy gap varies between 305 and 329 kJ mol-1 when a single thymine is treated at the quantum chemistry level, depending on its position in the model double-stranded 16-bp oligonucleotide. The energy of triplet state decreases up to 300 kJ mol-1 , due to polarization effects, when we consider coupled stacked nucleobases up to the inclusion of four nucleobases. Our value lies in good agreement with the energy inferred experimentally by Miranda and coworkers (270 kJ mol-1 ), and our theoretical exploration opens the door to investigations toward other more complex and biologically relevant environments, such as thymines embedded in nucleosome core particles. Our investigations also provide a reference for further studies using semi-empirical approaches such as density functional-based tight-binding, allowing to further rationalize sequence effects.

Molecular Mechanisms Associated with Clustered Lesion-Induced Impairment of 8-oxoG Recognition by the Human Glycosylase OGG1.

The 8-oxo-7,8-dihydroguanine, referred to as 8-oxoG, is a highly mutagenic DNA lesion that can provoke the appearance of mismatches if it escapes the DNA Damage Response. The specific recognition of its structural signature by the hOGG1 glycosylase is the first step along the Base Excision Repair pathway, which ensures the integrity of the genome by preventing the emergence of mutations. 8-oxoG formation, structural features, and repair have been matters of extensive research; more recently, this active field of research expended to the more complicated case of 8-oxoG within clustered lesions. Indeed, the presence of a second lesion within 1 or 2 helix turns can dramatically impact the repair yields of 8-oxoG by glycosylases. In this work, we use µs-range molecular dynamics simulations and machine-learning-based postanalysis to explore the molecular mechanisms associated with the recognition of 8-oxoG by hOGG1 when embedded in a multiple-lesion site with a mismatch in 5' or 3'. We delineate the stiffening of the DNA-protein interactions upon the presence of the mismatches, and rationalize the much lower repair yields reported with a 5' mismatch by describing the perturbation of 8-oxoG structural features upon addition of an adjacent lesion.

Influence of Divalent Cations in the Protein Crystallization Process Assisted by Lanthanide-Based Additives.

The use of lanthanide complexes as powerful auxiliaries for biocrystallography prompted us to systematically analyze the influence of the commercial crystallization kit composition on the efficiency of two lanthanide additives: [Eu(DPA)3]3- and Tb-Xo4. This study reveals that the tris(dipicolinate) complex presents a lower chemical stability and a strong tendency toward false positives, which are detrimental for its use in a high-throughput robotized crystallization platform. In particular, the crystal structures of (Mg(H2O)6)3[Eu(DPA)3]2·7H2O (1), {(Ca(H2O)4)3[Eu(DPA)3]2}n·10nH2O (2), and {Cu(DPA)(H2O)2}n (3), resulting from spontaneous crystallization in the presence of a divalent alkaline-earth cation and transmetalation, are reported. On the other hand, Tb-Xo4 is perfectly soluble in the crystallization media, stable in the presence of alkaline-earth dications, and slowly decomposes (within days) by transmetalation with transition metals. The original structure of [Tb4L4(H2O)4]Cl4·15H2O (4) is also described, where L represents a bis(pinacolato)triazacyclononane ligand. This paper also highlights a potential synergy of interactions between Tb-Xo4 and components of the crystallization mixtures, leading to the formation of complex adducts like {AdkA/Tb-Xo4/Mg2+/glycerol} in the protein binding sites. The observation of such multicomponent adducts illustrated the complexity and versatility of the supramolecular chemistry occurring at the surface of the proteins.