Probing interaction of a trilysine peptide with DNA underlying formation of guanine-lysine cross-links: insights from molecular dynamics.

DNA-protein cross-links constitute bulky DNA lesions that interfere with the cellular machinery. Amongst these stable covalently tethered adducts, the efficient nucleophilic addition of the free amino group of lysines onto the guanine radical cation has been evidenced. In vitro addition of a trilysine peptide onto a guanine radical cation generated in a TGT oligonucleotide is so efficient that competitive addition of a water molecule, giving rise to 8-oxo-7,8-dihydroguanine, is not observed. This suggests a spatial proximity between guanine and lysine for the stabilization of the prereactive complex. We report all-atom microsecond scale molecular dynamics simulations that probe the structure and interactions of the trilysine peptide (KKK) with two oligonucleotides. Our simulations reveal a strong, electrostatically driven yet dynamic interaction, spanning several association modes. Furthermore, the presence of neighbouring cytosines has been identified as a factor favoring KKK binding. Relying on ab initio molecular dynamics on a model system constituted of guanine and methylammonium, we also corroborate a mechanistic pathway involving fast deprotonation of the guanine radical cation followed by hydrogen transfer from ammonium leaving as a result a nitrogen reactive species that can subsequently cross-link with guanine. Our study sheds new light on a ubiquitous mechanism for DNA-protein cross-links also stressing out possible sequence dependences.

Molecular Dynamics Insights into Polyamine-DNA Binding Modes: Implications for Cross-Link Selectivity.

Biogenic polyamines, which play a role in DNA condensation and stabilization, are ubiquitous and are found at millimolar concentration in the nucleus of eukaryotic cells. The interaction modes of three polyamines-putrescine (Put), spermine (Spm), and spermidine (Spd)-with a self-complementary 16 base pair (bp) duplex, are investigated by all-atom explicit-solvent molecular dynamics. The length of the amine aliphatic chain leads to a change of the interaction mode from minor groove binding to major groove binding. Through all-atom dynamics, noncovalent interactions that stabilize the polyamine-DNA complex and prefigure the reactivity, leading to the low-barrier formation of deleterious DNA-polyamine cross-links, after one-electron oxidation of a guanine nucleobase, are unraveled. The binding strength is quantified from the obtained trajectories by molecular mechanics generalized Born surface area post-processing (MM-GBSA). The values of binding free energies provide the same affinity order, Put

Recognition of a tandem lesion by DNA bacterial formamidopyrimidine glycosylases explored combining molecular dynamics and machine learning.

The combination of several closely spaced DNA lesions, which can be induced by a single radical hit, constitutes a hallmark in the DNA damage landscape and radiation chemistry. The occurrence of such a tandem base lesion gives rise to a strong coupling with the double helix degrees of freedom and induces important structural deformations, in contrast to DNA strands containing a single oxidized nucleobase. Although such complex lesions are known to be refractory to repair by DNA glycosylases, there is still a lack of structural evidence to rationalize these phenomena. In this contribution, we explore, by numerical modeling and molecular simulations, the behavior of the bacterial glycosylase responsible for base excision repair (MutM), specialized in excising oxidatively-damaged defects such as 7,8-dihydro-8-oxoguanine (8-oxoG). The difference in lesion recognition between a simple damage and a tandem lesion featuring an additional abasic site is assessed at atomistic resolution owing to microsecond molecular dynamics simulations and machine learning postprocessing, allowing to extensively pinpoint crucial differences in the interaction patterns of the damaged bases. Our results reveal substantial changes in the interaction network surrounding the 8-oxoG upon addition of an adjacent abasic site, leading to the perturbation of the intercalation triad which is crucial for lesion recognition and processing. The recognition process might also be impacted by a more constrained MutM-DNA binding upon tandem damage, as shown by the machine learning post-processing. This work advocates for the use of such high throughput numerical simulations for exploring the complex combinatorial chemistry of tandem DNA lesions repair and more generally local multiple damaged sites of the utmost significance in radiation chemistry.

Effect of the Ligand Binding Strength on the Morphology of Functionalized Gold Nanoparticles.

Functionalized gold nanoparticles are investigated by density functional theory calculations in the context of cancer radiotherapy. Several typical experimental shapes, including nanostars, nanospheres, and nanorods, are modeled by optimizing Au clusters covered by organic monolayers composed of hydrated short-chain polyethylene glycol (PEG) ligands. The PEGylation stabilizes significantly the stellation of decahedral Au54 by deforming significantly its geometry at the spikes. The higher stability of the PEG molecules adsorbed on this stellated nanocluster with respect to the more spherical icosahedral Au55 and truncated octahedral Au79 leads to a larger energy cost to desorb them and thus a weaker propensity for the starred nanoparticle to exchange ligands with the cell membrane, in agreement with experiments. These results open interesting possibilities for advancing our understanding of the cellular uptake of gold nanoparticles.

A Water Solvation Shell Can Transform Gold Metastable Nanoparticles in the Fluxional Regime.

Solvated gold nanoparticles have been modeled in the fluxional regime by density functional theory including dispersion forces for an extensive set of conventional morphologies. The study of isolated adsorption of one water molecule shows that the most stable adsorption forms are similar (corners and edges) regardless of the nanoparticle shape and size, although the adsorption strength differs significantly (0.15 eV). When a complete and explicit water solvation shell interacts with gold nanoclusters, metastable in vacuum and presenting a predominance of (100) square facets (ino-decahedra Au55 and Au147), these nanoparticles are found unstable and transform into the closest morphologies exhibiting mainly (111) triangular facets and symmetries. The corresponding adsorption strength per water molecule becomes independent of shape and size and is enhanced by the formation of two hydrogen bonds on average. For applications in radiotherapy, this study suggests that the shapes of small gold nanoparticles should be homogenized by interacting with the biological environment.