Energetics of complex formation of the DNA hairpin structure d(GCGAAGC) with aromatic ligands

The energy contributions of various physical interactions to the total Gibbs energy of complex formation of the biologically important DNA hairpin d(GCGAAGC) with aromatic antitumor antibiotics daunomycin and novantron and the mutagens ethidium and proflavine have been calculated. It has been shown that the relatively small value of the total energy of binding of the ligands to the hairpin is the sum of components great in absolute value and different in sign. The contributions of van der Waals interactions and both intra- and intermolecular hydrogen bonds and bonds with aqueous environment have been studied. According to the calculations, the hydrophobic and van der Waals components are energetically favorable in complex formation of the ligands with the DNA pairpin d(GCGAAGC), whereas the electrostatic (with consideration of hydrogen bonds) and entropic components are unfavorable.     

Contribution of enthalpy to the energetics of complex formation of aromatic ligands with DNA

The energy contributions of electrostatic, van der Waals interactions, hydrogen bonds, and interactions of charge transfer type to the enthalpy of complex formation of the double-stand DNA with the antitumor antibiotics daunomycin, nogalamycin, and novantron, as well as the mutagens ethidium bromide and proflavine have been calculated. According to the calculations, the van der Waals component (except for nogalamycin) is energetically favorable during complex formation of the antibiotics with DNA, and the contributions of H bonds and electrostatic interactions are unfavorable, with the probability of charge transfer in the complexes being low. It has been shown that the relatively low value of the experimental enthalpy of binding is the sum of components greater in absolute value and different in the sign, which is the cause of large errors in estimating the total enthalpy of complex formation of aromatic ligands with DNA.     

Energy of ligand-RNA complex formation

This work presents an energetic analysis of complex formation of eleven ligands with different structure and charge with RNA aptamers. The most entire set of the components of the total Gibbs energy of complex formation has been first calculated using different physical factors: van der Waals, electrostatic and hydrophobic interactions, hydrogen bonds and specific factors being predominantly of entropic character. The calculated Gibbs energy is found to be in a good agreement with experimental data. Different energy components which stabilize and destabilize complexes are lined up according to the degree of importance. The results obtained provide an understanding of the role of different physical interactions in a ligand-RNA complex formation.     

Molecular Modeling-Based Energy Analysis of Dimeric Binding of Ligands to the Minor DNA Groove

Molecular-modeling methods have been used to perform energy analysis of dimeric complex formation between lexitropsins and double-stranded DNA. Stabilization of dimeric complexes by hydrophobic and van der Waals interactions has been demonstrated. Electrostatic interactions and the contributions of hydrogen bonds and changes in the number of degrees of freedom had a destabilizing influence. The energy of monomeric and dimeric binding has been compared.     

Solvation thermodynamics of biopolymers. I. Separation of the volume and surface interactions with estimates for proteins

The present paper is a systematic first approach to the problem of solvation thermodynamics of biomolecules. Most previous approaches have been only crude estimates of solvent contributions, and have simply assessed solvation free energy as proportional to surface areas. Here we estimate the various contributions and divide them into (a) hard-core interactions dependent upon the entire volume of solute and (b) the remainder of interactions manifested through surfaces, such as van der Waals, charge-charge, or hydrogen bonds. We have estimated the work to create a cavity with scaled-particle theory (SPT), the van der Waals interactions on the surface, and hydrogen bonds between the surface and the solvent. The conclusion here is that this latter term is the largest component of the solvation free energy of proteins. From estimates on nine diverse proteins, it is clear that the larger the protein, the more dominant is the hydrogen-bond term. In the next paper, we indicate that correlations between hydrogen-bonding groups on the surfaces could increase the magnitude of the hydrogen-bond contribution.     

Probing the electronic structures and properties of neutral and anionic ScSi n (0,−1) (n = 1–6) clusters using ccCA-TM and G4 theory

We apply molecular docking, molecular dynamics (MD) simulation, and binding free energy calculation to investigate and reveal the binding mechanism between five xanthine inhibitors and DPP-4. The electrostatic and van der Waals interactions of the five inhibitors with DPP-4 are analyzed and discussed. The computed binding free energies using MM-PBSA method are in qualitatively agreement with experimental inhibitory potency of five inhibitors. The hydrogen bonds of inhibitors with Ser630 and Asp663 can stabilize the inhibitors in binding sites. The van der Waals interactions, especially the key contacts with His740, Asn710, Trp629, and Tyr666 have larger contributions to the binding free energy and play important roles in distinguishing the variant bioactivity of five inhibitors.     

Dielectric screening effect of electronic polarization and intramolecular hydrogen bonding

Recent site‐resolved hydrogen exchange measurements have uncovered significant discrepancies between simulations and experimental data during protein folding, including the excessive intramolecular hydrogen bonds in simulations. This finding indicates a possibility that intramolecular charge–charge interactions have not included sufficient dielectric screening effect of the electronic polarization. Scaling down peptide atomic charges according to the optical dielectric constant is tested in this study. As a result, the number of intramolecular hydrogen bonds is lower than using unscaled atomic charges while reaching the same levels of helical contents or β‐hairpin backbone hydrogen bonds, because van der Waals interactions contribute substantially to peptide folding in water. Reducing intramolecular charge–charge interactions and hydrogen bonding increases conformational search efficiency. In particular, it reduces the equilibrium helical content in simulations using AMBER force field and the energy barrier in folding simulations using CHARMM force field.     

Study of beta-turns in globular proteins

The formation of beta-turns in globular proteins has been studied by the method of molecular mechanics. Statistical method of discriminant analysis was applied to calculate energy components and sequences of oligopeptide segments, and after this prediction of I type beta-turns has been drawn. The accuracy of true positive prediction is 65%. Components of conformational energy considerably affecting beta-turn formation were delineated. There are torsional energy, energy of hydrogen bonds, and van der Waals energy.     

Understanding beta-hairpin formation by molecular dynamics simulations of unfolding

We have studied the mechanism of formation of a 16-residue beta-hairpin from the protein GB1 using molecular dynamics simulations in an aqueous environment. The analysis of unfolding trajectories at high temperatures suggests a refolding pathway consisting of several transient intermediates. The changes in the interaction energies of residues are related with the structural changes during the unfolding of the hairpin. The electrostatic energies of the residues in the turn region are found to be responsible for the transition between the folded state and the hydrophobic core state. The van der Waals interaction energies of the residues in the hydrophobic core reflect the behavior of the radius of gyration of the core region. We have examined the opposing influences of the protein-protein (PP) energy, which favors the native state, and the protein-solvent (PS) energy, which favors unfolding, in the formation of the beta-hairpin structure. It is found that the behavior of the electrostatic components of PP and PS energies reflects the structural changes associated with the loss of backbone hydrogen bonding. Relative changes in the PP and PS van der Waals interactions are related with the disruption of the hydrophobic core of a protein. The results of the simulations support the hydrophobic collapse mechanism of beta-hairpin folding.     

Amino acid–base interactions: a three-dimensional analysis of protein–DNA interactions at an atomic level

To assess whether there are universal rules that govern amino acid–base recognition, we investigate hydrogen bonds, van der Waals contacts and water-mediated bonds in 129 protein–DNA complex structures. DNA–backbone interactions are the most numerous, providing stability rather than specificity. For base interactions, there are significant base–amino acid type correlations, which can be rationalised by considering the stereochemistry of protein side chains and the base edges exposed in the DNA structure. Nearly two-thirds of the direct read-out of DNA sequences involves complex networks of hydrogen bonds, which enhance specificity. Two-thirds of all protein–DNA interactions comprise van der Waals contacts, compared to about one-sixth each of hydrogen and water-mediated bonds. This highlights the central importance of these contacts for complex formation, which have previously been relegated to a secondary role. Although common, water-mediated bonds are usually non-specific, acting as space-fillers at the protein–DNA interface. In conclusion, the majority of amino acid–base interactions observed follow general principles that apply across all protein–DNA complexes, although there are individual exceptions. Therefore, we distinguish between interactions whose specificities are ‘universal’ and ‘context-dependent’. An interactive Web-based atlas of side chain–base contacts provides access to the collected data, including analyses and visualisation of the three-dimensional geometry of the interactions.     

Interaction of colchicine with human serum albumin investigated by spectroscopic methods

We investigated the interaction between colchicine and human serum albumin (HSA) by fluorescence and UV-vis absorption spectroscopy. In the mechanism discussion, it was proved that the fluorescence quenching of HSA by colchicine is a result of the formation of colchicines-HSA complex; van der Waals interactions and hydrogen bonds play a major role in stabilizing the complex. The modified Stern-Volmer quenching constant K(a) and corresponding thermodynamic parameters deltaH, deltaG, deltaS at different temperatures were calculated. The distance r between donor (Trp214) and acceptor (colchicine) was obtained according to fluorescence resonance energy transfer (FRET).     

Van der Waals forces. Special characteristics in lipid-water systems and a general method of calculation based on the Lifshitz theory

A practical method for examining and calculating van der Waals forces is derived from Lifshitz'' theory. Rather than treat the total van der Waals energy as a sum of pairwise interactions between atoms, the Lifshitz theory treats component materials as continua in which there are electromagnetic fluctuations at all frequencies over the entire body. It is necessary in principle to use total macroscopic dielectric data from component substances to analyze the permitted fluctuations; in practice it is possible to use only partial information to perform satisfactory calculations. The biologically interesting case of lipid-water systems is considered in detail for illustration. The method gives good agreement with measured van der Waals energy of interaction across a lipid film. It appears that fluctuations at infrared frequencies and microwave frequencies are very important although these are usually ignored in preference to UV contributions. “Retardation effects” are such as to damp out high frequency fluctuation contributions; if interaction specificity is due to UV spectra, this will be revealed only at interactions across <200 angstrom (A). Dependence of van der Waals forces on material electric properties is discussed in terms of illustrative numerical calculations.     

Thermodynamics of barnase unfolding.

The thermodynamics of barnase denaturation has been studied calorimetrically over a broad range of temperature and pH. It is shown that in acidic solutions the heat denaturation of barnase is well approximated by a 2-state transition. The heat denaturation of barnase proceeds with a significant increase of heat capacity, which determines the temperature dependencies of the enthalpy and entropy of its denaturation. The partial specific heat capacity of denatured barnase is very close to that expected for the completely unfolded protein. The specific denaturation enthalpy value extrapolated to 130 degrees C is also close to the value expected for the full unfolding. Therefore, the calorimetrically determined thermodynamic characteristics of barnase denaturation can be considered as characteristics of its complete unfolding and can be correlated with structural features--the number of hydrogen bonds, extent of van der Waals contacts, and the surface areas of polar and nonpolar groups. Using this information and thermodynamic information on transfer of protein groups into water, the contribution of various factors to the stabilization of the native structure of barnase has been estimated. The main contributors to the stabilization of the native state of barnase appear to be intramolecular hydrogen bonds. The contributions of van der Waals interactions between nonpolar groups and those of hydration effects of these groups are not as large if considered separately, but the combination of these 2 factors, known as hydrophobic interactions, is of the same order of magnitude as the contribution of hydrogen bonding.     

Protein docking using continuum electrostatics and geometric fit

The computer program DOT quickly finds low-energy docked structures for two proteins by performing a systematic search over six degrees of freedom. A novel feature of DOT is its energy function, which is the sum of both a Poisson-Boltzmann electrostatic energy and a van der Waals energy, each represented as a grid-based correlation function. DOT evaluates the energy of interaction for many orientations of the moving molecule and maintains separate lists scored by either the electrostatic energy, the van der Waals energy or the composite sum of both. The free energy is obtained by summing the Boltzmann factor over all rotations at each grid point. Three important findings are presented. First, for a wide variety of protein-protein interactions, the composite-energy function is shown to produce larger clusters of correct answers than found by scoring with either van der Waals energy (geometric fit) or electrostatic energy alone. Second, free-energy clusters are demonstrated to be indicators of binding sites. Third, the contributions of electrostatic and attractive van der Waals energies to the total energy term appropriately reflect the nature of the various types of protein-protein interactions studied.     

Investigation of Atomic Level Patterns in Protein—Small Ligand Interactions

Background

Shape complementarity and non-covalent interactions are believed to drive protein-ligand interaction. To date protein-protein, protein-DNA, and protein-RNA interactions were systematically investigated, which is in contrast to interactions with small ligands. We investigate the role of covalent and non-covalent bonds in protein-small ligand interactions using a comprehensive dataset of 2,320 complexes.

Methodology and Principal Findings

We show that protein-ligand interactions are governed by different forces for different ligand types, i.e., protein-organic compound interactions are governed by hydrogen bonds, van der Waals contacts, and covalent bonds; protein-metal ion interactions are dominated by electrostatic force and coordination bonds; protein-anion interactions are established with electrostatic force, hydrogen bonds, and van der Waals contacts; and protein-inorganic cluster interactions are driven by coordination bonds. We extracted several frequently occurring atomic-level patterns concerning these interactions. For instance, 73% of investigated covalent bonds were summarized with just three patterns in which bonds are formed between thiol of Cys and carbon or sulfur atoms of ligands, and nitrogen of Lys and carbon of ligands. Similar patterns were found for the coordination bonds. Hydrogen bonds occur in 67% of protein-organic compound complexes and 66% of them are formed between NH- group of protein residues and oxygen atom of ligands. We quantify relative abundance of specific interaction types and discuss their characteristic features. The extracted protein-organic compound patterns are shown to complement and improve a geometric approach for prediction of binding sites.

Conclusions and Significance

We show that for a given type (group) of ligands and type of the interaction force, majority of protein-ligand interactions are repetitive and could be summarized with several simple atomic-level patterns. We summarize and analyze 10 frequently occurring interaction patterns that cover 56% of all considered complexes and we show a practical application for the patterns that concerns interactions with organic compounds.     

Stereochemical complementarity of progesterone, RU486 and cavities between base pairs in partially unwound double stranded DNA assessed by computer modeling and energy calculations

Computer modeling was used to examine the relative fit of progesterone and RU486 in cavities constructed between base pairs in double stranded DNA. Progesterone was capable of forming two stereospecific hydrogen bonds between the carbonyl groups and protonated phosphate groups on adjacent strands. Favorable van der Waals and electrostatic energies were exhibited upon insertion of progesterone into DNA indicating an excellent fit. While RU486 could be accommodated between the base pairs and formed hydrogen bonds, there was a high van der Waals energy in the resulting complex. When the complexes were subjected to energy minimization, the conformation of the DNA was significantly altered in the RU486-DNA complex but not in the progesterone-DNA complex. No mechanistic interpretation of these results is proffered; however, such information may have evolutionary significance and could prove useful in designing new progesterone agonists and antagonists.     

Models of binding of 4'-nitrophenyl alpha-D-mannopyranoside to the lectin concanavalin A

Molecular models for the complex formed between the lectin concanavalin A (Con A) and the saccharide derivative 4'-nitrophenyl-alpha-D-mannopyranoside (alpha-PNM) are presented, combining evidence from 1H-n.m.r. measurements, semi-empirical energy calculations and interactive graphics modelling. The models are in good agreement with the experimental data. Close examination of the models suggests that hydrophobic interactions together with van der Waals interactions and hydrogen bonds contribute to the stability of the complexes. It appears that there is a limited number of possible modes of binding of alpha-PNM to Con A.     

The 2.3-A resolution structure of the maltose- or maltodextrin-binding protein, a primary receptor of bacterial active transport and chemotaxis

The three-dimensional structure of the maltose- or maltodextrin-binding protein (Mr = 40,622) with bound maltose has been obtained by crystallographic analysis at 2.8-A resolution. The structure, which has been partially refined at 2.3 A, is ellipsoidal with overall dimensions of 30 x 40 x 65 A and divided into two distinct globular domains by a deep groove. Although each domain is built from two peptide segments from the amino- and carboxyl-terminal halves, both domains exhibit similar supersecondary structure, consisting of a central beta-pleated sheet flanked on both sides with two or three parallel alpha-helices. The groove, which has a depth of 18 A and a base of about 9 x 18 A, contains the maltodextrin-binding site. We have previously observed the same general features in the well-refined structures of six other periplasmic receptors with specificities for L-arabinose, D-galactose/D-glucose, sulfate, phosphate, leucine/isoleucine/valine, and leucine. The bound maltose is buried in the groove and almost completely inaccessible to the bulk solvent. The groove is heavily populated by polar and aromatic groups many of which are involved in extensive hydrogen-bonding and van der Waals interactions with the maltose. All the disaccharide hydroxyl groups, which form a peripheral polar surface approximately in the plane of the sugar rings, are tied in a total of 11 direct hydrogen bonds with six charged side chains, one Trp side chain, and one peptide backbone NH, and five indirect hydrogen bonds via water molecules. The maltose is wedged between four aromatic side chains. The resulting stacking of these aromatic residues on the faces of the glucosyl units provides a majority of the van der Waals contacts in the complex. The nonreducing glucosyl unit of the maltose is involved in approximately twice as many hydrogen bonds and van der Waals contacts as the glucosyl unit at the reducing end. The binding protein-maltose complex shows the best example of the extensive use of polar and aromatic residues in binding oligosaccharides. The tertiary structure of the maltodextrin-binding protein, along with the results of genetic studies by a number of investigators, has also enabled us for the first time to map the different regions on the surface of the protein involved in the interactions with the membrane-bound protein components necessary for transport of and chemotaxis toward maltodextrins. These sites permit distinction of the "open cleft" (without bound sugar) and closed (with bound sugar) conformations of the binding protein by the chemotactic signal transducer with which the maltodextrin-binding protein interacts.(ABSTRACT TRUNCATED AT 250 WORDS)     

Novel homodimer model of the β-adrenergic receptor in complex with free fatty acids and cholesterol: first-principles calculation studies

We propose a theoretical novel homodimer model of the β- adrenergic receptor (βAR) in complex with a heterogeneous mixture of free fatty acids (FFAs) and cholesterol based on first-principles calculations. We used the density-functional-based tight binding with dispersion (DFTB-D) method, which accurately evaluates van der Waals interactions between FFAs and amino acid residues in the receptor. The calculations suggest that a stable homodimer of bAR can form a complex with FFAs and cholesterol by extensive van der Waals interactions in the cell membrane, and that the heterogeneous composition of the FFAs is important for the stability of the homodimer complex. The stable van der Waals interactions propagate from one of the bAR to the other through the cholesterol and FFAs in the homodimer complex. The energy propagation in the complex has the potential to enhance molecular signaling in adipocytes, because the stability of the complex can influence anti-adiposity effects after oral treatment of the FFA components.