On the Role of Flexibility in Protein–Ligand Interactions: the Example of p53 Tetramerization Domain

The recognition of protein surfaces by designed ligands has become an attractive approach in drug discovery. However, the variable nature and irregular behavior of protein surfaces defy this new area of research. The easy to understand “lock‐and‐key” model is far from being the ideal paradigm in biomolecular interactions and, hence, any new finding on how proteins and ligands behave in recognition events paves a step of the way. Herein, we illustrate a clear example on how an increase in flexibility of both protein and ligand can result in an increase in the stability of the macromolecular complex. The biophysical study of the interaction between a designed flexible tetraguanidinium‐calix[4]arene and the tetramerization domain of protein p53 (p53TD) and its natural mutant p53TD‐R337H shows how the floppy mutant domain interacts more tightly with the ligand than the well‐packed wild‐type protein. Moreover, the flexible calixarene ligand interacts with higher affinity to both wild‐type and mutated protein domains than a conformationally rigid calixarene analog previously reported. These findings underscore the crucial role of flexibility in molecular recognition processes, for both small ligands and large biomolecular surfaces.     

Molecular Tweezers: Concepts and Applications

Taken to the molecular level, the concept of “tweezers” opens a rich and fascinating field at the convergence of molecular recognition, biomimetic chemistry and nanomachines. Composed of a spacer bridging two interaction sites, the behaviour of molecular tweezers is strongly influenced by the flexibility of their spacer. Operating through an “induced‐fit” recognition mechanism, flexible molecular tweezers select the conformation(s) most appropriate for substrate binding. Their adaptability allows them to be used in a variety of binding modes and they have found applications in chirality signalling. Rigid spacers, on the contrary, display a limited number of binding states, which lead to selective and strong substrate binding following a “lock and key” model. Exquisite selectivity may be expressed with substrates as varied as C₆₀, nanotubes and natural cofactors, and applications to molecular electronics and enzyme inhibition are emerging. At the crossroad between flexible and rigid spacers, stimulus‐responsive molecular tweezers controlled by ionic, redox or light triggers belong to the realm of molecular machines, and, applied to molecular tweezing, open doors to the selective binding, transport and release of their cargo. Applications to controlled drug delivery are already appearing. The past 30 years have seen the birth of molecular tweezers; the next many years to come will surely see them blooming in exciting applications.     

Numerical interpretation of molecular surface field in dielectric modeling of solvation

Continuum solvent models, particularly those based on the Poisson‐Boltzmann equation (PBE), are widely used in the studies of biomolecular structures and functions. Existing PBE developments have been mainly focused on how to obtain more accurate and/or more efficient numerical potentials and energies. However to adopt the PBE models for molecular dynamics simulations, a difficulty is how to interpret dielectric boundary forces accurately and efficiently for robust dynamics simulations. This study documents the implementation and analysis of a range of standard fitting schemes, including both one‐sided and two‐sided methods with both first‐order and second‐order Taylor expansions, to calculate molecular surface electric fields to facilitate the numerical calculation of dielectric boundary forces. These efforts prompted us to develop an efficient approximated one‐dimensional method, which is to fit the surface field one dimension at a time, for biomolecular applications without much compromise in accuracy. We also developed a surface‐to‐atom force partition scheme given a level set representation of analytical molecular surfaces to facilitate their applications to molecular simulations. Testing of these fitting methods in the dielectric boundary force calculations shows that the second‐order methods, including the one‐dimensional method, consistently perform among the best in the molecular test cases. Finally, the timing analysis shows the approximated one‐dimensional method is far more efficient than standard second‐order methods in the PBE force calculations. © 2017 Wiley Periodicals, Inc.     

Molecular Recognition of Zwitterions with Artificial Receptors

Many biomolecules exist as internal ion pairs or zwitterions within a biologically relevant pH range. Despite their importance, the molecular recognition of this type of systems is specially challenging due to their strong solvation in aqueous media, and their trend to form folded or self‐assembled structures by pairing of charges of different sign. In this Minireview, we will discuss the molecular recognition of zwitterions using non‐natural, synthetic receptors. This contribution does not intend to make a full in‐depth revision of the existing research in the field, but a personal overview with selected representative examples from the recent literature.     

Modulation of Kinetic Pathways of Enzyme–Substrate Interaction in a Microfluidic Channel: Nanoscopic Water Dynamics as a Switch

Enzyme-mediated catalysis is attributed to enzyme–substrate interactions, with models such as “induced fit” and “conformational selection” emphasizing the role of protein conformational transitions. The dynamic nature of the protein structure, thus, plays a crucial role in molecular recognition and substrate binding. As large-scale protein motions are coupled to water motions, hydration dynamics play a key role in protein dynamics, and hence, in enzyme catalysis. Here, microfluidic techniques and time-dependent fluorescence Stokes shift (TDFSS) measurements are employed to elucidate the role of nanoscopic water dynamics in the interaction of an enzyme, α-Chymotrypsin (CHT), with a substrate, Ala-Ala-Phe-7-amido-4-methylcoumarin (AMC) in the cationic reverse micelles of benzylhexadecyldimethylammonium chloride (BHDC/benzene) and anionic reverse micelles of sodium bis(2-ethylhexyl)sulfosuccinate (AOT/benzene). The kinetic pathways unraveled from the microfluidic setup are consistent with the “conformational selection” fit for the interaction of CHT with AMC in the cationic reverse micelles, whereas an “induced fit” mechanism is indicated for the anionic reverse micelles. In the cationic reverse micelles of BHDC, faster hydration dynamics (≈550 ps) aid the pathway of “conformational selection”, whereas in the anionic reverse micelles of AOT, the significantly slower dynamics of hydration (≈1600 ps) facilitate an “induced fit” mechanism for the formation of the final enzyme–substrate complex. The role of water dynamics in dictating the mechanism of enzyme–substrate interaction becomes further manifest in the neutral reverse micelles of Brij-30 and Triton X-100. In the former, the faster water dynamics aid the “conformational selection” pathway, whereas the significantly slower dynamics of water molecules in the latter are conducive to the “induced fit” mechanism in the enzyme–substrate interaction. Thus, nanoscopic water dynamics act as a switch in modulating the pathway of recognition of an enzyme (CHT) by the substrate (AMC) in reverse micelles.     

Development of a Dual‐Functional Hydrogel Using RGD and Anti‐VEGF Aptamer

Synthetic molecular libraries hold great potential to advance the biomaterial development. However, little effort is made to integrate molecules with molecular recognition abilities selected from different libraries into a single biomolecular material. The purpose of this work is to incorporate peptides and nucleic acid aptamers into a porous hydrogel to develop a dual‐functional biomaterial. The data show that an anti‐integrin peptide can promote the attachment and growth of endothelial cells in a 3D porous poly(ethylene glycol) hydrogel and an antivascular endothelial growth factor aptamer can sequester and release VEGF of high bioactivity. Importantly, the dual‐functional porous hydrogel enhances the growth and survival of endothelial cells. This work demonstrates that molecules selected from different synthetic libraries can be integrated into one system for the development of novel biomaterials.     

Molecular mechanics models for tetracycline analogs

Tetracyclines (Tcs) are an important family of antibiotics that bind to the ribosome and several proteins. To model Tc interactions with protein and RNA, we have developed a molecular mechanics force field for 12 tetracyclines, consistent with the CHARMM force field. We considered each Tc variant in its zwitterionic tautomer, with and without a bound Mg(2+). We used structures from the Cambridge Crystallographic Data Base to identify the conformations likely to be present in solution and in biomolecular complexes. A conformational search by simulated annealing was undertaken, using the MM3 force field, for tetracycline, anhydrotetracycline, doxycycline, and tigecycline. Resulting, low-energy structures were optimized with an ab initio method. We found that Tc and its analogs all adopt an extended conformation in the zwitterionic tautomer and a twisted one in the neutral tautomer, and the zwitterionic-extended state is the most stable in solution. Intermolecular force field parameters were derived from a standard supermolecule approach: we considered the ab initio energies and geometries of a water molecule interacting with each Tc analog at several different positions. The final, rms deviation between the ab initio and force field energies, averaged over all forms, was 0.35 kcal/mol. Intramolecular parameters were adopted from either the standard CHARMM force field, the ab initio structure, or the earlier, plain Tc force field. The model reproduces the ab initio geometry and flexibility of each Tc. As tests, we describe MD and free energy simulations of a solvated complex between three Tcs and the Tet repressor protein.     

Comparison of structural,thermodynamic, kinetic and mass transport properties of Mg2+ ion models commonly used in biomolecular simulations

The prevalence of Mg²⁺ ions in biology and their essential role in nucleic acid structure and function has motivated the development of various Mg²⁺ ion models for use in molecular simulations. Currently, the most widely used models in biomolecular simulations represent a nonbonded metal ion as an ion‐centered point charge surrounded by a nonelectrostatic pairwise potential that takes into account dispersion interactions and exchange effects that give rise to the ion's excluded volume. One strategy toward developing improved models for biomolecular simulations is to first identify a Mg²⁺ model that is consistent with the simulation force fields that closely reproduces a range of properties in aqueous solution, and then, in a second step, balance the ion–water and ion–solute interactions by tuning parameters in a pairwise fashion where necessary. The present work addresses the first step in which we compare 17 different nonbonded single‐site Mg²⁺ ion models with respect to their ability to simultaneously reproduce structural, thermodynamic, kinetic and mass transport properties in aqueous solution. None of the models based on a 12‐6 nonelectrostatic nonbonded potential was able to reproduce the experimental radial distribution function, solvation free energy, exchange barrier and diffusion constant. The models based on a 12‐6‐4 potential offered improvement, and one model in particular, in conjunction with the SPC/E water model, performed exceptionally well for all properties. The results reported here establish useful benchmark calculations for Mg²⁺ ion models that provide insight into the origin of the behavior in aqueous solution, and may aid in the development of next‐generation models that target specific binding sites in biomolecules. © 2015 Wiley Periodicals, Inc.     

Potential of mean force of water–proton bath and molecular dynamic simulation of proteins at constant pH

An advanced implicit solvent model of water–proton bath for protein simulations at constant pH is presented. The implicit water–proton bath model approximates the potential of mean force of a protein in water solvent in a presence of hydrogen ions. Accurate and fast computational implementation of the implicit water–proton bath model is developed using the continuum electrostatic Poisson equation model for calculation of ionization equilibrium and the corrected MSR6 generalized Born model for calculation of the electrostatic atom–atom interactions and forces. Molecular dynamics (MD) method for protein simulation in the potential of mean force of water–proton bath is developed and tested on three proteins. The model allows to run MD simulations of proteins at constant pH, to calculate pH‐dependent properties and free energies of protein conformations. The obtained results indicate that the developed implicit model of water–proton bath provides an efficient way to study thermodynamics of biomolecular systems as a function of pH, pH‐dependent ionization‐conformation coupling, and proton transfer events. © 2012 Wiley Periodicals, Inc.     

Efficient quantification of the importance of contacts for the dynamical stability of proteins

Understanding the stability of the native state and the dynamics of a protein is of great importance for all areas of biomolecular design. The efficient estimation of the influence of individual contacts between amino acids in a protein structure is a first step in the reengineering of a particular protein for technological or pharmacological purposes. At the same time, the functional annotation of molecular evolution can be facilitated by such insight. Here, we use a recently suggested, information theoretical measure in biomolecular design - the Kullback-Leibler-divergence - to quantify and therefore rank residue-residue contacts within proteins according to their overall contribution to the molecular mechanics. We implement this protocol on the basis of a reduced molecular model, which allows us to use a well-known lemma of linear algebra to speed up the computation. The increase in computational performance is around 10(1)- to 10(4)-fold. We applied the method to two proteins to illustrate the protocol and its results. We found that our method can reliably identify key residues in the molecular mechanics and the protein fold in comparison to well-known properties in the serine protease inhibitor. We found significant correlations to experimental results, e.g., dissociation constants and Φ values.     

Enhanced stability of the model mini‐protein in amino acid ionic liquids and their aqueous solutions

Using molecular dynamics simulations, the structure of model mini‐protein was thoroughly characterized in the imidazolium‐based amino acid ionic liquids and their aqueous solutions. Complete substitution of water by organic cations and anions further results in hindered conformational flexibility of the mini‐protein. This observation suggests that amino acid‐based ionic liquids are able to defend proteins from thermally induced denaturation. We show by means of radial distributions that the mini‐protein is efficiently solvated by both solvents due to a good mutual miscibility. Amino acid‐based anions prevail in the first coordination sphere of positively charged sites of the mini‐protein whereas water molecules prevail in the first coordination sphere of negatively charged sites of the mini‐protein. © 2015 Wiley Periodicals, Inc.     

Shape‐Persistent and Adaptive Multivalency: Rigid Transgeden (TGD) and Flexible PAMAM Dendrimers for Heparin Binding

This study investigates transgeden (TGD) dendrimers (polyamidoamine (PAMAM)‐type dendrimers modified with rigid polyphenylenevinylene (PPV) cores) and compares their heparin‐binding ability with commercially available PAMAM dendrimers. Although the peripheral ligands are near‐identical between the two dendrimer families, their heparin binding is very different. At low generation (G1), TGD outperforms PAMAM, but at higher generation (G2 and G3), the PAMAMs are better. Heparin binding also depends strongly on the dendrimer/heparin ratio. We explain these effects using multiscale modelling. TGD dendrimers exhibit “shape‐persistent multivalency”; the rigidity means that small clusters of surface amines are locally well optimised for target binding, but it prevents the overall nanoscale structure from rearranging to maximise its contacts with a single heparin chain. Conversely, PAMAM dendrimers exhibit “adaptive multivalency”; the flexibility means individual surface ligands are not so well optimised locally to bind heparin chains, but the nanostructure can adapt more easily and maximise its binding contacts. As such, this study exemplifies important new paradigms in multivalent biomolecular recognition.     

Controlling Electrical Conductance through a π‐Conjugated Cruciform Molecule by Selective Anchoring to Gold Electrodes

Tuning charge transport at the single‐molecule level plays a crucial role in the construction of molecular electronic devices. Introduced herein is a promising and operationally simple approach to tune two distinct charge‐transport pathways through a cruciform molecule. Upon in situ cleavage of triisopropylsilyl groups, complete conversion from one junction type to another is achieved with a conductance increase by more than one order of magnitude, and it is consistent with predictions from ab initio transport calculations. Although molecules are well known to conduct through different orbitals (either HOMO or LUMO), the present study represents the first experimental realization of switching between HOMO‐ and LUMO‐dominated transport within the same molecule.     

Generative topographic mapping of binding pocket of β2 receptor and three‐way partial least squares modeling of inhibitory activities

The visualization and characterization of protein pockets is the starting point for many structure‐based drug design projects. The size and shape of protein pockets dictate 3D geometry of ligands that can strongly inhibit the following biological events. Thus, a minimal requirement for inhibition is that a molecule sterically binds the active site with some allowance for induced fit. Methods for direct display of active sites in a protein have become prevalent in recent years. In this study, a new mapping method, generative topographic mapping, is investigated to describe the 3D surface of protein pocket. The β2 receptor protein is used as a benchmark. By mapping the molecular surface points and assigning the associated molecular electrostatic potential (MEP) values, the original 3D structure of the active site is well reproduced by the 2D latent map in generative topographic mapping. The distributions of MEP values of two 2D latent maps derived from the inhibitor and the β2 receptor protein are well complemented. Using three‐way partial least squares modeling, a predictive model linking the inhibitory activity and their MEP values can be constructed, which was not feasible in the previous spherical self‐organizing map studies. The resulting regression coefficient matrix of the three‐way partial least squares model has many insights for understanding the structural requirements for β2 inhibitory activity. Copyright © 2014 John Wiley & Sons, Ltd.     

Balancing target flexibility and target denaturation in computational fragment‐based inhibitor discovery

Accounting for target flexibility and selecting “hot spots” most likely to be able to bind an inhibitor continue to be challenges in the field of structure‐based drug design, especially in the case of protein–protein interactions. Computational fragment‐based approaches using molecular dynamics (MD) simulations are a promising emerging technology having the potential to address both of these challenges. However, the optimal MD conditions permitting sufficient target flexibility while also avoiding fragment‐induced target denaturation remain ambiguous. Using one such technology (Site Identification by Ligand Competitive Saturation, SILCS), conditions were identified to either prevent denaturation or identify and exclude trajectories in which subtle but important denaturation was occurring. The target system used was the well‐characterized protein cytokine IL‐2, which is involved in a protein–protein interface and, in its unliganded crystallographic form, lacks surface pockets that can serve as small‐molecule binding sites. Nonetheless, small‐molecule inhibitors have previously been discovered that bind to two “cryptic” binding sites that emerge only in the presence of ligand binding, highlighting the important role of IL‐2 flexibility. Using the above conditions, SILCS with hydrophobic fragments was able to identify both sites based on favorable fragment binding while avoiding IL‐2 denaturation. An important additional finding was that acetonitrile, a water‐miscible fragment, fails to identify either site yet can induce target denaturation, highlighting the importance of fragment choice. © 2012 Wiley Periodicals, Inc.     

Live‐Cell‐Templated Dynamic Combinatorial Chemistry

Dynamic covalent chemistry combines in a single step the screening and synthesis of ligands for biomolecular recognition. In order to do that, a chemical entity is used as template within a dynamic combinatorial library of interconverting species, so that the stronger binders are amplified due to the efficient interaction with the target. Here we employed whole A549 living cells as template in a dynamic mixture of imines, for which amplification reflects the efficient and selective interaction with the corresponding extracellular matrix. The amplified polyamine showed strong interaction with the A549 extracellular matrix in on‐cell NMR experiments, while combination of NMR, SPR, and molecular dynamics simulations in model systems provided insights on the molecular recognition event. Notably, our work pioneers the use of whole living cells in dynamic combinatorial chemistry, which paves the way towards the discovery of new bioactive molecules in a more biorelevant environment.     

(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.     

Energy Landscape of Aptamer/Protein Complexes Studied by Single‐Molecule Force Spectroscopy

Aptamers are single‐stranded nucleic acid molecules selected in vitro to bind to a variety of target molecules. Aptamers bound to proteins are emerging as a new class of molecules that rival commonly used antibodies in both therapeutic and diagnostic applications. With the increasing application of aptamers as molecular probes for protein recognition, it is important to understand the molecular mechanism of aptamer–protein interaction. Recently, we developed a method of using atomic force microscopy (AFM) to study the single‐molecule rupture force of aptamer/protein complexes. In this work, we investigate further the unbinding dynamics of aptamer/protein complexes and their dissociation‐energy landscape by AFM. The dependence of single‐molecule force on the AFM loading rate was plotted for three aptamer/protein complexes and their dissociation rate constants, and other parameters characterizing their dissociation pathways were obtained. Furthermore, the single‐molecule force spectra of three aptamer/protein complexes were compared to those of the corresponding antibody/protein complexes in the same loading‐rate range. The results revealed two activation barriers and one intermediate state in the unbinding process of aptamer/protein complexes, which is different from the energy landscape of antibody/protein complexes. The results provide new information for the study of aptamer–protein interaction at the molecular level.     

Self‐Sorting of Two Hydrocarbon Receptors with One Carbonaceous Ligand

Non‐directional van der Waals forces in biological and synthetic supramolecular systems play important roles in molecular assembly, particularly in determining the distances of the interacting species. The van der Waals forces are normally used in combination with other directional forces and are considered to play a secondary role in achieving specificity and fidelity in molecular recognition. Using an ideal supramolecular system consisting solely of hydrogen and carbon atoms, we found that the van der Waals interactions enable the high‐fidelity sorting of two homomeric receptors during ligand‐induced assembly. The self‐sorting occurred in a narcissistic manner by repulsion of a competing diastereoisomeric receptor from the assembly. The structure–sorting relationship study with enantiomers further revealed the dominant role of the van der Waals forces in shape recognition for high‐fidelity self‐sorting.     

Covalent Assembly and Characterization of Nonsymmetrical Single‐Molecule Nodes

The covalent linking of molecular building blocks on surfaces enables the construction of specific molecular nanostructures of well‐defined shape. Molecular nodes linked to various entities play a key role in such networks, but represent a particular challenge because they require a well‐defined arrangement of different building blocks. Herein, we describe the construction of a chemically and geometrically well defined covalent architecture made of one central node and three molecular wires arranged in a nonsymmetrical way and thus encoding different conjugation pathways. Very different architectures of either very limited or rather extended size were obtained depending on the building blocks used for the covalent linking process on the Au(111) surface. Electrical measurements were carried out by pulling individual molecular nodes with the tip of a scanning tunneling microscope. The results of this challenging procedure indicate subtle differences if the nodes are contacted at inequivalent termini.