Investigation of the Stereochemical-Dependent DNA and RNA Binding of Arginine-Based Nucleopeptides

Nucleopeptides represent an intriguing class of nucleic acid analogues, in which nucleobases are placed in a peptide structure. The incorporation of Dand/or L-amino acids in nucleopeptide molecules allows the investigation of the role of backbone stereochemistry in determining the formation of DNA and RNA hybrids. Circular Dichroism (CD) spectroscopic studies indicated the nucleopeptide as having fully l-backbone configuration-formed stable hybrid complexes with RNA molecules. Molecular Dynamics (MD) simulations suggested a potential structure of the complex resulting from the interaction between the l-nucleopeptide and RNA strand. From this study, both the backbone (ionics and H-bonds) and nucleobases (pairing and π-stacking) of the chiral nucleopeptide appeared to be involved in the hybrid complex formation, highlighting the key role of the backbone stereochemistry in the formation of the nucleopeptide/RNA complexes.


Introduction
Nucleic acid-protein interactions are crucial in numerous cell functions, including transcription, translation and RNA maturation, which in turn modulate numerous physio-pathological processes [1,2].They regulate genetic information transfer by proper interaction between protein binding domains and specific regions of nucleic acids [3][4][5][6].The formation of highly stable and specific complexes depends on the strict recognition rules regulating the nucleobase association for nucleic acids and on a wide array of aminoacid-nucleobase interactions, including hydrophobic, ionic and hydrogen bonds for the nucleic acid-protein complexes.In the present day, the use of nucleobases constituting DNA and RNA aptamers is a well-established approach to interfere with both nucleic acid and protein targets for the diagnosis and the therapeutic treatment of numerous diseases, including cancer, inflammation and infections [7].In spite of their great potential, the full capitalization of these chemotypes as therapeutics is still limited by several drawbacks, such as low chemical variability, binding affinity and metabolic stability.In this regard, DNA aptamers containing pyrimidine bases endowed with amino acid residues (SOMAmer) [8,9] represent an enhancement in chemical diversity providing wide application as diagnostic and therapeutic tools.Alternatively, nucleobases have been added to peptides as recognition units by employing nucleobase amino acids (NBA) [10][11][12][13][14][15].To date, no study has been reported describing the use of NBA as elements to improve the selective protein recognition of target nucleic acids.Nevertheless, this basic concept has been proposed in several studies [16] involving peptide-based molecules, in order to make novel compounds that bind DNA and RNA specifically, such as Peptide Nucleic Acid (PNA) [17][18][19][20][21]. PNAs are nucleic acid analogues first described by Nielsen and coworkers, which have been largely applied for chemical biology investigations.They are able to efficiently mimic DNA molecules through a pseudopeptide backbone formed by achiral N-(2-aminoethyl) glycine (aeg) moieties bearing nucleobases linked via a methylene carbonyl bridge.Although PNAs selectively recognize DNA and RNA sequences affording very stable hybrids, the absence of chiral centers in its backbone hampers the folding in specific conformations, reducing its ability to discriminate highly structured DNA and RNA molecules.Additionally, PNAs exhibit numerous disadvantages, including low solubility in an aqueous medium, ambiguous directional selectivity of DNA/RNA binding, and low cell permeability.
To address some of the above-mentioned issues, several modifications have been performed on the original PNA structure [22][23][24].For instance, the incorporation of N-(2-aminoethyl)-D-Lysine residues led to the synthesis of cationic PNAs, which combined the nucleic acid recognition ability of PNA with the cell permeability of cationic cell penetrating peptides (CPP).These new PNA variants displayed advantages in terms of solubility and binding to DNA [23,24].Crystallographic studies demonstrated that the backbone chirality could conformationally modulate the PNA helical direction, thereby enhancing the recognition specificity of the target [25][26][27].In line with the results achieved by studying CPPs, the presence of the positive charge of Lys side-chains enhanced their solubility and especially their cell permeability, most likely through an endosomal pathway or direct cell membrane penetration [28].Similarly, new classes of chiral PNAs (α, and metry 2019, 11 FOR PEER REVIEW 2 ternatively, nucleobases have been added to peptides as recognition units by employing cleobase amino acids (NBA) [10][11][12][13][14][15].To date, no study has been reported describing the use of A as elements to improve the selective protein recognition of target nucleic acids.Nevertheless, s basic concept has been proposed in several studies [16] involving peptide-based molecules, in er to make novel compounds that bind DNA and RNA specifically, such as Peptide Nucleic Acid A) [17][18][19][20][21]. PNAs are nucleic acid analogues first described by Nielsen and coworkers, which ve been largely applied for chemical biology investigations.They are able to efficiently mimic DNA lecules through a pseudopeptide backbone formed by achiral N-(2-aminoethyl) glycine (aeg) ieties bearing nucleobases linked via a methylene carbonyl bridge.Although PNAs selectively ognize DNA and RNA sequences affording very stable hybrids, the absence of chiral centers in its ckbone hampers the folding in specific conformations, reducing its ability to discriminate highly uctured DNA and RNA molecules.Additionally, PNAs exhibit numerous disadvantages, luding low solubility in an aqueous medium, ambiguous directional selectivity of DNA/RNA ding, and low cell permeability.
To address some of the above-mentioned issues, several modifications have been performed on original PNA structure [22][23][24].For instance, the incorporation of N-(2-aminoethyl)-D-Lysine idues led to the synthesis of cationic PNAs, which combined the nucleic acid recognition ability of A with the cell permeability of cationic cell penetrating peptides (CPP).These new PNA variants played advantages in terms of solubility and binding to DNA [23,24].Crystallographic studies monstrated that the backbone chirality could conformationally modulate the PNA helical ection, thereby enhancing the recognition specificity of the target [25][26][27].In line with the results ieved by studying CPPs, the presence of the positive charge of Lys side-chains enhanced their ubility and especially their cell permeability, most likely through an endosomal pathway or direct l membrane penetration [28].Similarly, new classes of chiral PNAs (α, and ɣ-GPNA) have been ently designed by merging the recognition ability of PNA with the uptake features of cationic anidinium groups [29,30].These studies showed that GPNAs with α-D and/or ɣ-L-backbone nfiguration exhibited improved cellular uptake, increased binding affinity for DNA and RNA, and hanced thermal stability of the resulting hybrids.These favorable properties were ascribable to the nformational preorganization of chiral PNAs and the electrostatic interactions between cationic anidine and phosphate groups.In particular, they demonstrated that the backbone reochemistry contributed more with respect to the electrostatic forces in stabilizing the hybrid plex.
In this regard, nucleopeptides represent alternative chemotypes in which nucleobases are mbined with a peptide structure [31][32][33].This fascinating class of molecules are characterized by presence of amino acids-bearing nucleobases as side chains that promote nucleic acid recognition ough both base pairing and stacking interactions.In nucleopeptides, the nucleobasectionalized amino acids are alternated with underivatized amino acids at defined positions in the ptide backbone.Diederichsen et al. [34] and Mihara et al. [35,36,15] reported the synthesis of α and eptides with nucleobase-functionalized side chains and the use of L-α-amino-ɣ-nucleobasetyric acid to assemble RNA binding molecules and functional peptides, respectively.More ently, we have reported new synthetic strategies to obtain both homo-and hetero-sequences of ionic nucleopeptides that exhibited advantageous properties in terms of binding, specificity, cell rmeability and delivery for RNA, DNA and PNA molecules (Figure 1) [37,38].Although the rapeutic potential of nucleopeptides is less known with the respect to the parent PNA molecules, y are under investigation for their ability to interfere with the Reverse Transcriptase of the Human munodeficiency Virus (RT-HIV), which is fundamental for HIV infection [39].Additionally, cleopeptides have been successfully employed in combination with Doxorubicin, to improve its icacy in Multidrug Resistant Cancer (MRD) by sequestering adenosine triphosphate (ATP) [40].
-GPNA) have been recently designed by merging the recognition ability of PNA with the uptake features of cationic guanidinium groups [29,30].These studies showed that GPNAs with α-d and/or 9, 11 FOR PEER REVIEW 2 ly, nucleobases have been added to peptides as recognition units by employing amino acids (NBA) [10][11][12][13][14][15].To date, no study has been reported describing the use of ments to improve the selective protein recognition of target nucleic acids.Nevertheless, oncept has been proposed in several studies [16] involving peptide-based molecules, in ake novel compounds that bind DNA and RNA specifically, such as Peptide Nucleic Acid 21].PNAs are nucleic acid analogues first described by Nielsen and coworkers, which largely applied for chemical biology investigations.They are able to efficiently mimic DNA through a pseudopeptide backbone formed by achiral N-(2-aminoethyl) glycine (aeg) earing nucleobases linked via a methylene carbonyl bridge.Although PNAs selectively NA and RNA sequences affording very stable hybrids, the absence of chiral centers in its ampers the folding in specific conformations, reducing its ability to discriminate highly DNA and RNA molecules.Additionally, PNAs exhibit numerous disadvantages, low solubility in an aqueous medium, ambiguous directional selectivity of DNA/RNA d low cell permeability.dress some of the above-mentioned issues, several modifications have been performed on l PNA structure [22][23][24].For instance, the incorporation of N-(2-aminoethyl)-D-Lysine d to the synthesis of cationic PNAs, which combined the nucleic acid recognition ability of the cell permeability of cationic cell penetrating peptides (CPP).These new PNA variants advantages in terms of solubility and binding to DNA [23,24].Crystallographic studies ted that the backbone chirality could conformationally modulate the PNA helical hereby enhancing the recognition specificity of the target [25][26][27].In line with the results y studying CPPs, the presence of the positive charge of Lys side-chains enhanced their nd especially their cell permeability, most likely through an endosomal pathway or direct rane penetration [28].Similarly, new classes of chiral PNAs (α, and ɣ-GPNA) have been signed by merging the recognition ability of PNA with the uptake features of cationic m groups [29,30].These studies showed that GPNAs with α-D and/or ɣ-L-backbone on exhibited improved cellular uptake, increased binding affinity for DNA and RNA, and hermal stability of the resulting hybrids.These favorable properties were ascribable to the onal preorganization of chiral PNAs and the electrostatic interactions between cationic and phosphate groups.In particular, they demonstrated that the backbone istry contributed more with respect to the electrostatic forces in stabilizing the hybrid s regard, nucleopeptides represent alternative chemotypes in which nucleobases are with a peptide structure [31][32][33].This fascinating class of molecules are characterized by e of amino acids-bearing nucleobases as side chains that promote nucleic acid recognition oth base pairing and stacking interactions.In nucleopeptides, the nucleobasezed amino acids are alternated with underivatized amino acids at defined positions in the ckbone.Diederichsen et al. [34] and Mihara et al. [35,36,15] reported the synthesis of α and with nucleobase-functionalized side chains and the use of L-α-amino-ɣ-nucleobaseid to assemble RNA binding molecules and functional peptides, respectively.More e have reported new synthetic strategies to obtain both homo-and hetero-sequences of cleopeptides that exhibited advantageous properties in terms of binding, specificity, cell ty and delivery for RNA, DNA and PNA molecules (Figure 1) [37,38].Although the c potential of nucleopeptides is less known with the respect to the parent PNA molecules, der investigation for their ability to interfere with the Reverse Transcriptase of the Human ficiency Virus (RT-HIV), which is fundamental for HIV infection [39].Additionally, tides have been successfully employed in combination with Doxorubicin, to improve its Multidrug Resistant Cancer (MRD) by sequestering adenosine triphosphate (ATP) [40].
-l-backbone configuration exhibited improved cellular uptake, increased binding affinity for DNA and RNA, and enhanced thermal stability of the resulting hybrids.These favorable properties were ascribable to the conformational preorganization of chiral PNAs and the electrostatic interactions between cationic guanidine and phosphate groups.In particular, they demonstrated that the backbone stereochemistry contributed more with respect to the electrostatic forces in stabilizing the hybrid duplex.
In this regard, nucleopeptides represent alternative chemotypes in which nucleobases are combined with a peptide structure [31][32][33].This fascinating class of molecules are characterized by the presence of amino acids-bearing nucleobases as side chains that promote nucleic acid recognition through both base pairing and stacking interactions.In nucleopeptides, the nucleobase-functionalized amino acids are alternated with underivatized amino acids at defined positions in the peptide backbone.Diederichsen et al. [34] and Mihara et al. [15,35,36] reported the synthesis of α and β-peptides with nucleobase-functionalized side chains and the use of l-α-amino-Symmetry 2019, 11 FOR PEER REVIEW 2 Alternatively, nucleobases have been added to peptides as recognition units by employing nucleobase amino acids (NBA) [10][11][12][13][14][15].To date, no study has been reported describing the use of NBA as elements to improve the selective protein recognition of target nucleic acids.Nevertheless, this basic concept has been proposed in several studies [16] involving peptide-based molecules, in order to make novel compounds that bind DNA and RNA specifically, such as Peptide Nucleic Acid (PNA) [17][18][19][20][21]. PNAs are nucleic acid analogues first described by Nielsen and coworkers, which have been largely applied for chemical biology investigations.They are able to efficiently mimic DNA molecules through a pseudopeptide backbone formed by achiral N-(2-aminoethyl) glycine (aeg) moieties bearing nucleobases linked via a methylene carbonyl bridge.Although PNAs selectively recognize DNA and RNA sequences affording very stable hybrids, the absence of chiral centers in its backbone hampers the folding in specific conformations, reducing its ability to discriminate highly structured DNA and RNA molecules.Additionally, PNAs exhibit numerous disadvantages, including low solubility in an aqueous medium, ambiguous directional selectivity of DNA/RNA binding, and low cell permeability.
To address some of the above-mentioned issues, several modifications have been performed on the original PNA structure [22][23][24].For instance, the incorporation of N-(2-aminoethyl)-D-Lysine residues led to the synthesis of cationic PNAs, which combined the nucleic acid recognition ability of PNA with the cell permeability of cationic cell penetrating peptides (CPP).These new PNA variants displayed advantages in terms of solubility and binding to DNA [23,24].Crystallographic studies demonstrated that the backbone chirality could conformationally modulate the PNA helical direction, thereby enhancing the recognition specificity of the target [25][26][27].In line with the results achieved by studying CPPs, the presence of the positive charge of Lys side-chains enhanced their solubility and especially their cell permeability, most likely through an endosomal pathway or direct cell membrane penetration [28].Similarly, new classes of chiral PNAs (α, and ɣ-GPNA) have been recently designed by merging the recognition ability of PNA with the uptake features of cationic guanidinium groups [29,30].These studies showed that GPNAs with α-D and/or ɣ-L-backbone configuration exhibited improved cellular uptake, increased binding affinity for DNA and RNA, and enhanced thermal stability of the resulting hybrids.These favorable properties were ascribable to the conformational preorganization of chiral PNAs and the electrostatic interactions between cationic guanidine and phosphate groups.In particular, they demonstrated that the backbone stereochemistry contributed more with respect to the electrostatic forces in stabilizing the hybrid duplex.
In this regard, nucleopeptides represent alternative chemotypes in which nucleobases are combined with a peptide structure [31][32][33].This fascinating class of molecules are characterized by the presence of amino acids-bearing nucleobases as side chains that promote nucleic acid recognition through both base pairing and stacking interactions.In nucleopeptides, the nucleobasefunctionalized amino acids are alternated with underivatized amino acids at defined positions in the peptide backbone.Diederichsen et al. [34] and Mihara et al. [35,36,15] reported the synthesis of α and β-peptides with nucleobase-functionalized side chains and the use of L-α-amino-ɣ-nucleobasebutyric acid to assemble RNA binding molecules and functional peptides, respectively.More recently, we have reported new synthetic strategies to obtain both homo-and hetero-sequences of cationic nucleopeptides that exhibited advantageous properties in terms of binding, specificity, cell permeability and delivery for RNA, DNA and PNA molecules (Figure 1) [37,38].Although the therapeutic potential of nucleopeptides is less known with the respect to the parent PNA molecules, they are under investigation for their ability to interfere with the Reverse Transcriptase of the Human -nucleobase-butyric acid to assemble RNA binding molecules and functional peptides, respectively.More recently, we have reported new synthetic strategies to obtain both homo-and hetero-sequences of cationic nucleopeptides that exhibited advantageous properties in terms of binding, specificity, cell permeability and delivery for RNA, DNA and PNA molecules (Figure 1) [37,38].Although the therapeutic potential of nucleopeptides is less known with the respect to the parent PNA molecules, they are under investigation for their ability to interfere with the Reverse Transcriptase of the Human Immunodeficiency Virus (RT-HIV), which is fundamental for HIV infection [39].Additionally, nucleopeptides have been successfully employed in combination with Doxorubicin, to improve its efficacy in Multidrug Resistant Cancer (MRD) by sequestering adenosine triphosphate (ATP) [40].Herein, we describe the use of our synthetic approach to introduce specific structural modifications, such as the incorporation of D-and/or L-amino acids, in nucleopeptides with the aim to investigate the role of backbone stereochemistry in determining the formation of DNA and RNA hybrids.CD spectroscopy and UV-thermal stability investigations provided evidence of the impact of the backbone stereochemistry on the ability of nucleopeptides to hybridize DNA and RNA.Finally, the interaction of the best performing chiral nucleopeptide was studied by molecular dynamics (MD) simulations to achieve a new insight on the structure of the resulting complexes with RNA complementary strands.

Synthesis of Nucleopeptides np-1/4
Rink Amide AM-PS resin (60 mg, 0.033 mmol, 0.55 mmol/g) was suspended for 30 min in DCM/DMF 1:1 and then washed with DCM (3 × 1 min) and DMF (3 × 1 min).The Fmoc deprotection was carried out by treatment with a solution of 20% (v/v) piperidine in DMF (1 × 5 min; 1 × 25 min).The succeeding amide couplings were performed using 4 equivalents (according to the original loading of the resin) of the Fmoc-protected L or D-aminoacids.To a DMF (1.2 mL) solution of Fmoc-L (or D)-Arg(Pbf)-OH (84,5 mg, 0.13 mmol, 4 equivalents), HBTU (50 mg, 0.13 mmol, 4 equivalents) and HOBt (19 mg, 0.13 mmol, 4 equivalents) DIPEA (45 µ L, 0.26 mmol, 8 equivalents) was added and Herein, we describe the use of our synthetic approach to introduce specific structural modifications, such as the incorporation of d-and/or l-amino acids, in nucleopeptides with the aim to investigate the role of backbone stereochemistry in determining the formation of DNA and RNA hybrids.CD spectroscopy and UV-thermal stability investigations provided evidence of the impact of the backbone stereochemistry on the ability of nucleopeptides to hybridize DNA and RNA.Finally, the interaction of the best performing chiral nucleopeptide was studied by molecular dynamics (MD) simulations to achieve a new insight on the structure of the resulting complexes with RNA complementary strands.

Synthesis of Nucleopeptides np-1/4
Rink Amide AM-PS resin (60 mg, 0.033 mmol, 0.55 mmol/g) was suspended for 30 min in DCM/DMF 1:1 and then washed with DCM (3 × 1 min) and DMF (3 × 1 min).The Fmoc deprotection was carried out by treatment with a solution of 20% (v/v) piperidine in DMF (1 × 5 min; 1 × 25 min).The succeeding amide couplings were performed using 4 equivalents (according to the original loading of the resin) of the Fmoc-protected l or d-aminoacids.To a DMF (1.2 mL) solution of Fmoc-l (or d)-Arg(Pbf)-OH (84,5 mg, 0.13 mmol, 4 equivalents), HBTU (50 mg, 0.13 mmol, 4 equivalents) and HOBt (19 mg, 0.13 mmol, 4 equivalents) DIPEA (45 µL, 0.26 mmol, 8 equivalents) was added and the resulting mixture was added to the support.The mixture was shaken for 2 h at room temperature and then washed with DMF (3 × 1 min) and DCM (3 × 1 min).Next, the unreacted free amine groups were acetylated by treatment of the resin with 1 mL of DMF capping solution containing of Ac 2 O and of DIPEA (2 and 3 equivalents relative to the initial loading of the resin, respectively) and shaken for 5 min.The completion of the peptide sequences was achieved by repeating the cycles of the Fmoc deprotection and coupling reaction with Fmoc-l (or d)-Lys(Alloc)-OH and Fmoc-l (or d)-Arg(Pbf)-OH as above reported.These steps were monitored by Kaiser ninhydrine and TNBS tests and ascertained by quantitative UV spectroscopic measurements of Fmoc chromophore.After the last Fmoc removal, the treatment with a mixture of Ac 2 O (6.2 µL, 2 equiv) and DIPEA (17.3 µL, 3 equivalents) in DMF (1 mL) for 10 min, allowed N-terminal acetylation of the peptides.
The Alloc groups were removed by treating the supported full-protected peptide twice with a solution of Tetrakis(triphenylphosphine)palladium (0) (11 mg, 0.01 mmol, 5% mol relative to each Alloc group) and DMBA (123 mg, 0.0.79 mmol, 4 equivalents relative to each Alloc protective group) in dry DMF/DCM (1:1) for 1 h.The resulting solution was drained off and the resin was washed (3 × 1 min) with DMF and with DCM.The catalyst traces were eliminated by treating the support with a DMF solution of potassium N,N-diethylcarbamodithioate (0.06 M, 1.5 mL) for 1 h.To obtain the thyminyl-functionalized nucleopeptides, Thymin-1-yl-acetic acid (108,1 mg, 0.59 mmol, 3 equivalents for each amine group), HOAt (87 mg, 0.59 mmol) and PyAOP (308 mg, 0.59 mmol), DIPEA (205 µL, 1.08 mmol, 6 equivalents for each amine group) were dissolved in DMF (2 mL) and then added to resin for reacting over 8 h at room temperature.This latter step was repeated once.

Cleavage, HPLC Purification and Characterization of Nucleopeptides np-1/4
After washing with DMF (3 × 1 min), DCM (3 × 1 min), and Et 2 O (3 × 1 min) the supported nucleopeptides np-1/4 were then dried exhaustively.The np-1/4 were released from the solid support by 3 h treatment with TFA/TIS (95:5, 1.5 mL) at room temperature.Then, the resin was filtered off and the crude nucleopeptides were precipitated from cold Et 2 O (15 mL) and then centrifuged (6000 rpm × 15 min).After removal of supernatant, the precipitate was suspended once in Et 2 O, centrifuged and the supernatant removed.The resulting wet solid was dried under vacuum, dissolved in water/ACN (9:1) and purified by RP-HPLC (solvent A: water + 0.1% TFA; solvent B: acetonitrile + 0.1% TFA; from 10 to 70% of solvent B over 25 min, flow rate: 10 mL min −1 ).The collected fractions were evaporated from organic solvents under reduced pressure, frozen and then lyophilized.The final products were analyzed by analytical RP-HPLC (10-90% acetonitrile in water (0.1% TFA) over 20 min, flow rate of 1.0 mL/min) and characterized by MALDI-TOF spectrometry (Table 1).Yield, purity, retention times, and analytical data are reported in Table 1.

UV and CD procedures
UV measurements were performed by a spectrophotometer equipped with a Peltier temperature controller.Nucleopeptide stock solutions were prepared using doubly distilled water and RNAse free water.Their concentration was calculated by UV measurements at 260 nm at 90 • C by using the following molar extinction coefficient: 4 to the support.The mixture was shaken for 2 h at room temperature 1 min) and DCM (3 × 1 min).Next, the unreacted free amine groups he resin with 1 mL of DMF capping solution containing of Ac2O and lative to the initial loading of the resin, respectively) and shaken for ptide sequences was achieved by repeating the cycles of the Fmoc ion with Fmoc-L (or D)-Lys(Alloc)-OH and Fmoc-L (or D)-Arg(Pbf)steps were monitored by Kaiser ninhydrine and TNBS tests and spectroscopic measurements of Fmoc chromophore.After the last ith a mixture of Ac2O (6.2 μL, 2 equiv) and DIPEA (17.3 μL, 3 0 min, allowed N-terminal acetylation of the peptides.oved by treating the supported full-protected peptide twice with a osphine)palladium (0) (11  formed by a spectrophotometer equipped with a Peltier temperature solutions were prepared using doubly distilled water and RNAse as calculated by UV measurements at 260 nm at 90 °C by using the icient: Ɛ 260 = 13.7 mL/(μmol × cm) for A, and Ɛ 260 = 8.8 mL/(μmol ing curves were acquired in 1 cm path-length quartz cell, on a Jasco ometer equipped with a Peltier temperature programmer.Samples xes were formed by mixing equimolar concentrations of np-1/4 (10 260 = 13.7 mL/(µmol × cm) for A, and Symmetry 2019, 11 FOR PEER REVIEW 4 the resulting mixture was added to the support.The mixture was shaken for 2 h at room temperature and then washed with DMF (3 × 1 min) and DCM (3 × 1 min).Next, the unreacted free amine groups were acetylated by treatment of the resin with 1 mL of DMF capping solution containing of Ac2O and of DIPEA (2 and 3 equivalents relative to the initial loading of the resin, respectively) and shaken for 5 min.The completion of the peptide sequences was achieved by repeating the cycles of the Fmoc deprotection and coupling reaction with Fmoc-L (or D)-Lys(Alloc)-OH and Fmoc-L (or D)-Arg(Pbf)-OH as above reported.These steps were monitored by Kaiser ninhydrine and TNBS tests and ascertained by quantitative UV spectroscopic measurements of Fmoc chromophore.After the last Fmoc removal, the treatment with a mixture of Ac2O (6.2 μL, 2 equiv) and DIPEA (17.3 μL, 3 equivalents) in DMF (1 mL) for 10 min, allowed N-terminal acetylation of the peptides.The Alloc groups were removed by treating the supported full-protected peptide twice with a solution of Tetrakis(triphenylphosphine)palladium (0) (11 mg, 0.01 mmol, 5% mol relative to each Alloc group) and DMBA (123 mg, 0.0.79 mmol, 4 equivalents relative to each Alloc protective group) in dry DMF/DCM (1:1) for 1 h.The resulting solution was drained off and the resin was washed (3 × 1 min) with DMF and with DCM.The catalyst traces were eliminated by treating the support with a DMF solution of potassium N,N-diethylcarbamodithioate (0.06 M , 1.5 mL) for 1 h.To obtain the thyminyl-functionalized nucleopeptides, Thymin-1-yl-acetic acid (108,1 mg, 0.59 mmol, 3 equivalents for each amine group), HOAt (87 mg, 0.59 mmol) and PyAOP (308 mg, 0.59 mmol), DIPEA (205 µ L, 1.08 mmol, 6 equivalents for each amine group) were dissolved in DMF (2 mL) and then added to resin for reacting over 8 h at room temperature.This latter step was repeated once.

UV and CD procedures
UV measurements were performed by a spectrophotometer equipped with a Peltier temperature controller.Nucleopeptide stock solutions were prepared using doubly distilled water and RNAse free water.Their concentration was calculated by UV measurements at 260 nm at 90 °C by using the following molar extinction coefficient: Ɛ 260 = 13.7 mL/(μmol × cm) for A, and Ɛ 260 = 8.8 mL/(μmol × cm) for T. CD spectra and melting curves were acquired in 1 cm path-length quartz cell, on a Jasco J-810 Circular Dichroism Spectrometer equipped with a Peltier temperature programmer.Samples were lyophilized for 16 h.Duplexes were formed by mixing equimolar concentrations of np-1/4 (10 260 = 8.8 mL/(µmol × cm) for T. CD spectra and melting curves were acquired in 1 cm path-length quartz cell, on a Jasco J-810 Circular Dichroism Spectrometer equipped with a Peltier temperature programmer.Samples were lyophilized for 16 h.Duplexes were formed by mixing equimolar concentrations of np-1/4 (10 µM each base) and single stranded complementary polyA (M.W. average 100-500 kDa) and polydA (M.W. average 85-170 kDa) (10 µM each base) in a 10 mM phosphate buffer, 100 mM NaCl, pH 7.0 solution.The samples were annealed by incubating at 90 • C for 5 min, slowly cooling down at room temperature and then kept at 4 • C overnight.UV and CD melting was monitored at λ = 260 nm (for UV) and at at λ = 267 nm for the polyA containing hybrids and at λ = 247 nm for polydA containing hybrids, with a scan rate of 1 • C min −1 in the 5-85 • C temperature range.Melting temperatures (Tm) were extrapolated by the first derivatives of CD melting curves.CD spectra were performed in the 210-300 nm spectral range, keeping the temperature constant at 20 • C.

CD Measurements of Apparent Kd Determination of np-1/polyA and np1/polydA
The titration experiment was carried out by the addition of increasing concentrations of np-1 (0, 1.87, 3.76, 5.60, 7.47, 9.33, 11.20, 13.07, 14.93, 16.8, 18.67 10 µM each base) to a solution containing the single stranded complementary polyA (10 µM each base) and polydA (10 µM each base) in a 10 mM phosphate buffer, 100 mM NaCl, pH 7.0 solution.CD spectra were acquired at 20 • C, monitoring the signal decrements at λ = 260 nm for the polyA and at λ = 216 nm for polydA.The detected changes in CD intensity were plotted to generate saturation isotherm curves as a function of added np-1 amounts.Assuming the formation of 1:1 complex for each nucleobase, we calculated the apparent Kd values by using the following equation: where Θ is the CD intensity of the characteristic hybrid wavelength, Θmin and Θmax are the minimum and maximum value of CD intensity, respectively, and cA and cT are the Adenine and Thymine nucleobase concentrations, respectively.

Building of np-1/RNA Complex
Maestro [41,42] was used to build cp-1, the complex structure, using a double-stranded B-DNA hexamer of complementary A-T bases as a template.The phosphofuranose backbone chain of the Thymine-containing single-strand was removed and a nucleopeptide chain, consisting of 6 L-Arg interpolated with six L-Lys residues (np-1), was constructed in its place in an antiparallel fashion, while keeping the A-T base paring fixed.Finally, hydroxyl groups with R configuration were added to the 2' carbon atoms of the remaining furanose rings, and the C-and N-terminal ends of np-1 was capped with an amino and an acetyl group respectively.The initial structure of cp-1 was constructed base on the experimentally determined 2:1 stoichiometric interaction of two nucleopeptide molecules with a single 12 mer RNA oligonucleotide [37] To finally assemble the structure of cp-1 composed of a single-stranded 12 mer Adenine RNA oligonucleotide (A12 RNA ss) and two molecules of nucleopeptide np-1, the previously obtained cp-1 complex was duplicated and the RNA chains were linked covalently.The union was made by setting one copy of cp-1 on top of the other and connecting the 3' end of one RNA strands with the 5' end of the other through a newly added phosphate group.

Molecular Dynamics Simulation
Cp-1 complex was initially minimized with MacroModel using an OPLS 2005 force field in water and a Polak-Ribiere Conjugate Gradient minimization protocol of 100,000 maximum iterations and a convergence threshold of 0.001.A constrain in the base paring distances of A-T residues was used to fix the double-strand conformation during the process.When satisfactory minimized structures were obtained (i.e.no anomalous bond angles, amide dihedral angles, bond distances, etc.), molecular dynamic (MD) simulations were launched using Desmond as implemented in the Schrodinger suite [43][44][45].The model system was created using the System Builder tool of Desmond with an explicit TIP3P water solvent model.The positive charge of the complex was neutralized by 1 chloride ion, and additional ions up to a total concentration of 0.15 M KCl were added, all embedded into a rectangular box of 15 Å in each direction.Relaxation of the model system before performing the MD simulation was carried out with the default Desmond parameters, followed by a 200 ns unrestrained MD simulation with an NPT ensemble (300 K, 1.013 bar).An integration time-step of 2 fs was used for the multistep protocol with an interpolating function for electrostatic interactions (cut-off radius of 9.0 Å).Steady temperature and pressure values were maintained using the Nose-Hoover chain thermostat [45].

RMSD Clustering Analysis
The clustering analysis was carried out using the MD movie Tool in UCSF Chimera 1.10.1 (University of California, San Francisco, CA, USA) [45] as follows: the final MD simulation frames were exported with Maestro (Schrödinger, Cambridge, MA, USA) to separate pdb files and opened with UCSF Chimera MD movie Tool.A cluster analysis was next carried out, setting frame 1 as the starting one, and a step size of 150.The remaining options were set as default.The final cluster-representative structures given by Chimera were minimized with MacroModel (Schrödinger, Cambridge, MA, USA) with an OPLS 2005 force field in a Polak-Ribier Conjugate Gradient method for a maximum of 5000 iterations and a convergence threshold of 0.05 in water.

Synthesis
The synthesis of the nucleopeptides set was performed by a general synthetic path, that was developed by using Allyloxycarbonyl [46][47][48] (Alloc) as temporary protection strategy for the amino group connecting the nucleobase.The Alloc approach allowed using commercially available building blocks avoiding the solution phase synthetic steps required for the preparation of the N Alloc group) and DMBA (123 mg, 0.0.79 mmol, 4 equivalents relative to each Alloc protective group) in dry DMF/DCM (1:1) for 1 h.The resulting solution was drained off and the resin was washed (3 × 1 min) with DMF and with DCM.The catalyst traces were eliminated by treating the support with a DMF solution of potassium N,N-diethylcarbamodithioate (0.06 M , 1.5 mL) for 1 h.To obtain the thyminyl-functionalized nucleopeptides, Thymin-1-yl-acetic acid (108,1 mg, 0.59 mmol, 3 equivalents for each amine group), HOAt (87 mg, 0.59 mmol) and PyAOP (308 mg, 0.59 mmol), DIPEA (205 µ L, 1.08 mmol, 6 equivalents for each amine group) were dissolved in DMF (2 mL) and then added to resin for reacting over 8 h at room temperature.This latter step was repeated once.

UV and CD procedures
UV measurements were performed by a spectrophotometer equipped with a Peltier temperature controller.Nucleopeptide stock solutions were prepared using doubly distilled water and RNAse free water.Their concentration was calculated by UV measurements at 260 nm at 90 °C by using the following molar extinction coefficient: Ɛ 260 = 13.7 mL/(μmol × cm) for A, and Ɛ 260 = 8.8 mL/(μmol × cm) for T. CD spectra and melting curves were acquired in 1 cm path-length quartz cell, on a Jasco J-810 Circular Dichroism Spectrometer equipped with a Peltier temperature programmer.Samples were lyophilized for 16 h.Duplexes were formed by mixing equimolar concentrations of np-1/4 (10 -protected derivatives or nucleoamino acids needed for the subsequent assembly of nucleopeptide sequences on the solid support.The use of Alloc protecting group is widely employed to preserve amino acid side chains, including Lysine and Ornitine, those are required to be selectively released for direct functionalization on solid support (i.e.FITC, fatty acid etc) or lactam bridged cyclization [45,46].Interestingly, in spite of the application of an Alloc-protecting strategy for single amino acid being consolidated, few investigations have been reported for the protection of multiple sites.Thus, we aimed to explore its use also for oligomers containing numerous amino groups, those can be simultaneously released for further on-line derivatization, including the introduction of nucleobases in the case of nucleopeptides.As described in Scheme 1, the 12 mer amino acidic backbone of the nucleopeptides consisting of alternation of l-Arg and l-Lys (np-1) or d-Arg and l-Lys (np-2) or l-Arg and d-Lys (np-3) or d-Arg and d-Lys (np-4) was assembled using a Rink-amide-aminomethyl polystyrene (AM-PS) resin.After acetylation of the N-terminal function, all the Lys side chains were concurrently deprotected by a DMF/DCM solution of Pd(PPh 3 ) 4 and N,N-Dimethylbarbituric acid (DMBA) as scavenger.The treatment was reiterated three times to assure the complete removal of Alloc protecting groups as confirmed by HPLC-MS analyses.Next, the Thymine-1-acetic acids (T-CH 2 COOH) were placed by one-pot reaction by a 3 h reaction in the presence of HBTU/HOBt as activating/additive agents.All the nucleopeptides (np-1/4) were eventually released from the resin, purified by RP-HPLC and characterized by MALDI-TOF MS.

CD spectroscopic studies
CD studies were carried out with np-1/4 and with the hybrids obtained by combining the synthesized nucleopeptides np-1/4 with polyA and polydA ss targets.All the CD analyses were performed using 10 μM concentration of np-1/4 in 10 mM phosphate buffer (pH = 7.0) at 20 °C.As depicted in Figure 2, all the single stranded nucleopeptides (np-1/4) exhibited the distinctive CD spectra of linear cationic peptides in an unordered state endowed with a single band below 200 nm.

CD Spectroscopic Studies
CD studies were carried out with np-1/4 and with the hybrids obtained by combining the synthesized nucleopeptides np-1/4 with polyA and polydA ss targets.All the CD analyses were performed using 10 µM concentration of np-1/4 in 10 mM phosphate buffer (pH = 7.0) at 20 • C. As depicted in Figure 2, all the single stranded nucleopeptides (np-1/4) exhibited the distinctive CD spectra of linear cationic peptides in an unordered state endowed with a single band below 200 nm.

CD spectroscopic studies
CD studies were carried out with np-1/4 and with the hybrids obtained by combining the synthesized nucleopeptides np-1/4 with polyA and polydA ss targets.All the CD analyses were performed using 10 μM concentration of np-1/4 in 10 mM phosphate buffer (pH = 7.0) at 20 °C.As depicted in Figure 2, all the single stranded nucleopeptides (np-1/4) exhibited the distinctive CD spectra of linear cationic peptides in an unordered state endowed with a single band below 200 nm.As expected, enantiomeric nucleopeptides containing l-Arg/l-Lys and d-Arg/d-Lys (np-1 and -4) and d-Arg/l-Lys and l-Arg/d-Lys (np-2 and -3) exhibited mirror CD spectra.Next, we investigated the spectroscopic behavior of the hybrids composed of nucleopeptides np-1/4 with PolyA and PolydA (1:1 ratio in nucleobases).Long oligomers, such as polyA and polydA, were chosen for the CD experiments to amplify the resulting binding effects in CD spectra and to take advantage from "the flanking base effect", that allows the nucleobases next to 5' and 3' termini to stabilize the conformation of inner bases even if they do not contribute to the base pairing.Compared to polyA and polydA ss, the CD spectra of np-1/4/polyA (Figure 3a) and np-1/4/polydA (Figure 3b) hybrids exhibited hypochromic effects that could be ascribable to the cationic character of the side chain of Arginine residues.
Symmetry 2019, 11 FOR PEER REVIEW 8 As expected, enantiomeric nucleopeptides containing L-Arg/L-Lys and D-Arg/D-Lys (np-1 and -4) and D-Arg/L-Lys and L-Arg/D-Lys (np-2 and -3) exhibited mirror CD spectra.Next, we investigated the spectroscopic behavior of the hybrids composed of nucleopeptides np-1/4 with PolyA and PolydA (1:1 ratio in nucleobases).Long oligomers, such as polyA and polydA, were chosen for the CD experiments to amplify the resulting binding effects in CD spectra and to take advantage from "the flanking base effect", that allows the nucleobases next to 5' and 3' termini to stabilize the conformation of inner bases even if they do not contribute to the base pairing.Compared to polyA and polydA ss, the CD spectra of np-1/4/polyA (Figure 3a) and np-1/4/polydA (Figure 3b) hybrids exhibited hypochromic effects that could be ascribable to the cationic character of the side chain of Arginine residues.To assess the interaction between nucleopeptides and nucleic acids, CD binding experiments were performed in tandem cells, collecting ''sum'' CD spectra relative to the separated np-1/4 and ss polyA and polydA, and ''mix'' CD spectra, acquired after cell mixing and annealing of np-1/4 with polyA and polydA ss.The differences observed comparing the ''sum'' and ''mix'' CD spectra strongly suggested the occurrence of the hybridization event in both cases (Figure 4a and 4b).In particular, we observed dramatic changes in the CD spectra of hybrids that could be attributed to the loss of nucleic acids helicity resulting from non-specific electrostatic interactions with cationic nucleopeptides.Intriguingly, the strongest hypochromic effect was observed for hybrids containing full L-nucleopeptide, suggesting that besides the electrostatic contribution, the backbone stereochemistry plays a main role in determining the conformation assumed by the hybrid duplexes.
The thermal stability (melting temperature, Tm) of hybrids involving nucleopeptides np-1/4 and both polyA and polydA complements was determined by CD and UV experiments, and compared with Tm of the hexathymine oligodeoxynucleotide (T6)/PolyA and PolydA and ss.(Table 2).To assess the interaction between nucleopeptides and nucleic acids, CD binding experiments were performed in tandem cells, collecting "sum" CD spectra relative to the separated np-1/4 and ss polyA and polydA, and "mix" CD spectra, acquired after cell mixing and annealing of np-1/4 with polyA and polydA ss.The differences observed comparing the "sum" and "mix" CD spectra strongly suggested the occurrence of the hybridization event in both cases (Figure 4a,b).
Symmetry 2019, 11 FOR PEER REVIEW 8 As expected, enantiomeric nucleopeptides containing L-Arg/L-Lys and D-Arg/D-Lys (np-1 and -4) and D-Arg/L-Lys and L-Arg/D-Lys (np-2 and -3) exhibited mirror CD spectra.Next, we investigated the spectroscopic behavior of the hybrids composed of nucleopeptides np-1/4 with PolyA and PolydA (1:1 ratio in nucleobases).Long oligomers, such as polyA and polydA, were chosen for the CD experiments to amplify the resulting binding effects in CD spectra and to take advantage from "the flanking base effect", that allows the nucleobases next to 5' and 3' termini to stabilize the conformation of inner bases even if they do not contribute to the base pairing.Compared to polyA and polydA ss, the CD spectra of np-1/4/polyA (Figure 3a) and np-1/4/polydA (Figure 3b) hybrids exhibited hypochromic effects that could be ascribable to the cationic character of the side chain of Arginine residues.To assess the interaction between nucleopeptides and nucleic acids, CD binding experiments were performed in tandem cells, collecting ''sum'' CD spectra relative to the separated np-1/4 and ss polyA and polydA, and ''mix'' CD spectra, acquired after cell mixing and annealing of np-1/4 with polyA and polydA ss.The differences observed comparing the ''sum'' and ''mix'' CD spectra strongly suggested the occurrence of the hybridization event in both cases (Figure 4a and 4b).In particular, we observed dramatic changes in the CD spectra of hybrids that could be attributed to the loss of nucleic acids helicity resulting from non-specific electrostatic interactions with cationic nucleopeptides.Intriguingly, the strongest hypochromic effect was observed for hybrids containing full L-nucleopeptide, suggesting that besides the electrostatic contribution, the backbone stereochemistry plays a main role in determining the conformation assumed by the hybrid duplexes.
The thermal stability (melting temperature, Tm) of hybrids involving nucleopeptides np-1/4 and both polyA and polydA complements was determined by CD and UV experiments, and compared with Tm of the hexathymine oligodeoxynucleotide (T6)/PolyA and PolydA and ss.(Table 2).In particular, we observed dramatic changes in the CD spectra of hybrids that could be attributed to the loss of nucleic acids helicity resulting from non-specific electrostatic interactions with cationic nucleopeptides.Intriguingly, the strongest hypochromic effect was observed for hybrids containing full L-nucleopeptide, suggesting that besides the electrostatic contribution, the backbone stereochemistry plays a main role in determining the conformation assumed by the hybrid duplexes.
The thermal stability (melting temperature, Tm) of hybrids involving nucleopeptides np-1/4 and both polyA and polydA complements was determined by CD and UV experiments, and compared with Tm of the hexathymine oligodeoxynucleotide (T6)/PolyA and PolydA and ss.(Table 2).
Notably, the highest detected thermal stability was observed for np-1/polyA hybrid (Tm = 45 • C), revealing that np-1 represented the most stabilizing RNA binder even if compared to np-3/4/polyA and polydA hybrids.Moreover, the melting profile appeared as a sigmoid shape for np-1/polyA when compared to melting profiles of the any other hybrids (Figure 5).In particular, the melting profiles of polyA hybrids in the presence of l-Arg/d-Lys and d-Arg/d-Lys containing nucleopeptides (np-2 and -3) appeared very complex, suggesting the existence of more than two species and a multistate equilibrium (Figure 5).3'-TTTTTT-5' Not Detected Not Detected Notably, the highest detected thermal stability was observed for np-1/polyA hybrid (Tm = 45 °C), revealing that np-1 represented the most stabilizing RNA binder even if compared to np-3/4/polyA and polydA hybrids.Moreover, the melting profile appeared as a sigmoid shape for np-1/polyA when compared to melting profiles of the any other hybrids (Figure 5).In particular, the melting profiles of polyA hybrids in the presence of L-Arg/D-Lys and D-Arg/D-Lys containing nucleopeptides (np-2 and -3) appeared very complex, suggesting the existence of more than two species and a multistate equilibrium (Figure 5).The CD melting profiles suggested that the dissociation of np-1/polyA and np-1/polydA was a two-state process.Thus, we were able to estimate the apparent Kd values of np-1 to polyA and polydA by performing CD titration experiment with increasing concentrations of np-1 in the presence of a constant concentration of nucleic acids.As depicted in Figure 6, the apparent Kd values further corroborated the results achieved by us and other groups [37,39], indicating that cationic nucleopeptides are endowed with higher affinity for RNA (Figure 6a) with the respect to DNA molecules (Figure 6b).The CD melting profiles suggested that the dissociation of np-1/polyA and np-1/polydA was a two-state process.Thus, we were able to estimate the apparent Kd values of np-1 to polyA and polydA by performing CD titration experiment with increasing concentrations of np-1 in the presence of a constant concentration of nucleic acids.As depicted in Figure 6, the apparent Kd values further corroborated the results achieved by us and other groups [37,39], indicating that cationic nucleopeptides are endowed with higher affinity for RNA (Figure 6a) with the respect to DNA molecules (Figure 6b).

Molecular Dynamics Simulations
To gain structural insights into the binding capabilities of chiral nucleopeptide np-1 toward RNA, we performed molecular dynamics (MD) simulations of the single nucleopeptide consisting of interpolated L-Arg and L-Lys residues covalently linked to Thymine nucleobases.In this case, np-1 contained six L-Arg and six L-Lys residues with the Thymine nucleobases linked to the ε-NH2-group of the Lysine side chain through an acetyl group (Figure 6).This polymer was chosen as a model for the MD simulations due to its interesting spectroscopic behavior toward RNA ss.The interactions of the molecular complex cp1 (Figure 7) made up by a single-stranded 12 mer Adenine RNA oligonucleotide (A12 RNA ss) and two nucleopeptide np-1 molecules were assessed through MD simulations.The initial simulation structure was constructed based on the experimentally determined 2:1 stoichiometric interaction of np-1 with a complementary singlestranded 12 mer Adenine RNA oligonucleotide, and assuming an initial recognition event between the three macromolecules guided by the base-pairing interactions of the three species [40] (Figure 7).
After running the 200 ns long MD simulation (see Materials and Methods section), we analyzed the conformational stability of the simulated complex cp-1 in terms of the root-mean-square deviation (RMSD) of all the non-hydrogen atoms in the complex (Figure 8).Overall, the analysis of the MD trajectories demonstrated that after an initial equilibration period the complex reaches a fairly stable conformation with an 8 Å RMSD deviation from the initial input structure.This increase in RMSD is followed by a moderate second change, after approximately 40 ns, to reach a second conformation with an 11 Å RMSD deviation from the initial structure.This last conformation demonstrates to be fairly stable all through the remaining part of the simulation (160 ns).

Molecular Dynamics Simulations
To gain structural insights into the binding capabilities of chiral nucleopeptide np-1 toward RNA, we performed molecular dynamics (MD) simulations of the single nucleopeptide consisting of interpolated l-Arg and l-Lys residues covalently linked to Thymine nucleobases.In this case, np-1 contained six l-Arg and six l-Lys residues with the Thymine nucleobases linked to the ε-NH 2 -group of the Lysine side chain through an acetyl group (Figure 6).This polymer was chosen as a model for the MD simulations due to its interesting spectroscopic behavior toward RNA ss.
The interactions of the molecular complex cp1 (Figure 7) made up by a single-stranded 12 mer Adenine RNA oligonucleotide (A 12 RNA ss) and two nucleopeptide np-1 molecules were assessed through MD simulations.The initial simulation structure was constructed based on the experimentally determined 2:1 stoichiometric interaction of np-1 with a complementary single-stranded 12 mer Adenine RNA oligonucleotide, and assuming an initial recognition event between the three macromolecules guided by the base-pairing interactions of the three species [40] (Figure 7).

Molecular Dynamics Simulations
To gain structural insights into the binding capabilities of chiral nucleopeptide np-1 toward RNA, we performed molecular dynamics (MD) simulations of the single nucleopeptide consisting of interpolated L-Arg and L-Lys residues covalently linked to Thymine nucleobases.In this case, np-1 contained six L-Arg and six L-Lys residues with the Thymine nucleobases linked to the ε-NH2-group of the Lysine side chain through an acetyl group (Figure 6).This polymer was chosen as a model for the MD simulations due to its interesting spectroscopic behavior toward RNA ss.The interactions of the molecular complex cp1 (Figure 7) made up by a single-stranded 12 mer Adenine RNA oligonucleotide (A12 RNA ss) and two nucleopeptide np-1 molecules were assessed through MD simulations.The initial simulation structure was constructed based on the experimentally determined 2:1 stoichiometric interaction of np-1 with a complementary singlestranded 12 mer Adenine RNA oligonucleotide, and assuming an initial recognition event between the three macromolecules guided by the base-pairing interactions of the three species [40] (Figure 7).
After running the 200 ns long MD simulation (see Materials and Methods section), we analyzed the conformational stability of the simulated complex cp-1 in terms of the root-mean-square deviation (RMSD) of all the non-hydrogen atoms in the complex (Figure 8).Overall, the analysis of the MD trajectories demonstrated that after an initial equilibration period the complex reaches a fairly stable conformation with an 8 Å RMSD deviation from the initial input structure.This increase in RMSD is followed by a moderate second change, after approximately 40 ns, to reach a second conformation with an 11 Å RMSD deviation from the initial structure.This last conformation demonstrates to be fairly stable all through the remaining part of the simulation (160 ns).After running the 200 ns long MD simulation (see Materials and Methods section), we analyzed the conformational stability of the simulated complex cp-1 in terms of the root-mean-square deviation (RMSD) of all the non-hydrogen atoms in the complex (Figure 8).Overall, the analysis of the MD trajectories demonstrated that after an initial equilibration period the complex reaches a fairly stable conformation with an 8 Å RMSD deviation from the initial input structure.This increase in RMSD is followed by a moderate second change, after approximately 40 ns, to reach a second conformation with an 11 Å RMSD deviation from the initial structure.This last conformation demonstrates to be fairly stable all through the remaining part of the simulation (160 ns).To better determine the conformational behavior of the simulated complex, we conducted an RMSD-based clustering analysis of the MD 1002 trajectory frames obtained after the MD simulation.A total of three clusters, were obtained accounting for the most representative structures in the MD simulation of cp-1.Of the three clusters, the last one stands for the highest number of frames and includes the structures calculated for the last 160 ns of the simulations.The most representative structure of this cluster is depicted in Figure 9a, corresponding to a wrapped conformation of the complex.Next, we analyzed the structure depicted in Figure 9a to gather some information about the way np-1 interacts with the RNA counterpart.As can be observed in Figure 9a, the base-pairing of both To better determine the conformational behavior of the simulated complex, we conducted an RMSD-based clustering analysis of the MD 1002 trajectory frames obtained after the MD simulation.A total of three clusters, were obtained accounting for the most representative structures in the MD simulation of cp-1.Of the three clusters, the last one stands for the highest number of frames and includes the structures calculated for the last 160 ns of the simulations.The most representative structure of this cluster is depicted in Figure 9a, corresponding to a wrapped conformation of the complex.To better determine the conformational behavior of the simulated complex, we conducted an RMSD-based clustering analysis of the MD 1002 trajectory frames obtained after the MD simulation.A total of three clusters, were obtained accounting for the most representative structures in the MD simulation of cp-1.Of the three clusters, the last one stands for the highest number of frames and includes the structures calculated for the last 160 ns of the simulations.The most representative structure of this cluster is depicted in Figure 9a, corresponding to a wrapped conformation of the complex.Next, we analyzed the structure depicted in Figure 9a to gather some information about the way np-1 interacts with the RNA counterpart.As can be observed in Figure 9a, the base-pairing of both Next, we analyzed the structure depicted in Figure 9a to gather some information about the way np-1 interacts with the RNA counterpart.As can be observed in Figure 9a, the base-pairing of both np-1 molecules with the RNA chain is almost completely lost, with other chemical interactions taking place in the complex.More importantly, the RNA chain has lost its helical conformation and is adopting an almost linear disposition as can be seen in a zenithal view of the complex (Figure 9b).Responsible for this loss of helicity are the two np-1 molecules.Of these, the np-1 molecule (np-1-top) originally attached to the 5' end of the complex, uses its three N-terminal Arginines (R1, R2, and R3) to form stable ionic interactions with the RNA phosphate backbone at the 3' end.Moreover, np-1-top also forms H-bonds between its first and second N-terminal Thymine rings (T1 and T2) with the two middle RNA Adenine (A6 and A7), helping to keep the RNA chain in a near-linear conformation (Figure 10a).Similarly, important is the H-bond formed by the amide group in Thymine T3 arm with Adenine A7, which helps to hold the RNA chain before the turn that takes place between Adenine A7 and A8.The ionic interactions, however, are not the only chemical interactions between np-1-top and RNA.The complex also shows π-stacking interactions among C-terminal Thymine rings (T5 and T6) and Adenine A10 and A12.Both T5 and T6 are sandwiched between A10 and A12, holding the top half of the RNA chain with a successive six nucleobases π-stacking.
Symmetry 2019, 11 FOR PEER REVIEW 12 np-1 molecules with the RNA chain is almost completely lost, with other chemical interactions taking place in the complex.More importantly, the RNA chain has lost its helical conformation and is adopting an almost linear disposition as can be seen in a zenithal view of the complex (Figure 9b).Responsible for this loss of helicity are the two np-1 molecules.Of these, the np-1 molecule (np-1top) originally attached to the 5' end of the complex, uses its three N-terminal Arginines (R1, R2, and R3) to form stable ionic interactions with the RNA phosphate backbone at the 3' end.Moreover, np-1-top also forms H-bonds between its first and second N-terminal Thymine rings (T1 and T2) with the two middle RNA Adenine (A6 and A7), helping to keep the RNA chain in a near-linear conformation (Figure 10a).Similarly, important is the H-bond formed by the amide group in Thymine T3 arm with Adenine A7, which helps to hold the RNA chain before the turn that takes place between Adenine A7 and A8.The ionic interactions, however, are not the only chemical interactions between np-1-top and RNA.The complex also shows π-stacking interactions among Cterminal Thymine rings (T5 and T6) and Adenine A10 and A12.Both T5 and T6 are sandwiched between A10 and A12, holding the top half of the RNA chain with a successive six nucleobases πstacking.For the second np-1 molecule (np-1-btm) (Figure 10b), only the two C-terminal Arginines R5 and R6 are forming stable ionic interactions with the RNA phosphate backbone in the middle of the chain.The rest of the np-1-btm molecule, however, seems to be interacting exclusively at the 3' end of the RNA, where the majority of the Thymine rings are grouped in a cluster with Adenine (A2, A3, and A4) by π-stacking interactions, participating in the stabilization of the wrapped RNA conformation.For the second np-1 molecule (np-1-btm) (Figure 10b), only the two C-terminal Arginines R5 and R6 are forming stable ionic interactions with the RNA phosphate backbone in the middle of the chain.The rest of the np-1-btm molecule, however, seems to be interacting exclusively at the 3' end of the RNA, where the majority of the Thymine rings are grouped in a cluster with Adenine (A2, A3, and A4) by π-stacking interactions, participating in the stabilization of the wrapped RNA conformation.
To account for the importance of the ionic interactions described above, we plotted its distances during the complete MD simulation (Figure 11).The results show that all the interactions are highly conserved during the most part of the simulation.For np-1-top, after an initial equilibration period of approximately 40 ns, a fair stabilization of the H-bond distances occurs for Arginines R1, R2 and R3 (Figure 11a), while for Thymines T1, T2, and T3 the H-bond distances (two for T1, one for T2 and one for T3) are quite stable from the beginning of the simulation, particularly those between T1 and A6, and T2 and A7 (Figure 11b).As for np-1-btm, only Arginine R5 is able to form stable and significant ionic interactions with the phosphate backbone after the initial 40 ns equilibration period, contributing to the stabilization of the extended RNA chain (Figure 11c).These data remark the importance of the Arginine chains in the nucleopeptide np-1 to induce and maintain an extended conformation of the RNA chain, as well as the role of the nucleobase groups, not only for the initial recognition of the RNA chain but also to stabilize such an interaction with RNA.To account for the importance of the ionic interactions described above, we plotted its distances during the complete MD simulation (Figure 11).The results show that all the interactions are highly conserved during the most part of the simulation.For np-1-top, after an initial equilibration period of approximately 40 ns, a fair stabilization of the H-bond distances occurs for Arginines R1, R2 and R3 (Figure 11a), while for Thymines T1, T2, and T3 the H-bond distances (two for T1, one for T2 and one for T3) are quite stable from the beginning of the simulation, particularly those between T1 and A6, and T2 and A7 (Figure 11b).As for np-1-btm, only Arginine R5 is able to form stable and significant ionic interactions with the phosphate backbone after the initial 40 ns equilibration period, contributing to the stabilization of the extended RNA chain (Figure 11c).These data remark the importance of the Arginine chains in the nucleopeptide np-1 to induce and maintain an extended conformation of the RNA chain, as well as the role of the nucleobase groups, not only for the initial recognition of the RNA chain but also to stabilize such an interaction with RNA.

Conclusions
In this study, we report the synthesis of Arginine-based chiral nucleopeptides in which stereochemically-defined nucleobase and underivatized amino acids were alternated in the peptide backbone.With the purpose to investigate the role of backbone stereochemistry in determining the formation of DNA and RNA hybrids, we performed spectroscopic studies.CD experiments highlighted that the nucleopeptide having a fully l-backbone configuration and consisting of interpolated l-Arginine and l-Lysine residues covalently linked to Thymine nucleobases, formed the most intriguing and stable hybrid complexes with RNA molecules.Molecular Dynamics (MD) simulations provided insight into potential structure of the resulting complex of the l-nucleopeptide with RNA complementary strands, pointing out detailed chemical interactions involving both the l-backbone (ionics and H-bonds) and nucleobases (pairing and π-stacking) of the chiral nucleopeptide.In conclusion, we provided the first evidence of the effect of the backbone stereochemistry on the chemical interactions involved in the formation of the nucleopeptide/RNA complexes.These findings added a new piece to the complex nucleopeptide/nucleic interaction puzzle that could contribute to the future design of peptide-based molecules acting as antisense agents or diagnostic probes.

Figure 7 .
Figure 7. Schematic representation of complex cp-1 submitted to MD.

Figure 7 .
Figure 7. Schematic representation of complex cp-1 submitted to MD.

Figure 7 .
Figure 7. Schematic representation of complex cp-1 submitted to MD.

Figure 8 .
Figure 8. RMSD plot of MD trajectories from their respective initial frame in complex cp-1.

Figure 9 .
Figure 9. (a) Representative structure of the most populated cluster in cp-1; (b) Simplified representation of a zenithal view of cp-1.

Figure 8 .
Figure 8. RMSD plot of MD trajectories from their respective initial frame in complex cp-1.

Symmetry 2019, 11 FOR PEER REVIEW 11 Figure 8 .
Figure 8. RMSD plot of MD trajectories from their respective initial frame in complex cp-1.

Figure 9 .
Figure 9. (a) Representative structure of the most populated cluster in cp-1; (b) Simplified representation of a zenithal view of cp-1.

Figure 9 .
Figure 9. (a) Representative structure of the most populated cluster in cp-1; (b) Simplified representation of a zenithal view of cp-1.

Figure 10 .
Figure 10.(a) Simplified representation of the most important interactions between np-1top and RNA; (b) Simplified representation of the most important interactions between np-1-btm and RNA.RNA chain in orange, np-1-top in cyan blue, np-1-btm in pink, and ionic interactions in purple.

Figure 10 .
Figure 10.(a) Simplified representation of the most important interactions between np-1-top and RNA; (b) Simplified representation of the most important interactions between np-1-btm and RNA.RNA chain in orange, np-1-top in cyan blue, np-1-btm in pink, and ionic interactions in purple.

Table 1 .
Yield, purity, retention times, and analytical data of the described nucleopeptides.

Table 1 .
ion times, and analytical data of the described nucleopeptides.

Table 1 .
Yield, purity, retention times, and analytical data of the described nucleopeptides.

Table 1 .
Yield, purity, retention times, and analytical data of the described nucleopeptides.