Molecular dynamics and Raman optical activity spectra reveal nucleotide conformation ratios in solution †

Nucleotide conformational flexibility aﬀects their biological functions. Although the spectroscopy of Raman optical activity (ROA) is well suited to structural analyses in aqueous solutions, the link between the spectral shape and the nucleotide geometry is not fully understood. We recorded the Raman and ROA spectra of model nucleotides (rAMP, rGMP, rCMP, and dTMP) and interpreted them on the basis of molecular dynamics (MD) combined with density functional theory (DFT). The relation between the sugar puckering, base conformation and spectral intensities is discussed. Hydrogen bonds between the sugar’s C3 0 hydroxyl and the phosphate groups were found to be important for the sugar puckering. The simulated spectra correlated well with the experimental data and provided an understanding of the dependence of the spectral shapes on conformational dynamics. Most of the strongest spectral bands could be assigned to vibrational molecular motions. Decomposition of the experimental spectra into calculated subspectra based on arbitrary maps of free energies provided experimental conformer populations, which could be used to verify and improve the MD predictions. The analyses indicate some flaws of common MD force fields, such as being unable to describe the fine conformer distribution. Also the accuracy of conformer populations obtained from the spectroscopic data depends on the simulations, improvement of which is desirable for gaining a more detailed insight in the future. Improvement of the spectroscopic and computational methodology for nucleotides also provides opportunities for its application to larger nucleic acids.


Introduction
2][3] Their structural variations affect the conformation and impact nucleotide recognition in living cells, ligand binding, and enzymatic activity.For example, adenosine diphosphate and triphosphate are involved in muscle work and act as an energy pool in many enzymatic processes.Remdesivir, a drug recommended for the treatment of the COVID-19 flu, is metabolized to an adenosine analogue as the active form. 4yclic adenosine monophosphate serves as a messenger as it controls and mediates activities of peptide hormones.Other nucleotides and oligonucleotides have been widely applied as therapeutic agents, where their biological activity could be manipulated by their chemical modifications. 5part from the more usual NMR, 6 solutions' properties of nucleotides can be conveniently studied using chiroptical spectroscopies, where the spectra are very sensitive to geometry variations.As a standard in the past, electronic circular dichroism (ECD) provides useful information, although the bands corresponding to electronic transitions are not well resolved and are difficult to interpret. 7The ECD signal also primarily arises from the nitrogenous base which does not reflect much conformation of the ribose or deoxyribose residues. 7Vibrational optical activity (VOA) techniques including vibrational circular dichroism (VCD) and Raman optical activity (ROA) in principle provide more numerous and more resolved bands than ECD.These come from all molecular parts and are better linked to a local molecular structure. 8owever, although VCD is very sensitive to the structure of oligonucleotides [9][10][11] and nucleic acids, [12][13][14] the signal of mononucleotides is rather weak. 15The aqueous environment and strong water absorption also cause problems for VCD.ROA pursued in the current study is measurable more easily, and a wider region of wavenumbers includes more bands providing information on the conformation of the glycosidic bond, sugar, and the phosphate group.Unlike the one-photon techniques (ECD and VCD), ROA spectroscopy detects a small difference in the scattering of the left and right circularly polarized light.][18] For some time, interpretation of ROA experiments relied on characteristic spectral markers of similar compounds 19,20 or statistical methods. 213][24] However, as for oligopeptides and proteins, 25,26 the computations are complicated by nucleotides' flexibility and polarity.In addition, while protein ROA essentially stems from the peptide backbone, 27 more diverse chemical groups are involved in nucleotides.Spectral intensities of aromatic base, sugar, and phosphate may be calculated with different precision and the resultant spectrum is distorted.Below, we need to address this problem specifically for the PQO stretching frequency.For the solvent, simplified dielectric solvent models 28,29 provide only limited accuracy 30,31 and computations explicitly involving the solvent-solute interactions are preferred.For the description of the molecular dynamics, ab initio molecular dynamics and time-dependent approaches 32,33 may be technically superior, but they are currently computationally too demanding for our systems.
In the present study, we are therefore using a combination of (classical) molecular dynamics (MD) and density functional theory (DFT), where solute-solvent clusters are extracted from MD trajectories.Using DFT, after a minimal optimization of geometries, spectral intensities are calculated and many clusters are averaged. 31,34The spectral analysis can be significantly enhanced by the decomposition of the experimental spectra into the simulated conformer sub-spectra. 35,36Instead of fitting the conformer populations directly, we use an intermediate fit using an arbitrary energy surface 37 that makes it possible to process a large number of clusters (B1000) at once.The information about conformer populations can also be used to verify and improve MD force fields, although their development goes beyond the topic of this work.Overall, in spite of the approximations involved, the MD/ DFT simulations provide an in-depth understanding of the observed Raman and ROA spectral features and they enable us to link spectral intensities to molecular vibrations and thus make the VOA spectroscopy even more useful in understanding the behavior of nucleotides and nucleic acids.

Spectral measurement
Deoxythymidine 5 0 -monophosphate disodium salt (dTMP, Z95%), cytidine 5 0 -monophosphate disodium salt (rCMP, Z99%), adenosine 5 0 -monophosphate (rAMP, Z97%) and guanosine 5 0 -monophosphate disodium salt (rGMP, Z99%) (Fig. 1) were purchased from Sigma-Aldrich and dissolved in deionized water to a concentration of 0.05 mol L À1 .This concentration, relatively low for ROA spectroscopy, was used to avoid nucleotide stacking.Solutions' pH was adjusted to 8 with NaOH to ensure fully deprotonated phosphate groups.Raman and ROA spectra were recorded at 25 1C on a spectrometer located at the Palacky Úniversity Olomouc. 38The experiments were performed using scattered circular polarization (SCP) modulation and backscattering geometry (1801).The spectra were recorded in a fused silica cell of 4 mm optical path length (3 mm width), 70 mL sample volume, in a range of 50-4500 cm À1 , using 532 nm excitation wavelength, 0.5 s illumination time cycles, and 8 cm À1 resolution.Laser power for the sample was 800 mW and the spectra were recorded for about 30 h.From the Raman spectra, the solvent (H 2 O) signal was subtracted, and minor baseline corrections were made.

Molecular dynamics (MD)
All MD simulations were carried out within the Amber18 program. 39For rGMP, four different force fields were tested: GAFF, 40 GAFF2, 41 RNA.OL3, 42 and DNA.Bsc1. 43GAFF2 was then used as a default.The TIP3P force field 44 was applied for water.Atomic partial charges were obtained using a standard electrostatic field fitting procedure, 45 based on ab initio calculations performed at the B3LYP 46 /6-31++G(d,p)/CPCM 47 level of theory.For all quantum chemistry, we used the Gaussian 16 program. 48he nucleotides, sodium counterions, and water were kept in 40 3 Å 3 cubic boxes, an NVT dynamics was run at the temperature of 300 K, using 1 fs integration time step, Langevin thermostat with 2 ps À1 collision frequency, particle mesh Ewald method, and 8 Å Lennard-Jones cutoff.After 1 ns equilibration, production runs at least 7 ns long followed, the trajectories were analyzed using the cpptraj module of Amber18, and intramolecular H-bonds were counted using 3 Å distance between heavy atoms and 1351 angle cutoffs.In addition to unrestricted MD, dependencies of the free energy on two torsion angles were obtained using the weighted histogram analysis method (WHAM). 49In this case, angles w and g were incremented by 361, harmonic force constants of 3 kcal mol À1 rad 2 were applied, and 1 ns trajectories were investigated.
Fig. 1 Investigated nucleotides, characteristic torsion angles g, w, and t 0 À t 4 are indicated for dTMP or rCMP, and sugar puckering was defined in a standard way 1 as P = (t 2 + t 4 Àt 1 À t 3 )/(3.078Â t 0 ) with the amplitude of t m = t 1 /cos(P).
This journal is © the Owner Societies 2023

Simulated Raman and ROA spectra
The spectra were computed for MD snapshot geometries by default for every 10 ps.From the free dynamics, the spectra of 100-700 snapshots were averaged as specified below, waters farther than 3.6 Å from the nucleotides were deleted, and the remaining clusters partially optimized in normal mode coordinates 50,51 at the B3PW91 52 /6-31++G**/CPCM level.At the same level, harmonic Raman and ROA intensities were calculated, with polarizability derivatives related to the water molecules set to zero.The B3PW91 functional provided reasonable results for ROA earlier 27,53 and in trial computations gave almost same results as the more common B3LYP functional (Fig. S1, ESI †).The intensities were convoluted with Lorentzian bands of 10 cm À1 full width at half-height.The calculated frequency of the phosphate symmetric stretching band was significantly lower than observed.Therefore, the PQO bond force constants were scaled to match the experiment, using a protocol previously used for the CQO bonds. 54For comparison with the experiment, resultant frequencies were scaled by 0.98.

Spectral fitting and decomposition algorithm
The frequency scaling enabled us to directly decompose an experimental spectrum S into subspectra S i of many snapshots, where n is the frequency.In the decomposition, the coefficients c i were expressed using arbitrary free energy maps as described in ref. 37, and current implementation is schematically overviewed in Fig. 2. The w and g angles (Fig. 1) served as the arbitrary variables, and we set the coefficients to c i = exp(ÀDG i / kT), where kT is the Boltzmann temperature quantum, and we used m = 5.Then we minimized the error as a function of the coefficients d ab , and the integration ran from 200 to 1800 cm À1 .Integrating the resultant probability p(w,g) = exp(ÀDG/kT) over the relevant coordinates, we could determine the relative populations of conformers I-VI defined below.The algorithm described above requires geometries with evenly spaced coordinates. 37Differently populated parts of the (w,g) coordinate plane would contribute differently to the error.The MD snapshots, however, are concentrated around energy-favored coordinates.Therefore, 400 intermediate (arbitrary) spectra (S i ) were calculated and used in the decomposition instead.Each intermediate spectrum corresponded to a point (w i ,g i ) on a grid in the plane (angle increment d = 181, 20 Â 20 points).If some MD geometries in the (w i À d/2,w i + d/2) Â (g i -d/2,g i + d/2) square were available, the arbitrary spectrum was calculated as an average of these clusters.Otherwise, MD geometry and spectrum closest to (w i ,g i ) were taken.
Fig. 2 The process of fitting an experimental spectrum with MD snapshot spectra, for two-dimensional arbitrary maps.
For the decomposition, MD snapshots (typically more than 1000) were taken from the WHAM computations covering a large portion of the conformational space.For the fitting, explicit water molecules were not included in computations of the spectra, which cause only minor differences in Raman and ROA spectral intensities (Fig. S2, ESI †).Raman bands that could be unambiguously assigned were shifted to the experimental positions, and frequencies between them interpolated.Note that although for an ideal case, true maps of free energies are obtained using this procedure, 37 given the approximation used we take them as intermediate constructs of the fitting procedure corresponding to the exact dependencies only approximately.

Effect of the force fields on rGMP geometry and spectra
The rGMP conformation was rather sensitive to the force field variant used; distributions of sugar puckering, w angle, and puckering amplitude obtained with the GAFF, GAFF2, RNA.OL3 and DNA.Bsc1 force field are plotted in Fig. 3. GAFF and GAFF2 both favor the anti conformer (w B 1801), in line with NMR observations. 55However, they give different sugar puckerings, where for GAFF, all values are rather equally possible, while GAFF2 prefers the C3 0 -endo pucker with P B 01.The RNA.OL3 and DNA.Bsc1 results are similar, and both force fields favor the syn conformer (w B 601) and C2 0 -endo pucker (P B 1801).
As discussed elsewhere 56 the w angle affects the occurrences of intramolecular hydrogen bonds.For example, for GAFF, the C'OHÁ Á ÁOP hydrogen bond was present in 25% of the snapshots analyzed, while GAFF2 gave 70%.The puckering amplitude (t m ) occurs most often around 331 for GAFF; for the other three force fields, the amplitude is higher (B391) and its distribution is narrower.In Table S1 (ESI †) are provided some hints of the differences between the four force fields; supposedly variations of the torsion angle parameters affect the conformer ratios most strongly.
The geometry variations are reflected in the average spectra calculated using the four force fields.In the Raman (Fig. 4, right) spectra, the variations are rather subtle.For example, the experimentally observable doublet 1324/1367 cm À1 is best reproduced with the GAFF2 force field, although the calculated splitting is still limited.ROA spectra (Fig. 4, left) differ much more and no force field gives satisfactory agreement with the experiment.GAFF and GAFF2 give well the 1340 cm À1 experimental negative band and the trends in the signal below 300 cm À1 ; the AE868/891 cm À1 bands are reasonably obtained only by RNA.Ol3, for the AE564/588 cm À1 couplet, DNA.Bsc1 works best, etc.The reasonable agreement of the computed Raman spectra and a worse one for ROA suggest that no force field gives very realistic conformer distribution of rGMP.We choose the GAFF2 force field as an advanced alternative of GAFF as a default.The RNA.Ol3 and DNA.Bsc1 force fields might be biased towards canonical nucleic acid geometries.

Nucleotide conformers
The pyrimidine nucleotides (dTMP, rCMP) behave slightly differently from those derived from purine (rGMP, rAMP).The calculated dependencies of the free energy on the g and w angles are plotted in Fig. 5, left, and the parameters for selected energy minima are listed in Table S2 (ESI †).The w torsion angle determines the position of the base to the sugar and most often adopts two values, syn (w B 801) and anti (w B 1801); the angle g adopts about three values (g + /g À /t B 601/À601/AE 1801). 57,58erhaps contra-intuitively, the smaller pyrimidine residues in dTMP and rCMP make larger areas energetically forbidden (e.g., with E 4 10 kcal mol À1 ) in these plots.For dTMP, the most favored conformer II is anti/t, for rCMP, it is anti/g À (I), both in agreement with the crystal structure. 59The anti/g À most populated conformer I of rGMP is consistent with ref. 55, and the preference of rAMP for anti/t (I) corresponds to the X-ray. 59evertheless, the simulations strongly suggest that for all nucleotides, multiple conformers/conformer classes are present in solutions.
An average distribution of the sodium ions around the molecule and the radial distribution function related to sodium and phosphorus as obtained by MD is for rGMP documented in Fig. 5, right.The probability to encounter a sodium molecule in the vicinity of the phosphate group is higher than elsewhere, and some ions are momentarily bonded to it.However, such vicinity involves relatively small volume, and most sodium atoms freely float in the solution.In our (40 Å) 3 box, only about 6% of the sodium atoms stayed closer than 4 Å to phosphorus.This justifies the approach used in the calculation of the nucleotide spectra where sodium atoms were included in MD and were then only included implicitly in DFT, as a part of the dielectric environment.

Intramolecular hydrogen bonds
The nucleotide conformation is also given by a fine equilibrium between intra and inter-molecular hydrogen bonds.The occurrence of snapshots where internal hydrogen bonds are formed as obtained by MD is for GAFF2 summarized in Table 1.The H-bonds somewhat correlate with the conformation of the sugar, which is shown for rGMP and rCMP in Fig. 6 where Hbond distances and sugar puckering as functions of time are plotted.Time dependencies of other coordinates for all four nucleotides are shown in Fig. S3 (ESI †).The C1 0 -exo (P B 1101) conformation of rCMP appears during the C3 0 OÁ Á ÁOP H-bonding, but the correlation is not particularly strong.It is stronger for rGMP, where the C3 0 OÁ Á ÁOP interaction (small d 1 ) stabilizes the sugar in the C2 0 -endo (PB1801) form; at the same time, H-bonding between C2 0 OH and the base is broken (increase of d 2 ).
A very frequent (74%, Table 1) H-bonding between the PO and C3 0 hydroxyl groups of dTMP results in the ''east type'' sugar conformations with P B 1001. 1 In rAMP, the C2 0 OHÁ Á ÁN3 H-bond occurs in 9% of the snapshots, and the overall low percentage of intramolecular hydrogen bonds makes possible an even distribution of the C2 0 -endo and C3 0 -endo conformations.Other relations among the principle coordinates (pseudorotation angle P, angles w, g, b) can be seen in conformer populations recalculated to the free energy maps in Fig. S4 (ESI †).As in Fig. 5 for w and g, the maps also confirm that the nucleotides have many degrees of freedom, but the coordinates are correlated, and a limited number of distinct conformers exist in solutions.

Spectral fitting and conformer ratios
The experimental Raman and ROA spectra of all four nucleotides are compared to the simulated ones in Fig. 7.A plain average from the GAFF2 MD snapshots is shown along with the fit where the calculated frequencies were scaled.Assignments of the main Raman and ROA bands are summarized in Table S3 (ESI †).For rGMP, the calculated spectra from 700 snapshots (Fig. 7) are reasonably close to those obtained with 100 snapshots (Fig. 4), which shows that little improvement would be achieved with longer MD/more geometries.This is also indicated by Fig. S5 (ESI †), top, where for rGMP, spectra generated from 7 ns, and 1 ms MD runs differ by less than B15% in maximal intensities.For rAMP, similarly as for rGMP, 100  ) in rGMP stabilizes C2 0 -endo (P B 1801), and for rCMP, the correlation is not so tight.
This journal is © the Owner Societies 2023 snapshots averaging produced quite similar spectra as for the 700 snapshots (Fig. S5, ESI †, bottom).The plain snapshot averaging (blue in Fig. 7) wellreproduces the main Raman features and also Raman spectra differences among the four nucleotides.This is also reflected in the high values of the similarity indices 60 (40.85), which are equal to one for identical spectra.A band-to-band assignment is nevertheless problematic in the region around 1300 cm À1 with a large number of vibrational modes, most frequently related to CH bending (Table S3, ESI †).As expected for the vibrational optical activity, 61 the similarity indices are much lower for ROA (0.42-0.60), although the main ROA spectral features are reproduced by the computations as well.As discussed in detail for rGMP (Fig. 4), inconsistencies occur in particular below 1300 cm À1 ; for dTMP, for example, ROA below 150 cm À1 is negative, while the simulation gives a bisignate pattern, and the experimental positive band at 1228 cm À1 is not reproduced, etc.Overall, however, the accuracy is reasonable and corresponds to that achieved for similarly complex (flexible and polar) molecules. 34,62Interesting is also the large lowfrequency ROA, close to the 50 cm À1 instrumental limit, where the parent Raman spectra are virtually featureless and start to diverge towards the Rayleigh line.Similar ROA bands were previously observed for sugars, 63 peptides 64 and chiral liquids, 65 and seem to provide interesting, although still relatively unexplored, possibility to expand ROA structural studies.The fitted ROA spectra (red in Fig. 7) are closer to the experiment than the plain average.This is not always true for the Raman spectra, indicating some benefits of the explicit water model used for the MD averaged (blue) spectra.The conformer populations obtained from MD and from the fit of the Raman and ROA spectra with calculated subspectra are listed in Table 2.As discussed before, Raman spectra are rather independent of the conformation, and conformer populations derived from them may be less reliable than for ROA. 36,63The limited sensitivity of Raman intensities can be also seen in Fig. S6 (ESI †), where calculated spectra of individual conformers are compared.Nevertheless, the Raman results significantly differ from the ROA populations only for rAMP.For dTMP, there is a large discrepancy between the MD and ROA I : II conformer ratios (8 : 92 vs. 54 : 46), suggesting that the GAFF2 force field may not be appropriate for this molecule.For rCMP, the experimental conformer ratios almost copy the MD results, within the expected accuracy of the fitting/decomposition process (B10-20%). 35,36,66A large difference also appears for conformers of rGMP, e.g., 68% of I with WHAM, while both Raman and ROA give 6%.For rAMP, the ROA and MD conformer populations are closer.
Although the arbitrary dependencies of the free energies on the (g,w) angles obtained during the spectral fitting are hampered by the computational errors, some similarities with those obtained by WHAM indicate that they are partially realistic (Fig. S7, ESI †).The free energy maps obtained from ROA spectra are closer to the MD/WHAM results, confirming that the conformational predictions based on ROA are more reliable than for Raman spectra.As follows also from Table 2, in some cases, the MD and spectra-derived conformer populations differ a lot, and the energy minima obtained by the two approaches do not have the same coordinates.Overall, the spectral analysis in light of the theoretical simulations suggests rather inadequate conformer distribution obtained from MD, and although the accuracy of the numbers obtained from the fitting is also limited, it provides enormous potential to study potential energy surfaces of molecules in solutions.

Spectra and local molecular structure
Similarly as for nucleosides, 16 the nucleotide ROA spectrum can be roughly divided in regions sensitive to coupling of the sugar ring and the base (B1200-1600 cm À1 ), sugar puckering (B900-1200 cm À1 ), and sugar and base deformations (o900 cm À1 ).The (PO 3 ) 2À group stretching contributes to the signal around 1150 cm À1 , with weak ROA.Conformation of the glycosidic bond itself is reflected by many ROA bands; its effect on the strongest bands for rGMP in the 1300-1500 cm À1 region is shown in Fig. 8.The comparison with the experiment suggests that the syn conformer with a positive signal around 1380 cm À1 is significantly present in the sample, but not dominant, in accordance with the MD and spectra analyses (Table 2).Similar comparison can be done in other spectral regions and for other nucleotides (Fig. S8, ESI †), which on one side documents the enormous sensitivity of ROA to the conformation and on the other one also shows the limits of the contemporary methodology.In particular, better precision of the simulation would be desirable to increase the effective spatial resolution of this kind of spectroscopy.
To understand the nucleotide spectral origin more deeply, in a computer experiment, the intensity tensors (polarizability derivatives) 67 for different molecular parts were set to zero.Thus approximate contributions of nucleotides' chromophores and interactions to Raman and ROA signals could be estimated for rGMP, as shown in Fig. 9. From the left side of the figure, we see that in the entire spectral region, the signal of the aromatic base dominates in the Raman spectrum, except for the phosphate band at 1007 cm À1 .In spite of its inherent achirality (planarity), the base also provides large ROA contributions.However, the ROA signal becomes much bigger when interaction between the chromophores is enabled by including polarizabilities from the base and sugar, or from the phosphate and sugar (right side of the figure).In particular, the basesugar interaction gives rise to strong positive and negative bands around 1350 cm À1 .These bands are caused by the base and sugar CH bending modes, and their coupling also explains the ROA sensitivity to the variation of the w angle (Fig. 8).
We did not attempt to analyze the lowest-frequency (o100 cm À1 ) ROA signal, because of its complexity and extensive computer demands needed. 65Nevertheless, the experimental spectra plotted in the full range (Fig. S9, ESI †) also show in this region significant differences between the nucleotides; this signal is particularly strong for rAMP, stronger than in any other wavenumber region.Also the CH stretching Raman and ROA bands (B2900 cm À1 ) provide so far limited information, as they are weak and prone to error in experiment, and eventual simulation would require rather expensive computations involving anharmonic corrections. 38As explained before, 63 ROA in the OH/NH stretching region (43000 cm À1 ) cannot be measured reliably, due to the large Raman scattering of water.Maybe in the future, it can provide some information about the molecular hydration shell.

Conclusions
We measured and analyzed the Raman and ROA spectra of four common nucleotides to investigate the link of the spectral shapes to the geometry, and to probe the quality of the molecular dynamics force fields and spectral simulation techniques.The results appeared helpful in the effort to better understand the nucleotide conformational behavior and to develop the spectroscopy of vibrational optical activity as a universal tool in molecular structural studies.
An analysis of molecular dynamics trajectories revealed that sugar puckering was strongly dependent on intramolecular H-bonding, which differed in different nucleotides.In dTMP, H-bonding between the PO and C3 0 hydroxyl groups stabilized the sugar-pucker in C1 0 -exo conformation.The H-bonding was less extended in rAMP where it allowed for an equilibrium of the C2 0 -endo and C3 0 -endo sugar conformation.
The conformer distribution obtained from molecular dynamics combined with DFT provided a reasonable basis to interpret the experimental spectra.Raman and ROA bands were assigned, which was useful for monitoring the conformation of the base and sugar puckering.The phosphate group gave a rather weak ROA signal.The decomposition of the experimental spectra into the theoretical ones indicated further potential of the optical spectroscopy, where in principle the true molecular potential energy surface can be obtained.However, the procedure currently has a drawback of the limited precision of the molecular dynamics and DFT computations.Nevertheless, Raman optical activity already appears as a useful tool for the study of the solution structure of nucleic acid components.

Fig. 3
Fig. 3 Calculated probabilities p of puckering P, torsion angle w, and amplitude t m for rGMPs and GAFF, GAFF2, RNA.OL3 and DNA.Bsc1 force fields.

Fig. 4
Fig. 4 Calculated and experimental Raman (I R + I L ) and ROA (I R À I L ) spectra of rGMP.Calculated averages from 100 snapshots obtained with four force fields (GAFF, GAFF2, RNA.OL3 and DNA.Bsc1) are shown, with explicit water molecules, and frequencies are scaled by 0.98.

Fig. 5 (
Fig. 5 (left) Dependence of the free energy on the torsion angles w and g (in degrees) as obtained from the metadynamics, (right) 3-D distribution of the sodium counterions around rGMP and the corresponding radial distribution function related to P and Na atoms.The purple surface encloses the region with elevated Na + concentration.

Fig. 7
Fig. 7 Raman (I R + I L ) and ROA (I R À I L ) spectra of studied nucleotides, plain average from 700 MD snapshots (blue), a fit (red), and experiment (black).The calculated spectra were normalized to the strongest Raman band, and similarity indices 60 between scaled calculated and experimental spectra, Ð S cal S exp dn= Ð S cal 2 dn Ð S exp 2 dn À Á 1=2 , are indicated for the integrals within 200-1900 cm À1 .

Table 1
Percentages of MD snapshots with detected intramolecular hydrogen bonds

Table 2
Conformer populations obtained from WHAM (%), and decomposition of Raman and ROA spectra, see also Fig.5and TableS2(ESI) for the conformer definition Fig. 8 Calculated ROA spectra of the four rGMP conformers (MD average, frequency multiplied by 0.98) and the experiment.This journal is © the Owner Societies 2023