The systematic approach to describing conformational rearrangements in G-quadruplexes.

Conformational changes in DNA G-quadruplex (GQ)-forming regions affect genome function and, thus, compose an interesting research topic. Computer modelling may yield insight into quadruplex folding and rearrangement, particularly molecular dynamics simulations. Here, we show that specific parameters, which are distinct from those commonly used in DNA conformational analyses, must be introduced for adequate interpretation and, most importantly, convenient visual representation of the quadruplex modelling results. We report a set of parameters that comprehensively and systematically describe GQ geometry in dynamics. The parameters include those related to quartet planarity, quadruplex twist, and quartet stacking; they are used to quantitatively characterise various types of quadruplexes and rearrangements, such as quartet distortion/disruption or deviation/bulging of a single nucleotide from the quartet plane. Our approach to describing conformational changes in quadruplexes using the new parameters is exemplified by telomeric quadruplex rearrangement, and the benefits of applying this approach to analyse other structures are discussed.

Introduction

Much research in the fields of molecular biology and biophysical chemistry is devoted to noncanonical nucleic acid structures. Conformational rearrangements in G-quadruplexes (GQs), which are planar arrays of guanine tetrads paired by Hoogsteen bonding, are particularly important because they may account for multiple transcription/translation/recombination and splicing-related events (Bochman, Paeschke, & Zakian, 2012; Burge, Parkinson, Hazel, Todd, & Neidle, 2006; Maizels, 2006; Murat & Balasubramanian, 2014). Molecular modelling is commonly used to complement NMR and XRD data on known GQ structures or to predict rearrangements in GQs upon interaction with various ligands (Haider & Neidle, 2010; Šponer, Cang, & Cheatham, 2012). A number of software packages and web servers, including VMD (Humphrey, Dalke, & Schulten, 1996), CURVES (Blanchet, Pasi, Zakrzewska, & Lavery, 2011), 3DNA (Colasanti, Lu, & Olson, 2013) and some others, allow for visualisation of noncanonical nucleic acid structures and the assessment of basic structural parameters. However, specific parameters are necessary to quantitatively describe GQ geometry. Ideally, such parameters should aid in characterising various types of GQ topologies in dynamics. Basic parameters, such as backbone torsion angles, glycosidic torsion angles and sugar pucker, are informative and applicable to flexibility characterisation (Cang, Sponer, & Cheatham, 2011; Li, Luo, Xue, & Li, 2010) but not GQ characterisation. Similarly, the traditional root mean square deviations (RMSD) assessment cannot be used to comprehensively characterise GQ conformational fluctuations as a function of time. Several researchers have attempted to introduce GQ-specific parameters, but a systematic study has not been reported.

Lech et al. described stacked G-tetrad geometry using two specific parameters in a study on the influence of stacking on GQ stability: relative rotation angle and the distance between parallel planar tetrads (separation distance). For opposite-polarity-stacked tetrads, the angle between stacked tetrads was defined as the rotation of one tetrad in the direction of its hydrogen bond polarity with respect to another tetrad through the central axis that is perpendicular to the tetrad plane. For same-polarity-stacked tetrads, the angle was arbitrarily defined as the rotation of the tetrad containing the 5′-end guanine in the direction of its hydrogen bond polarity (Lech, Heddi, & Phan, 2013). Both the rotation angle (also referred to as the “twist angle”) and separation distance are key characteristics of GQs, and their fluctuations indicate conformational changes in GQs. For example, twist angle time plots have been used to monitor dynamics of the thrombin-binding GQ aptamer (TBA) (Reshetnikov, Golovin, Spiridonova, Kopylov, & Šponer, 2010) (the twist angle was determined as the angle between the vectors that join adjacent guanine C1′ atoms in each of the two TBA tetrads). In another recent TBA modelling study (Reshetnikov, Golovin, & Kopylov, 2010), a more in-depth analysis of GQ dynamics was performed, and two original parameters were reported. The distance between the G-quartet guanine centres of mass (COMs) was used to monitor “in plane” GQ motions (fluctuations in a quartet plane), and the distance between the two tetragon COMs (the outer tetragon was formed by the N9 atoms, and the inner tetragon formed by the O6 atoms) was used to monitor “out of plane” GQ motions (the tendency of a quartet to form a tepee structure). An alternative approach to assessing “out of plane” GQ motions was proposed by Zhu, Xiao, and Liang (2013). They demonstrated G-tetrad deviation from the plane conformation in telomeric GQs using RMSD for the axis perpendicular to the tetrad plane. Additional specific parameters have been introduced to describe conformational changes in GQ loops (Zhu et al., 2013) and ion capture (Reshetnikov et al., 2011), but those types of changes are outside the scope of this article.

The above studies were not aimed at comprehensive GQ characterisation but were focused on different individual aspects of GQ conformational diversity and dynamics. We analysed and reconsidered the approaches used in those studies to propose an original set of parameters for monitoring GQ quartet planarity and integrity. Relative orientation of GQ quartets can be characterised using known parameters (Lech et al., 2013), i.e. twist angle and separation distance. Fluctuations of the new parameters indicate both “in plane” and “out of plane” GQ motions and may be used to monitor GQ distortion/disruption or single-nucleotide bulging. The latter type of GQ rearrangement has not been quantitatively analysed so far, but it is important given the recently reported new types of quadruplexes (Mukundan & Phan, 2013).

Herein, we use a model of the human telomeric quadruplex (telGQ) to illustrate the advantages of our new parameters in describing various types of motions because telGQ is a well-characterised structure, and its folding as well as unfolding pathways have been investigated in detail (Gray, Trent, & Chaires, 2014; Lane, Chaires, Gray, & Trent, 2008; Stadlbauer, Trantírek, Cheatham, Koča, & Šponer, 2014). Human telomere sequences can adopt different topologies, including an anti-parallel basket structure with one diagonal and two lateral loops (Wang & Patel, 1993), a parallel structure with all strand-reversal loops (Parkinson, Lee, & Neidle, 2002), [3+1] hybrid structures with one strand-reversal and two lateral loops (Luu, Phan, Kuryavyi, Lacroix, & Patel, 2006; Phan, Luu, & Patel, 2006) and certain higher-order structures (Petraccone et al., 2011). In a K+ solution, natural human telomere fragments comprising four G tracks and one- or two-nucleotide flanks are slowly exchanging conformers of the hybrid GQ, “hybrid-1” and “hybrid-2.” The conformers differ in the loop arrangement (Phan et al., 2006), and their ratio depends on the flanks; d[TAGGG(TTAGGG)3] forms up to 70% of hybrid-1 (PDB entry 2JSM), and d[TAGGG(TTAGGG)3TT] forms up to 70% of hybrid-2 (PDB entry 2JSL). In d[TAGGG(TTA GGG)3T], the two conformers form at comparable proportions (Phan, Kuryavyi, Luu, & Patel, 2007).

We also discuss a TBA analogue (the thrombin-binding aptamer with a non-natural (triazole) internucleotide linkage) to demonstrate that the new parameters also apply to modified and artificial GQ structures.

We report explicit solvent and implicit solvent molecular dynamics simulations (MD) for the TBA analogue and the telomeric GQ, respectively. The MD trajectories are relatively short because we did not mean to ascertain the conformations of the already known TBA and telGQ structures – we rather use them as model GQs. The systematic approach described herein and the new parameters can be employed for interpretation and convenient visual representation of any simulations of GQ dynamics.

Materials and methods

We modelled a telGQ hybrid-1 (d[TAGGG(TTAGGG)3]; PDB: 2JSM) in implicit solvent and a TBA analogue with a triazole internucleotide linkage (d[GGTTGGTT*TGGTTGG], where * is a modified fragment; for schematic representation, see Supplementary material) in explicit solvent (T = 300 K; the negative charges were neutralised by adding Na+ ions, and K+ ions were used to stabilise the triazole–TBA GQ structure) with the Amber 10 suite. Two models of the triazole–TBA complex with thrombin were obtained using the NMR-based and XRD-based models of the TBA/thrombin complex (PDB entries 1HAO and 1HUT, respectively). For more information and the details of the explicit solvent MD simulation, see (Varizhuk et al., 2013). Implicit solvent MD simulation of telGQ was performed as follows.

The conditions described in (Varizhuk et al., 2013) did not yield substantial conformational rearrangements in the telomeric GQ, and therefore, we increased the temperature to facilitate the rearrangements (T = 400 K). Implicit solvent in generalised Born/solvent-accessible surface area (GB/SA) formalism (Bashford & Case, 2000) was employed using the OBC application developed by Onufriev, Bashford, and Case (2000) in the presence of 0.1 M NaCl. The DNA interatomic interactions were calculated using a parmbsc0 force field. Prior to MD simulation, the model structure was minimised by 250 steps of steepest descent followed by 250 steps of conjugate gradient. Next, gradual heating to 400 K over 20 ps was performed. The SHAKE algorithm (Ryckaert, Ciccotti, & Berendsen, 1977) was applied to constrain the hydrogen atom bonds, which allowed us to use a 2 fs step. The dielectric constants 1 (interior) and 80 (exterior) were employed in the GB-MD simulations. The production phase of the GB-MD simulations continued until the quadruplex core was disrupted (approximately 34 ns). To control the temperature, we used the Langevin thermostat with the collision frequency 5 ps−1. Snapshot visualisation and conformational rearrangement analyses were performed using VMD program. Snapshots were collected every 0.1 ns. The new GQ-related parameters were calculated using programs written with the internal VMD atom selection language. To conveniently represent the parameter plots, short-term fluctuations were smoothed by the moving average method (the span for the moving average was five).

Results and discussion

GQ-related parameters: general description

We developed a set of parameters that facilitate comprehensive quantitative characterisation of potentially any GQ core in dynamics. A schematic representation of the hybrid-1 telGQ core is shown in Figure 1(A) as an illustration of our numeration system. We used core-related numeration, which is independent from loop arrangement and does not coincide with numeration along the chain, but is universal and convenient for indexing parameters. Quartet numeration begins at the 5′-terminus. The guanines in the 1st quartet are numbered in the direction of the hydrogen bonding polarity, from the donor to the acceptor, and from the 5′-terminus. In the case of telGQ, the direction is counterclockwise from the perspective of the 1st quartet. The guanines in other quartets are numbered in accordance with quartet 1 (e.g. the guanine stacked with Gua1 1 (the first guanine in quartet 1) will also be the 1st in quartet 2 (Gua1 2) and so forth). Here and below, the lower index of a guanine residue, particular atom or relative parameter indicates the guanine’s position in a quartet, while the upper index is a quartet number. For example, Gua2 3 is the second Gua in quartet 3 and N92 3 is N9 in Gua2 3.

Figure 1.

GQ-related parameters and the basic types of quartet bending.

First, we suggest using a gyration radius (the known parameter) to illustrate major structural distortions as an alternative or in addition to the commonly used parameter RMSD; the gyration radius can be estimated for the entire GQ or a particular quartet. A dramatic increase in this parameter indicates GQ/quartet collapse. Thus, a time plot of the gyration radius provides a general approximation of the MD trajectory.

Next, for a more detailed analysis, we suggest 6 basic subsets of parameters that indicate relative quartet positions (Figure 1(A); parameters 1 and 2) and motions in each particular quartet (Figure 1(B); parameters 3–6).

The distance between the quartet COMs, which is also referred to as the “separation distance” (Lech et al., 2013). It is usually equal to 3.5 Å unless some major GQ distortions occur.

The twist angle (Θ).

The twist angle for quartets α and β is the angle between the vectors N9 -N9 and N9 -N9 irrespectively of the quartet polarity. Each twist angle is the average of four values (i = 1–3 is the Gua number in a quartet; the fourth pair is N91–N94).

The twist angle depends significantly on the type of stacking interactions (Syn vs. Anti Gua conformations) and thus takes different values for quartet pairs with similar and opposite polarities (Lech et al., 2013). In our numeration system, it is typically close to 10° (Anti-Syn stacking) or 40° (Syn-Anti stacking) in the case of the opposite quartet polarities and approximately 25–30° in the case of similar polarities (Anti-Anti). The comparison of the separation distances and twist angles for several well-known GQ structures, available in PDB, is presented in the Supplementary material, Table S1.

The dihedral angle of four N1 atoms.

This parameter is basically similar to that proposed by Mashimo, Yagi, Sannohe, Rajendran, and Sugiyama (2010) (i.e. the dihedral angle of four O6 atoms);

The distances between the quartet COM and COMs of particular guanines or the COM(GQ)–COM(Gua) distances;

The angles between the normals to the Gua planes and the axis of a quartet () or the angle between the normals to Gua (delta ). The axis of a quartet is the vector that joins its COM with the COM of a neighbouring quartet. It does not coincide with a normal to a quartet. In the case of telGQ, the axis of quartet 1 joins the COMs of quartets 1 and 2. The axis of quartet 2 coincides with the axis of quartet 3 and joins their COMs.

The parameters and delta indicate the relative orientations of all guanines in a quartet plane and, therefore, are more informative than previously reported planarity-related parameters, i.e. COM(N9-tetragon)–COM(O6-tetragon) distances (Reshetnikov, Golovin, Spiridonova et al., 2010) and RMSDs in the axis perpendicular to the quartet plane (Zhu et al., 2013).

The angles between Gua axes (Ω). The Gua axis is the vector that joins the middle of the C4–C5 interval with the C8 atom. Ω indicates relative Gua rotation and is normally equal to ~90° for neighbouring guanines and ~180° for non-neighbouring guanines.

We additionally propose specific characteristics of particular types of quartet plane bending. Thus far, only “tepee-like” bending has been described in detail (Reshetnikov, Golovin, & Kopylov, 2010). The parameters used by Reshetnikov et al. (i.e. the COM(N9-tetragon)–COM(O6-tetragon) distances) appear reasonable and convenient for monitoring this type of quartet plane deformation. Herein, we describe “lengthwise” (Figure 1(C)) and “diagonal” (Figure 1(D)) quartet bending.

“Lengthwise” bending implies that two pairs of neighbouring guanines are approximately planar. For instance, in Figure 1(C), Gua is approximately coplanar with Gua, and Gua is approximately coplanar with Gua. The vectors N9–C2 and N9–C8 lie in the nominal GuaGua plane. The vectors N9–C2 and N9–C8 lie in the nominal GuaGua plane. Alternative “nominal planes” could be constructed using relative Gua COMs and the quartet COM. However, that approach would produce understated values of the angle between the nominal planes. “Diagonal” quartet bending implies that two non-neighbouring guanines are coplanar, but both or one of the other two guanines deviate. For instance, in Figure 1(D), Gua is coplanar with Gua; Gua/Gua or both deviate. The vectors that join the COMs of Gua, Gua and Gua lie in the GuaGuaGua nominal plane. The vectors that join the COMs of Gua, Gua and Gua lie in the GuaGuaGua nominal plane. The angle () between the nominal planes formed by Gua pairs or triads indicates “lengthwise” or “diagonal” quartet bending, respectively. In the latter case, this parameter is somewhat similar to the N1 dihedral angle.

Quantitative evaluation of conformational rearrangements in GQs

On the basis of the above parameters, we propose the following systematic approach to describing GQ dynamics with a focus on conformational changes in the GQ core.

(1) General characterisation of major structural changes (parameters: RMSD and gyration radii (R g) of the quartets and the whole GQ).

where is the radius vector of atom i, m is its mass, N is a number of atoms in the quartet (for R g quartet) or in the whole GQ (for R g GQ), and is the radius vector of the quartet COM (for R g quartet) or GQ COM (for R g GQ)

(2) Analysis of relative quartet positions (parameters: the separation distance (COM distance) and the twist angle Θ).

where α and β are quartet numbers, and is the vector that joins N9 atoms of Gua and Gua in the respective quartet.

(3) Analysis of motions within a quartet (parameters: the distances between the Gua COMs COM(Gua)–COM(GQ) or COM(Gua)–COM(quartet), the angles between the Gua planes ( or delta ), the dihedral angle of the N1 atoms and the angles between Gua axes).

The dihedral angle (Ξ) of the N1 atoms:

where is the vector that joins N1 atoms of Gua and Gua.

The angle between the normal to Gua ( and the axis :

where is the vector that joins COM of quartet α with the COM of the neighbouring quartet β, is the vector that joins N9 and N2 atoms in Gua , and joins N9 and O6 atoms in Gua .

The angle between the normals to Gua and Gua:

The angle between the axes of Gua and Gua:

where is the vector that joins the middle of the C4–C5 interval with the C8 atom in Gua, and joins the middle of the C4–C5 interval with the C8 atom in Gua.

(4) Analysis of quartet bending type (parameters: the angles () between nominal planes formed by Gua pairs and triads).

(4.1.) “Lengthwise” bending (the nominal planes are GuaGua and GuaGua).

where and are normals to the nominal planes GuaGua and GuaGua, respectively; is the vector that joins N9 and C2 atoms in Gua; joins N9 and C2 atoms in Gua; joins N9 with C8, and joins N9 with C8.

(4.2.) “Diagonal” bending (the nominal planes are GuaGuaGua and GuaGuaGua).

where and are normals to the nominal planes GuaGuaGua and GuaGuaGua, respectively; is the vector that joins COMs of Gua and Gua; is the vector that joins COMs of Gua and Gua, etc.

By comparing time plots for the above parameters, one can distinguish between particular types of GQ motions or estimate their relative contributions to GQ distortion. For instance, significant fluctuations in Ω and/or Gua COM distance plots, but not in and N1 dihedral angle plots, indicate in-plane quartet motions. If the Gua COM distance, Ω and the N1 dihedral angle fluctuate, while remains constant, the Gua deviates from but remains parallel to the quartet plane. N1 dihedral angle and/or fluctuations without significant COM fluctuations indicate quartet bending without disruption.

We illustrate this approach by the analysis of telGQ dynamics.

Application of the new parameters for the analysis of GQ MD trajectory

A 34 ns MD simulation was performed for the hybrid-1 telGQ (the schematic representation is given in Figure 1(A)). Despite the relatively low energetic barrier between the two hybrid telGQ conformations, we did not observe hybrid-2 formation, mainly because the simulation time was limited. NMR data indicate that interconversion is slow (Dai, Carver, Punchihewa, Jones, & Yang, 2007). Nevertheless, we observed certain conformational changes that indicate the beginning of major rearrangements (Figure 2).

Figure 2.

TelGQ MD simulation snapshots.

At 0 ns (Figure 2(A)), the quartets were planar, and each quartet was stabilised by 8 H-bonds. During the first 20 ns, we only observed minor fluctuations, such as moderate quartet bending and reversible H-bonding disruption. Major quartet distortions were observed at approximately 21–22 ns and 27–30 ns. At 21 ns, fluctuations in quartet 1 yielded significant deviation and, finally, bulging of Gua1 1 (Figure 2(B)). Theoretically, this could indicate the beginning of detachment of the 5′-terminal G-track from the GQ core, which is considered the primary pathway for hybrid-1 → hybrid-2 interconversion via the intermediate triplex structure (Mashimo et al., 2010). However, in our simulation, Gua1 1 bulging was reversible, and at 22 ns, quartet 1 recovered its planarity (Figure 2(C)); however, its H-bonds were substantially rearranged (Figure 2(D)). The NH(G1 1)–O6(G2 1) bond was broken, and the NH2(G2 1)–O6(G3 1) bond formed instead of the NH2(G2 1)–N7(G3 1) bond. Thereafter, Gua2 1 lost its H-bonds with the neighbouring guanines, which led to Gua2 1 bulging from the quartet at 27–28 ns (Figure 2(E)). At 29 ns, quartet 1 was fully disintegrated, and the remaining two quartets underwent bending (Figure 2(F)). At 30 ns, telGQ almost returned to its 28-ns conformation; however, thereafter, telGQ fully collapsed due to increasing fluctuation.

The telGQ MD trajectory was described quantitatively using the proposed algorithm.

Time plots of the gyration radii, which provide a general estimation of telGQ integrity; separation distances and twist angles, which indicate relative quartet positions, and the N1 dihedral angels, which allow for the rough assessment of quartet planarity, are shown in Figure 3.

Figure 3.

Basic GQ-related parameters and their fluctuations during a telGQ MD simulation.

The plots of COM(Gua)–COM(quartet) distances and angles, which illustrate the changes of quartet integrity and planarity, respectively, are shown in Figure 4. (For the COM(Gua)–COM(GQ) plots and delta plots, see Supplementary material, Figure S1). Figure 5 illustrates “diagonal” and “lengthwise” bending of telGQ quartets.

Figure 4.

Parameters that describe quartet integrity (A–C) and planarity (D–F) in telGQ.

Figure 5.

The angles between nominal planes formed by the Gua pairs and triads that characterise telGQ quartet bending.

Moderate gyration radii fluctuations at approximately 21–22 ns (Figure 3(A)) indicate reversible quartet distortions. Major fluctuations at 27–28 ns suggest substantial conformational changes that lead to a dramatic increase in the gyration radii after 30 ns, which indicates GQ collapse. The separation distance fluctuations at 21 ns (Figure 3(C), quartets 1 and 2) might appear to indicate mutual quartet divergence; however, as we show below, this is actually not the case. Quartet 1 does not deviate from quartet 2; only particular guanines move out of the quartet plane, and the plane bends significantly, which is evident from the N1 dihedral angle plot (Figure 3(D)).

The COM distance plots, N1 dihedral angle plots and plots in Figures 3 and 4 clearly illustrate the following conformational changes in telGQ: Gua1 1 deviation from quartet 1 at approximately 21 ns; Gua2 1 deviation from quartet 1 at approximately 27 ns and subsequent collapse of quartet 1 with simultaneous bending in quartets 2 and 3.

The disrupted integrity of quartet 1 at 21–22 ns and after 27 ns is evident in Figure 4(A) and (D). The maximum COM(Gua)–COM(quartet) (Figure 4(A)) and COM(Gua)–COM(GQ) (Figure S1(A)) distances were observed for Gua1 1 (21–22 ns) and Gua2 1 (27–28 ns). The absence of substantial fluctuations at 21–22 ns in the quartet 1 plots (Figure 4(D)) indicates that Gua1 1 deviates from but remains parallel to the quartet plane (which is also evident in the 22-ns MD snapshot, Figure 2(B)). After 27 ns, all fluctuations in quartet 1 increased (Figure 4(D)), and the maximum amplitude was observed for Gua2 1. In the delta plots (Figure S1(D)), maximum amplitudes were observed for the pairs Gua1 1–Gua2 1, Gua2 1–Gua3 1 and Gua2 1–Gua4 1. The latter observations indicate that Gua2 1 not only deviates from but also rotates out of the quartet plane (Figure 2(F)). The fluctuations in quartets 2 and 3 are relatively small compared with quartet 1 until 27 ns. Figure 4(B) and (C) illustrates the integrity of quartets 2 and 3, and their planarity is most evident in Figure 3(D). Figure 4(E) and (F) illustrates a gradual increase in the relative Gua plane inclination. Apparently, quartet 2 is the least flexible.

Relative rotation of the guanines (Ω deviation from the normal values) was insignificant in all quartets until 27 ns, when major distortions began (Supplementary material, Figure S2). Mixed “lengthwise”/“diagonal” bending was of quartets 1 and 2 (especially quartet 1) after 27 ns is evident from Figure 5(A) and (B).

The new set of GQ core-related parameters is clearly applicable not only to natural DNA/RNA GQs but also their analogues with artificial loops (e.g. various chemically modified GQ aptamers and their complexes with proteins or small molecule ligands). In our recent study, we described the conformation and dynamics of a modified thrombin-binding aptamer, triazole–TBA, in complex with thrombin. Triazole–TBA is an active thrombin inhibitor with enhanced resistance to nuclease hydrolysis due to a non-natural triazole internucleotide linkage in a central loop of the two-tetrad antiparallel GQ structure (Sequence: d[GGTTGGTT*TGGTTGG], * is a triazole internucleotide linkage).

Molecular dynamics simulations were used to predict the most likely triazole–TBA/thrombin-binding mode. Two models of the aptamer/protein complex were considered: model 1, in which the aptamer binds the thrombin exosite-1 through the TT loops, and model 2, in which the modified (TT*T) loop participates in the binding. The triazole–TBA models 1 and 2 were constructed using NMR and XRD structures of the unmodified TBA, PDB:1HAO and PDB:1HUT, respectively. MD simulation results suggest that model 1 is valid; the GQ was stable and did not dissociate from the protein throughout the 35 ns dynamics simulation, but in model 2, we observed significant deformation of the GQ core. Details on the triazole–TBA interaction with thrombin are described in (Varizhuk et al., 2013). Here, we focused on GQ conformational changes in models 1 and 2.

Figure S3 (Supplementary material) shows MD simulation snapshots that illustrate the major deformations of triazole–TBA GQ in complex with thrombin; the new set of parameters was used to monitor, illustrate and analyse those deformations. In general, triazole–TBA model 1 underwent only minor changes (i.e. insignificant “lengthwise” bending in quartet 2 at approximately 9.5 ns (Figure S3(B)) and in quartet 1 at approximately 30 ns (Figure S3(C))). Alternatively, in model 2, the quartet planarity was severely distorted (Figures S3(F–H) and S4). Major fluctuations began at approximately 27 ns (Figures S2(F), S4 and S5(B) and (D)). The and delta plots clearly show that Gua2 1 and Gua3 1 turned away/rotated from quartet 1 (Figure S6(E) and (G)), and Gua3 2 deviated from quartet 2 (Figure S6(F) and (H)). Profound “lengthwise” bending of quartets 1 and 2 was observed at 33 ns (Figure S7(C) and (D)); at the end of the simulation, most binding contacts with the protein were lost.

Conclusion

We described a systematic approach to discerning conformational changes in GQs based on a consecutive characterisation of GQ core integrity in general, the relative orientation of G-quartets and the motions in each quartet. New quartet integrity- and quartet planarity-related parameters were developed. Taken together, these characterisations will aid in discriminating between particular types of quartet deformation (e.g. bulging of a single guanine, inclination of a single-guanine, tepee- or wedge-like bending of a quartet.). The advantages of the new parameters were demonstrated using hybrid-1 telGQ and TBA, but they may be applied to any GQ structure.

Supplementary material

Supplementary material dealing with the comparison of relative quartet positions in different GQ types (Table S1), quartet integrity and polarity in telGQ (Figure S1), relative Gua rotation in telGQ (Figure S2), MD simulation for triazole-TBA (Figure S3) and quartet integrity/planarity in triazole-TBA (Figures S4-S7) in addition to detailed instructions for using the VMD scripts is available online at: http://dx.doi.org/10.1080/07391102.2015.1055303.

VMD scripts for calculating all GQ-related parameters are available from the authors at: http://niifhm.ru/wp-content/uploads/2015/04/VMD_script_for_describing_2_tetrad_GQs.txt; http://niifhm.ru/wp-content/uploads/2015/04/VMD_script_for_describing_3_tetrad_GQs.txt; http://niifhm.ru/wp-content/uploads/2015/04/VMD_script_for_describing_4_tetrad_GQs.txt

Disclosure statement

No potential conflict of interest was reported by the authors.

Click here for additional data file.

Click here for additional data file.