Anchorage‐dependent binding of integrin I‐domain to adhesion ligands

The dynamic interactions between leukocyte integrin receptors and ligands in the vascular endothelium, extracellular matrix, or invading pathogens result in leukocyte adhesion, extravasation, and phagocytosis. This work examined the mechanical strength of the connection between iC3b, a complement component that stimulates phagocytosis, and the ligand‐binding domain, the I‐domain, of integrin α_Mβ₂. Single‐molecule force measurements of α_M I‐domain–iC3b complexes were conducted by atomic force microscope. Strikingly, depending on loading rates, immobilization of the I‐domain via its C‐terminus resulted in a 1.3‐fold to 1.5‐fold increase in unbinding force compared with I‐domains immobilized via the N‐terminus. The force spectra (unbinding force versus loading rate) of the I‐domain–iC3b complexes revealed that the enhanced mechanical strength is due to a 2.4‐fold increase in the lifetime of the I‐domain–iC3b bond. Given the structural and functional similarity of all integrin I‐domains, our result supports the existing allosteric regulatory model by which the ligand binding strength of integrin can be increased rapidly when a force is allowed to stretch the C‐terminus of the I‐domain. This type of mechanism may account for the rapid ligand affinity adjustment during leukocyte migration. Copyright © 2015 John Wiley & Sons, Ltd.     

Landscape of protein–small ligand binding modes

Elucidating the mechanisms of specific small‐molecule (ligand) recognition by proteins is a long‐standing conundrum. While the structures of these molecules, proteins and ligands, have been extensively studied, protein–ligand interactions, or binding modes, have not been comprehensively analyzed. Although methods for assessing similarities of binding site structures have been extensively developed, the methods for the computational treatment of binding modes have not been well established. Here, we developed a computational method for encoding the information about binding modes as graphs, and assessing their similarities. An all‐against‐all comparison of 20,040 protein–ligand complexes provided the landscape of the protein–ligand binding modes and its relationships with protein‐ and chemical spaces. While similar proteins in the same SCOP Family tend to bind relatively similar ligands with similar binding modes, the correlation between ligand and binding similarities was not very high (R² = 0.443). We found many pairs with novel relationships, in which two evolutionally distant proteins recognize dissimilar ligands by similar binding modes (757,474 pairs out of 200,790,780 pairs were categorized into this relationship, in our dataset). In addition, there were an abundance of pairs of homologous proteins binding to similar ligands with different binding modes (68,217 pairs). Our results showed that many interesting relationships between protein–ligand complexes are still hidden in the structure database, and our new method for assessing binding mode similarities is effective to find them.     

Love–Hate ligands for high resolution analysis of strain in ultra-stable protein/small molecule interaction

The pathway of ligand dissociation and how binding sites respond to force are not well understood for any macromolecule. Force effects on biological receptors have been studied through simulation or force spectroscopy, but not by high resolution structural experiments. To investigate this challenge, we took advantage of the extreme stability of the streptavidin–biotin interaction, a paradigm for understanding non-covalent binding as well as a ubiquitous research tool. We synthesized a series of biotin-conjugates having an unchanged strong-binding biotin moiety, along with pincer-like arms designed to clash with the protein surface: ‘Love–Hate ligands’. The Love–Hate ligands contained various 2,6-di-ortho aryl groups, installed using Suzuki coupling as the last synthetic step, making the steric repulsion highly modular. We determined binding affinity, as well as solving 1.1–1.6 Å resolution crystal structures of streptavidin bound to Love–Hate ligands. Striking distortion of streptavidin’s binding contacts was found for these complexes. Hydrogen bonds to biotin’s ureido and thiophene rings were preserved for all the ligands, but biotin’s valeryl tail was distorted from the classic conformation. Streptavidin’s L3/4 loop, normally forming multiple energetically-important hydrogen bonds to biotin, was forced away by clashes with Love–Hate ligands, but Ser45 from L3/4 could adapt to hydrogen-bond to a different part of the ligand. This approach of preparing conflicted ligands represents a direct way to visualize strained biological interactions and test protein plasticity.     

Molecular recognition force spectroscopy study of the dynamic interaction between aptamer GBI‐10 and extracellular matrix protein tenascin‐C on human glioblastoma cell

Molecular recognition force spectroscopy (MR‐FS) was applied to investigate the dynamic interaction between aptamer GBI‐10 and tenascin‐C (TN‐C) on human glioblastoma cell surface at single‐molecule level. The unbinding force between aptamer GBI‐10 and TN‐C was 39 pN at the loading rate of 0.3 nN sec^?1. A series of kinetic parameters concerning interaction process such as the unbinding force f_u, the association rate constant k_on, dissociation rate constant at zero force k_off, and dissociation constant K_D for aptamer GBI‐10/TN‐C complexes were acquired. In addition, the interaction of aptamer GBI‐10 with TN‐C depended on the presence of Mg²⁺. This work demonstrates that MR‐FS can be used as an attractive tool for exploring the interaction forces and dynamic process of aptamer and ligand at the single‐molecule level. As a future perspective, MR‐FS may be used as a potential diagnostic and therapeutic tool by combining with other techniques. Copyright © 2012 John Wiley & Sons, Ltd.     

The relationship between ligand-binding thermodynamics and protein-ligand interaction forces measured by atomic force microscopy.

The interaction forces between biotin and a set of streptavidin site-directed mutants with altered biotin-binding equilibrium and activation thermodynamics have been measured by atomic force microscopy. The AFM technique readily discriminates differences in interaction force between the site-directed (Trp to Phe or Ala) mutants. The interaction force is poorly correlated with both the equilibrium free energy of biotin binding and the activation free energy barrier to dissociation of the biotin-streptavidin complex. The interaction force is generally well correlated with the equilibrium biotin-binding enthalpy as well as the enthalpic activation barrier, but in the one mutant where these two parameters are altered in opposite directions, the interaction force is clearly correlated with the activation enthalpy of dissociation. These results suggest that the AFM force measurements directly probe the enthalpic activation barrier to ligand dissociation.     

Energy landscape of streptavidin-biotin complexes measured by atomic force microscopy

The dissociation of ligand and receptor involves multiple transitions between intermediate states formed during the unbinding process. In this paper, we explored the energy landscape of the streptavidin-biotin interaction by using the atomic force microscope (AFM) to measure the unbinding dynamics of individual ligand-receptor complexes. The rupture force of the streptavidin-biotin bond increased more than 2-fold over a range of loading rates between 100 and 5000 pN/s. Moreover, the force measurements showed two regimes of loading in the streptavidin-biotin force spectrum, revealing the presence of two activation barriers in the unbinding process. Parallel experiments carried out with a streptavidin mutant (W120F) were used to investigate the molecular determinants of the activation barriers. From these experiments, we attributed the outer activation barrier in the energy landscape to the molecular interaction of the '3-4' loop of streptavidin that closes behind biotin.     

Ligand binding and self‐association cooperativity of β‐lactoglobulin

Unlike most small globular proteins, lipocalins lack a compact hydrophobic core. Instead, they present a large central cavity that functions as the primary binding site for hydrophobic molecules. Not surprisingly, these proteins typically exhibit complex structural dynamics in solution, which is intricately modified by intermolecular recognition events. Although many lipocalins are monomeric, an increasing number of them have been proven to form oligomers. The coupling effects between self‐association and ligand binding in these proteins are largely unknown. To address this issue, we have calorimetrically characterized the recognition of dodecyl sulfate by bovine β‐lactoglobulin, which forms weak homodimers at neutral pH. A thermodynamic analysis based on coupled‐equilibria revealed that dimerization exerts disparate effects on the ligand‐binding capacity of β‐lactoglobulin. Protein dimerization decreases ligand affinity (or, reciprocally, ligand binding promotes dimer dissociation). The two subunits in the dimer exhibit a positive, entropically driven cooperativity. To investigate the structural determinants of the interaction, the crystal structure of β‐lactoglobulin bound to dodecyl sulfate was solved at 1.64 Å resolution. Copyright © 2013 John Wiley & Sons, Ltd.     

Molecular recognition of Fc‐specific ligands binding onto the consensus binding site of IgG: insights from molecular simulation

Immunoglobulin G (IgG) plays an important role in clinical diagnosis and therapeutics. Meanwhile, the consensus binding site (CBS) on the Fc domain of IgG is responsible for ligand recognition, especially for Fc‐specific ligands. In this study, molecular simulation methods were used to investigate molecular interactions between the CBS of the Fc domain and seven natural Fc‐specific ligands. The analysis on the binding energy of the Fc–ligand complex indicated that hydrophobic interactions provide the main driving force for the Fc–ligand binding processes. The hot spots on the ligands and Fc were identified with the computational alanine scanning approach. It was found that the residues of tryptophan and tyrosine on the ligands have significant contributions for the Fc–ligand binding, while Met252, Ile253, Asn434, His435, and Tyr436 are the key residues of Fc. Moreover, two binding modes based on tryptophan or tyrosine were summarized and constructed according to the pairwise interaction analysis. Guidelines for the rational design of CBS‐specific ligands with high affinity and specificity were proposed. Copyright © 2014 John Wiley & Sons, Ltd.     

Conformational diversity of ligands bound to proteins

The phenomenon of molecular recognition, which underpins almost all biological processes, is dynamic, complex and subtle. Establishing an interaction between a pair of molecules involves mutual structural rearrangements guided by a highly convoluted energy landscape, the accurate mapping of which continues to elude us. Increased understanding of the degree to which the conformational space of a ligand is restricted upon binding may have important implications for docking studies, structure refinement and for function prediction methods based on geometrical comparisons of ligands or their binding sites. Here, we present an analysis of the conformational variability exhibited by three of the most ubiquitous biological ligands in nature, ATP, NAD and FAD. First, we demonstrate qualitatively that these ligands bind to proteins in widely varying conformations, including several cases in which parts of the molecule assume energetically unfavourable orientations. Next, by comparing the distribution of bound ligand shapes with the set of all possible molecular conformations, we provide a quantitative assessment of previous observations that ligands tend to unfold when binding to proteins. We show that, while extended forms of ligands are indeed common in ligand-protein structures, instances of ligands in almost maximally compact arrangements can also be found. Thirdly, we compare the conformational variation in two sets of ligand molecules, those bound to homologous proteins, and those bound to unrelated proteins. Although most superfamilies bind ligands in a fairly conserved manner, we find several cases in which significant variation in ligand configuration is observed.     

Binding specificity and the ligand dissociation process in the E. coli biotin holoenzyme synthetase

The binding of the Escherichia coli biotin holoenzyme synthetase to the two ligands, biotin and bio-5'-AMP, is coupled to disorder-to-order transitions in the protein. In the structure of the biotin complex, a "glycine-rich" loop that is disordered in the apo-enzyme is folded over the ligand. Mutations in three residues in this loop result in significant changes in the affinity of the enzyme for both biotin and bio-5'-AMP. The kinetic basis of these losses in the affinity resides primarily in changes in the unimolecular rates of dissociation of the complexes. In this work, isothermal titration calorimetry has been employed to examine the detailed thermodynamics of binding of three loop mutants to biotin and bio-5'-AMP. The energetic features of dissociation of the protein*ligand complexes also have been probed by measuring the temperature dependencies of the unimolecular dissociation rates. Analysis of the data using the Eyring formalism yielded entropic and enthalpic contributions to the energetic barrier to dissociation. The thermodynamic results coupled with the known structures of the apo-enzyme and biotin complex have been used to formulate a model for progression from the ground-state complex to the transition state in biotin dissociation. In this model, the transition-state is characterized by both partial disruption of noncovalent bonds and acquisition of some of the disorder that characterizes the glycine-rich loop in the absence of ligand.     

Characterization of enhanced monovalent and bivalent thrombin DNA aptamer binding using single molecule force spectroscopy

Thrombin aptamer binding strength and stability is dependent on sterical parameters when used for atomic force microscopy sensing applications. Sterical improvements on the linker chemistry were developed for high-affinity binding. For this we applied single molecule force spectroscopy using two enhanced biotinylated thrombin aptamers, BFF and BFA immobilized on the atomic force microscopy tip via streptavidin. BFF is a dimer composed of two single-stranded aptamers (aptabody) connected to each other by a complementary sequence close to the biotinylated end. In contrast, BFA consists of a single DNA strand and a complementary strand in the supporting biotinylated part. By varying the pulling velocity in force-distance cycles the formed thrombin-aptamer complexes were ruptured at different force loadings allowing determination of the energy landscape. As a result, BFA aptamer showed a higher binding force at the investigated loading rates and a significantly lower dissociation rate constant, k_off, compared to BFF. Moreover, the potential of the aptabody BFF to form a bivalent complex could clearly be demonstrated.     

Comparison of human and mouse E-selectin binding to Sialyl-Lewis<Superscript>x</Superscript>

Background

During inflammation, leukocytes are captured by the selectin family of adhesion receptors lining blood vessels to facilitate exit from the bloodstream. E-selectin is upregulated on stimulated endothelial cells and binds to several ligands on the surface of leukocytes. Selectin:ligand interactions are mediated in part by the interaction between the lectin domain and Sialyl-Lewis x (sLe^x), a tetrasaccharide common to selectin ligands. There is a high degree of homology between selectins of various species: about 72 and 60 % in the lectin and EGF domains, respectively. In this study, molecular dynamics, docking, and steered molecular dynamics simulations were used to compare the binding and dissociation mechanisms of sLe^x with mouse and human E-selectin. First, a mouse E-selectin homology model was generated using the human E-selectin crystal structure as a template.

Results

Mouse E-selectin was found to have a greater interdomain angle, which has been previously shown to correlate with stronger binding among selectins. sLe^x was docked onto human and mouse E-selectin, and the mouse complex was found to have a higher free energy of binding and a lower dissociation constant, suggesting stronger binding. The mouse complex had higher flexibility in a few key residues. Finally, steered molecular dynamics was used to dissociate the complexes at force loading rates of 2000–5000 pm/ps². The mouse complex took longer to dissociate at every force loading rate and the difference was statistically significant at 3000 pm/ps². When sLe^x-coated microspheres were perfused through microtubes coated with human or mouse E-selectin, the particles rolled more slowly on mouse E-selectin.

Conclusions

Both molecular dynamics simulations and microsphere adhesion experiments show that mouse E-selectin protein binds more strongly to sialyl Lewis x ligand than human E-selectin. This difference was explained by a greater interdomain angle for mouse E-selectin, and greater flexibility in key residues. Future work could introduce similar amino acid substitutions into the human E-selectin sequence to further modulate adhesion behavior.

     

A modular design of low‐background bioassays based on a high‐affinity molecular pair barstar:barnase

High‐affinity molecular pairs provide a convenient and flexible modular base for the design of molecular probes and protein/antigen assays. Specificity and sensitivity performance indicators of a bioassay critically depend on the dissociation constant (K_D) of the molecular pair, with avidin:biotin being the state‐of‐the‐art molecular pair (K_D ～ 1 fM) used almost universally for applications in the fields of nanotechnology and proteomics. In this paper, we present an alternative high‐affinity protein pair, barstar:barnase (K_D ～ 10 fM), which addresses several shortfalls of the avidin:biotin system, including non‐negligible background due to the non‐specific binding. A quantitative assessment of the non?specific binding carried out using a model assay revealed inherent irreproducibility of the [strept]avidin:biotin‐based assays, attributed to the avidin binding to solid phases, endogenous biotin molecules and serum proteins. On the other hand, the model assays assembled via a barstar:barnase protein linker proved to be immune to such non‐specific binding, showing good prospects for high‐sensitivity rare biomolecular event nanoproteomic assays.     

Molecular dynamics simulation of the induced‐fit binding process of DNA aptamer and L‐argininamide

Aptamers are rare functional nucleic acids with binding affinity to and specificity for target ligands. Recent experiments have lead to the proposal of an induced‐fit binding mechanism for L ‐argininamide (Arm) and its binding aptamer. However, at the molecular level, this mechanism between the aptamer and its coupled ligand is still poorly understood. The present study used explicit solvent molecular dynamics (MD) simulations to examine the critical bases involved in aptamer‐Arm binding and the induced‐fit binding process at atomic resolution. The simulation results revealed that the Watson‐Crick pair (G10‐C16), C9, A12, and C17 bases play important roles in aptamer‐Arm binding, and that binding of Arm results in an aptamer conformation optimized through a general induced‐fit process. In an aqueous solution, the mechanism has the following characteristic stages: (a) adsorption stage, the Arm anchors to the binding site of aptamer with strong electrostatic interaction; (b) binding stage, the Arm fits into the binding site of aptamer by hydrogen‐bond formation; and (c) complex stabilization stage, the hydrogen bonding and electrostatic interactions cooperatively stabilize the complex structure. This study provides dynamics information on the aptamer‐ligand induced‐fit binding mechanism. The critical bases in aptamer‐ligand binding may provide a guideline in aptamer design for molecular recognition engineering.     

Different combinations of atomic interactions predict protein‐small molecule and protein‐DNA/RNA affinities with similar accuracy

Interactions between proteins and other molecules play essential roles in all biological processes. Although it is widely held that a protein's ligand specificity is determined primarily by its three‐dimensional structure, the general principles by which structure determines ligand binding remain poorly understood. Here we use statistical analyses of a large number of protein?ligand complexes with associated binding‐affinity measurements to quantitatively characterize how combinations of atomic interactions contribute to ligand affinity. We find that there are significant differences in how atomic interactions determine ligand affinity for proteins that bind small chemical ligands, those that bind DNA/RNA and those that interact with other proteins. Although protein‐small molecule and protein‐DNA/RNA binding affinities can be accurately predicted from structural data, models predicting one type of interaction perform poorly on the others. Additionally, the particular combinations of atomic interactions required to predict binding affinity differed between small‐molecule and DNA/RNA data sets, consistent with the conclusion that the structural bases determining ligand affinity differ among interaction types. In contrast to what we observed for small‐molecule and DNA/RNA interactions, no statistical models were capable of predicting protein?protein affinity with >60% correlation. We demonstrate the potential usefulness of protein‐DNA/RNA binding prediction as a possible tool for high‐throughput virtual screening to guide laboratory investigations, suggesting that quantitative characterization of diverse molecular interactions may have practical applications as well as fundamentally advancing our understanding of how molecular structure translates into function. Proteins 2015; 83:2100–2114. © 2015 The Authors. Proteins: Structure, Function, and Bioinformatics Published by Wiley Periodicals, Inc.     

A ruthenium dimer complex with a flexible linker slowly threads between DNA bases in two distinct steps

Several multi-component DNA intercalating small molecules have been designed around ruthenium-based intercalating monomers to optimize DNA binding properties for therapeutic use. Here we probe the DNA binding ligand [μ-C4(cpdppz)₂(phen)₄Ru₂]⁴⁺, which consists of two Ru(phen)₂dppz²⁺ moieties joined by a flexible linker. To quantify ligand binding, double-stranded DNA is stretched with optical tweezers and exposed to ligand under constant applied force. In contrast to other bis-intercalators, we find that ligand association is described by a two-step process, which consists of fast bimolecular intercalation of the first dppz moiety followed by ∼10-fold slower intercalation of the second dppz moiety. The second step is rate-limited by the requirement for a DNA-ligand conformational change that allows the flexible linker to pass through the DNA duplex. Based on our measured force-dependent binding rates and ligand-induced DNA elongation measurements, we are able to map out the energy landscape and structural dynamics for both ligand binding steps. In addition, we find that at zero force the overall binding process involves fast association (∼10 s), slow dissociation (∼300 s), and very high affinity (K_d ∼10 nM). The methodology developed in this work will be useful for studying the mechanism of DNA binding by other multi-step intercalating ligands and proteins.     

AFM force spectroscopy reveals how subtle structural differences affect the interaction strength between Candida albicans and DC‐SIGN

The fungus Candida albicans is the most common cause of mycotic infections in immunocompromised hosts. Little is known about the initial interactions between Candida and immune cell receptors, such as the C‐type lectin dendritic cell‐specific intracellular cell adhesion molecule‐3 (ICAM‐3)‐grabbing non‐integrin (DC‐SIGN), because a detailed characterization at the structural level is lacking. DC‐SIGN recognizes specific Candida‐associated molecular patterns, that is, mannan structures present in the cell wall of Candida. The molecular recognition mechanism is however poorly understood. We postulated that small differences in mannan‐branching may result in considerable differences in the binding affinity. Here, we exploit atomic force microscope‐based dynamic force spectroscopy with single Candida cells to gain better insight in the carbohydrate recognition capacity of DC‐SIGN. We demonstrate that slight differences in the N‐mannan structure of Candida, that is, the absence or presence of a phosphomannan side chain, results in differences in the recognition by DC‐SIGN as follows: (i) it contributes to the compliance of the outer cell wall of Candida, and (ii) its presence results in a higher binding energy of 1.6 k_BT. The single‐bond affinity of tetrameric DC‐SIGN for wild‐type C. albicans is ~10.7 k_BT and a dissociation constant k_D of 23 μM, which is relatively strong compared with other carbohydrate–protein interactions described in the literature. In conclusion, this study shows that DC‐SIGN specifically recognizes mannan patterns on C. albicans with high affinity. Knowledge on the binding pocket of DC‐SIGN and its pathogenic ligands will lead to a better understanding of how fungal‐associated carbohydrate structures are recognized by receptors of the immune system and can ultimately contribute to the development of new anti‐fungal drugs. Copyright © 2015 John Wiley & Sons, Ltd.     

Structural studies of bovine,equine, and leporine serum albumin complexes with naproxen

Serum albumin, a protein naturally abundant in blood plasma, shows remarkable ligand binding properties of numerous endogenous and exogenous compounds. Most of serum albumin binding sites are able to interact with more than one class of ligands. Determining the protein‐ligand interactions among mammalian serum albumins is essential for understanding the complexity of this transporter. We present three crystal structures of serum albumins in complexes with naproxen (NPS): bovine (BSA‐NPS), equine (ESA‐NPS), and leporine (LSA‐NPS) determined to 2.58 Å (C2), 2.42 Å (P6₁), and 2.73 Å (P2₁2₁2₁) resolutions, respectively. A comparison of the structurally investigated complexes with the analogous complex of human serum albumin (HSA‐NPS) revealed surprising differences in the number and distribution of naproxen binding sites. Bovine and leporine serum albumins possess three NPS binding sites, but ESA has only two. All three complexes of albumins studied here have two common naproxen locations, but BSA and LSA differ in the third NPS binding site. None of these binding sites coincides with the naproxen location in the HSA‐NPS complex, which was obtained in the presence of other ligands besides naproxen. Even small differences in sequences of serum albumins from various species, especially in the area of the binding pockets, influence the affinity and the binding mode of naproxen to this transport protein. Proteins 2014; 82:2199–2208. © 2014 Wiley Periodicals, Inc.     

Topology and thermodynamics of gaseous ligands diffusion paths in human neuroglobin

The physiological role of recently discovered human neuroglobin (Ngb) is still unknown. Sound hypothesis says that it protects brain during hypoxia. In this paper the advanced potential of mean force by implicit ligand sampling (PMF/ILS) method is used to study the free energy landscape of Ngb for O₂, NO and CO ligands. The multiple diffusion paths are discovered and four ligand binding cavities are determined. The data show that certain regions are easily accessible by O₂ and NO but are protected from CO. Free energy landscapes provide realistic data for stochastic models of ligand diffusion in proteins.