Ultrafast conformation-dependent charge transfer in N, N, N', N'-tetramethyl-1,3-propanediamine: Effect of flexible carbon skeleton on electron lone pair interactions

Flexible three-carbon skeleton makes N, N, N', N'-tetramethyl-1,3-propanediamine (TMPDA) an important diamine system to investigate the conformation-dependent electron lone pair interactions and charge delocalization. The charge transfer process linked to structural motions of the three-carbon skeleton has been monitored in real time by the Rydberg electron binding energy (BE) spectra of TMPDA coupled with quantum chemical calculations. Optical excitation to the 3p state with a 200 nm pump pulse initially generated a localized charge on one of the two nitrogen atoms that may partially transfer to the other one. Rapid internal conversion (IC) from the 3p to 3s state occurred within 430 fs, resulting in an initial charge delocalized 3s_h/3s_l population ratio of 23.6 %/76.4 %. A final 3s_h/3s_l (51.9 %/48.1 %) equilibrium proceeded within about 2.64 ps. The 3s_h (TTTT+, GG'TG+ and G'GG'G+) and 3s_l (GG'GG'+ and GG'G'G+) (see text for structure definitions) are identified as the extended and folded conformers, respectively. Two types of electron lone pair interactions, i.e., through-space interaction (TSI) and through-bond interaction (TBI), are found to coexist in TMPDA to drive charge transfer. The GG'GG'+ and GG'G'G+ structures exhibit TSI, while the TTTT+ structure shows TBI. The GG'TG+ and G'GG'G+ structures exhibit both TSI and TBI. Flexible three-carbon skeleton provide more opportunities for the two N-electron lone pairs to overlap in space (i.e., TSI), making TMPDA to be favorable for the most stably folded conformation.

What tunes guanine ionization potential in a nucleosome? An all-in-one systematic QM/MM assessment

Guanine radical cations are precursors to oxidatively induced DNA lesions, and the determination of oxidative DNA hot spots beyond oligonucleotides remains a current challenge. In order to rationalize the finetuned ionization properties of the ∼60 guanines in a nucleosome core particle, we report a robust molecular dynamics-then-FO-DFTB/MM (fragment-orbital tight-binding density functional theory/molecular mechanics) simulation protocol spanning 20 μs. Our work allows us to identify several factors governing guanine ionization potential and map oxidative hotspots. Our results highlight the predominant role of the proximity of positively charged histone residues in the modulation of the guanine ionization potential up to 0.6 eV. Consequently, fast, long-range hole transfer in nucleosomal DNA could be tuned by the proximity of histone tails, which differs, from a biological point of view, on the chromatin state.

Contribution of oxidation reactions to photo-induced damage to cellular DNA

This review article is aimed at providing updated information on the contribution of immediate and delayed oxidative reactions to the photo-induced damage to cellular DNA/skin under exposure to UVB/UVA radiations and visible light. Low-intensity UVC and UVB radiations that operate predominantly through direct excitation of the nucleobases are very poor oxidizing agents giving rise to very low amounts of 8-oxo-7,8-dihydroguanine and DNA strand breaks with respect to the overwhelming bipyrimidine dimeric photoproducts. The importance of these two classes of oxidatively generated damage to DNA significantly increases together with a smaller contribution of oxidized pyrimidine bases upon UVA irradiation. This is rationalized in terms of sensitized photooxidation reactions predominantly mediated by singlet oxygen together with a small contribution of hydroxyl radical that appear to also be implicated in the photodynamic effects of the blue light component of visible light. Chemiexcitation-mediated formation of "dark" cyclobutane pyrimidine dimers in UVA-irradiated melanocytes is a recent major discovery that implicates in the initial stage, a delayed generation of reactive oxygen and nitrogen species giving rise to triplet excited carbonyl intermediate and possibly singlet oxygen. High-intensity UVC nanosecond laser radiation constitutes a suitable source of light to generate pyrimidine and purine radical cations in cellular DNA via efficient biphotonic ionization.

Synthesis of Original Cyclic Dinucleotide Analogues Using the Sulfo-click Reaction

The stimulator of interferon genes (STING) protein plays a crucial role in the activation of the innate immune response. Activation of STING is initiated by cyclic dinucleotides (CDNs) which prompted the community to synthesize structural analogues to enhance their biological properties. We present here the synthesis and biological evaluation of four novel CDN analogues composed of an N-acylsulfonamide linkage. These CDNs were obtained in high overall yields via the sulfo-click reaction as a key step.

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.

Comprehensive analysis unveils altered binding kinetics of 5-/6-methylCytosine/adenine modifications in R2R3-DNA system

Recent studies have shown that DNA methylation is an important epigenetic marker. Two prominent forms are methylation of the C5 position of cytosine and methylation of the C6 position of adenine. Given the vital significance of DNA methylation, investigating the mechanisms that influence protein binding remains a compelling pursuit. This study used molecular dynamics simulations to investigate the binding patterns of R2R3 protein and four differentially methylated DNAs. The alanine scanning combined with interaction entropy method was used to identify key residues that respond to different methylation patterns. The order of protein binding ability to DNA is as follows: unmethylated DNA > A11 methylation (5'-A6mAC-3') (6m2A system) > A10 methylation (5'-6mAAC-3') (6m1A system) > both A10 and A11 methylation (5'-6mA6mAC-3') (6mAA system) > C12 methylation (5'-AA5mC-3') (5mC system). All methylation systems lead to the sixth α helix (H6) (residues D105 to L116) moving away from the binding interface, and in the 5mC and 6m1A systems, the third α helix (H3) (residues G54 to L65) exhibits a similar trend. When the positively charged amino acids in H3 and H6 move away from the binding interface, their electrostatic and van der Waals interactions with the negatively charged DNA are weakened. Structural changes induced by methylation contributed to the destabilization of the hydrogen bond network near the original binding site, except for the 6m2A system. Moreover, there is a positive correlation between the number of methylated sites and the probability of distorting the DNA structure. Our study explores how different methylation patterns affect binding and structural adaptability, and have implications for drug discovery and understanding diseases related to abnormal methylation.

Quantitative investigation of the effects of DNA modifications and protein mutations on MeCP2-MBD-DNA interactions

DNA methylation as an important epigenetic marker, has gained attention for the significance of three oxidative modifications (hydroxymethyl-C (hmC), formyl-C (fC), and carboxyl-C (caC)). Mutations occurring in the methyl-CpG-binding domain (MBD) of MeCP2 result in Rett. However, uncertainties persist regarding DNA modification and MBD mutation-induced interaction changes. Here, molecular dynamics simulations were used to investigate the underlying mechanisms behind changes due to different modifications of DNA and MBD mutations. Alanine scanning combined with the interaction entropy method was employed to accurately evaluate the binding free energy. The results show that MBD has the strongest binding ability for mCDNA, followed by caC, hmC, and fCDNA, with the weakest binding ability observed for CDNA. Further analysis revealed that mC modification induces DNA bending, causing residues R91 and R162 closer to the DNA. This proximity enhances van der Waals and electrostatic interactions. Conversely, the caC/hmC and fC modifications lead to two loop regions (near K112 and K130) closer to DNA, respectively. Furthermore, DNA modifications promote the formation of stable hydrogen bond networks, however mutations in the MBD significantly reduce the binding free energy. This study provides detailed insight into the effects of DNA modifications and MBD mutations on binding ability. It emphasizes the necessity for research and development of targeted Rett compounds that induce conformational compatibility between MBD and DNA, enhancing the stability and strength of their interactions.

Effectiveness of a novel gene nanotherapy based on putrescine for cancer treatment

Gene therapy has long been proposed for cancer treatment. However, the use of therapeutic nucleic acids presents several limitations such as enzymatic degradation, rapid clearance, and poor cellular uptake and efficiency. In this work we propose the use of putrescine, a precursor for higher polyamine biosynthesis for the preparation of cationic nanosystems for cancer gene therapy. We have formulated and characterized putrescine-sphingomyelin nanosystems (PSN) and studied their endocytic pathway and intracellular trafficking in cancer cells. After loading a plasmid DNA (pDNA) encoding the apoptotic Fas Ligand (FasL), we proved their therapeutic activity by measuring the cell death rate after treatment of MDA-MB-231 cells. We have also used xenografted zebrafish embryos as a first in vivo approach to demonstrate the efficacy of the proposed PSN-pDNA formulation in a more complex model. Finally, intratumoral and intraperitoneal administration to mice-bearing MDA-MB-231 xenografts resulted in a significant decrease in tumour cell growth, highlighting the potential of the developed gene therapy nanoformulation for the treatment of triple negative breast cancer.

Capturing Protein-Nucleic Acid Interactions by High-Intensity Laser-Induced Covalent Crosslinking

Interactions of DNA with structural proteins such as histones, regulatory proteins, and enzymes play a crucial role in major cellular processes such as transcription, replication and repair. The in vivo mapping and characterization of the binding sites of the involved biomolecules are of primary importance for a better understanding of genomic deployment that is implicated in tissue and developmental stage-specific gene expression regulation. The most powerful and commonly used approach to date is immunoprecipitation of chemically cross-linked chromatin (XChIP) coupled with sequencing analysis (ChIP-seq). While the resolution and the sensitivity of the high-throughput sequencing techniques have been constantly improved little progress has been achieved in the crosslinking step. Because of its low efficiency the use of the conventional UVC lamps remains very limited while the formaldehyde method was established as the "gold standard" crosslinking agent. Efficient biphotonic crosslinking of directly interacting nucleic acid-protein complexes by a single short UV laser pulse has been introduced as an innovative technique for overcoming limitations of conventionally used chemical and photochemical approaches. In this survey, the main available methods including the laser approach are critically reviewed for their ability to generate DNA-protein crosslinks in vitro model systems and cells.

Molecular investigation of the dual inhibition mechanism for targeted P53 regulator MDM2/MDMX inhibitors

Inhibitors that competitively bind MDM2/MDMX can block the inhibition of P53 by MDM2/MDMX and restart its tumor-suppressive effect. Molecular studies targeting MDM2/MDMX inhibitors have always been a hot topic in anticancer drug design. Although numerous inhibitors have been designed previously against MDM2/MDMX, their dual inhibition efficacy has not been demonstrated, and few studies assessed the general causes affecting the dual inhibition of MDM2/MDMX by these inhibitors. Here, molecular dynamics simulations and alanine scanning combined with the interaction entropy method were employed to precisely investigate whether 16 inhibitors could dually inhibit MDM2/MDMX and the similarities and differences in the interaction modes. Thereby addressing the key residue sites affecting dual inhibition. Residues L54/M53, I61/60, M62/61, Y67/66, and V93/92 of MDM2/MDMX, which are in corresponding positions in both protein structures, provide significant conditions for these inhibitors to bind to MDM2/MDMX tightly. In addition, most of these inhibitors prefer to bind MDM2 than MDMX, and residues H96 and I99 in MDM2 are attractive targets for inhibitors, resulting in inhibitors binding to MDM2/MDMX with different affinity. These key residues should be considered in the development of dual inhibitors. For these 16 inhibitors, most have dual inhibitory potential for MDM2/MDMX based on the binding affinity of the complexes. Still, it is questionable whether they can exert excellent dual inhibition considering the assessment of the hot-spots. At least their binding affinity for MDMX is not superior to that for MDM2 due to the difference in energy of the van der Waals interactions at the key sites. Furthermore, based on the analysis of three representative inhibitors (TUZ/HRH and HRQ with different binding preferences for MDM2/MDMX), 3-chloropyridine in TUZ leads to the differential binding affinity between the inhibitor and MDM2/MDMX. It readily forms hydrophobic interactions with the surrounding residues H96 and I99. But this phenomenon does not occur in the TUZ-MDMX system, implying the critical role of residues H96/P95 and I99/L98. And the completely different binding mechanism of HRQ binding to MDM2/MDMX explains its inability to inhibit MDM2 well. Thus, we are cautious about its dual inhibitory ability. Besides, HRH is more prone to strong van der Waals interactions with MDM2 than MDMX whereas its 2-chlorofluorobenzene is detrimental to this. We hope that these findings will provide reliable molecular insights for the screening and optimization of targeting MDM2/MDMX dual inhibitors.

Spatial and Temporal Resolution of the Oxygen-Independent Photoinduced DNA Interstrand Cross-Linking by a Nitroimidazole Derivative

DNA damage is ubiquitous in nature and is at the basis of emergent treatments such as photodynamic therapy, which is based on the activation of highly oxidative reactive oxygen species by photosensitizing O₂. However, hypoxia observed in solid tumors imposes the necessity to devise oxygen-independent modes of action able to induce DNA damage under a low oxygen concentration. The complexity of these DNA damage mechanisms in realistic environments grows exponentially when taking into account light absorption and subsequent excited-state population, photochemical and (photo)-redox reactions, the multiple species involved in different electronic states, noncovalent interactions, multiple reaction steps, and the large number of DNA reactive sites. This work tackles all the intricate reactivity of a photosensitizer based on a nitroimidazole derivative reacting toward DNA in solution under UV light exposition. This is performed through a combination of ground- and excited-state quantum chemistry, classical molecular dynamics, and hybrid QM/MM simulations to rationalize in detail the formation of DNA interstrand cross-links (ICLs) exerted by the noncanonical noncovalent photosensitizer. Unprecedented spatial and temporal resolution of these phenomena is achieved, revealing that the ICL is sequence-specific and that the fastest reactions take place at AT, GC, and GT steps involving either the opposite nucleobases or adjacent Watson-Crick base pairs. The N7 and O6 positions of guanine, the N7 and N3 sites of adenine, the N4 position of cytosine, and the O2 atom of thymine are deemed as the most nucleophile sites and are positively identified to participate in the ICL productions. This work provides a multiscale computational protocol to study DNA reactivity with noncovalent photosensitizers, and contributes to the understanding of therapies based on photoinduced DNA damage at molecular and electronic levels. In addition, we believe the depth understanding of these processes should assist the design of new photosensitizers considering their molecular size, electronic properties, and the observed regioselectivity toward nucleic acids.

Origin of the tight binding mode to ACE2 triggered by multi-point mutations in the omicron variant: a dynamic insight

The continuous spread of the newly emerged SARS-CoV-2 Omicron variant (B.1.1.529) has become an important reason for the surge in COVID-19 infections. Its numerous mutated residues containing key sites on the receptor-binding domain (RBD) undoubtedly pose new challenges for epidemic control. Although the preventive measures are becoming more sophisticated, the effects of mutations on the binding of the virus to the receptor protein remain to be elucidated. Here, we used molecular dynamics (MD) simulations to investigate the differences in the binding mode between the Omicron variant and the angiotensin-converting enzyme 2 (ACE2) compared to the wild-type strain (WT). Multi-point mutations in the Omicron variant RBD could cause the conformation shift in the large Loop (where T478K and E484A are located), which makes it easier to wrap the N-terminal helix of ACE2 and form tighter contacts. The stronger electrostatic interaction was the main reason for its enhanced binding affinity as compared to WT. This was due to the large number of positively charged patches (N440K, T478K, Q493R, Q498R, and Y505H) formed by the substitution of neutral amino acids at multiple sites. The appearance of these highly polar hydrophilic amino acids may cause local perturbations and affect the electrostatic complementarity of RBD with the ACE2, and further mediate conformational changes. Thus, a more extensive interaction network was found in the mutation system and the complex interaction cluster was formed near E37@ACE2, which was essential for the stable binding of the two. In addition, we speculated that these mutations may affect the electrostatic complementarity with the four potential antibodies to reduce the sensitivity of the virus to antibodies. This study reveals the key details of the Omicron variant binding to ACE2 and provides important theoretical views for the enhanced infectivity of this variant. We hope that these observations can provide timely molecular insights for responding to the Omicron variant pandemic.

Caging Polycations: Effect of Increasing Confinement on the Modes of Interaction of Spermidine3+ With DNA Double Helices

Polyamines have important roles in the modulation of the cellular function and are ubiquitous in cells. The polyamines putrescine²⁺, spermidine³⁺, and spermine⁴⁺ represent the most abundant organic counterions of the negatively charged DNA in the cellular nucleus. These polyamines are known to stabilize the DNA structure and, depending on their concentration and additional salt composition, to induce DNA aggregation, which is often referred to as condensation. However, the modes of interactions of these elongated polycations with DNA and how they promote condensation are still not clear. In the present work, atomistic molecular dynamics (MD) computer simulations of two DNA fragments surrounded by spermidine³⁺ (Spd³⁺) cations were performed to study the structuring of Spd³⁺ “caged” between DNA molecules. Microsecond time scale simulations, in which the parallel DNA fragments were constrained at three different separations, but allowed to rotate axially and move naturally, provided information on the conformations and relative orientations of surrounding Spm³⁺ cations as a function of DNA-DNA separation. Novel geometric criteria allowed for the classification of DNA-Spd³⁺ interaction modes, with special attention given to Spd³⁺ conformational changes in the space between the two DNA molecules (caged Spd³⁺). This work shows how changes in the accessible space, or confinement, around DNA affect DNA-Spd³⁺ interactions, information fundamental to understanding the interactions between DNA and its counterions in environments where DNA is compacted, e.g. in the cellular nucleus.

Insights from computational analysis: how does the SARS-CoV-2 Delta (B.1.617.2) variant hijack ACE2 more effectively?

The SARS-CoV-2 Delta (B.1.617.2) variant was identified in India in October 2020, and it has quickly become the mainstream strain with strong toxicity and spread, posing great challenges to epidemic...

Interstrand Crosslinking Involving Guanine: A New Major UVC Laser-Induced Biphotonic Oxidatively Generated DNA Damage

Several classes of oxidatively generated DNA damage including oxidized purine and pyrimidine bases, interstrand base crosslinks and DNA-protein crosslinks have been previously shown to be generated in both isolated DNA and cellular DNA upon exposure to either 266 nm laser irradiation or one-electron oxidants. In this study, we provide evidence that biphotonic ionization of guanine bases by UVC laser irradiation of double-stranded deoxyoligonucleotides in aerated aqueous solutions induces the formation of interstrand cross-links (ICLs). This is supported by various experiments including sequencing gel analyses of formed photoproducts and effects of UVC laser intensity on their formation. This constitutes a novel example of the diversity of reactions of guanine radical cation that can be generated by various one-electron oxidants including UVC laser biphotonic ionization, direct effect of ionization radiation and type I photosensitizers. However, the exact structure of the interstrand base adducts that is a challenging analytical issue remains to be further established. Examples of relevant biochemical/structural applications of biphotonic induction of ICLs in DNA samples by high-intensity UVC laser pulses are provided.

Targeting mechanism for SARS-CoV-2 in silico: interaction and key groups of TMPRSS2 toward four potential drugs

The global dissemination of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) has seriously endangered human health. The number of confirmed cases is still increasing; however, treatment options are limited. Transmembrane protease serine 2 (TMPRSS2), as a key protease that primes the binding of SARS-CoV-2 spike protein and angiotensin-converting enzyme 2 (ACE2), has become an attractive target and received widespread attention. Thus, four potential drugs (bromhexine, camostat, gabexate, and nafamostat) were used to explore the mechanism of binding with TMPRSS2 in this work. A 65 ns molecular dynamics simulation was performed three times for each drug-TMPRSS2 system for reliable energy calculation and conformational analysis, of which the simulations of nafamostat-TMPRSS2 systems were further extended to 150 ns three times due to the discovery of two binding modes. Through the results of calculating binding free energy by nine methods, the binding affinity of camostat, gabexate, and nafamostat to TMPRSS2 showed great advantages compared with bromhexine, where the nafamostat was surprisingly found to present two reasonable binding conformations (forward and reverse directions). Two negatively charged amino acids (Asp435 and Glu299) can clamp the two positively charged groups (guanidinium group and amidinium group) in either forward or reverse fashion, and the forward one is more stable than the reverse. In addition, compared with gabexate, the dimethylamino group in camostat forms more van der Waals interactions with surrounding hot-spots His296 and Val280, resulting in a stronger affinity to TMPRSS2. For bromhexine, multiple binding sites are displayed in the binding pocket due to its small molecular structure, and van der Waals interactions play the dominant role in the binding process. In particular, six typical hot-spots were identified in the last three serine protease inhibitor systems, i.e., Asp435, Ser436, Gln438, Trp461, Ser463, and Gly464. The guanidinium groups of the drugs have powerful interactions with adjacent residues due to the formation of more hydrogen bonds, suggesting that this may be the critical site for drug design against TMPRSS2. This work provides valuable molecular insight into these four drug-TMPRSS2 binding mechanisms and is helpful for designing and screening drugs targeting TMPRSS2.

Photosensitization Reactions of Biomolecules: Definition, Targets and Mechanisms

Photosensitization reactions have been demonstrated to be largely responsible for the deleterious biological effects of UV and visible radiation, as well as for the curative actions of photomedicine. A large number of endogenous and exogenous photosensitizers, biological targets and mechanisms have been reported in the past few decades. Evolving from the original definitions of the type I and type II photosensitized oxidations, we now provide physicochemical frameworks, classifications and key examples of these mechanisms in order to organize, interpret and understand the vast information available in the literature and the new reports, which are in vigorous growth. This review surveys in an extended manner all identified photosensitization mechanisms of the major biomolecule groups such as nucleic acids, proteins, lipids bridging the gap with the subsequent biological processes. Also described are the effects of photosensitization in cells in which UVA and UVB irradiation triggers enzyme activation with the subsequent delayed generation of superoxide anion radical and nitric oxide. Definitions of photosensitized reactions are identified in biomolecules with key insights into cells and tissues.

Inhibition mechanism and hot-spot prediction of nine potential drugs for SARS-CoV-2 Mpro by large-scale molecular dynamic simulations combined with accurate binding free energy calculations

Coronavirus disease 2019 (COVID-19), which is caused by a new coronavirus known as severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), is spreading around the world. However, a universally effective treatment regimen has not been developed to date. The main protease (M^pro), a key enzyme of SARS-CoV-2, plays a crucial role in the replication and transcription of this virus in cells and has become the ideal target for rational antiviral drug design. In this study, we performed molecular dynamics simulations three times for these complexes of M^pro (monomeric and dimeric) and nine potential drugs that have a certain effect on the treatment of COVID-19 to explore their binding mechanism. In addition, a total of 12 methods for calculating binding free energy were employed to determine the optimal drug. Ritonavir, Arbidol, and Chloroquine consistently showed an outstanding binding ability to monomeric M^pro under various methods. Ritonavir, Arbidol, and Saquinavir presented the best performance when binding to a dimer, which was independent of the protonated state of Hie41 (protonated at N_ε) and Hid41 (protonated at N_δ), and these findings suggest that Chloroquine may not effectively inhibit the activity of dimeric M^proin vivo. Furthermore, three common hot-spot residues of Met165, Hie41, and Gln189 of monomeric M^pro systems dominated the binding of Ritonavir, Arbidol, and Chloroquine. In dimeric M^pro, Gln189, Met165, and Met49 contributed significantly to binding with Ritonavir, Arbidol, and Saquinavir; therefore, Gln189 and Met165 might serve as the focus in the discovery and development of anti-COVID-19 drugs. In addition, the van der Waals interaction played a significant role in the binding process, and the benzene ring of the drugs showed an apparent inhibitory effect on the normal function of M^pro. The binding cavity had great flexibility to accommodate these different drugs. The results would be notably helpful for enabling a detailed understanding of the binding mechanisms for these important drug-M^pro interactions and provide valuable guidance for the design of potent inhibitors.

Molecular Dynamics Approach for Capturing Calixarene-Protein Interactions: The Case of Cytochrome C

Functionalized supramolecular cages are of growing importance in biology and biochemistry. They have recently been proposed as efficient auxiliaries to obtain high-resolution cocrystallized proteins. Here, we propose a molecular dynamics investigation of the supramolecular association of sulfonated calix-[8]-arenes to cytochrome c starting from initially distant proteins and ligands. We characterize two main binding sites for the sulfonated calixarene on the cytochrome c surface which are in perfect agreement with the previous experiments with regard to the structure (comparison with the X-ray structure PDB 6GD8) and the binding free energies [comparison between the molecular mechanics Poisson-Boltzmann surface area analysis and the isothermal titration calorimetry measurements]. The per-residue decomposition of the interaction energies reveals the detailed picture of this electrostatically driven association and notably the role of arginine R13 as a bridging residue between the two main anchoring sites. In addition, the analysis of the residue behavior by means of a supervised machine learning protocol unveils the formation of a hydrogen bond network far from the binding sites, increasing the rigidity of the protein. This study paves the way toward an automated procedure to predict the supramolecular protein-cage association, with the possibility of a computational screening of new promising derivatives for controlled protein assembly and protein surface recognition processes.

Exogenous Polyamines Only Indirectly Induce Stress Tolerance in Wheat Growing in Hydroponic Culture under Polyethylene Glycol-Induced Osmotic Stress

The aim of the present study was to evaluate the effect of osmotic stress caused by polyethylene glycol (PEG) 6000 in hydroponic culture on wheat seedlings of drought-resistant Chinese Spring (CS) and drought-susceptible SQ1 cultivar, and to examine the alleviative role of exogenous polyamines (PAs) applied to the medium. The assessment was based on physiological (chlorophyll a fluorescence kinetics, chlorophyll and water content) as well as biochemical (content of carbohydrates, phenols, proline, salicylic and abscisic acid, activity of low molecular weight antioxidants) parameters, measured after supplementation with PAs (putrescine, spermidine and spermine) on the 3rd, 5th and 7th day of the treatment. The results indicate that PAs ameliorate the effects of stress, indirectly and conditionally inducing stress tolerance of wheat seedlings. In contrast to the susceptible SQ1, the resistant CS cultivar activated its protective mechanisms, adjusting the degree of their activation to the level of the stress, depending on the genetic resources of the plant. Increased accumulation of antioxidants in the resistant CS in response to stress after the application of PAs confirms the hypothesis that PAs are involved in the signaling pathway determining the antioxidative response and the tolerance of wheat plants to drought stress.

Phos-tag-based micropipette-tip method for analysis of phosphomonoester-type impurities in synthetic oligonucleotides

Various chromatographic techniques, combined with mass spectrometry, have been developed for the analysis of impurities in oligonucleotide drugs, but those methods have generally been less focused on possible phosphomonoester-type compounds. Here, we introduce a simple method for separating terminally phosphorylated impurities from parent oligonucleotides by using a phosphate-affinity micropipette tip (Phos-tag tip). All steps for the phosphate-affinity separation (binding, washing, and elution) are conducted in aqueous buffers at neutral pH. The entire separation protocol requires less than 30 min per sample. In practical examples, we demonstrated that phosphorylated impurities in natural-type and chemically modified oligonucleotides can be efficiently separated by the Phos-tag tip method and subsequently characterized by using ion-pairing reversed-phase liquid chromatography mass spectrometry (IP-RPLC-MS). Thus, a combination of the Phos-tag tip method and IP-RPLC-MS is useful for characterizing and identifying phosphomonoester-type impurities in oligonucleotide drugs.

Pattern preferences of DNA nucleotide motifs by polyamines putrescine2+, spermidine3+ and spermine4

The interactions of natural polyamines (putrescine2+, spermidine3+ and spermine4+) with DNA double helix are studied to characterize their nucleotide sequence pattern preference. Atomistic Molecular Dynamics simulations have been carried out for three systems consisting of the same DNA fragment d(CGCGAATTCGCGAATTCGCG) with different polyamines. The results show that polyamine molecules are localized with well-recognized patterns along the double helix with different residence times. We observed a clear hierarchy in the residence times of the polyamines, with the longest residence time (ca 100ns) in the minor groove. The analysis of the sequence dependence shows that polyamine molecules prefer the A-tract regions of the minor groove - in its narrowest part. The preferable localization of putrescine2+, spermidine3+ and spermine4+ in the minor groove with A-tract motifs is correlated with modulation of the groove width by a specific nucleotide sequences. We did develop a theoretical model pointing to the electrostatic interactions as the main driving force in this phenomenon, making it even more prominent for polyamines with higher charges. The results of the study explain the specificity of polyamine interactions with A-tract region of the DNA double helix which is also observed in experiments.

The effect of spermidine on guanine decomposition via photoinduced electron transfer in DNA

The addition of spermidine caused the attenuation of guanine decomposition via photoinduced electron transfer in pyrene-modified DNA, and higher added concentrations of spermidine resulted in the promotion of decomposition in condensed DNA.

Reaction Kinetics, Product Branching, and Potential Energy Surfaces of 1O2-Induced 9-Methylguanine-Lysine Cross-Linking: A Combined Mass Spectrometry, Spectroscopy, and Computational Study

We report a kinetics and mechanistic study on the ¹O₂ oxidation of 9-methylguanine (9MG) and the cross-linking of the oxidized intermediate 2-amino-9-methyl-9H-purine-6,8-dione (9MOG^OX) with N^α-acetyl-lysine-methyl ester (abbreviated as LysNH₂) in aqueous solutions of different pH. Experimental measurements include the determination of product branching ratios and reaction kinetics using mass spectrometry and absorption spectroscopy, and the characterization of product structures by employing collision-induced dissociation. Strong pH dependence was revealed for both 9MG oxidation and the addition of nucleophiles (water and LysNH₂) at the C5 position of 9MOG^OX. The ¹O₂ oxidation rate constant of 9MG was determined to be 3.6 × 10⁷ M^-1·s^-1 at pH 10.0 and 0.3 × 10⁷ M^-1·s^-1 at pH 7.0, both of which were measured in the presence of 15 mM LysNH₂. The ωB97XD density functional theory coupled with various basis sets and the SMD implicit solvation model was used to explore the reaction potential energy surfaces for the ¹O₂ oxidation of 9MG and the formation of C5-water and C5-LysNH₂ adducts of 9MOG^OX. Computational results have shed light on reaction pathways and product structures for the different ionization states of the reactants. The present work has confirmed that the initial ¹O₂ addition represents the rate-limiting step for the oxidative transformations of 9MG. All of the downstream steps are exothermic with respect to the starting reactants. The C5-cross-linking of 9MOG^OX with LysNH₂ significantly suppressed the formation of spiroiminodihydantoin (9MSp) resulting from the C5-water addition. The latter became dominant only at the low concentration (∼1 mM) of LysNH₂.

Specific effects of antitumor active norspermidine on the structure and function of DNA

We compared the effects of trivalent polyamines, spermidine (SPD) and norspermidine (NSPD), a chemical homologue of SPD, on the structure of DNA and gene expression. The chemical structures of SPD and NSPD are different only with the number of methylene groups between amine groups, [N-3-N-4-N] and [N-3-N-3-N], respectively. SPD plays vital roles in cell function and survival, including in mammals. On the other hand, NSPD has antitumor activity and is found in some species of plants, bacteria and algae, but not in humans. We found that both polyamines exhibit biphasic effect; enhancement and inhibition on in vitro gene expression, where SPD shows definitely higher potency in enhancement but NSPD causes stronger inhibition. Based on the results of AFM (atomic force microscopy) observations together with single DNA measurements with fluorescence microscopy, it becomes clear that SPD tends to align DNA orientation, whereas NSPD induces shrinkage with a greater potency. The measurement of binding equilibrium by NMR indicates that NSPD shows 4-5 times higher affinity to DNA than SPD. Our theoretical study with Monte Carlo simulation provides the insights into the underlying mechanism of the specific effect of NSPD on DNA.

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.

Counterion-Dependent Mechanisms of DNA Origami Nanostructure Stabilization Revealed by Atomistic Molecular Simulation

The DNA origami technique has proven to have tremendous potential for therapeutic and diagnostic applications like drug delivery, but the relatively low concentrations of cations in physiological fluids cause destabilization and degradation of DNA origami constructs preventing in vivo applications. To reveal the mechanisms behind DNA origami stabilization by cations, we performed atomistic molecular dynamics simulations of a DNA origami rectangle in aqueous solvent with varying concentrations of magnesium and sodium as well as polyamines like oligolysine and spermine. We explored the binding of these ions to DNA origami in detail and found that the mechanism of stabilization differs between ion types considerably. While sodium binds weakly and quickly exchanges with the solvent, magnesium and spermine bind close to the origami with spermine also located in between helices, stabilizing the crossovers characteristic for DNA origami and reducing repulsion of parallel helices. In contrast, oligolysine of length ten prevents helix repulsion by binding to adjacent helices with its flexible side chains, spanning the gap between the helices. Shorter oligolysine molecules with four subunits are weak stabilizers as they lack both the ability to connect helices and to prevent helix repulsion. This work thus shows how the binding modes of ions influence the stabilization of DNA origami nanostructures on a molecular level.

Tuning Protein Frameworks via Auxiliary Supramolecular Interactions

Protein crystals with their precise, periodic array of functional building blocks have potential applications in biomaterials, sensing, and catalysis. This paper describes how a highly porous crystalline framework of a cationic redox protein and an anionic macrocycle can be modulated by a small cationic effector. Ternary composites of protein (∼13 kDa), calix[8]arene (∼1.5 kDa), and effector (∼0.2 kDa) formed distinct crystalline architectures, dependent on the effector concentration and the crystallization technique. A combination of X-ray crystallography and density functional theory (DFT) calculations was used to decipher the framework variations, which appear to be dependent on a calixarene conformation change mediated by the effector. This "switch" calixarene was observed in three states, each of which is associated with a different interaction network. Two structures obtained by co-crystallization with the effector contained an additional protein "pillar", resulting in framework duplication and decreased porosity. These results suggest how protein assembly can be engineered by supramolecular host-guest interactions.

Cellular and Animal Model Studies on the Growth Inhibitory Effects of Polyamine Analogues on Breast Cancer

Polyamine levels are elevated in breast tumors compared to those of adjacent normal tissues. The female sex hormone, estrogen is implicated in the origin and progression of breast cancer. Estrogens stimulate and antiestrogens suppress the expression of polyamine biosynthetic enzyme, ornithine decarboxylate (ODC). Using several bis(ethyl)spermine analogues, we found that these analogues inhibited the proliferation of estrogen receptor-positive and estrogen receptor negative breast cancer cells in culture. There was structure-activity relationship in the efficacy of these compounds in suppressing cell growth. The activity of ODC was inhibited by these compounds, whereas the activity of the catabolizing enzyme, spermidine/spermine N¹-acetyl transferase (SSAT) was increased by 6-fold by bis(ethyl)norspermine in MCF-7 cells. In a transgenic mouse model of breast cancer, bis(ethyl)norspermine reduced the formation and growth of spontaneous mammary tumor. Recent studies indicate that induction of polyamine catabolic enzymes SSAT and spermine oxidase (SMO) play key roles in the anti-proliferative and apoptotic effects of polyamine analogues and their combinations with chemotherapeutic agents such as 5-fluorouracil (5-FU) and paclitaxel. Thus, polyamine catabolic enzymes might be important therapeutic targets and markers of sensitivity in utilizing polyamine analogues in combination with other therapeutic agents.