Protein environment affects the water-tryptophan binding mode. MD, QM/MM, and NMR studies of engrailed homeodomain mutants.

Water molecules can interact with aromatic moieties using either their O-H bonds or their lone-pairs of electrons. In proteins, water-π interactions have been reported to occur with tryptophan and histidine residues, and dynamic exchange between O-Hπ hydrogen bonding and lone-pairπ interactions was suggested to take place, based on ab initio calculations. Here we used classical and QM/MM molecular dynamics simulations, complemented with an NMR study, to examine a specific water-indole interaction observed in the engrailed homeodomain and in its mutants. Our simulations indicate that the binding mode between water and indole can adapt to the potential created by the surrounding amino acids (and by the residues at the DNA surface in protein-DNA complexes), and support the model of dynamic switching between the O-Hπ hydrogen bonding and lone-pairπ binding modes.

Water-Tryptophan Interactions: Lone-pair⋠⋠⋠π or O-H⋠⋠⋠π? Molecular Dynamics Simulations of ß-Galactosidase Suggest that Both Modes Can Co-exist.

In proteins, the indole side chain of tryptophan can interact with water molecules either in-plane, forming hydrogen bonds, or out-of-plane, with the water molecule contacting the aromatic π face. The latter interaction can be either of the lone pair⋅⋅⋅π (lp⋅⋅⋅π) type or corresponds to the O-H⋅⋅⋅π binding mode, an ambiguity that X-ray structures usually do not resolve. Here, we report molecular dynamics (MD) simulations of a solvated ß-galactosidase monomer, which illustrate how a water molecule located at the π face of an indole side chain of tryptophan can adapt to the position of proximate residues and "select" its binding mode. In one such site, the water molecule is predicted to rapidly oscillate between the O-H⋅⋅⋅π and lp⋅⋅⋅π binding modes, thus gaining entropic advantage. Our MD simulations provide support for the role of lp⋅⋅⋅π interactions in the stabilization of protein structures.

Lone pair-π interactions in biological systems: occurrence, function, and physical origin.

Lone pair-π interactions are now recognized as a supramolecular bond whose existence in biological systems is documented by a growing number of examples. They are commonly attributed to electrostatic forces. This review attempts to highlight some recent discoveries evidencing the important role which lone pair-π interactions, and anion-π interactions in particular, play in stabilizing the structure and affecting the function of biomolecules. Special attention is paid to studies exploring the physical origin of these at first glance counterintuitive interactions between a lone pair of electrons of one residue and the π-cloud of another. Recent theoretical work went beyond the popular electrostatic model and inquired the extent to which orbital interactions have to be taken into account. In at least one biologically relevant case-that of anion-flavin interactions-a substantial charge-transfer component has been shown to operate.

Anion-π Interactions in Flavoproteins Involve a Substantial Charge-Transfer Component.

Anion-π interactions have been shown to stabilize flavoproteins and to regulate the redox potential of the flavin cofactor. They are commonly attributed to electrostatic forces. Herein we show that anion-flavin interactions can have a substantial charge-transfer component. Our conclusion emanates from a multi-approach theoretical analysis and is backed by previously reported observations of absorption bands, originating from charge transfer between oxidized flavin and proximate cysteine thiolate groups. This partial covalency of anion-flavin contacts renders classical simulations of flavoproteins questionable.

Lone-pair-π interactions: analysis of the physical origin and biological implications.

Lone-pair-π (lp-π) interactions have been suggested to stabilize DNA and protein structures, and to participate in the formation of DNA-protein complexes. To elucidate their physical origin, we have carried out a theoretical multi-approach analysis of two biologically relevant model systems, water-indole and water-uracil complexes, which we compared with the structurally similar chloride-tetracyanobenzene (TCB) complex previously shown to contain a strong charge-transfer (CT) binding component. We demonstrate that the CT component in lp-π interactions between water and indole/uracil is significantly smaller than that stabilizing the Cl(-)-TCB reference system. The strong lp(Cl(-))-π(TCB) orbital interaction is characterized by a small energy gap and an efficient lp-π* overlap. In contrast, in lp-π interactions between water and indole or uracil, the corresponding energy gap is larger and the overlap less efficient. As a result, water-uracil and water-indole interactions are weak forces composed by smaller contributions from all energy components: electrostatics, polarization, dispersion, and charge transfer. In addition, indole exhibits a negative electrostatic potential at its π-face, making lp-π interactions less favorable than O-Hπ hydrogen bonding. Consequently, some of the water-tryptophan contacts observed in X-ray structures of proteins and previously interpreted as lp-π interactions [Luisi, et al., Proteins, 2004, 57, 1-8], might in fact arise from O-Hπ hydrogen bonding.

On the non-classical contribution in lone-pair-π interaction: IQA perspective.

In the present work the nature of lone-pair-π interactions between water molecules and a number of π-rings with different substituents/hetero-atoms in the light of quantum chemical topology approaches is studied. The Quantum Theory of Atoms in Molecules (QTAIM) and Interacting Quantum Atoms (IQA) were employed for distinguishing the role of heteroatoms and electron withdrawing substituents in the complex formation between water and π-rings. Our IQA study identified three classes of water-π complexes on the basis of the relative role of electrostatics (classical) and exchange-correlation (non-classical) factors in the interaction energy between the oxygen of water (the lone-pair donor) and the sp(2) atoms of the π-ring, i.e. the primary lp-π interaction. Considering both the primary and secondary (the rest of interatomic interactions except Owater-π-ring atoms) interactions demonstrates that the exchange-correlation is the dominant contributor to the binding energy. This proves a non-negligible contribution of non-classical factors in the stabilization of the lone-pair-π complexes. However, in spite of a relatively large contribution of the exchange-correlation, this part of the interaction energy is virtually counterbalanced by the deformation energy, i.e. the increase in atomic kinetic energy upon complexation. This finding clarifies why water-π interactions can be modelled by simple electrostatics without the need to invoke quantum effects.

Quantum chemical topology study of the water-platinum(II) interaction.

The "inverse hydration" of neutral complexes of Pt(II) by an axial water molecule, whose one OH-bond is oriented toward Pt, has been the subject of recent works, theoretical as well as experimental. To study the influence of the ligands on this non-conventional H-bond, we extend here our previous energy calculations, using the second-order Moeller-Plesset perturbation theory (MP2) method together with the Dolg-Pélissier pseudopotential for platinum, to various neutral complexes including the well-known chemotherapeutic agent "cisplatin". The stabilization energy, depending on the nature and the configuration of platinum ligands, is dominated by the same important dispersive component, for all the investigated complexes. For a further characterization of this particular H-bond, we used the atoms in molecules theory (AIM) and the topological analysis of the electron localization function (ELF). The charge transfer occurring from the complex to the water molecule and the Laplacian of the density at the bond critical point between water and Pt are identified as interesting AIM descriptors of this non-conventional H-bond. Beyond this AIM analysis, we show that the polarization of the ELF bonding O-H basin involved in the non-conventional H-bond is enhanced during the approach of the water molecule to the Pt complexes. When the water medium, treated in an implicit solvation model, is taken into account, the interaction energies become independent on the nature and configuration of platinum ligands. However, the topological descriptors remain qualitatively unchanged.

Cisplatin GG-crosslinks within single-stranded DNA: origin of the preference for left-handed helicity.

Molecular dynamics (MD) simulations of the single-stranded DNA trinucleotide TG*G*, with the G* guanines crosslinked by the antitumor drug cisplatin, were performed with explicit representation of the water as solvent. The purpose of the simulations was to explain previous NMR observations indicating that in single-stranded cisplatin-DNA adducts, the crosslinked guanines adopt a left-handed helical orientation, whereas in duplexes, the orientation is right-handed. The analysis of the MD trajectory of TG*G* has ascribed a crucial role to hydrogen-bonding (direct or through-water) interactions of the 5'-oriented NH(3) ligand of platinum with acceptor groups at the 5'-side of the crosslink, namely the TpG* phosphate and the terminal 5'-OH group. These interactions bring about some strain into the trinucleotide which is slightly but significantly (1-1.5 kcal.mol(-1)) higher for the right-handed orientation than for the left-handed one. During the unconstrained, 3 ns long MD simulation, left-handed conformations were ~15 times more abundant than the right-handed ones. This sampling difference agrees roughly with the calculated energy difference in strain energy. Overall, these results show that the Pt-GG crosslink within single-stranded DNA is malleable and can access different conformations at a moderate energy cost. This malleability could be of importance in interactions between the platinated DNA and cellular proteins, in which the DNA is locally unwound.

Platinum-DNA interstrand crosslinks: molecular determinants of bending and unwinding of the double helix.

Platinum diamine complexes are able to crosslink the guanines of d(GC)(2) dinucleotides within double-stranded DNA. The interstrand crosslink thus formed causes a bend of the double helix toward the minor groove and the helical sense changes locally to left-handed, resulting in a considerable unwinding. The bend and unwinding angles have been shown to depend on the platinum ligands. Here, we have used molecular dynamics simulations to investigate the DNA 20-mer d(C(1)T(2)C(3)T(4)C(5)C(6)T(7)T(8)G*(9)C(10)T(11)C(12)T(13)C(14)C(15)T(16)T(17)C(18)T(19)C(20))-d(G(21)A(22)G(23)A(24)A(25)G(26)G(27)A(28)G(29)A(30)G*(31)C(32)A(33)A(34)G(35)G(36)A(37)G(38)A(39)G(40)) with the G* guanines crosslinked by cis-Pt(NH(3))(2)(2+), Pt(R,R-DACH)(2+), or Pt(S,S-DACH)(2+). Previous investigations on cisplatin interstrand adducts indicated that the structure is similar in solid state and in solution; thus, we used the reported X-ray structure of a cisplatin adduct as a starting model. Replacing in the MD-relaxed model for the DNA duplex crosslinked with cis-Pt(NH(3))(2)(2+) the two NH(3) platinum ligands by R,R-DACH or S,S-DACH led to clashes between the DACH residue and the deoxyribose of C(12). Confrontation of MD-derived models with gel shift measurements suggested that these clashes are avoided differently in the adducts of Pt(R,R-DACH)(2+)versus Pt(S,S-DACH)(2+). The R,R-isomer avoids the clash by untwisting the T(11)/A(30)-C(12)/G(29) step, thus increasing the global unwinding. In contrast, the S,S-isomer modifies the shift and slide parameters of this step, which dislocates the helical axis and enhances the bend angle. The clash that leads to the differentiation of the structures as a function of the diamine ligand is related to a hydrogen bond between the platinum complex and the T(11) base and could be characteristic of interstrand crosslinks at d(pyG*Cpy)-d(puG*Cpu) sequences.

Evaluation of dissociation constants from competition binding experiments based on the relative binding ratio.

Methods probing protein-DNA associations include direct binding titrations and competition binding experiments. For the latter, we present here a simple procedure allowing the quantitative evaluation of dissociation constants. We show that the ratio between the fraction of a DNA probe bound to protein in the absence of competitor and that in the presence of competitor is, at large competitor concentrations, a linear function of the competitor concentration, and we derive equations allowing the dissociation constant for the protein-competitor complex to be evaluated from the slope. We show further that a self-competition experiment, where the DNA probe and competitor are chemically the same species, can be used as a complement to a direct titration to determine the fraction of protein that is correctly folded for specific DNA binding. Thus, such a combination of direct and self-competition titration can be used as a check of the conformational purity of DNA binding proteins.

Cisplatin adducts on a GGG sequence within a DNA duplex studied by NMR spectroscopy and molecular dynamics simulations.

The antitumor drug cisplatin(cis-[PtCl2(NH3)2]) reacts with cellular DNA to form GG intrastrand adducts between adjacent guanines as predominant lesions. GGG sites have been shown to be hotspots of platination. To study the structural perturbation induced by binding of cisplatin to two adjacent guanines of a GGG trinucleotide,we examined here the decanucleotide duplex d[(G1C2C3G*4 G*5 G6T7-C8G9C10).d(G11C12G13A14C15C16C17G18-G19C20)] (dsCG*G*G) intrastrand cross-linked at the G* guanines by cis-{Pt(NH3)2}2+ using NMR spectroscopy and molecular dynamics (MD) simulations.The NMR spectra of dsCG*G*G were found to be similar to those of previously characterized DNA duplexes cross-linked by cisplatin at apyG*G*X site (py=pyrimidine; X=C,T, A). This similarity of NMR spectra indicates that the base at the 3'-side of the G*G*-Pt cross-link does not affect the structure to a large extent. An unprecedented reversible isomerization between the duplex dsCG*G*G (bearing a G*4 G*5 -Pt chelate) and duplex dsGG*G*T (bearing a G*5 G*6 -Pt chelate)was observed, which yielded a 40:60 equilibrium between the two intrastrand GG-Pt cross-links. No formation of interstrand cross-links was observed.NMR spectroscopic data of dsCG*G*G indicated that the deoxyribose of the 5'-G* adopts an N-type conformation, and the cytidines C3, C15,and C16 have average phase angles intermediate between S and N. The NMR spectroscopic chemical shifts of dsGG*G*T showed some fundamental differences to those of pyG*G*-platinum adducts but were in agreement with the NMR spectra reported previously for the DNA duplexes crosslinked at an AG*G*C sequence by cisplatin or oxaliplatin. The presence of apurine instead of a pyrimidine at the 5'-side of the G*G* cross-link seems therefore to affect the structure of the XG* step significantly.

Synthetic route to dinuclear platinum(II) complexes [{trans-PtCl(NH3)2}2(mu-L)]2+ (L = aliphatic or heterocyclic diamine) as potential antitumor agents, exploiting the mutual activation of hydroxido ligands and ammonium groups.

A simple and efficient method for the synthesis of potentially antitumor-active dinuclear platinum complexes of the general formula [{trans-PtCl(NH3)2}2(mu-L)]((n+2)+) (L = aliphatic or heterocyclic diamine; n = charge of L) is presented. The procedure is based on the mutual in situ activation of trans-[PtCl(OH)(NH3)2] and the linker L in the form of a diammonium salt. This synthetic pathway yielded the Farrell compound [{trans-PtCl(NH3)2}2{mu-NH2(CH2)6NH2}]Cl2 (BBR3005) in quantitative yield. Using the same procedure, we prepared the new pyrazolate-bridged compound [{trans-PtCl(NH3)2}2(mu-pz)]Cl, determined its X-ray structure, and tested its cytotoxicity against three wild-type and one cisplatin-resistant cell lines.

LEF-1 recognition of platinated GG sequences within double-stranded DNA. Influence of flanking bases.

The lymphoid enhancer-binding factor 1 (LEF-1) recognizes a double-stranded 9 base-pairs (bp) long motif in DNA which is significantly bent upon binding. This bend is centered at two destacked adenines whose geometry closely resembles that of two adjacent guanines crosslinked by the antitumor drug cisplatin. It has been proposed that cisplatin-GG crosslinks could hijack high mobility group (HMG) box containing transcription factors such as LEF-1. In order to examine such a possibility, we used electrophoretic mobility shift assays to determine the affinity of the HMG box of LEF-1 for a series of 25 oligonucleotides containing a central GG sequence, free or site-specifically modified by cisplatin. The binding affinity of the GG-platinated oligonucleotides was 3-6-fold higher than that determined for the corresponding unplatinated oligonucleotides, however, the binding to all cisplatin-modified oligonucleotides was at least 1 order of magnitude weaker than that to the 25 bp oligonucleotide containing the recognition 9 bp motif. The binding affinity was dependent on the nature of bases flanking the cisplatin-crosslinked G(*)G(*) dinucleotide, the AG(*)G(*)T sequence displaying the strongest affinity and CG(*)G(*)T showing the strongest binding enhancement upon platination. In contrast, modification of the AGGT sequence with the third-generation platinum antitumor drug oxaliplatin did not enhance the affinity significantly. These results suggest that the cisplatin-caused bending of DNA does produce a target for LEF-1 binding, however, the cisplatinated DNA does not appear to be a strong competitor for the LEF-1 recognition sequence.

Study of intramolecular competition between carboxylate and phosphonate for Pt(II) with the aid of a novel tridentate carboxylato-thioether-phosphonato ligand.

The tridentate dianionic ligand 2-[2'-(hydroxyisopropoxyphosphoryl)phenylsulfanyl]benzoate (L(2-)) reacts with cis-[Pt(NH(3))(2)(H(2)O)(2)](2+) to form an S,O-chelate in which the O-coordinated group is either carboxylate or phosphonate, depending on the degree of protonation of the complex. Carboxylate appears to be the stronger ligand, and the stoichiometric reaction between cis-[Pt(NH(3))(2)(H(2)O)(2)](2+) and L(2-) yields the neutral species [Pt(L)(NH(3))(2)], with L bound by sulfanyl and carboxylate groups, both in solution and in the solid state. Upon protonation of [Pt(L)(NH(3))(2)], the stronger basicity of the carboxylate causes the Pt coordination to switch from carboxylate to phosphonate, and the uncoordinated carboxylate group becomes protonated. In methanolic solution, the first-order kinetics of this rearrangement could be observed by (31)P NMR spectroscopy. Both complexes-the carboxylate-bound neutral complex [Pt(L)(NH(3))(2)].H(2)O (triclinic, P1 (no. 2), a=9.529(6), b=9.766(6), c=12.299(7) angstroms, alpha=106.91(2), beta=101.71(2), gamma=102.05(2) degrees, Z=2) and the perchlorate salt of the phosphonate-bound complex [Pt(LH)(NH(3))(2)]ClO(4).H(2)O (monoclinic, P2(1)/c (no. 14), a=12.095(2), b=14.046(2), c=14.448(2) angstroms, beta=95.55(2) degrees, Z=4)-were characterized by X-ray crystallography.

A pyrazolato-bridged dinuclear platinum(II) complex induces only minor distortions upon DNA-binding.

The cytotoxic, pyrazolato-bridged dinuclear platinum(II) complex [(cis-{Pt(NH3)2})2(mu-OH)(mu-pz)]2+ (pz=pyrazolate) has been found to cross-link two adjacent guanines of a double-stranded DNA decamer without destabilizing the duplex and without changing the directionality of the helix axis. A 1H NMR study of the oligonucleotide d(CTCTG*G*TCTC)-d(GAGACCAGAG), cross-linked at the two G* guanines by [(cis-{Pt(NH3)2})2(mu-pz)]3+, and molecular dynamics simulations of the explicitly solvated duplex were performed to characterize the structural details of the adduct. The dinuclear platinum cross-link unwinds the helix by approximately 15 degrees , that is, to a similar extent as the widely used antitumor drug cisplatin, but, in contrast to the latter, induces no significant bend in the helix axis. The Watson-Crick base-pairing remains intact, and the melting temperature of the duplex is unaffected by the cross-link. The helical twist is considerably reduced between the two platinated bases, as becomes manifest in an unusually short sequential H1'-H1' distance. This unwinding also affects the sugar ring of the guanosine in the 3'-position to the cross-link, which presents an N<-->S equilibrium. This is the first cytotoxic platinum complex that has been successfully designed by envisioning the structural consequences of its binding to DNA.

Fixing the conformations of diamineplatinum(II)-GpG chelates: NMR and CD signatures of individual rotamers.

The bulky, asymmetric analog of the antitumor drug cisplatin, [PtCl(2)(tmen)] (tmen = N,N,N'-trimethylethylenediamine), was used to produce crosslinks with the dinucleotide d(GpG), modeling the most frequent lesions that cisplatin and its analogs cause to DNA. The ligand tmen was chosen because it is expected to constrain the guanine cis to the NMe(2) group in the adduct [Pt(tmen){d(GpG)}](+) to an orientation perpendicular to the coordination plane and to stabilize the other guanine in an oblique orientation, thus maintaining a head-to-head geometry typical of cisplatin-d(GpG) crosslinks within single- and double-stranded DNA. Of the four possible combinations of tmen chirality (R or S symmetry of the coordinated NHMe group) and crosslink direction (5'-G bound cis to the secondary or the tertiary amino group of tmen), two isomers were preponderantly formed, [Pt(R-tmen){d(GpG)}](+) with 5'-G bound cis to NMe(2) and [Pt(S-tmen){d(GpG)}](+) with 5'-G bound cis to NHMe. The former was shown to have a right-handed R2 orientation of guanines similar to that found in duplex DNA, whereas the latter had a left-handed L1 orientation that modeled cisplatin-d(GpG) adducts within single-stranded DNA. The R2 rotamer was found to be in an equilibrium (as observed using EXSY spectroscopy) with a minor fraction (< or =4%) of a Delta-HT rotamer related to R2 by rotation of the 3'-G about the Pt-N7 bond. The major rotamers R2 and L1 were isolated using reverse-phase HPLC, and their NMR and CD signatures were compared to those of the corresponding rotamers of the less hindered adduct [Pt(dmen)(GpG)](+) (dmen = N,N-dimethylethylenediamine). From this and other comparisons with previously reported platinum dinucleotide complexes, and from molecular modeling, it could be concluded that both steric repulsion between guanine and substituents of the cis amino group and N-H...O6 hydrogen bonding are significant effects favoring the oblique orientation of one guanine base typical of the HH rotamers of [Pt(diamine){d(GpG)}](+) and [Pt(diamine)(GpG)](+) complexes.

Recognition complex between the HMG domain of LEF-1 and its cognate DNA studied by molecular dynamics simulations with explicit solvation.

Molecular dynamics simulations of the complex formed between the HMG box of the lymphoid enhancer-binding factor (LEF-1) and its cognate DNA duplex were carried out with explicit inclusion of water. The simulation started with an NMR-based model (pdb code 2LEF) and the dynamics was pursued for 10 nanoseconds without constraints. It revealed that water intervenes in many ionic/polar interactions, establishing in particular local equilibria between direct and water-mediated hydrogen bonds, and thus increasing the entropy of the complex. Quite unexpectedly, the simulation indicated that a binding pocket for a specific water molecule may be reversibly formed at the apex of the bend induced in the DNA helix by LEF-1 binding, where a methionine side chain intercalates between two destacked adenines. We observed that the specific water molecule can temporarily replace the intercalated S-CH(3) group, acting as a sort of "extension" of the side chain. The residence time of this water molecule was about 3.5 ns. Simulations of the cognate DNA alone showed that this sequence has no intrinsic tendency to bend; therefore, the bending occurs solely as a consequence of the recognition, following the "induced-fit" mechanism.

Influence of dT20 and [d(AT)10]2 on cisplatin hydrolysis studied by two-dimensional [1H,15N] HMQC NMR spectroscopy.

The influence of the presence of DNA on the kinetics of cisplatin (cis-[PtCl2(NH3)2]) aquation (replacement of Cl- by H2O) and anation (replacement of H2O by Cl-) involved in the hydrolysis of cisplatin have been determined by two-dimensional [1H,15N] HMQC NMR spectroscopy. Single-stranded dT20 and double-stranded [d(AT)10]2 oligonucleotides were used as DNA models, avoiding guanines which are known to react rapidly with aquated cisplatin forms. Reactions starting from cis-[PtCl2(15NH3)2], or from a stoichiometric mixture of cis-[Pt(15NH3)2(H2O)2]2+ and Cl- (all 0.5 mM Pt(II); in ionic strength, adjusted to 0.095 M or 0.011 M with NaClO4, pH between 3.0 and 4.0) were followed in an NMR tube in both the absence and presence of 0.7 mM dT20 or [d(AT)10]2. In the presence of dT20, we observed a slight and ionic-strength-independent decrease (15-20 %) of the first aquation rate constant, and a more significant decrease of the second anation rate constant. The latter was more important at low ionic strength, and can be explained by efficient condensation of cis-[Pt(15NH3)2(H2O)2]2+ on the surface of single-stranded DNA, in a region depleted of chloride anions. At low ionic strength, we observed an additional set of [1H,15N] HMQC spectral signals indicative of an asymmetric species of PtN2O2 coordination, and we assigned them to phosphate-bound monoadducts of cis-[Pt(15NH3)2(H2O)2]2+. Double-stranded [d(AT)10]2 slowed down the first aquation step also by approximately 15 %; however, we could not determine the influence on the second hydrolysis step because of a significant background reaction with cis-[Pt(NH3)2(H2O)2]2+.