Carbohydrate chain of ganglioside GM1 as a ligand: identification of the binding strategies of three 15 mer peptides and their divergence from the binding modes of growth-regulatory galectin-1 and cholera toxin

The branched pentasaccharide chain of ganglioside GM1 is a prominent cell surface ligand, for example, for cholera toxin or tumor growth-regulatory homodimeric galectins. This activity profile via protein recognition prompted us to examine the binding properties of peptides with this specificity. Our study provides insights into the mechanism of molecular interaction of this thus far unexplored size limit of the protein part. We used three pentadecapeptides in a combined approach of mass spectrometry, NMR spectroscopy and molecular modelling to analyze the ligand binding in solution. Availability of charged and hydrophobic functionalities affected the intramolecular flexibility of the peptides differently. Backfolding led to restrictions in two cases; the flexibility was not reduced significantly by association of the ligand in its energetically privileged conformations. Major contributions to the interaction energy arise from the sialic acid moiety contacting Arg/Lys residues and the N-terminal charge. Considerable involvement of stacking between the monovalent ligand and aromatic rings could not be detected. This carbohydrate binding strategy is similar to how an adenoviral fiber knob targets sialylated glycans. Rational manipulation for an affinity enhancement can now be directed to reduce the flexibility, exploit the potential for stacking and acquire the cross-linking capacity of the natural lectins by peptide attachment to a suitable scaffold.     

Identification and Characterization of a Single High‐Affinity Fatty Acid Binding Site in Human Serum Albumin

A single high‐affinity fatty acid binding site in the important human transport protein serum albumin (HSA) is identified and characterized using an NBD (7‐nitrobenz‐2‐oxa‐1,3‐diazol‐4‐yl)‐C₁₂ fatty acid. This ligand exhibits a 1:1 binding stoichiometry in its HSA complex with high site‐specificity. The complex dissociation constant is determined by titration experiments as well as radioactive equilibrium dialysis. Competition experiments with the known HSA‐binding drugs warfarin and ibuprofen confirm the new binding site to be different from Sudlow‐sites I and II. These binding studies are extended to other albumin binders and fatty acid derivatives. Furthermore an X‐ray crystal structure allows locating the binding site in HSA subdomain IIA. The knowledge about this novel HSA site will be important for drug depot development and for understanding drug‐protein interaction, which are important prerequisites for modulation of drug pharmacokinetics.     

Influence of hydrogen bonds between sidechains and cooperativity on self‐assembly of the cyclopeptides

The self‐assembly of four cyclic D,L‐octapeptides, [‐(D‐Ala‐Gln)₄‐], [‐(D‐Val‐Gln)₄‐], [‐(D‐Leu‐Gln)₄‐], and [‐(D‐Phe‐Gln)₄‐], was investigated on the theory level in detail. Based on these cyclic peptides, which contain L‐Gln residues and possess C₄ symmetry, a series of oligomers were constructed according to different stacking modes as well as interaction patterns. We employed the semiempirical molecular orbital method AM1 to optimize the structures of all the oligomers, some of which were further studied using density functional method B3PW91/6‐31G to calculate the interaction energies. The studies indicate that when these cyclopeptides aggregate to form oligomers, or even nanotubes, four more hydrogen bonds could form between the sidechains of L‐Gln residues in addition to eight hydrogen bonds formed between the backbones of adjacent two cyclic peptides, a result that would clearly affect the self‐assembling process of cyclic peptides. The main effects can be summarized as follows. First, the dimers of these cyclic peptides with C₄ symmetry are more stable than those with D₄ symmetry due to their additional H‐bonds between Gln sidechains. Second, for the self‐assembly of the cyclopeptides, there is a competition between parallel and antiparallel stacking modes in lower oligomers such as dimers. However, with an increasing degree of oligomerization, energetically there is an increased possibility for the cyclic peptides to take the parallel stacking mode in assembly. Finally, the synergetic effect of weak interactions is the fundamental driving force for cyclic peptides to form stable nanotubes. © 2005 Wiley Periodicals, Inc. Int J Quantum Chem, 2005     

Including ligand‐induced protein flexibility into protein tunnel prediction

In proteins with buried active sites, understanding how ligands migrate through the tunnels that connect the exterior of the protein to the active site can shed light on substrate specificity and enzyme function. A growing body of evidence highlights the importance of protein flexibility in the binding site on ligand binding; however, the influence of protein flexibility throughout the body of the protein during ligand entry and egress is much less characterized. We have developed a novel tunnel prediction and evaluation method named IterTunnel, which includes the influence of ligand‐induced protein flexibility, guarantees ligand egress, and provides detailed free energy information as the ligand proceeds along the egress route. IterTunnel combines geometric tunnel prediction with steered molecular dynamics in an iterative process to identify tunnels that open as a result of ligand migration and calculates the potential of mean force of ligand egress through a given tunnel. Applying this new method to cytochrome P450 2B6, we demonstrate the influence of protein flexibility on the shape and accessibility of tunnels. More importantly, we demonstrate that the ligand itself, while traversing through a tunnel, can reshape tunnels due to its interaction with the protein. This process results in the exposure of new tunnels and the closure of preexisting tunnels as the ligand migrates from the active site. © 2014 Wiley Periodicals, Inc.     

Site‐Specific Investigation of the Steady‐State Kinetics and Dynamics of the Multistep Binding of Bile Acid Molecules to a Lipid Carrier Protein

The investigation of multi‐site ligand–protein binding and multi‐step mechanisms is highly demanding. In this work, advanced NMR methodologies such as 2D ¹H–¹⁵N line‐shape analysis, which allows a reliable investigation of ligand binding occurring on micro‐ to millisecond timescales, have been extended to model a two‐step binding mechanism. The molecular recognition and complex uptake mechanism of two bile salt molecules by lipid carriers is an interesting example that shows that protein dynamics has the potential to modulate the macromolecule–ligand encounter. Kinetic analysis supports a conformational selection model as the initial recognition process in which the dynamics observed in the apo form is essential for ligand uptake, leading to conformations with improved access to the binding cavity. Subsequent multi‐step events could be modelled, for several residues, with a two‐step binding mechanism. The protein in the ligand‐bound state still exhibits a conformational rearrangement that occurs on a very slow timescale, as observed for other proteins of the family. A global mechanism suggesting how bile acids access the macromolecular cavity is thus proposed.     

High Relaxivity Supramolecular Adducts Between Human‐Liver Fatty‐Acid‐Binding Protein and Amphiphilic GdIII Complexes: Structural Basis for the Design of Intracellular Targeting MRI Probes

Gadolinium complexes linked to an apolar fragment are known to be efficiently internalized into various cell types, including hepatocytes. Two lipid‐functionalized gadolinium chelates have been investigated for the targeting of the human liver fatty acid binding protein (hL‐FABP) as a means of increasing the sensitivity and specificity of intracellular‐directed MRI probes. hL‐FABP, the most abundant cytosolic lipid binding protein in hepatocytes, displays the ability to interact with multiple ligands involved in lipid signaling and is believed to be an obligate carrier to escort lipidic drugs across the cell. The interaction modes of a fatty acid and a bile acid based gadolinium complex with hL‐FABP have been characterized by relaxometric and NMR experiments in solution with close‐to‐physiological protein concentrations. We have introduced the analysis of paramagnetic‐induced protein NMR signal intensity changes as a quantitative tool for the determination of binding stoichiometry and of precise metal‐ion‐center positioning in protein–ligand supramolecular adducts. A few additional NMR‐derived restraints were then sufficient to locate the ligand molecules in the protein binding sites by using a rapid data‐driven docking method. Relaxometric and ¹³C NMR competition experiments with oleate and the gadolinium complexes revealed the formation of heterotypic adducts, which indicates that the amphiphilic compounds may co‐exist in the protein cavity with physiological ligands. The differences in adduct formation between fatty acid and bile acid based complexes provide the basis for an improved molecular design of intracellular targeted probes.     

An Organometallic Compound which Exhibits a DNA Topology‐Dependent One‐Stranded Intercalation Mode

Understanding how small molecules interact with DNA is essential since it underlies a multitude of pathological conditions and therapeutic interventions. Many different intercalator compounds have been studied because of their activity as mutagens or drugs, but little is known regarding their interaction with nucleosomes, the protein‐packaged form of DNA in cells. Here, using crystallographic methods and molecular dynamics simulations, we discovered that adducts formed by [(η⁶‐THA)Ru(ethylenediamine)Cl][PF₆] (THA=5,8,9,10‐tetrahydroanthracene; RAED‐THA‐Cl[PF₆]) in the nucleosome comprise a novel one‐stranded intercalation and DNA distortion mode. Conversely, the THA group in fact remains solvent exposed and does not disrupt base stacking in RAED‐THA adducts on B‐form DNA. This newly observed DNA binding mode and topology dependence may actually be prevalent and should be considered when studying covalently binding intercalating compounds.     

Synthetic dsDNA‐Binding Peptides Using Natural Compounds as Model

We have developed a series of short DNA‐binding peptides containing newly synthesized, unnatural as well as natural amino acid building blocks. By a combinatorial‐library approach, oligopeptides were developed with moderate dsDNA‐binding affinities. Two strategies were used to further enhance the binding affinity of the lead peptides: Ac‐Arg‐Ual‐Sar‐Chi‐Chi‐Chi‐Arg‐NH₂ and Ac‐Arg‐Cbg‐Cha‐Chi‐Chi‐Tal‐Arg‐NH₂. Site‐selective amino acid substitutions increased the binding affinities up to 2 × 10^?5 M . Further enhancement of the binding affinities could be achieved by coupling of an acridine intercalating unit, using linker arms of different length and flexibility. With the introduction of a new lysine‐based acridine unit, different types of oligopeptide–acridine conjugates were designed using known dsDNA‐binding ligands as model compounds. The binding capacities of these new oligopeptide–acridine conjugates have been investigated by a fluorescent intercalator (ethidium bromide) displacement (FID) assay. With the synthesis of the dipeptide–acridine conjugates, binding affinities in the low micromolar range were obtained (6.4 × 10^?6 M ), which is similar to the binding strength of the well‐known DNA binder Hoechst 33258.     

The Binding Mode of a Tau Peptide with Tubulin

The microtubule‐associated protein Tau promotes the polymerization of tubulin and modulates the function of microtubules. As a consequence of the dynamic nature of the Tau–tubulin interaction, the structural basis of this complex has remained largely elusive. By using NMR methods optimized for ligand–receptor interactions in combination with site‐directed mutagenesis we demonstrate that the flanking domain downstream of the four microtubule‐binding repeats of Tau binds competitively to a site on the α‐tubulin surface. The binding process is complex, involves partial coupling of different interacting regions, and is modulated by phosphorylation at Y394 and S396. This study strengthens the hypothesis of an intimate relationship between Tau phosphorylation and tubulin binding and highlights the power of the INPHARMA NMR method to characterize the interaction of peptides derived from intrinsically disordered proteins with their molecular partners.     

Polycyclic Azoniahetarenes: Assessing the Binding Parameters of Complexes between Unsubstituted Ligands and G‐Quadruplex DNA

Polycyclic azoniahetarenes were employed to determine the effect of the structure of unsubstituted polyaromatic ligands on their quadruplex‐DNA binding properties. The interactions of three isomeric diazoniadibenzo[b,k]chrysenes ( 4 a – c ), diazoniapentaphene ( 5 ), diazoniaanthra[1,2‐a]anthracene ( 6 ), and tetraazoniapentapheno[6,7‐h]pentaphene ( 3 ) with quadruplex DNA were examined by DNA melting studies (FRET melting) and fluorimetric titrations. In general, penta‐ and hexacyclic azoniahetarenes bind to quadruplex DNA (K_b≈10⁶ M ^?1) even in the absence of additional functional side chains. The binding modes of 4 a – c and 3 were studied in more detail by ligand displacement experiments, isothermal titration calorimetry, and CD and NMR spectroscopy. All experimental data indicate that terminal π stacking of the diazoniachrysenes to the quadruplex is the major binding mode; however, because of different electron distributions of the π systems of each isomer, these ligands align differently in the binding site to achieve ideal binding interactions. It is proposed that tetraazonia ligand 3 binds to the quadruplex by terminal stacking with a small portion of its π system, whereas a significant part of the bulky ligand most likely points outside the quadruplex structure, and is thus partially placed in the grooves. Notably, 3 and the known tetracationic porphyrin TMPyP4 exhibit almost the same binding properties towards quadruplex DNA, with 3 being more selective for quadruplex than for duplex DNA. Overall, studies on azonia‐type hetarenes enable understanding of some parameters that govern the quadruplex‐binding properties of parent ligand systems. Since unsubstituted ligands were employed in this study, complementary and cooperative effects of additional substituents, which may interfere with the ligand properties, were eliminated.     

Synthesis,Coordination, and Spectroscopic Properties of 2,2′‐Bipyridyldimesitylpalladium(II) and 6‐Mesityl‐2,2′‐bipyridyldimesitylpalladium(II)

The synthetic route to the dimesitylpalladium(II) complex [(bpy)PdMes₂] ( 1 ) (Mes = mesityl = 2,4,6‐trimethyl phenyl) does not only give the desired compound but also the 6‐mesityl‐2,2′bipyridyldimesitylpalladium [(6‐Mes‐bpy)PdMes₂] ( 2 ) complex and the free ligand 6,6′‐dimesityl‐2,2′‐bipyridine in reasonable yields. Single crystals of 2 were examined by X‐Ray diffraction. The compound reveals a sterically crowded molecular structure. An intramolecular π‐stacking interaction was found between the mesityl substituent on the bipyridine ligand and the adjacent mesityl ligand. The electrochemical behaviour of 1 and 2 together with a related compound was examined at various temperatures showing two reversible reduction reactions and reversible one‐electron oxidation steps at low temperatures. The latter are assigned to Pd^II/Pd^III couples.     

(Thermo)dynamic Role of Receptor Flexibility,Entropy, and Motional Correlation in Protein–Ligand Binding

The binding of 2‐amino‐5‐methylthiazole to the W191G cavity mutant of cytochrome c peroxidase is an ideal test case to investigate the entropic contribution to the binding free energy due to changes in receptor flexibility. The dynamic and thermodynamic role of receptor flexibility are studied by 50 ns‐long explicit‐solvent molecular dynamics simulations of three separate receptor ensembles: W191G binding a K⁺ ion, W191G–2a5mt complex with a closed 190–195 gating loop, and apo with an open loop. We employ a method recently proposed to estimate accurate absolute single‐molecule configurational entropies and their differences for systems undergoing conformational transitions. We find that receptor flexibility plays a generally underestimated role in protein–ligand binding (thermo)dynamics and that changes of receptor motional correlation determine such large entropy contributions.     

The self‐assembled of cyclic D,L‐α‐peptide systems: Insights into the structure and energetics

Cyclic D,L ‐α‐peptides are able to self‐assemble to nanotubes, although the inherent reason of the stability of this kind of nanotube as well as the intrinsic driving force of self‐assembly of the cyclic D ,L ‐α‐peptides still remain elusive. In this work, using several computational approaches, we investigated the structural and energy characteristics of a series of cyclo[(‐L ‐Phe‐D ‐Ala‐)₄] and cyclo[(‐L ‐Ala‐D ‐Ala‐)₄] oligomers. The results reveal that the thermodynamic stability, cooperativity, and self‐assembly patterns of cyclic D ,L ‐α‐peptide nanotubes are mainly determined by the interactions between cross‐strand side chains instead of those between backbones. For cyclo[(‐L ‐Phe‐D ‐Ala‐)₄] oligomers, the steric interaction between cross‐strand side chains, especially the electrostatic repulsion between the phenyls in Phe residues, brings anticooperative effect into parallel stacking mode, which is responsible for the preference of self‐assembling nanotube in antiparallel vs. parallel stacking orientation. Based on our results, a novel self‐assembling mechanism is put forward—it is the L ‐L antiparallel dimer of cyclo[(‐L ‐Phe‐D ‐Ala‐)₄], instead of the commonly presumed monomer, that acts as the basic building block in self assembly. It explains why these cyclic peptides uniquely self‐assemble to form antiparallel nanotubes. © 2009 Wiley Periodicals, Inc. Int J Quantum Chem, 2010     

Mechanistic Studies on the Stereoselectivity of the Serotonin 5‐HT1A Receptor

G‐protein‐coupled receptors (GPCRs) are involved in a wide range of physiological processes, and they have attracted considerable attention as important targets for developing new medicines. A central and largely unresolved question in drug discovery, which is especially relevant to GPCRs, concerns ligand selectivity: Why do certain molecules act as activators (agonists) whereas others, with nearly identical structures, act as blockers (antagonists) of GPCRs? To address this question, we employed all‐atom, long‐timescale molecular dynamics simulations to investigate how two diastereomers (epimers) of dihydrofuroaporphine bind to the serotonin 5‐HT_1A receptor and exert opposite effects. By using molecular interaction fingerprints, we discovered that the agonist could mobilize nearby amino acid residues to act as molecular switches for the formation of a continuous water channel. In contrast, the antagonist epimer remained firmly stabilized in the binding pocket.     

Complex Formation of the Tetracycline‐Binding Aptamer Investigated by Specific Cross‐Relaxation under DNP

While dynamic nuclear polarization (DNP) under magic‐angle spinning (MAS) is generally a powerful method capable of greatly enhancing the sensitivity of solid‐state NMR spectroscopy, hyperpolarization also gives rise to peculiar spin dynamics. Here, we elucidate how specific cross‐relaxation enhancement by active motions under DNP (SCREAM‐DNP) can be utilized to selectively obtain MAS‐NMR spectra of an RNA aptamer in a tightly bound complex with a methyl‐bearing ligand (tetracycline) due to the effective CH₃‐reorientation at an optimized sample temperature of approximately 160 K. SCREAM‐DNP can spectrally isolate the complex from non‐bound species in an RNA mixture. This selectivity allows for a competition assay between the aptamer and a mutant with compromised binding affinity. Variations in molecular structure and methyl dynamics, as observed by SCREAM‐DNP, between free tetracycline and RNA‐bound tetracycline are discussed.     

NMR Structure of a Triangulenium‐Based Long‐Lived Fluorescence Probe Bound to a G‐Quadruplex

An NMR structural study of the interaction between a small‐molecule optical probe (DAOTA‐M2) and a G‐quadruplex from the promoter region of the c‐myc oncogene revealed that they interact at 1:2 binding stoichiometry. NMR‐restrained structural calculations show that binding of DAOTA‐M2 occurs mainly through π–π stacking between the polyaromatic core of the ligand and guanine residues of the outer G‐quartets. Interestingly, the binding affinities of DAOTA‐M2 differ by a factor of two for the outer G‐quartets of the unimolecular parallel G‐quadruplex under study. Unrestrained MD calculations indicate that DAOTA‐M2 displays significant dynamic behavior when stacked on a G‐quartet plane. These studies provide molecular guidelines for the design of triangulenium derivatives that can be used as optical probes for G‐quadruplexes.     

Twist‐dependent stacking energy of base‐pair steps in B‐DNA geometry: A density functional theory approach

Stacking energy of all the 10 unique DNA base‐pair steps (bp step) are calculated using density functional theory within the ultrasoft pseudopotential plane wave method and local density approximation for the exchange‐correlation functional. We have studied the dependence of stacking energy on twist angle, an aspect found difficult to explain using classical theory. We have found that the twist angle for different bp steps at stacking energy minimum matches extremely well with the values of average twist obtained from B‐DNA crystal structure data. This indicates that the use of a proper quantum chemical method to calculate the π‐π electronic interactions may explain stacking energy without incorporating hydrophobic interaction through solvent or effect of backbone through pseudobond. From the twist angle‐dependent stacking energy profile, we have also generated the probability distributions of twist for all the bp steps and calculated the variance of the distribution. Our calculated variances show similar trend to that of the experimental data for which sufficient numbers of data are available. The TA, AT, and CG doublets show large variances among the 10 possible bp steps, indicating their maximum flexibility. This might be the case of unusual deformation observed at the TATA‐box while binding to TBP protein. © 2008 Wiley Periodicals, Inc. Int J Quantum Chem, 2008     

Organic Ligand Binding by a Hydrophobic Cavity in a Designed Tetrameric Coiled‐Coil Protein

The design and characterization of a hydrophobic cavity in de novo designed proteins provides a wide range of information about the functions of de novo proteins. We designed a de novo tetrameric coiled‐coil protein with a hydrophobic pocketlike cavity. Tetrameric coiled coils with hydrophobic cavities have previously been reported. By replacing one Leu residue at the a position with Ala, hydrophobic cavities that did not flatten out due to loose peptide chains were reliably created. To perform a detailed examination of the ligand‐binding characteristics of the cavities, we originally designed two other coiled‐coil proteins: AM2, with eight Ala substitutions at the adjacent a and d positions at the center of a bundled structure, and AM2W, with one Trp and seven Ala substitutions at the same positions. To increase the association of the helical peptides, each helical peptide was connected with flexible linkers, which resulted in a single peptide chain. These proteins exhibited CD spectra corresponding to superhelical structures, despite weakened hydrophobic packing. AM2W exhibited binding affinity for size‐complementary organic compounds. The dissociation constants, K_d, of AM2W were 220 nM for adamantane, 81 μM for 1‐adamantanol, and 294 μM for 1‐adamantaneacetic acid, as measured by fluorescence titration analyses. Although it was contrary to expectations, AM2 did not exhibit any binding affinity, probably due to structural defects around the designed hydrophobic cavity. Interestingly, AM2W exhibited incremental structure stability through ligand binding. Plugging of structural defects with organic ligands would be expected to facilitate protein folding.     

Electronic Olfactory Sensor Based on A. mellifera Odorant‐Binding Protein 14 on a Reduced Graphene Oxide Field‐Effect Transistor

An olfactory biosensor based on a reduced graphene oxide (rGO) field‐effect transistor (FET), functionalized by the odorant‐binding protein 14 (OBP14) from the honey bee (Apis mellifera) has been designed for the in situ and real‐time monitoring of a broad spectrum of odorants in aqueous solutions known to be attractants for bees. The electrical measurements of the binding of all tested odorants are shown to follow the Langmuir model for ligand–receptor interactions. The results demonstrate that OBP14 is able to bind odorants even after immobilization on rGO and can discriminate between ligands binding within a range of dissociation constants from K_d=4 μM to K_d=3.3 mM . The strongest ligands, such as homovanillic acid, eugenol, and methyl vanillate all contain a hydroxy group which is apparently important for the strong interaction with the protein.