The covariation between TpA deficiency, CpG deficiency, and G+C content of human isochores is due to a mathematical artifact

CpG and TpA dinucleotides are underrepresented in the human genome. The CpG deficiency is due to the high mutation rate from C to T in methylated CpG's. The TpA suppression was thought to reflect a counterselection against TpA's destabilizing effect in RNA. Unexpectedly, the TpA and CpG deficiencies vary according to the G+C contents of sequences. It has been proposed that the variation in CpG suppression was correlated with a particular chromatin organization in G+C-rich isochores. Here, we present an improved model of dinucleotide evolution accounting for the overlap between successive dinucleotides. We show that an increased mutation rate from CpG to TpG or CpA induces both an apparent TpA deficiency and a correlation between CpG and TpA deficiencies and G+C content. Moreover, this model shows that the ratio of observed over expected CpG frequency underestimates the real CpG deficiency in G+C-rich sequences. The predictions of our model fit well with observed frequencies in human genomic data. This study suggests that previously published selectionist interpretations of patterns of dinucleotide frequencies should be taken with caution. Moreover, we propose new criteria to identify unmethylated CpG islands taking into account this bias in the measure of CpG depletion.     

NMR investigation of Hoogsteen base pairing in quinoxaline antibiotic--DNA complexes: comparison of 2:1 echinomycin, triostin A and [N-MeCys3,N-MeCys7] TANDEM complexes with DNA oligonucleotides.

Hoogsteen base pairs have been demonstrated to occur in base pairs adjacent to the CpG binding sites in complexes of triostin A and echinomycin with a variety of DNA oligonucleotides. To understand the relationship of these unusual base pairs to the sequence specificity of these quinoxaline antibiotics, the conformation of the base pairs flanking the YpR binding sites of the 2:1 drug-DNA complexes of triostin A with [d(ACGTACGT)]2 and of the TpA specific [N-MeCys3, N-MeCys7] TANDEM with [d(ATACGTAT)]2 have been studied by 1H NMR spectroscopy. In both the 2:1 triostin A-DNA complex and the 2:1 [N-MeCys3, N-MeCys7] TANDEM-DNA complex, the terminal A.T base pairs are Hoogsteen base paired with the 5' adenine in the syn conformation. This indicates that both TpA specific and CpG specific quinoxaline antibiotics are capable of inducing Hoogsteen base pairs in DNA. However, in both 2:1 complexes, Hoogsteen base pairing is limited to the terminal base pairs. In the 2:1 triostin A complex, the internal adenines are anti and in the 2:1 [N-MeCys3, N-MeCys7] TANDEM-DNA complex, the internal guanines are anti regardless of pH, which indicates that the central base pairs of both complexes form Watson-Crick base pairs. This indicates that the sequence dependent nature of Hoogsteen base pairing is the same in TpA specific and CpG specific quinoxaline antibiotic-DNA complexes. We have calculated a low resolution three-dimensional structure of the 2triostin A-[d(ACGTACGT)]2 complex and compared it with other CpG specific quinoxaline antibiotic-DNA complexes. The role of stacking in the formation of Hoogsteen base pairs in these complexes is discussed.     

CpG and TpA frequencies in the plant system.

Higher plant nuclear sequences reveal avoidance of CpG and TpA doublets. Chloroplast sequences avoid the TpA doublet in all codon positions. The chloroplast genome is not methylated but codon positions II-III and untranslated regions avoid CpG. The mitochondrial genome, also unmethylated, avoids CpG in all codon positions. We therefore deduce that methylation is not sufficient to explain CpG avoidance in the higher plant systems. Other factors must be taken into account such as amino acid composition, codon choices and perhaps stability of the DNA helix.     

Solution structure of 2-(pyrido[1,2-e]purin-4-yl)amino-ethanol intercalated in the DNA duplex d(CGATCG)2

The solution structure of the complex formed between d(CGATCG)(2) and 2-(pyrido[1,2-e]purin-4-yl)amino-ethanol, a new antitumor drug under design, has been resolved using NMR spectroscopy and restrained molecular dynamic simulations. The drug molecule intercalates between each of the CpG dinucleotide steps with its side chain lying in the minor groove. Analysis of NMR data establishes a weak stacking interaction between the intercalated ligand and the DNA bases; however, the drug/DNA affinity is enhanced by a hydrogen bond between the hydroxyl group of the end of the intercalant side chain and the amide group of guanine G6. Unrestrained molecular dynamic simulations performed in a water box confirm the stability of the intercalation model. The structure of the intercalated complex enables insight into the structure-activity relationship, allowing rationalization of the design of new antineoplasic agents.     

DNA minor groove recognition of a non-self-complementary AT-rich sequence by a tris-benzimidazole ligand.

The crystal structure of the non-self-complementary dodecamer DNA duplex formed by d(CG[5BrC]ATAT-TTGCG) and d(CGCAAATATGCG) has been solved to 2.3 A resolution, together with that of its complex with the tris-benzimidazole minor groove binding ligand TRIBIZ. The inclusion of a bromine atom on one strand in each structure enabled the possibility of disorder to be discounted. The native structure has an exceptional narrow minor groove, of 2.5-2.6 A in the central part of the A/T region, which is increased in width by approximately 0.8 A on drug binding. The ligand molecule binds in the central part of the sequence. The benzimidazole subunits of the ligand participate in six bifurcated hydrogen bonds with A:T base pair edges, three to each DNA strand. The presence of a pair of C-H...O hydrogen bonds has been deduced from the close proximity of the pyrrolidine group of the ligand to the TpA step in the sequence.     

Crystal structure of the topoisomerase II poison 9-amino-[N-(2-dimethylamino)ethyl]acridine-4-carboxamide bound to the DNA hexanucleotide d(CGTACG)2.

The structure of the complex formed between d(CGTACG)(2) and the antitumor agent 9-amino-[N-(2-dimethylamino)ethyl]acridine-4-carboxamide has been solved to a resolution of 1.6 A using X-ray crystallography. The complex crystallized in space group P6(4) with unit cell dimensions a = b = 30.2 A and c = 39.7 A, alpha = beta = 90 degrees, gamma = 120 degrees. The asymmetric unit contains a single strand of DNA, 1. 5 drug molecules, and 29 water molecules. The final structure has an overall R factor of 19.3%. A drug molecule intercalates between each of the CpG dinucleotide steps with its side chain lying in the major groove, and the protonated dimethylamino group partially occupies positions close to ( approximately 3.0 A) the N7 and O6 atoms of guanine G2. A water molecule forms bridging hydrogen bonds between the 4-carboxamide NH and the phosphate group of the same guanine. Sugar rings adopt the C2'-endo conformation except for cytosine C1 which moves to C3'-endo, thereby preventing steric collision between its C2' methylene group and the intercalated acridine ring. The intercalation cavity is opened by rotations of the main chain torsion angles alpha and gamma at guanines G2 and G6. Intercalation perturbs helix winding throughout the hexanucleotide compared to B-DNA, steps 1 and 2 being unwound by 8 degrees and 12 degrees, respectively, whereas the central TpA step is overwound by 17 degrees. An additional drug molecule, lying with the 2-fold axis in the plane of the acridine ring, is located at the end of each DNA helix, linking it to the next duplex to form a continuously stacked structure. The protonated N,N-dimethylamino group of this "end-stacked" drug hydrogen bonds to the N7 atom of guanine G6. In both drug molecules, the 4-carboxamide group is internally hydrogen bonded to the protonated N-10 atom of the acridine ring. The structure of the intercalated complex enables a rationalization of the known structure-activity relationships for inhibition of topoisomerase II activity, cytotoxicity, and DNA-binding kinetics for 9-aminoacridine-4-carboxamides.     

A molecular understanding of mitoxantrone-DNA adduct formation: effect of cytosine methylation and flanking sequences

When mitoxantrone is activated by formaldehyde it can form adducts with DNA. These occur preferentially at CpG and CpA sequences and are enhanced 2-3-fold at methylated CpG sequences compared with non-methylated sites. We sought to understand the molecular factors involved in enhanced adduct formation at these methylated sites. This required, first, clarification of factors that contributed to the formation of adducts at CpG sites. For this purpose mass spectrometry of an oligonucleotide duplex (containing a single CpG adduct site) was used to confirm the presence of an additional carbon atom (derived from formaldehyde) on the drug-DNA complex. The effect of 3'-flanking sequences was revealed by electrophoretic analysis of oligonucleotide-drug adducts, and the preferred adduct-forming site was identified as 5'-CGG-3'. Radiolabeled studies of drug-DNA adducts confirmed that the site of attachment involved the exocyclic amino of guanine. Molecular modeling analysis of the relative stability of the intercalated form of mitoxantrone was consistent with observed adduct-forming potential of CG sites with varying flanking sequences. The known preference for adduct formation at methylated CG sites was confirmed by energetics calculations and shown to be due to a shift of equilibrium of the intercalated form of the drug from the major groove (at CG sites) to the minor groove (at methylated CG sites). This increases the relative amount of drug that is located adjacent to the N-2 exocyclic amino of guanine in the minor groove, where covalent linkage is facilitated. These results account for the enhanced covalent binding of mitoxantrone to methylated CG sequences and provide a molecular model of the interactions.     

Actinomycin D facilitates transition of AT domains in molecules of sequence (AT)nAGCT(AT)n to a DNAse I detectable alternating structure.

The interaction of actinomycin D with (AT)nAGCT(AT)n (where n = 2, 3, or 4) was investigated using a combination of imino proton NMR and DNAse I digestion. The stoichiometry of the interaction appears to be one:one with the actinomycin chromophore intercalated between the two GC base pairs. This binding event facilitates the conversion of the flanking repetitive AT regions to an alternating conformation characterized by induced sensitivity of the ApT sequences to attack by DNAse I. The neighboring TpA sequences do not exhibit rate changes as a function of binding of the drug. The potential relevance of such ligand induced DNA structural alterations is discussed.     

Flexibility of the B-DNA backbone: effects of local and neighbouring sequences on pyrimidine-purine steps.

The structurally correlated dihedral angles epsilon and zeta are known for their large variability within the B-DNA backbone. We have used molecular modelling to study both energetic and mechanical features of these variations which can produce BI/BII transitions. Calculations were carried out on DNA oligomers containing either YpR or RpY dinucleotides steps within various sequence environments. The results indicate that CpA and CpG steps favour the BI/BII transition more than TpA or any RpY step. The stacking energy and its intra- and inter-strand components explain these effects. Analysis of neighbouring base pairs reveals that BI/BII transitions of CpG and CpA are easiest within (Y)n(R)n sequences. These can also induce a large vibrational amplitude for TpA steps within the BI conformation.     

First-neighbor specificities of actinomycin-DNA bindings by circular dichroism.

The circular dichroism spectra of eleven double-stranded DNAs, five natural with known nearest neighbor frequencies and six synthetic polydimers and polytrimers, were measured from 210 to 310 nm in the absence and presence of increasing amounts of actinomycin up to saturation. Based on the fact that the circular dichroism of nucleic acids is a nearest-neighbor frequency-dependent property, matrix analysis of the problem revealed which neighbor sets were perturbed by actinomycin, presumably by intercalation of the planar moiety of the molecule. The intercalation sites can be separated into three families. The first-neighbor units GpC and CpG are very favorable binding sites for actinomycin. ApG, CpC, ApC, TpC, and TpG appear to be less attractive sites, while ApT, TpA, and ApA are unfavorable sites.     

Association of anthracyclines and synthetic hexanucleotides. Structural factors influencing sequence specificity

The equilibrium and kinetic aspects of the interaction between four anthracyclines and two synthetic self-complementary hexanucleotides was investigated by fluorescence detection. Two of the studied anthracyclines are widely used antitumor drugs: doxorubicin (1, formerly adriamycin) and daunorubicin (2, formerly daunomycin). The other two, 9-deoxydoxorubicin (3) and 3'-deamino-3'-hydroxy-4'-epidoxorubicin (4), are doxorubicin analogues with modifications of the chemical groups that have been proposed as responsible for sequence specificity (Chen, K.-X., Gresh, N. and Pullman, B. (1985). J. Biomol. Struct. Dyn. 3, 445-466). One of the oligonucleotides, d(CGTACG), is identical to that used in the high resolution x-ray structure determination of the daunorubicin intercalative complex (Wang, A. H.-J., Ughetto, G., Quigley, G. J. & Rich, A. (1987). Biochemistry 26, 1152-1163). Binding to this hexanucleotide is compared with intercalation into the d(CGCGCG) duplex, revealing sequence preferences of the four anthracyclines. Taking into account the anthracycline aggregation and the dissociation of the hexanucleotide double standard form, results can be interpreted with a model that assumes complete fluorescence quenching at intercalative sites containing the CG base pair, and a large residual fluorescence after intercalation within the TpA fragment. All four anthracyclines show preferential intercalation at sites near the ends of both hexanucleotide duplexes, partly as a result of positive cooperativity in the formation of di-intercalated species at these sites. Within the limits of experimental error, complete site specificity for the CpG fragment is found in the intercalation of 1 and 2 into d(CGTACG) duplex, whereas analogues 3 and 4 give increasing evidence of intercalation at other sites including the fluorescence-preserving TpA fragment. Site specificity is less pronounced in the association with d(CGCGCG), when cooperativity is taken into account. Kinetic data corroborate the results of equilibrium studies and are interpreted with a mechanism that includes formation of an intermediate bound species followed by drug redistribution to preferential sites. Finally, from a comparison of pertinent site binding constants, approximate free energy contributions to sequence specific DNA interaction, due to C9-OH on the aglycone and -NH3+ on daunosamine, are estimated not to exceed 2 kcal/mol.     

DNA bis-intercalation: Theoretical calculation of binding curves

A statistical mechanical calculation of the binding properties of DNA bis-intercalators is presented, based on the sequence-generating function method of Lifson. The effects of binding by intercalation of one or both chromophores of a bifunctional intercalating agent are examined. The secular equation for a general model that includes the effects of neighbor (nearest and non-nearest) exclusion and/or cooperativity in the binding of both singly and doubly intercalated ligands is derived. Numerical results for binding curves are presented for a more restricted model in which each type of bound ligand rigorously excludes its nearest neighbor and the total number of sites covered by a doubly intercalated ligand is variable. At low values of free ligand concentration bis-intercalation dominates the binding process, while at high value of free ligand concentration, intercalation of only one chromophore per ligand becomes significant due to the unavailability of contiguous free sites required for bis-intercalation. Also, depending on the binding parameters, the free energy of the system can be lowered by a loss of doubly intercalated ligands in favor of singly intercalated ones. Corresponding to this transition in binding mode, the average number of sites occupied by a bound ligand decreases from that characteristic of bis-intercalation to that characteristic of mono-intercalation as free ligand concentration increases. An analysis of Scatchard plots describing bis-intercalation is presented.     

DNA sequence recognition by under-methylated analogues of triostin A.

Two new analogues of TANDEM (des-N-tetramethyl triostin A) have been synthesised in an effort to elucidate the molecular basis of DNA nucleotide sequence recognition in this series of compounds. Their binding preferences have been investigated by DNAase I footprinting and differential inhibition of restriction nuclease attack. The presence of a single N-methyl group on only one valine residue (in [N-MeVal4] TANDEM) abolishes the ability to recognise DNA, presumably because this antibiotic analogue has suffered an unfavourable conformational change in the depsipeptide ring. A bis-methylated analogue, [N-MeCys3, N-MeCys7]TANDEM, was found to interact quite strongly with DNA and afforded binding sites, rich in AT residues, identical to those of TANDEM. Footprinting with various DNA fragments of known sequence showed that this analogue recognises sequences containing the dinucleotide TpA, although we cannot exclude the possibility that it binds to ApT as well. [N-MeCys3, N-MeCys7]TANDEM inhibits cutting by RsaI, a restriction enzyme that recognises GTAC but not by Sau3AI which recognises GATC. This provides further supportive evidence that the ligand (and, by extension, TANDEM itself) prefers binding to sequences containing the dinucleotide step TpA.     

Interaction of Hoechst 33258 and echinomycin with nucleosomal DNA fragments containing isolated ligand binding sites

We have examined the interaction of Hoechst 33258 and echinomycin with nucleosomal DNA fragments which contain isolated ligand binding sites. A 145 base pair fragment was prepared on the basis of the sequence of tyrT DNA, which contained no CpG or (A/T)(4) binding sites for these ligands. Isolated binding sites were introduced into this fragment at discrete locations where the minor groove is known to face toward or away from the protein core when reconstituted onto nucleosome core particles. The interaction of ligands with target sites on these nucleosomal DNA fragments was assessed by DNase I footprinting. We find that Hoechst 33258 can bind to single nucleosomal sites which face both toward and away from the protein core, without affecting the nucleosome structure. Hoechst binding is also observed on nucleosomal fragments which contain two or more drug binding sites, though in these cases the footprints are accompanied by the presence of new cleavage products in positions which suggest that the ligand has caused a proportion of the DNA molecules to adopt a new rotational positioning on the protein surface. Hoechst 33258 does not affect nucleosome reconstitution with any of these fragments. In contrast, the bifunctional intercalating antibiotic echinomycin is not able to bind to single nucleosomal CpG sites. Echinomycin footprints are observed on nucleosomal fragments containing two or more CpG sites, but there are no changes in the cleavage patterns in the remainder of the fragment. Echinomycin abolishes nucleosome reconstitution when included in the reconstitution mixture.     

Probing the conformations of eight cloned DNA dodecamers; CGCGAATTCGCG, CGCGTTAACGCG, CGCGTATACGCG, CGCGATATCGCG, CGCAAATTTGCG, CGCTTTAAAGCG, CGCGGATCCGCG and CGCGGTACCGCG.

The self complementary DNA dodecamers d(CGCGAATTCGCG), d(CGCGTTAACGCG), d(CGCGTATACGCG), d(CGCGATATCGCG), d(CGCAAATTTGCG), d(CGCTTTAAAGCG), d(CGCGGATCCGCG) and d(CGCGGTACCGCG) have been cloned into the Smal site of plasmid pUC19. Radiolabelled polylinker fragments containing these inserts have been digested with nucleases and chemical agents, probing the structure of the central AT base pairs. The sequences AATT and AAATTT are relatively resistant to digestion by DNase I, micrococcal nuclease and hydroxyl radicals, consistent with the suggestion that they possess a narrow minor groove. Nuclease digestion of TTAA is much more even, and comparable to that at mixed sequence DNA. TpA steps in ATAT, TATA and GTAC are cut less well by DNAse I than in TTAA. DNasel cleavage of surrounding bases, especially CpG is strongly influenced by the nature of the central sequence.     

The effects of actinomycin on the structure of dAn.dTn and (dA-dT)n regions surrounding its GC binding site. A footprinting study.

The effect of actinomycin on the structure of DNA fragments containing the sequences (AT)5GC(AT)5, (TA)5GC(TA)5, A9GCT9, and T9GCA9, cloned into the SmaI site of pUC19, has been studied by footprinting analysis using a variety of probes known to be sensitive to DNA structure. In each case clear footprints are found around the central GC sites. DNase I cleavage of fragments containing alternating AT shows much greater cutting at ApT than TpA; in the presence of actinomycin, although this preference is retained, there is a large increase in the cutting efficiency at the closest TpA steps. DNase I cleavage in homopolymeric regions of A and T, which is normally very poor, is greatly enhanced by drug binding. With T9GCA9 the enhancements are propagated in both directions, whereas changes are only found to the 5'-side of the GC site in A9GCT9. The results are confirmed by similar experiments with micrococcal nuclease and DNase II. Small increases in sensitivity to diethylpyrocarbonate are found at adenines proximal to GC. Experiments performed at 4 degrees C suggest that conformational changes are a necessary consequence of drug binding.     

Hydration of nucleic acid fragments: comparison of theory and experiment for high-resolution crystal structures of RNA, DNA, and DNA-drug complexes.

A computationally efficient method to describe the organization of water around solvated biomolecules is presented. It is based on a statistical mechanical expression for the water-density distribution in terms of particle correlation functions. The method is applied to analyze the hydration of small nucleic acid molecules in the crystal environment, for which high-resolution x-ray crystal structures have been reported. Results for RNA [r(ApU).r(ApU)] and DNA [d(CpG).d(CpG) in Z form and with parallel strand orientation] and for DNA-drug complexes [d(CpG).d(CpG) with the drug proflavine intercalated] are described. A detailed comparison of theoretical and experimental data shows positional agreement for the experimentally observed water sites. The presented method can be used for refinement of the water structure in x-ray crystallography, hydration analysis of nuclear magnetic resonance structures, and theoretical modeling of biological macromolecules such as molecular docking studies. The speed of the computations allows hydration analyses of molecules of almost arbitrary size (tRNA, protein-nucleic acid complexes, etc.) in the crystal environment and in aqueous solution.     

CpG Dinucleotide Frequencies Reveal the Role of Host Methylation Capabilities in Parvovirus Evolution

Parvoviruses are rapidly evolving viruses that infect a wide range of hosts, including vertebrates and invertebrates. Extensive methylation of the parvovirus genome has been recently demonstrated. A global pattern of methylation of CpG dinucleotides is seen in vertebrate genomes, compared to “fractional” methylation patterns in invertebrate genomes. It remains unknown if the loss of CpG dinucleotides occurs in all viruses of a given DNA virus family that infect host species spanning across vertebrates and invertebrates. We investigated the link between the extent of CpG dinucleotide depletion among autonomous parvoviruses and the evolutionary lineage of the infected host. We demonstrate major differences in the relative abundance of CpG dinucleotides among autonomous parvoviruses which share similar genome organization and common ancestry, depending on the infected host species. Parvoviruses infecting vertebrate hosts had significantly lower relative abundance of CpG dinucleotides than parvoviruses infecting invertebrate hosts. The strong correlation of CpG dinucleotide depletion with the gain in TpG/CpA dinucleotides and the loss of TpA dinucleotides among parvoviruses suggests a major role for CpG methylation in the evolution of parvoviruses. Our data present evidence that links the relative abundance of CpG dinucleotides in parvoviruses to the methylation capabilities of the infected host. In sum, our findings support a novel perspective of host-driven evolution among autonomous parvoviruses.     

Interaction of minor groove binding ligands with long AT tracts.

We have used quantitative DNase I footprinting to examine the ability of distamycin and Hoechst 33258 to discriminate between different arrangements of AT residues, using synthetic DNA fragments containing multiple blocks of (A/T)6or (A/T)10in identical sequence environments. Previous studies have shown that these ligands bind less well to (A/T)4sites containing TpA steps. We find that in (A/T)6tracts distamycin shows little discrimination between the various sites, binding approximately 2-fold stronger to TAATTA than (TA)3, T3A3and GAATTC. In contrast, Hoechst 33258 binds approximately 20-fold more tightly to GAATTC and TAATTA than T3A3and (TA)3. Hydroxyl radical footprinting reveals that both ligands bind in similar locations at the centre of each AT tract. At (A/T)10sites distamycin binds with similar affinity to T5A5, (TA)5and AATT, though bands in the centre of (TA)5are protected at approximately 50-fold lower concentration than those towards the edges. Hoechst 33258 shows a similar pattern of preference, with strong binding to AATT, T5A5and the centre of (TA)5. Hydroxyl radical footprinting reveals that at low concentrations both ligands bind at the centre of (TA)5and A5T5, while at higher concentrations ligand molecules bind to each end of the (A/T)10tracts. At T5A5two ligand molecules bind at either end of the site, even at the lowest ligand concentration, consistent with the suggestion that these compounds avoid the TpA step. Similar DNase I footprinting experiments with a DNA fragment containing T n (n = 3-6) tracts reveals that both ligands bind in the order T3< T4 << T5 = T6.     

Role of stacking interactions in the binding sequence preferences of DNA bis-intercalators: insight from thermodynamic integration free energy simulations

The major structural determinant of the preference to bind to CpG binding sites on DNA exhibited by the natural quinoxaline bis-intercalators echinomycin and triostin A, or the quinoline echinomycin derivative, 2QN, is the 2-amino group of guanine (G). However, relocation of this group by means of introduction into the DNA molecule of the 2-aminoadenine (=2,6-diaminopurine, D) base in place of adenine (A) has been shown to lead to a drastic redistribution of binding sites, together with ultratight binding of 2QN to the sequence DTDT. Also, the demethylated triostin analogs, TANDEM and CysMeTANDEM, which bind with high affinity to TpA steps in natural DNA, bind much less tightly to CpI steps, despite the fact that both adenosine and the hypoxanthine-containing nucleoside, inosine (I), provide the same hydrogen bonding possibilities in the minor groove. To study both the increased binding affinity of 2QN for DTDT relative to GCGC sites and the remarkable loss of binding energy between CysMeTANDEM and ICIC compared with ATAT, a series of thermodynamic integration free energy simulations involving conversions between DNA base pairs have been performed. Our results demonstrate that the electrostatic component of the stacking interactions between the heteroaromatic rings of these compounds and the bases that make up the intercalation sites plays a very important role in the modulation of their binding affinities.