Thermodynamics aspects of interactions between small ligands and DNA

DNA is a molecular target for many anticancer and antiviral drugs. Therefore, a clear understanding of the interaction of small molecules with DNA is important in the rational design of ligands that can bind to DNA with high affinity and selectivity. There are several methods to investigate interactions between drug and DNA. Some of them measures changing into DNA structures, such as lengthening and untwisting of helix of DNA. Other techniques measure changing in drug environment. With the increasing availability of sensitive microcalorimeters, particular interest has arisen in the thermodynamics of drug-DNA interaction. Using such methods permit direct determination of enthalpy changes associated with reactions. One experiment permits to obtain also binding constant, hence an almost complete thermodynamic profile can be established. This profile offers key insights into the molecular forces that drive complex formation and permit to estimate which kind of interaction are responsible of forming these complexes.     

Thermodynamics of water mediating protein-ligand interactions in cytochrome P450cam: a molecular dynamics study.

Ordered water molecules are observed by crystallography and nuclear magnetic resonance to mediate protein-ligand interactions. Here, we examine the energetics of hydrating cavities formed at protein-ligand interfaces using molecular dynamics simulations. The free energies of hydrating two cavities in the active site of two liganded complexes of cytochrome P450cam were calculated by multiconfigurational thermodynamic integration. The complex of cytochrome P450cam with 2-phenyl-imidazole contains a crystallographically well defined water molecule mediating hydrogen bonds between the protein and the inhibitor. We calculate that this water molecule is stabilized by a binding free energy of -11.6 +/- kJ/mol. The complex of cytochrome P450cam with its natural substrate, camphor, contains a cavity that is empty in the crystal structure although a water molecule in it could make a hydrogen bond to camphor. Here, solvation of this cavity is calculated to be unfavorable by +15.8 +/- 5.0 kJ/mol. The molecular dynamics simulations can thus distinguish a hydrated interfacial cavity from an empty one. They also provide support for the notion that protein-ligand complexes can accommodate empty interfacial cavities and that such cavities are likely to be unhydrated unless more than one hydrogen bond can be made to a water molecule in the cavity.     

Free energy calculations of protein-ligand interactions

In the calculation of free energies of binding for protein-ligand complexes, we distinguish endpoint methods, methods involving alchemical modifications and methods that physically displace the ligand from the protein. Most methodological advances seem to come from a clever combination of multiple existing methods to enhance the sampling or to utilize specific advantages of various approaches. The coupling parameters common in thermodynamic integration and in Hamiltonian replica exchange are for instance combined to yield replica exchange thermodynamic integration. As new methods mostly aim to improve efficiency or to attain more complete sampling, there are good prospects to understand and tackle the sampling problem better and to shift the focus towards the scoring problem in the context of more robust and accurate force fields.     

Understanding protein-ligand interactions: the price of protein flexibility

In order to design selective, high-affinity ligands to a target protein, it is advantageous to understand the structural determinants for protein-ligand complex formation at the atomic level. In a model system, we have successively mapped the factor Xa binding site onto trypsin, showing that certain mutations influence both protein structure and inhibitor specificity. Our previous studies have shown that introduction of the 172SSFI175 sequence of factor Xa into rat or bovine trypsin results in the destabilisation of the intermediate helix with burial of Phe174 (the down conformation). Surface exposure of the latter residue (the up conformation) is critical for the correct formation of the aromatic box found in factor Xa-ligand complexes. In the present study, we investigate the influence of aromatic residues in position 174. Replacement with the bulky tryptophan (SSWI) shows reduced affinity for benzamidine-based inhibitors (1) and (4), whereas removal of the side-chain (alanine, SSAI) or exchange with a hydrophilic residue (arginine, SSRI) leads to a significant loss in affinity for all inhibitors studied. The variants could be crystallised in the presence of different inhibitors in multiple crystal forms. Structural characterisation of the variants revealed three different conformations of the intermediate helix and 175 loop in SSAI (down, up and super-up), as well as a complete disorder of this region in one crystal form of SSRI, suggesting that the compromised affinity of these variants is related to conformational flexibility. The influence of Glu217, peripheral to the ligand-binding site in factor Xa, was investigated. Introduction of Glu217 into trypsin variants containing the SSFI sequence exhibited enhanced affinity for the factor Xa ligands (2) and (3). The crystal structures of these variants also exhibited the down and super-up conformations, the latter of which could be converted to up upon soaking and binding of inhibitor (2). The improved affinity of the Glu217-containing variants appears to be due to a shift towards the up conformation. Thus, the reduction in affinity caused by conformational variability of the protein target can be partially or wholly offset by compensatory binding to the up conformation. The insights provided by these studies will be helpful in improving our understanding of ligand binding for the drug design process.     

Analysis of protein-ligand interactions by fluorescence polarization

Quantification of the associations between biomolecules is required both to predict and understand the interactions that underpin all biological activity. Fluorescence polarization (FP) provides a nondisruptive means of measuring the association of a fluorescent ligand with a larger molecule. We describe an FP assay in which binding of fluorescein-labeled inositol 1,4,5-trisphosphate (IP(3)) to N-terminal fragments of IP(3) receptors can be characterized at different temperatures and in competition with other ligands. The assay allows the standard Gibbs free energy (ΔG°), enthalpy (ΔH°) and entropy (ΔS°) changes of ligand binding to be determined. The method is applicable to any purified ligand-binding site for which an appropriate fluorescent ligand is available. FP can be used to measure low-affinity interactions in real time without the use of radioactive materials, it is nondestructive and, with appropriate care, it can resolve ΔH° and ΔS°. The first part of the protocol, protein preparation, may take several weeks, whereas the FP measurements, once they have been optimized, would normally take 1-6 h.     

Nicotiana benthamiana: its history and future as a model for plant-pathogen interactions

Nicotiana benthamiana is the most widely used experimental host in plant virology, due mainly to the large number of diverse plant viruses that can successfully infect it. Additionally, N. benthamiana is susceptible to a wide variety of other plant-pathogenic agents (such as bacteria, oomycetes, fungi, and so on), making this species a cornerstone of host-pathogen research, particularly in the context of innate immunity and defense signaling. Moreover, because it can be genetically transformed and regenerated with good efficiency and is amenable to facile methods for virus-induced gene silencing or transient protein expression, N. benthamiana is rapidly gaining popularity in plant biology, particularly in studies requiring protein localization, interaction, or plant-based systems for protein expression and purification. Paradoxically, despite being an indispensable research model, little is known about the origins, genetic variation, or ecology of the N. benthamiana accessions currently used by the research community. In addition to addressing these latter topics, the purpose of this review is to provide information regarding sources for tools and reagents that can be used to support research in N. benthamiana. Finally, we propose that N. benthamiana is well situated to become a premier plant cell biology model, particularly for the virology community, who as a group were the first to recognize the potential of this unique Australian native.     

Screening of protein-ligand interactions by affinity chromatography

This paper examines affinity chromatography (AC) as an alternative tool for the determination of protein-ligand interactions for the particular case in which the ligand is the same protein. The methodology is less labor-intensive and more sample-efficient than traditional methods used to measure the second virial coefficient (B(22)), a parameter commonly used to evaluate protein-protein interactions. The chromatographic capacity factor (k') was studied for lysozyme and equine serum albumin for a wide range of experimental solution conditions such as crystallizing agent concentration, protein concentration and pH. Parallel experiments using AC to determine k' and static light scattering (SLS) to determine B(22) showed that the two parameters were highly correlated. Two different column volumes ( approximately 1 and approximately 0.1 mL) were tested and gave essentially the same values for k', showing the feasibility of miniaturization.     

Extent of enthalpy-entropy compensation in protein-ligand interactions

The extent of enthalpy-entropy compensation in protein-ligand interactions has long been disputed because negatively correlated enthalpy (ΔH) and entropy (TΔS) changes can arise from constraints imposed by experimental and analytical procedures as well as through a physical compensation mechanism. To distinguish these possibilities, we have created quantitative models of the effects of experimental constraints on isothermal titration calorimetry (ITC) measurements. These constraints are found to obscure any compensation that may be present in common data representations and regression analyses (e.g., in ΔH vs. -TΔS plots). However, transforming the thermodynamic data into ΔΔ-plots of the differences between all pairs of ligands that bind each protein diminishes the influence of experimental constraints and representational bias. Statistical analysis of data from 32 diverse proteins shows a significant and widespread tendency to compensation. ΔΔH versus ΔΔG plots reveal a wide variation in the extent of compensation for different ligand modifications. While strong compensation (ΔΔH and -TΔΔS opposed and differing by < 20% in magnitude) is observed for 22% of modifications (twice that expected without compensation), 15% of modifications result in reinforcement (ΔΔH and -TΔΔS of the same sign). Because both enthalpy and entropy changes arise from changes to the distribution of energy states on binding, there is a general theoretical expectation of compensated behavior. However, prior theoretical studies have focussed on explaining a stronger tendency to compensation than actually found here. These results, showing strong but imperfect compensation, will act as a benchmark for future theoretical models of the thermodynamic consequences of ligand modification.     

General and targeted statistical potentials for protein-ligand interactions

We present a novel atom-atom potential derived from a database of protein-ligand complexes. First, we clarify the similarities and differences between two statistical potentials described in the literature, PMF and Drugscore. We highlight shortcomings caused by an important factor unaccounted for in their reference states, and describe a new potential, which we name the Astex Statistical Potential (ASP). ASP's reference state considers the difference in exposure of protein atom types towards ligand binding sites. We show that this new potential predicts binding affinities with an accuracy similar to that of Goldscore and Chemscore. We investigate the influence of the choice of reference state by constructing two additional statistical potentials that differ from ASP only in this respect. The reference states in these two potentials are defined along the lines of Drugscore and PMF. In docking experiments, the potential using the new reference state proposed for ASP gives better success rates than when these literature reference states were used; a success rate similar to the established scoring functions Goldscore and Chemscore is achieved with ASP. This is the case both for a large, general validation set of protein-ligand structures and for small test sets of actives against four pharmaceutically relevant targets. Virtual screening experiments for these targets show less discrimination between the different reference states in terms of enrichment. In addition, we describe how statistical potentials can be used in the construction of targeted scoring functions. Examples are given for cdk2, using four different targeted scoring functions, biased towards increasingly large target-specific databases. Using these targeted scoring functions, docking success rates as well as enrichments are significantly better than for the general ASP scoring function. Results improve with the number of structures used in the construction of the target scoring functions, thus illustrating that these targeted ASP potentials can be continuously improved as new structural data become available.     

Mesenchymal-epithelial interactions: past, present, and future

Abstract This review summarizes the history of research on mesenchymal–epithelial interactions in prostatic development from the first studies in 1970 to the present. From this study we have learned that prostatic development requires a reciprocal interaction between epithelium and mesenchyme in which urogenital sinus mesenchyme induces and patterns epithelial development and differentiation, while developing prostatic epithelium induces and patterns mesenchymal differentiation into smooth muscle and other resident cell types in the stroma. Prostatic development requires androgen action mediated by the androgen receptor (AR). Through analysis of tissue recombinants composed of wild-type and AR-null epithelium and mesenchyme, we have learned that many "androgenic effects" on prostatic epithelium do not require epithelial AR, but instead are elicited by the paracrine action of AR-positive mesenchyme. Present and future studies reviewed in this issue deal with the molecular mechanisms in this developmental communication between epithelium and mesenchyme.     

Protein-protein and protein-ligand interactions studied by electrospray-ionization mass spectrometry

Preservation of non-covalent interactions in biopolymer mass spectrometry offers new approaches to binding analysis. Recent work from our laboratory is reviewed here and discussed with reference to recent literature in the field. Three issues are considered in particular: hydrophobically stabilized complexes, pH-dependent transitions, and linked protein-ligand and protein-protein binding equilibria.     

Thermodynamics of drug-DNA interactions

Many anticancer, antibiotic, and antiviral drugs exert their primary biological effects by reversibly interacting with nucleic acids. Therefore, these biomolecules represent a major target in drug development strategies designed to produce next generation therapeutics for diseases such as cancer. In order to improve the clinical efficacy of existing drugs and also to design new ones it is necessary to understand the molecular basis of drug-DNA interactions in structural, thermodynamic, and kinetic detail. The past decade has witnessed an increase in the number of rigorous biophysical studies of drug-DNA systems and considerable knowledge has been gained in the energetics of these binding reactions. This is, in part, due to the increased availability of high-sensitivity calorimetric techniques, which have allowed the thermodynamics of drug-DNA interactions to be probed directly and accurately. The focus of this article is to review thermodynamic approaches to examining drug-DNA recognition. Specifically, an overview of a recently developed method of analysis that dissects the binding free energy of these reactions into five component terms is presented. The results of applying this analysis to the DNA binding interactions of both minor groove drugs and intercalators are discussed. The solvent water plays a key role in nucleic acid structure and consequently in the binding of ligands to these biomolecules. Any rational approach to DNA-targeted drug design requires an understanding of how water participates in recognition and binding events. Recent studies examining hydration changes that accompany DNA binding by intercalators will be reviewed. Finally some aspects of cooperativity in drug-DNA interactions are described and the importance of considering cooperative effects when examining these reactions is highlighted.     

Thermodynamics of lipid-peptide interactions

This review is focused on peptide molecules which exhibit a limited solubility in the aqueous phase and bind to the lipid membrane from the aqueous medium. Surface adsorption, membrane insertion, and specific binding are usually accompanied by changes in the heat content of the system and can be measured conveniently with isothermal titration calorimetry, avoiding the necessity of peptide labeling. The driving forces for peptide adsorption and binding are hydrophobicity, electrostatics, and hydrogen bonding. An exclusively hydrophobic interaction is exemplified by the immunosuppressant drug cyclosporine A. Its insertion into the membrane can be described by a simple partition equilibrium X(b)=K(0)C(eq). If peptide and membrane are both charged, electrostatic interactions are dominant leading to nonlinear binding curves. The concentration of the peptide near the membrane interface can then be much larger than its bulk concentration. Electrostatic effects must be accounted for by means of the Gouy-Chapman theory before conventional binding models can be applied. A small number of peptides and proteins bind with very high affinity to a specific lipid species only. This is illustrated for the lantibiotic cinnamycin (Ro 09-0198) which forms a 1:1 complex with phosphatidyethanolamine with a binding constant of 10(8) M(-1). Membrane adsorption and insertion can be accompanied by conformational transitions facilitated, in part, by hydrogen bonding mechanisms. The two membrane-induced conformational changes to be discussed are the random coil-to-alpha-helix transition of amphipathic peptides and the random coil-to-beta-structure transition of Alzheimer peptides.     

ConsDock: A new program for the consensus analysis of protein-ligand interactions

Protein-based virtual screening of chemical libraries is a powerful technique for identifying new molecules that may interact with a macromolecular target of interest. Because of docking and scoring limitations, it is more difficult to apply as a lead optimization method because it requires that the docking/scoring tool is able to propose as few solutions as possible and all of them with a very good accuracy for both the protein-bound orientation and the conformation of the ligand. In the present study, we present a consensus docking approach (ConsDock) that takes advantage of three widely used docking tools (Dock, FlexX, and Gold). The consensus analysis of all possible poses generated by several docking tools is performed sequentially in four steps: (i) hierarchical clustering of all poses generated by a docking tool into families represented by a leader; (ii) definition of all consensus pairs from leaders generated by different docking programs; (iii) clustering of consensus pairs into classes, represented by a mean structure; and (iv) ranking the different means starting from the most populated class of consensus pairs. When applied to a test set of 100 protein-ligand complexes from the Protein Data Bank, ConsDock significantly outperforms single docking with respect to the docking accuracy of the top-ranked pose. In 60% of the cases investigated here, ConsDock was able to rank as top solution a pose within 2 A RMSD of the X-ray structure. It can be applied as a postprocessing filter to either single- or multiple-docking programs to prioritize three-dimensional guided lead optimization from the most likely docking solution.