A Comprehensive Analysis of the Structural Recognition between KCTD Proteins and Cullin 3.

KCTD ((K)potassium Channel Tetramerization Domain-containing) proteins constitute an emerging class of proteins involved in fundamental physio-pathological processes. In these proteins, the BTB domain, which represents the defining element of the family, may have the dual role of promoting oligomerization and favoring functionally important partnerships with different interactors. Here, by exploiting the potential of recently developed methodologies for protein structure prediction, we report a comprehensive analysis of the interactions of all KCTD proteins with their most common partner Cullin 3 (Cul3). The data here presented demonstrate the impressive ability of this approach to discriminate between KCTDs that interact with Cul3 and those that do not. Indeed, reliable and stable models of the complexes were only obtained for the 15 members of the family that are known to interact with Cul3. The generation of three-dimensional models for all KCTD-Cul3 complexes provides interesting clues on the determinants of the structural basis of this partnership as clear structural differences emerged between KCTDs that bind or do not bind Cul3. Finally, the availability of accurate three-dimensional models for KCTD-Cul3 interactions may be valuable for the ad hoc design and development of compounds targeting specific KCTDs that are involved in several common diseases.

Computing the electronic circular dichroism spectrum of DNA quadruple helices of different topology: A critical test for a generalized excitonic model based on a fragment diabatization.

In this study, we exploit a recently developed fragment diabatization-based excitonic model, FrDEx, to simulate the electronic circular dichroism (ECD) spectra of three guanine-rich DNA sequences arranged in guanine quadruple helices with different topologies: thrombin binding aptamer (antiparallel), c-Myc promoter (parallel), and human telomeric sequence (3+1 hybrid). Starting from time-dependent density functional theory (TD-DFT) calculations with the M052X functional, we apply our protocol to parameterize the FrDEX Hamiltonian, which accounts for electron density overlap and includes both the coupling with charge transfer transitions and the effect of the surrounding bases on the local excitation of each chromophore. The TD-DFT/M052X spectral shapes are in good agreement with the experimental ones, the main source of discrepancy being related to the intrinsic error on the computed transition energies of guanine monomer. FrDEx spectra are fairly close to the reference TD-DFT ones, allowing a significant advance with respect to a more standard excitonic Hamiltonian. We also show that the ECD spectra are sensitive to the inclusion of the inner K + cation in the calculation.

Unveiling CD59-Antibody Interactions to Design Paratope-Mimicking Peptides for Complement Modulation.

CD59 is an abundant immuno-regulatory human protein that protects cells from damage by inhibiting the complement system. CD59 inhibits the assembly of the Membrane Attack Complex (MAC), the bactericidal pore-forming toxin of the innate immune system. In addition, several pathogenic viruses, including HIV-1, escape complement-mediated virolysis by incorporating this complement inhibitor in their own viral envelope. This makes human pathogenic viruses, such as HIV-1, not neutralised by the complement in human fluids. CD59 is also overexpressed in several cancer cells to resist the complement attack. Consistent with its importance as a therapeutical target, CD59-targeting antibodies have been proven to be successful in hindering HIV-1 growth and counteracting the effect of complement inhibition by specific cancer cells. In this work, we make use of bioinformatics and computational tools to identify CD59 interactions with blocking antibodies and to describe molecular details of the paratope-epitope interface. Based on this information, we design and produce paratope-mimicking bicyclic peptides able to target CD59. Our results set the basis for the development of antibody-mimicking small molecules targeting CD59 with potential therapeutic interest as complement activators.

Alphafold Predictions Provide Insights into the Structural Features of the Functional Oligomers of All Members of the KCTD Family.

Oligomerization endows proteins with some key properties such as extra-stabilization, long-range allosteric regulation(s), and partnerships not accessible to their monomeric counterparts. How oligomerization is achieved and preserved during evolution is a subject of remarkable scientific relevance. By exploiting the abilities of the machine-learning algorithms implemented in AlphaFold (AF) in predicting protein structures, herein, we report a comprehensive analysis of the structural states of functional oligomers of all members of the KCTD protein family. Interestingly, our approach led to the identification of reliable three-dimensional models for the pentameric states of KCNRG, KCTD6, KCTD4, KCTD7, KCTD9, and KCTD14 and possibly for KCTD11 and KCTD21 that are involved in key biological processes and that were previously uncharacterized from a structural point of view. Although for most of these proteins, the CTD domains lack any sequence similarity, they share some important structural features, such as a propeller-like structure with a central cavity delimited by five exposed and regular ß-strands. Moreover, the structure of the related proteins KCTD7 and KCTD14, although pentameric, appears to be characterized by a different organization of the CTD region, with the five chains forming a circle-like structure with a large cavity. Our predictions also suggest that other members of the family, such as KCTD10, KCTD13, and TNFAIP1, present a strong propensity to assume dimeric states. Although the structures of the functional oligomers reported herein represent models that require additional validations, they provide a consistent and global view of KCTD protein oligomerization.

Structure-Based Development of SARS-CoV-2 Spike Interactors.

Coronaviruses, including SARS-CoV-2 (the etiological agent of the current COVID-19 pandemic), rely on the surface spike glycoprotein to access the host cells, mainly through the interaction of their receptor-binding domain (RBD) with the human angiotensin-converting enzyme 2 (ACE2). Therefore, molecular entities able to interfere with the binding of the SARS-CoV-2 spike protein to ACE2 have great potential to inhibit viral entry. Starting from the available structural data on the interaction between SARS-CoV-2 spike protein and the host ACE2 receptor, we engineered a set of soluble and stable spike interactors, here denoted as S-plugs. Starting from the prototype S-plug, we adopted a computational approach by combining stability prediction, associated to single-point mutations, with molecular dynamics to enhance both S-plug thermostability and binding affinity to the spike protein. The best developed molecule, S-plug3, possesses a highly stable α-helical con-formation (with melting temperature Tm of 54 °C) and can interact with the spike RBD and S1 domains with similar low nanomolar affinities. Importantly, S-plug3 exposes the spike RBD to almost the same interface as the human ACE2 receptor, aimed at the recognition of all ACE2-accessing coronaviruses. Consistently, S-plug3 preserves a low nanomolar dissociation constant with the delta B.1.617.2 variant of SARS-CoV-2 spike protein (KD = 29.2 ± 0.6 nM). Taken together, we provide valid starting data for the development of therapeutical and diagnostic tools against coronaviruses accessing through ACE2.

Electronic Circular Dichroism Spectra of DNA Quadruple Helices Studied by Molecular Dynamics Simulations and Excitonic Calculations including Charge Transfer States.

We here investigate the Electronic Circular Dichroism (ECD) Spectra of two representative Guanine-rich sequences folded in a Quadruple helix (GQ), by using a recently developed fragment diabatisation based excitonic model (FrDEx). FrDEx can include charge transfer (CT) excited states and consider the effect of the surrounding monomers on the local excitations (LEs). When applied to different structures generated by molecular dynamics simulations on a fragment of the human telomeric sequence (Tel21/22), FrDEx provides spectra fully consistent with the experimental one and in good agreement with that provided by quantum mechanical (QM) method used for its parametrization, i.e., TD-M05-2X. We show that the ECD spectrum is moderately sensitive to the conformation adopted by the bases of the loops and more significantly to the thermal fluctuations of the Guanine tetrads. In particular, we show how changes in the overlap of the tetrads modulate the intensity of the ECD signal. We illustrate how this correlates with changes in the character of the excitonic states at the bottom of the La and Lb bands, with larger LE and CT involvement of bases that are more closely stacked. As an additional test, we utilised FrDEx to compute the ECD spectrum of the monomeric and dimeric forms of a GQ forming sequence T30695 (5'TGGGTGGGTGGGTGGG3'), i.e., a system containing up to 24 Guanine bases, and demonstrated the satisfactory reproduction of the experimental and QM reference results. This study provides new insights on the effects modulating the ECD spectra of GQs and, more generally, further validates FrDEx as an effective tool to predict and assign the spectra of closely stacked multichromophore systems.

Studying the excited electronic states of guanine rich DNA quadruplexes by quantum mechanical methods: main achievements and perspectives.

The main insights into the photoactivated dynamics of guanine quadruplexes (G4s) recently provided by quantum mechanical computations are concisely reviewed here. The experimental steady state absorption and circular dichroism spectra of different topologies can be reproduced and assigned. After light absorption from excited states delocalized over multiple bases, the most important decay pathways involve localization of the excitation over a single base or on two stacked guanines, excimers with different degrees of charge transfer character. Two main photochemical reactions are discussed. One involves the photodimerization of two stacked guanine bases on the 'neutral' excimer path. The other, ionization of guanine, which triggers deprotonation of the resulting cation to form (G-H2)˙ and (G-H1)˙ radicals. Both the static and dynamical properties of G4 excited states are ruled by their topology and modulated by the inner coordinated metal ions.

Multiple Electronic and Structural Factors Control Cyclobutane Pyrimidine Dimer and 6-4 Thymine-Thymine Photodimerization in a DNA Duplex.

The T-T photodimerization paths leading to the formation of cyclobutane pyrimidine dimer (CPD) and 6-4 pyrimidine pyrimidone (64-PP), the two main DNA photolesions, have been resolved for a T-T step in a DNA duplex by two complementary state-of-the-art quantum mechanical approaches: QM(CASPT2//CASSCF)/MM and TD-DFT/PCM. Based on the analysis of several different representative structures, we define a new-ensemble of cooperating geometrical and electronic factors (besides the distance between the reacting bonds) ruling T-T photodimerization in DNA. CPD is formed by a barrierless path on an exciton state delocalized over the two bases. Large interbase stacking and shift values, together with a small pseudorotation phase angle for T at the 3'-end, favor this reaction. The oxetane intermediate, leading to a 64-PP adduct, is formed on a singlet T→T charge-transfer state and is favored by a large interbase angle and slide values. A small energy barrier (<0.3 eV) is associated to this path, likely contributing to the smaller quantum yield observed for this process. Eventually, a clear directionality is always shown by the electronic excitation characterizing the singlet photoactive state driving the photodimerization process: an exciton that is more localized on T3 and a 5'-T→3'-T charge transfer for CPD and oxetane formation, respectively, thus calling for specific electronic constraints.

Molecular recognition of Cullin3 by KCTDs: insights from experimental and computational investigations.

Recent investigations have highlighted a key role of the proteins of the KCTD (K-potassium channel tetramerization domain containing proteins) family in several fundamental biological processes. Despite the growing importance of KCTDs, our current understanding of their biophysical and structural properties is very limited. Biochemical characterizations of these proteins have shown that most of them act as substrate adaptor in E3 ligases during protein ubiquitination. Here we present a characterization of the KCTD5-Cullin3 interactions which are mediated by the KCTD5 BTB domain. Isothermal titration calorimetry experiments reveal that KCTD5 avidly binds the Cullin3 (Cul3). The complex presents a 5:5 stoichiometry and a dissociation constant of 59 nM. Molecular modeling and molecular dynamics simulations clearly indicate that the two proteins form a stable (KCTD5-Cul3)(5) pinwheel-shaped heterodecamer in which two distinct KCTD5 subunits cooperate in the binding of each cullin chain. Molecular dynamics simulations indicate that different types of interactions contribute to the stability of the assembly. Interestingly, residues involved in Cul3 recognitions are conserved in the KCTD5 orthologs and paralogs implicated in important biological processes. These residues are also rather well preserved in most of the other KCTD proteins. By using molecular modeling techniques, the entire ubiquitination system including the E3 ligase, the E2 conjugating enzyme and ubiquitin was generated. The analysis of the molecular architecture of this complex machinery provides insights into the ubiquitination processes which involve E3 ligases with a high structural complexity.

The determinants of bond angle variability in protein/peptide backbones: A comprehensive statistical/quantum mechanics analysis.

The elucidation of the mutual influence between peptide bond geometry and local conformation has important implications for protein structure refinement, validation, and prediction. To gain insights into the structural determinants and the energetic contributions associated with protein/peptide backbone plasticity, we here report an extensive analysis of the variability of the peptide bond angles by combining statistical analyses of protein structures and quantum mechanics calculations on small model peptide systems. Our analyses demonstrate that all the backbone bond angles strongly depend on the peptide conformation and unveil the existence of regular trends as function of ψ and/or φ. The excellent agreement of the quantum mechanics calculations with the statistical surveys of protein structures validates the computational scheme here employed and demonstrates that the valence geometry of protein/peptide backbone is primarily dictated by local interactions. Notably, for the first time we show that the position of the H(α) hydrogen atom, which is an important parameter in NMR structural studies, is also dependent on the local conformation. Most of the trends observed may be satisfactorily explained by invoking steric repulsive interactions; in some specific cases the valence bond variability is also influenced by hydrogen-bond like interactions. Moreover, we can provide a reliable estimate of the energies involved in the interplay between geometry and conformations.

Bond distances in polypeptide backbones depend on the local conformation.

By combining quantum-mechanical analysis of small model peptides and statistical surveys of high-resolution protein structures, a systematic conformational dependence of bond lengths in polypeptide backbones has been unveiled which involves both the peptide bond (C-O and C-N) and those bonds centred on the C(α) atom. All of these bond lengths indeed display a systematic variability in the ψ angle according to both calculations and surveys of protein structures. The overall agreement between the computed and the statistical data suggests that these trends are essentially driven by local effects. The dependence of C(α) distances on ψ is governed by interactions between the σ system of the C(α) moiety and the C-O π system of the peptide bond. Maximum and minimum values for each bond distance are found for conformations with the specific bond perpendicular and parallel to the adjacent CONH peptide plane, respectively. On the other hand, the variability of the C-O and C-N distances is related to the strength of the interactions between the lone pair of the N atom and the C-O π* system, which is modulated by the ψ angle. The C-O and C-N distances are related but their trends are not strictly connected to peptide-bond planarity, although a correlation amongst all of these parameters is expected on the basis of the classical resonance model.

Effect of C5-methylation of cytosine on the photoreactivity of DNA: a joint experimental and computational study of TCG trinucleotides.

DNA methylation, occurring at the 5 position of cytosine, is a natural process associated with mutational hotspots in skin tumors. By combining experimental techniques (optical spectroscopy, HPLC coupled to mass spectrometry) with theoretical methods (molecular dynamics, DFT/TD-DFT calculations in solution), we study trinucleotides with key sequences (TCG/T5mCG) in the UV-induced DNA damage. We show how the extra methyl, affecting the conformational equilibria and, hence, the electronic excited states, increases the quantum yield for the formation of cyclobutane dimers while reducing that of (6-4) adducts.

A comprehensive analysis of SARS-CoV-2 missense mutations indicates that all possible amino acid replacements in the viral proteins occurred within the first two-and-a-half years of the pandemic.

The surveillance of COVID-19 pandemic has led to the determination of millions of genome sequences of the SARS-CoV-2 virus, with the accumulation of a wealth of information never collected before for an infectious disease. Exploring the information retrieved from the GISAID database reporting at that time >13 million genome sequences, we classified the 141,639 unique missense mutations detected in the first two-and-a-half years (up to October 2022) of the pandemic. Notably, our analysis indicates that 98.2 % of all possible conservative amino acid replacements occurred. Even non-conservative mutations were highly represented (73.9 %). For a significant number of residues (3 %), all possible replacements with the other nineteen amino acids have been observed. These observations strongly indicate that, in this time interval, the virus explored all possible alternatives in terms of missense mutations for all sites of its polypeptide chain and that those that are not observed severely affect SARS-CoV-2 integrity. The implications of the present findings go well beyond the structural biology of SARS-CoV-2 as the huge amount of information here collected and classified may be valuable for the elucidation of the sequence-structure-function relationships in proteins.

The photoactivated dynamics of dGpdC and dCpdG sequences in DNA: a comprehensive quantum mechanical study.

Study of alternating DNA GC sequences by different time-resolved spectroscopies has provided fundamental information on the interaction between UV light and DNA, a process of great biological importance. Multiple decay paths have been identified, but their interplay is still poorly understood. Here, we characterize the photophysics of GC-DNA by integrating different computational approaches, to study molecular models including up to 6 bases described at a full quantum mechanical level. Quantum dynamical simulations, exploiting a nonadiabatic linear vibronic coupling (LVC) model, coupled with molecular dynamics sampling of the initial structures of a (GC)5 DNA duplex, provide new insights into the photophysics in the sub-picosecond time-regime. They indicate a substantial population transfer, within 50 fs, from the spectroscopic states towards G → C charge transfer states involving two stacked bases (CTintra), thus explaining the ultrafast disappearance of fluorescence. This picture is consistent with that provided by quantum mechanical geometry optimizations, using time dependent-density functional theory and a polarizable continuum model, which we use to parametrize the LVC model and to map the main excited state deactivation pathways. For the first time, the infrared and excited state absorption signatures of the various states along these pathways are comprehensively mapped. The computational models suggest that the main deactivation pathways, which, according to experiment, lead to ground state recovery on the 10-50 ps time scale, involve CTintra followed by interstrand proton transfer from the neutral G to C-. Our calculations indicate that CTintra is populated to a larger extent and more rapidly in GC than in CG steps and suggest the likely involvement of monomer-like and interstrand charge transfer decay routes for isolated and less stacked CG steps. These findings underscore the importance of the DNA sequence and thermal fluctuations for the dynamics. They will also aid the interpretation of experimental results on other sequences.

Biochemical and structural characterization of recombinant short-chain NAD(H)-dependent dehydrogenase/reductase from Sulfolobus acidocaldarius highly enantioselective on diaryl diketone benzil.

The gene encoding a novel alcohol dehydrogenase that belongs to the short-chain dehydrogenases/reductases superfamily was identified in the aerobic thermoacidophilic crenarchaeon Sulfolobus acidocaldarius strain DSM 639. The saadh2 gene was heterologously overexpressed in Escherichia coli, and the resulting protein (SaADH2) was purified to homogeneity and both biochemically and structurally characterized. The crystal structure of the SaADH2 NADH-bound form reveals that the enzyme is a tetramer consisting of identical 27,024-Da subunits, each composed of 255 amino acids. The enzyme has remarkable thermophilicity and thermal stability, displaying activity at temperatures up to 80 °C and a 30-min half-inactivation temperature of ∼88 °C. It also shows good tolerance to common organic solvents and a strict requirement for NAD(H) as the coenzyme. SaADH2 displays a preference for the reduction of alicyclic, bicyclic and aromatic ketones and α-ketoesters, but is poorly active on aliphatic, cyclic and aromatic alcohols, showing no activity on aldehydes. Interestingly, the enzyme catalyses the asymmetric reduction of benzil to (R)-benzoin with both excellent conversion (98 %) and optical purity (98 %) by way of an efficient in situ NADH-recycling system involving a second thermophilic ADH. The crystal structure of the binary complex SaADH2-NADH, determined at 1.75 Å resolution, reveals details of the active site providing hints on the structural basis of the enzyme enantioselectivity.

Crystallographic and spectroscopic characterizations of Sulfolobus solfataricus TrxA1 provide insights into the determinants of thioredoxin fold stability.

Structural characterizations of thioredoxins (Trxs) are important for their involvement in severe pathologies and for their stable scaffold. Here we report a combined structural and spectroscopic characterization of a Trx isolated from the hyperthermophilic archaeon Sulfolobus solfataricus (SsTrxA1). Thermal denaturation unveils that SsTrxA1 is endowed with a remarkable stability in the explored temperature range 50-105°C. The structure of the oxidized form of SsTrxA1 determined at 1.9Å resolution presents a number of peculiar features. Although the protein was crystallized in a slightly acid medium (pH 6.5) as many as ten intramolecular/intermolecular carboxyl-carboxylate interactions involving glutamic and aspartic acid side chains are found in three independent SsTrxA1 molecules present in the asymmetric unit. Surprisingly for a hyperthermostable protein, the structure of SsTrxA1 is characterized by the presence (a) of a very limited number of intramolecular salt bridges and (b) of a cavity nearby Cys52, a residue that is frequently a phenylananine in other members of the family. Chemical denaturation investigations carried out on SsTrxA1 and SsTrxA2 show that both proteins present a significant stability against guanidine hydrochloride, thus indicating that ionic interactions play a minor role in their stabilization. Compared to Trxs from mesophilic sources, SsTrxA1 displays a longer α-helix 1 and a shorter loop connecting this α-helix with ß-strand 2. As these features are shared with Trxs isolated from thermophilic sources, the shortening of this loop may be a general strategy adopted to stabilize this fold. This feature may be exploited for the design of hyperthermostable Trx scaffolds.

Local Backbone Geometry Plays a Critical Role in Determining Conformational Preferences of Amino Acid Residues in Proteins.

The definition of the structural basis of the conformational preferences of the genetically encoded amino acid residues is an important yet unresolved issue of structural biology. In order to gain insights into this intricate topic, we here determined and compared the amino acid propensity scales for different (φ, ψ) regions of the Ramachandran plot and for different secondary structure elements. These propensities were calculated using the Chou-Fasman approach on a database of non-redundant protein chains retrieved from the Protein Data Bank. Similarities between propensity scales were evaluated by linear regression analyses. One of the most striking and unexpected findings is that distant regions of the Ramachandran plot may exhibit significantly similar propensity scales. On the other hand, contiguous regions of the Ramachandran plot may present anticorrelated propensities. In order to provide an interpretative background to these results, we evaluated the role that the local variability of protein backbone geometry plays in this context. Our analysis indicates that (dis)similarities of propensity scales between different regions of the Ramachandran plot are coupled with (dis)similarities in the local geometry. The concept that similarities of the propensity scales are dictated by the similarity of the NCαC angle and not necessarily by the similarity of the (φ, ψ) conformation may have far-reaching implications in the field.

Structural and Binding Properties of the Active Cell Wall Hydrolase RipA from M. tuberculosis, a Promising Biosensing Molecule for Early Warning Bacterial Detection.

BACKGROUND: Peptidoglycan is an essential component of the cell wall in all bacteria. In particular, the cell walls of Gram-positive bacteria are composed mostly of a thick layer of peptidoglycan. Its accessibility has important implications for their sensing in whole bacterial detection methodologies. Indeed, there is an urgent demand for rapid tests which can identify whole bacteria, e.g., directly at the point of care. OBJECTIVE: The aim of this work is to explore the suitability of RipA, a key cell division protein of M. tuberculosis, for whole cell biosensing of Gram-positive bacteria. METHODS: We here conducted Molecular Dynamics (MD) studies aimed at the understanding of the structural and dynamic features of active RipA and at the design of a suitable bioreceptor. Based on these studies, we engineered a RipA variant for covalent oriented immobilisation on golden surfaces and are able to bind peptidoglycan, albeit without degrading it. Surface Plasmon Resonance (SPR) was employed to check the ability of functionalized golden chips to recognize whole bacteria. RESULTS: MD analyses elucidated the structural details of the active form of RipA and suggested that this enzyme, once inactivated, presents a rigid and well-exposed peptidoglycan recognition cleft. We engineered RipA for proper oriented immobilisation on golden chips for SPR studies. Results show that once chemically coupled to a golden chip, the developed RipA-based bioreceptor is able to detect B. subtilis, used as a model in a concentration-dependent mode. CONCLUSION: Results highlight the potential of the engineered molecule to be integrated in the development of early warning biosensors for Gram-positive contamination in clinical diagnosis or food-borne infections.

Effect of the Thermal Fluctuations of the Photophysics of GC and CG DNA Steps: A Computational Dynamical Study.

Here we refine and assess two computational procedures aimed to include the effect of thermal fluctuations on the electronic spectra and the ultrafast excited state dynamics of multichromophore systems, focusing on DNA duplexes. Our approach is based on a fragment diabatization procedure that, from a given Quantum Mechanical (QM) reference method, can provide the parameters (energy and coupling) of the reference diabatic states on the basis of the isolated fragments, either for a purely electronic excitonic Hamiltonian (FrDEx) or a linear vibronic coupling Hamiltonian (FrD-LVC). After having defined the most cost-effective procedure for DNA duplexes on two smaller fragments, FrDEx is used to simulate the absorption and Electronic Circular Dichroism (ECD) spectra of (GC)5 sequences, including the coupling with the Charge Transfer (CT) states, on a number of structures extracted from classical Molecular Dynamics (MD) simulations. The computed spectra are close to the reference TD-DFT calculations and fully consistent with the experimental ones. We then couple MD simulations and FrD-LVC to simulate the interplay between local excitations and CT transitions, both intrastrand and interstrand, in GC and CG steps when included in a oligoGC or in oligoAT DNA sequence. We predict that for both sequences a substantial part of the photoexcited population on G and C decays, within 50-100 fs, to the corresponding intrastrand CT states. This transfer is more effective for GC steps that, on average, are more closely stacked than CG ones.

Early steps of oxidative damage in DNA quadruplexes are position-dependent: Quantum mechanical and molecular dynamics analysis of human telomeric sequence containing ionized guanine.

Guanine radical cation (G•+) is a key intermediate in many oxidative processes occurring in nucleic acids. Here, by combining mixed Quantum Mechanical/Molecular Mechanics calculations and Molecular Dynamics (MD) simulations, we study how the structural behaviour of a tract GGG(TTAGGG)3 (hereafter Tel21) of the human telomeric sequence, folded in an antiparallel quadruple helix, changes when one of the G bases is ionized to G•+ (Tel21+). Once assessed that the electron-hole is localized on a single G, we perform MD simulations of twelve Tel21+ systems, differing in the position of G•+ in the sequence. When G•+ is located in the tetrad adjacent to the diagonal loop, we observe substantial structural rearrangements, which can decrease the electrostatic repulsion with the inner Na+ ions and increase the solvent exposed surface of G•+. Analysis of solvation patterns of G•+ provides new insights on the main reactions of G•+, i.e. the deprotonation at two different sites and hydration at the C8 atom, the first steps of the processes producing 8oxo-Guanine. We suggest the main structural determinants of the relative reactivity of each position and our conclusions, consistent with the available experimental trends, can help rationalizing the reactivity of other G-quadruplex topologies.