Theoretical study of binding affinity for diamidine with DNA

Diamidine molecules, which have been recognized as the powerful gene drug candidates over the past decades, can bind in the DNA minor groove, inhibit the duplication of morbid sequences, and fight against a number of human and animal diseases. In this paper, on the basis of the binding models of a series of diamidines with DNA, the important influencing factors for the binding affinity of diamidines with DNA were systematically analyzed. The obtained results demonstrated that the curvature, length, distal group, and heteroaromatic ring of diamidine are four important factors, which could influence their binding affinities. Specifically, the better the curvature of the diamidine fits DNA minor groove, the higher the binding affinity is; increasing the molecular length within a certain range can make the binding affinity higher; changing the distal group of diamidine from amidino to imidazole or pyrimidine is favorable for improving the corresponding binding affinity; and the introduction of central heteroaromatic rings of diamidine molecules influences their binding affinities. One diamidine (named as DB103d) with ideal DNA binding affinity validates the four important factors proposed in the present work. The results obtained in this work might be helpful for the design of new efficient diamidine-based drug candidates.

     

In situ analysis of cisplatin binding to DNA: the effects of physiological ionic conditions

Platinum-based anti-cancer drugs form a major family of cancer chemotherapeutic agents. Cisplatin, the first member of the family, remains a potent anti-cancer drug and exhibits its clinical effect by inducing local DNA kinks and subsequently interfering with DNA metabolism. Although its mechanism is reasonably well understood, effects of intracellular ions on cisplatin activity are left to be elucidated because cisplatin binding to DNA, thus its drug efficacy, is modified by various ions. One such issue is the effect of carbonate ions: cisplatin binding to DNA is suppressed under physiological carbonate conditions. Here, we examined the role of common cellular ions (carbonate and chloride) by measuring cisplatin binding in relevant physiological buffers via a DNA micromanipulation technique. Using two orthogonal single-molecule methods, we succeeded in detecting hidden monofunctional adducts (kink-free, presumably clinically inactive form) and clearly showed that the major effect of carbonates was to form such adducts and to prevent them from converting to bifunctional adducts (kinked, clinically active). The chloride-rich environment also led to the formation of monofunctional adducts. Our approach is widely applicable to the study of the transient behaviours of various drugs and proteins that bind to DNA in different modes depending on various physical and chemical factors such as tension, torsion, ligands, and ions.     

Theoretical study of copper binding to GHK peptide

We report ligand field molecular mechanics, density functional theory and semi-empirical studies on the binding of Cu(II) to GlyHisLys (GHK) peptide. Following exhaustive conformational searching using molecular mechanics, we show that relative energy and geometry of conformations are in good agreement between GFN2-xTB semi-empirical and B3LYP-D DFT levels. Conventional molecular dynamics simulation of Cu-GHK shows the stability of the copper-peptide binding over 100 ps trajectory. Four equatorial bonds in 3N1O coordination remain stable throughout simulation, while a fifth in apical position from C-terminal carboxylate is more fluxional. We also show that the automated conformer and rotamer search algorithm CREST is able to correctly predict the metal binding position from a starting point consisting of separated peptide, copper and water.     

Theoretical study of the interaction pattern and the binding affinity between procaine and DNA bases

The interacting patterns and mechanism of the binding affinity between the local anaesthetic procaine and four DNA bases (adenine, cytosine, guanine and thymine) in neutral form have been investigated in gas phase using the Austin Model l and density functional methods. The results show that the complexes are mainly stabilized by the H-bonding interactions. The bond critical point properties of the optimized complexes were analyzed by using the atoms in molecules theory with DFT method and the results show that the presence of the C?H···O or C?H···N hydrogen bonding. The natural bond orbital analysis was performed to quantitatively evaluate the hydrogen bonding interaction. The interacting energy shows that the binding of procaine with guanine is the most strong, whereas its binding to cytosine exhibits relatively weaker stability. The strength order of the relevant transferred charge between procaine and DNA base with natural population analysis are consist with the HOMO–LUMO gap results for each complex. And the order is accord with the relevant electrochemical experimental results.     

Theoretical prediction of the elastic contribution to steric stabilization

Summary It is shown how it is possible to use the Flory network theory to develop simple analytic formulae for the elastic contribution (first proposed byJäckel) to the steric stabilization by macromolecules. Constant segment density and gaussian segment density distribution models are discussed. It seems likely that the elastic contribution is negligible in the interpenetration domain but becomes important as significant compression occurs. Incipient flocculation is inferred to occur in a pseudo-secondary minimum, the depth of which is almost independent of any elastic contribution.

---

Zusammenfassung Es wird gezeigt, daß die Netzwerktheorie vonFlory benutzt werden kann, einfache analytische Ausdrücke zu entwickeln für den elastischen Anteil der sterischen Stabilisierung durch Makromoleküle. Konstante Segmentdichte und Gaußverteilung der Segmentdichte werden vorausgesetzt. Es ist wahrscheinlich, daß der elastische Anteil im Bereich der wechselseitigen Durchdringung vernachlässigt werden kann, aber wirksam wird bei merklicher Kompression. Die Anfangsstadien der Flockung hängen mit dem pseudosekundären Minimum zusammen, dessen Tiefe nahezu unabhängig ist von jedem elastischen Beitrag.

---

---

With 7 figures     

吡啶-BH~3相互作用复合物的理论研究

对吡啶-BH~3复合物分别用MP2/6-31+G^*和B3LYP/6-31+G^*进行理论计算以预测该复合物的构型及解离能，得到四种构型，在MP2优化构型基础上作CCSD/6-31+G^*单点能量计算以验证MP2与B3LYP结果的可靠性，然后用B3LYP作振动频率分析，计算了各构型的垂直电离势，最后用更大基组作单点能量计算和自然键轨道(NBO)分析。结果表明，N-B直接相连的构型最稳定，其解离能为141.50kJ/mol,MP2和B3LYP对N-H接近的构型结果相关较大，另外两种构型稳定性介于二者之间，解离能分别为15.18kJ/mol,14.06kJ/mol(MP2/6-31+G^*)。     

Protein-DNA interaction: impedance study of MutS binding to a DNA mismatch

MutS binding to a double-stranded DNA containing a single nucleotide mismatch can be conveniently monitored by impedance spectroscopy and represents the first step in developing an electrochemical binding essay for single nucleotide mismatch detection.     

A model study of the laser-induced stabilization of a collision complex

The possibility of stabilizing a collision complex by stimulated emmission with coherent radiation is explored for a model. three-atom system. The lifetimes of several resonance states of this system were obtained previously and are used in a simple two-level model to describe the interaction of radiation with a resonance state and a true bound state. For the case where the energy difference between the two levels is time independent, simple, analytical expressions are used to obtain the time-dependent transition probability. Numerical solutions are obtained for the more realistic case where this energy difference is time dependent. For both cases, a substantial transition probability for stabilization is calculated for moderate laser intensities (e.g. 8.1 kW/cm² for an assumed dipole matrix element of 0.01 D).     

Theoretical study of cooperativity in multivalent polymers for colloidal stabilization

Multivalent polymers, i.e., copolymers with multiple binding sites, have been proposed recently for stabilization of fusogentic liposomes and other liposomal colloids useful for drug delivery. The performance of such polymers critically depends on their molecular architecture, in particular the strength and frequency of surface anchoring sites along the backbone of a highly soluble polymer. In this work, we investigate the adsorption and surface forces due to multivalent polymers based on coarse-grained polymer models. We find that for W-type polymers that form dangling tails when all anchoring segments are attached to a surface, increasing the chain length at fixed polymer composition leads to a stronger repulsive barrier in the polymer-mediated surface forces thereby increasing the ability of the polymer to stabilize colloidal particles. This prediction conforms to an earlier experiment indicating that increasing the number of hydrophobic anchors along poly(ethylene glycol) polymers results in the cooperative behavior for both surface adsorption and steric stabilization. For M-type multivalent polymers that have weakly anchoring sites placed at the ends, however, addition of binding sites at fixed polymer composition could lead to negative cooperativity, i.e., the more binding sites, the less the amount of adsorption or the weaker the ability of surface protection. The theory also predicts that polymers with two anchoring sites (e.g., telechelic copolymers) are most efficient for colloidal stabilization.     

Relationship of solution and protein-bound structures of DNA duplexes with the major intrastrand cross-link lesions formed on cisplatin binding to DNA.

DNA bases in the three-base-pair (3bp) region of duplexes with the two major lesions of cisplatin (cis-PtCl(2)(NH(3))(2)) with DNA, namely d(XGG) and d(XAG) ( = N7-platinated base), differ in their relative positions by as much as approximately 3.5 A in structures in the literature. Such large differences impede drug design and assessments of the effects of protein binding on DNA structure. One recent and several past structures based on NMR-restrained molecular dynamics (RMD) differ significantly from the reported X-ray structure of an HMG-bound XGG 16-mer DNA duplex (Ohndorf, U.-M.; Rould, M. A.; He, Q.; Pabo, C. O.; Lippard, S. J. Nature 1999, 399, 708). This 16-mer structure has several significant novel and unique features (e.g., a bp step with large positive shift and slide). Hypothesizing that novel structural features in the XGG or XAG region of duplexes elude discovery by NMR methods (especially because of the flexible nature of the 3bp region), we studied an oligomer with only G.C bp's in the XGGY site by NMR methods for the first time. This 9-mer gave a 5'-G N1H signal with a normal shift and intensity and showed clear NOE cross-peaks to C NHb and NHe. We assigned for the first time (13)C NMR signals of a duplex with a GG lesion. These data, by adding NMR-based criteria to those inherent in NOESY and COSY data, have more specifically defined the structural features that should be present in an acceptable model. In particular, our data indicated that the sugar of the X residue has an N pucker and that the GG cross-link should have a structure similar to the original X-ray structure of cis-Pt(NH(3))(2)(d(pGpG)) (Sherman S. E.; Gibson, D.; Wang, A. H.-J.; Lippard, S. J. J. Am. Chem. Soc. 1988, 110, 7368). With these restrictions added to NOE restraints, an acceptable model was obtained only when we started our modeling with the 16-mer structural features. The new X-ray/NMR-based model accounted for the NOESY data better than NOE-based models, was very similar in structure to the 16-mer, and differed from solely NOE-based models. We conclude that all XGG and XAG (X = C or T) duplexes undoubtedly have structures similar to those of the 16-mer and our model. Thus, protein binding does not change greatly the structure of the 3bp region. The structure of this region can now be used in understanding structure-activity relationships needed in the design of new carrier ligands for improving Pt anticancer drug activity.     

Slowing of cisplatin aquation in the presence of DNA but not in the presence of phosphate: improved understanding of sequence selectivity and the roles of monoaquated and diaquated species in the binding of cisplatin to DNA

1H-15N HSQC NMR spectroscopy is used to study the aquation reactions of cisplatin in 9 mM NaClO4 and 9 mM phosphate (pH 6) solutions at 298 K. For the first time in a single reaction and, therefore, under a single set of reaction conditions, the amounts of all species formed are followed and the rates of aquation, diaquation, and related anation processes are determined in both media. Binding of phosphate to aquated Pt species is observed, but the initial rate of aquation is not affected by the presence of 9 mM phosphate. The reaction between cisplatin and the 14-base-pair self-complementary oligonucleotide 5'-d(AATTGGTACCAATT)-3', having a GpG intrastrand binding site, is investigated. Various kinetic models for this reaction are evaluated and the most appropriate found to be that with a reversible aquation step and a single binding site for the self-complementary duplex. The rate constant for aquation is (1.62 +/- 0.02) x 10(-5) s-1, with the anation rate constant fixed at 4.6 x 10(-3) M-1 s-1, the value obtained from the aquation studies. The rate constants for monofunctional binding of cis-[PtCl(15NH3)2-(OH2)]+ to the sequence were 0.48 +/- 0.19 and 0.16 +/- 0.06 M-1 s-1 for the 3'- and 5'-guanine bases, respectively. Closure rate constants for the monofunctional adducts are (2.55 +/- 0.07) x 10(-5) and (0.171 +/- 0.011) x 10(-5) s-1, for the 3'- and 5'-guanines, respectively. The presence of DNA slows the aquation of cisplatin by 30-40% compared to that observed in 9 mM NaClO4 or 9 mM phosphate, and there is some evidence that the degree of slowing is sequence dependent. The possibility that cis-[Pt(OH)(NH3)2(OH2)]+ contributes to the binding of cisplatin to DNA is investigated, and it is found that about 1% followed this route, the majority of the binding occurring via the monoaquated species cis-[PtCl(NH3)2(OH2)]+. Comparison of the rates of disappearance of cisplatin in reactions at single defined GpG, ApG, GpA, GpTpG and 1,2-interstrand GG binding sites shows that the adduct profile is determined at the level of monofunctional adduct formation.     

Theoretical study of the amino acid interactions in the binding region of the trypsin–pancreatic trypsin inhibitor complex

The meaningful interactions in the contact region of the trypsin–pancreatic trypsin inhibitor complex have been evaluated using free energy simulation methods and appropriate thermodynamic cycles. Consequently, mutations on a few selected residues were performed to destroy the specificity of these interactions preserving the three-dimensional structure of the original species; therefore, the original and mutated residues involved had to be similar from a topological point of view but with opposite chemical properties. The thermodynamic perturbation, conventional thermodynamic integration, thermodynamic integration with end-point perturbational correction and thermodynamic integration conjugated with the extrapolation method of Brooks formalisms have been used in these calculations. The results obtained, with these four alternative approaches, are in reasonable, mutual agreement and reveal that the strength of a specific interaction increases with the charge separation between the amino acid residues involved, and with the hydrophilicity of the surrounding environment.     

The effect of DNA binding on initial 8-methoxypsoralen photochemistry.

Abstract—Flash photolysis experiments on the complex of 8-methoxypsoralen with calf thymus DNA have shown binding promotes the anaerobic quenching of the 8-MOP triplet state and inhibits the accessibility to oxygen. The results indicate dark psoralen-DNA interactions suppress the singlet oxygen generation observed in prior work.     

Thermodynamics of folding, stabilization, and binding in an engineered protein--protein complex

We analyzed the thermodynamics of a complex protein-protein binding interaction using the (engineered) Z(SPA)(-)(1) affibody and it's Z domain binding partner as a model. Free Z(SPA)(-)(1) exists in an equilibrium between a molten-globule-like (MG) state and a completely unfolded state, wheras a well-ordered structure is observed in the Z:Z(SPA)(-)(1) complex. The thermodynamics of the MG state unfolding equilibrium can be separated from the thermodynamics of binding and stabilization by combined analysis of isothermal titration calorimetry data and a separate van't Hoff analysis of thermal unfolding. We find that (i) the unfolding equilibrium of free Z(SPA)(-)(1) has only a small influence on effective binding affinity, that (ii) the Z:Z(SPA)(-)(1) interface is inconspicuous and structure-based energetics calculations suggest that it should be capable of supporting strong binding, but that (iii) the conformational stabilization of the MG state to a well-ordered structure in the Z:Z(SPA)(-)(1) complex is associated with a large change in conformational entropy that opposes binding.     

An electrochemical study of safranin O binding to DNA at the surface

The interaction of safranin O (SO) with double-stranded calf thymus DNA was investigated electrochemically, using a DNA-modified glassy carbon (GC-DNA) electrode. The results were compared with those obtained using a bare GC electrode. The formal potential of SO was more negative when using the GC-DNA electrode, although the rate of heterogeneous electron transfer was not altered. The GC-DNA electrode enabled preconcentration of the SO on the electrode surface, despite the fact that the mass transfer effects in the thin DNA layer adsorbed on the surface were still observed. The diffusion coefficient of SO and the binding ratio for the oxidized and reduced forms of the bound species were obtained. A binding isotherm for SO at the GC-DNA electrode was plotted from coulometric titrations, giving a binding constant of 5.8×10⁴ L mol⁻¹.