Configurational entropy and cooperativity between ligand binding and dimerization in glycopeptide antibiotics

Crystal structure of vancomycin bound to the resistance determinant D-alanine-D-serine

Vancomycin is a glycopeptide antibiotic that for decades has been a mainstay of treatment for persistent bacterial infections. However, the spread of antibiotic resistance threatens its continued utility. In particular, vancomycin-resistant enterococci (VRE) have become a pressing clinical challenge. Vancomycin acts by binding and sequestering the intermediate Lipid II in cell-wall biosynthesis, specifically recognizing a D-alanine-D-alanine dipeptide motif within the Lipid II molecule. VRE achieve resistance by remodeling this motif to either D-alanine-D-lactate or D-alanine-D-serine; the former substitution essentially abolishes recognition by vancomycin of Lipid II, whereas the latter reduces the affinity of the antibiotic by roughly one order of magnitude. The complex of vancomycin bound to D-alanine-D-serine has been crystallized, and its 1.20 Å X-ray crystal structure is presented here. This structure reveals that the D-alanine-D-serine ligand is bound in essentially the same position and same pose as the native D-alanine-D-alanine ligand. The serine-containing ligand appears to be slightly too large to be comfortably accommodated in this way, suggesting one possible contribution to the reduced binding affinity. In addition, two flexible hydroxyl groups - one from the serine side chain of the ligand, and the other from a glucose sugar on the antibiotic - are locked into single conformations in the complex, which is likely to contribute an unfavorable entropic component to the recognition of the serine-containing ligand.

Pursuit of next-generation glycopeptides: a journey with vancomycin

Vancomycin, a blockbuster antibiotic of the glycopeptide class, has been a life-saving therapeutic against multidrug-resistant Gram-positive infections. The emergence of glycopeptide resistance has however enunciated the need to develop credible alternatives with potent activity against vancomycin-resistant bacteria. Medicinal chemistry has responded to this challenge through various strategies, one of them being the development of semisynthetic analogues. Many groups, including ours, have been contributing towards the development of semisynthetic vancomycin analogues to tackle vancomycin-resistant bacteria. In this feature article, we have discussed our research contribution to the field of glycopeptides, which includes our strategies and designs of vancomycin analogues incorporating multimodal mechanisms of action. The strategies discussed here, such as conferring membrane activity, enhanced binding to target, multivalency, etc. involve semisynthetic modifications to vancomycin at the carboxy terminal and the amino group of the vancosamine sugar of vancomycin, to develop novel analogues. These analogues have demonstrated their superior efficacy in tackling the inherited forms of vancomycin resistance in Gram-positive and Gram-negative bacteria, including highly drug-resistant strains. More importantly, these analogues also possess the ability to tackle various non-inherited forms of bacterial resistance, such as metabolically dormant stationary-phase and persister cells, bacterial biofilms, and intracellular pathogens. Our derivatives also display superior pharmacokinetics, and less propensity for resistance development, owing to their different modes of action. Through this feature article, we present to the reader a concise picture of the multitude of approaches that can be used to tackle different types of resistance through semisynthetic modifications to vancomycin. We have also highlighted the challenges and lacunae in the field, and potential directions which future research can explore.

Hydrophobic residues of melittin mediate its binding to αA-crystallin

The molecular chaperone αA-crystallin, mainly localized in the human ocular lens, is believed to protect the lens from opacification and cataract, by suppressing the aggregation of the other lens proteins. The present study provides structural and thermodynamic insights into the ability of human αA-crystallin (HAA) to bind to its partially unfolded clients in the lens, using a small peptide, melittin from bee venom, as a model client. We characterized the thermodynamic parameters of the binding process between melittin and HAA through isothermal titration calorimetry (ITC), and found the binding to be endothermic and entropy-driven. We identified the amino acids in melittin important for binding to HAA by saturation-transfer difference (STD) nuclear magnetic resonance (NMR) experiments, and analysis of NMR line broadening upon titration of melittin with HAA. Our results suggest that hydrophobic residues Ile17 and Ile20 on the C-terminal region of melittin are in close contact with HAA in the melittin-HAA complex. Information obtained from NMR experiments was used to generate structural models of the melittin-HAA complex by molecular docking with high-ambiguity driven docking (HADDOCK). Structural models of the melittin-HAA complex reveal important principles underlying the interaction of HAA with its clients.

Ligand-regulated oligomerisation of allosterically interacting proteins

The binding of ligands to distinct sites at proteins or at protein clusters is often cooperative or anti-cooperative due to allosteric signalling between those sites. The allostery is usually attributed to a configurational change of the proteins from a relaxed to a configurationally different tense state. Alternatively, as originally proposed by Cooper and Dryden, a tense state may be achieved by merely restricting the thermal vibrations of the protein around its mean configuration. In this work, we provide theoretical tools to investigate fluctuation allostery using cooling and titration experiments in which ligands regulate dimerisation, or ring or chain formation. We discuss in detail how ligands may regulate the supramolecular (co)polymerisation of liganded and unliganded proteins.

Allosteric communication pathways and thermal rectification in PDZ-2 protein: a computational study

Allosteric communication in proteins is a fundamental and yet unresolved problem of structural biochemistry. Previous findings, from computational biology ( Ota, N.; Agard, D. A. J. Mol. Biol. 2005 , 351 , 345 - 354 ), have proposed that heat diffuses in a protein through cognate protein allosteric pathways. This work studied heat diffusion in the well-known PDZ-2 protein, and confirmed that this protein has two cognate allosteric pathways and that heat flows preferentially through these. Also, a new property was also observed for protein structures: heat diffuses asymmetrically through the structures. The underling structure of this asymmetrical heat flow was a normal length hydrogen bond (∼2.85 Å) that acted as a thermal rectifier. In contrast, thermal rectification was compromised in short hydrogen bonds (∼2.60 Å), giving rise to symmetrical thermal diffusion. Asymmetrical heat diffusion was due, on a higher scale, to the local, structural organization of residues that, in turn, was also mediated by hydrogen bonds. This asymmetrical/symmetrical energy flow may be relevant for allosteric signal communication directionality in proteins and for the control of heat flow in materials science.

An alternative interpretation of the ultracold methylhydroxycarbene rearrangement mechanism: cooperative effects

Recent studies have reported surprising results related to the rearrangement of carbenes under ultracold conditions, through quantum tunnelling. Here, we demonstrate that a rearrangement of methylhydroxycarbene is possible through cooperative effects.

Update on carbohydrate-containing antibacterial agents

BACKGROUND

Since the first known use of antibiotics > 2,500 years ago, a research field with immense importance for the welfare of mankind has been developed. After a decrease in interest in this topic by the end of the 20th century the occurrence of (poly-)resistant strains of bacteria induced a revival of antibiotics research. Health systems have been seeking viable and reliable solutions to this dangerous and expansive threat.

OBJECTIVE

This review will focus on carbohydrate-containing antibiotics and will give an outline of recently published novel isolated, semisynthetic as well as synthetic structures, their mechanism of action, if known, and the strategies for the design of compounds with potential by improved antibacterial properties.

METHODS

The literature between 2000 and 2008 was screened with main focus on recent examples of novel structures and strategies for the lead finding of exclusively antibacterial agents.

RESULTS/CONCLUSION

With the explanation of the role of the carbohydrate moieties in the respective antibacterial agents together with better synthetic strategies in carbohydrate chemistry as well as improvements in assay development for high throughput screening methods, carbohydrate-containing antibiotics can be used for the finding of potential drug leads that contribute to the fight against infections and diseases caused by (resistant) bacterial pathogens.

The allosteric modulation of thyroxine-binding globulin affinity is entropy driven

BACKGROUND

Thyroxine-binding globulin (TBG) is a non-inhibitory member of the serpin family of proteins whose main structural element is the reactive center loop (RCL), that, upon cleavage by proteases, is inserted into the protein core adopting a β-strand conformation (stressed to relaxed transition, S-to-R). After S-to-R transition thyroxine (T4) affinity decreases. However, crystallographic studies in the presence or absence of the hormone in different states are unable to show significant differences in the structure and interactions of the binding site. Experimental results also suggest the existence of several S states (differing in the number of inserted RCL residues), associated with a differential affinity.

METHODS

To shed light into the molecular basis that regulates T4 affinity according to the degree of RCL insertion in TBG, we performed extended molecular dynamics simulations combined with several thermodynamic analysis of the T4 binding to TBG in three different S states, and in the R state.

RESULTS

Our results show that, despite T4 binding in the protein by similar interactions in all states, a good correlation between the degree of RCL insertion and the binding affinity, driven by a change in TBG conformational entropy, was observed.

CONCLUSION

TBG allosteric regulation is entropy driven. The presence of multiple S states may allow more efficient T4 release due to protease activity.

GENERAL SIGNIFICANCE

The presented results are clear examples of how computer simulation methods can reveal the thermodynamic basis of allosteric effects, and provide a general framework for understanding serpin allosteric affinity regulation.

Vancomycin: ligand recognition, dimerization and super-complex formation

The antibiotic vancomycin targets lipid II, blocking cell wall synthesis in Gram-positive bacteria. Despite extensive study, questions remain regarding how it recognizes its primary ligand and what is the most biologically relevant form of vancomycin. In this study, molecular dynamics simulation techniques have been used to examine the process of ligand binding and dimerization of vancomycin. Starting from one or more vancomycin monomers in solution, together with different peptide ligands derived from lipid II, the simulations predict the structures of the ligated monomeric and dimeric complexes to within 0.1 nm rmsd of the structures determined experimentally. The simulations reproduce the conformation transitions observed by NMR and suggest that proposed differences between the crystal structure and the solution structure are an artifact of the way the NMR data has been interpreted in terms of a structural model. The spontaneous formation of both back-to-back and face-to-face dimers was observed in the simulations. This has allowed a detailed analysis of the origin of the cooperatively between ligand binding and dimerization and suggests that the formation of face-to-face dimers could be functionally significant. The work also highlights the possible role of structural water in stabilizing the vancomycin ligand complex and its role in the manifestation of vancomycin resistance.

Develop and test a solvent accessible surface area-based model in conformational entropy calculations

It is of great interest in modern drug design to accurately calculate the free energies of protein-ligand or nucleic acid-ligand binding. MM-PBSA (molecular mechanics Poisson-Boltzmann surface area) and MM-GBSA (molecular mechanics generalized Born surface area) have gained popularity in this field. For both methods, the conformational entropy, which is usually calculated through normal-mode analysis (NMA), is needed to calculate the absolute binding free energies. Unfortunately, NMA is computationally demanding and becomes a bottleneck of the MM-PB/GBSA-NMA methods. In this work, we have developed a fast approach to estimate the conformational entropy based upon solvent accessible surface area calculations. In our approach, the conformational entropy of a molecule, S, can be obtained by summing up the contributions of all atoms, no matter they are buried or exposed. Each atom has two types of surface areas, solvent accessible surface area (SAS) and buried SAS (BSAS). The two types of surface areas are weighted to estimate the contribution of an atom to S. Atoms having the same atom type share the same weight and a general parameter k is applied to balance the contributions of the two types of surface areas. This entropy model was parametrized using a large set of small molecules for which their conformational entropies were calculated at the B3LYP/6-31G* level taking the solvent effect into account. The weighted solvent accessible surface area (WSAS) model was extensively evaluated in three tests. For convenience, TS values, the product of temperature T and conformational entropy S, were calculated in those tests. T was always set to 298.15 K through the text. First of all, good correlations were achieved between WSAS TS and NMA TS for 44 protein or nucleic acid systems sampled with molecular dynamics simulations (10 snapshots were collected for postentropy calculations): the mean correlation coefficient squares (R²) was 0.56. As to the 20 complexes, the TS changes upon binding; TΔS values were also calculated, and the mean R² was 0.67 between NMA and WSAS. In the second test, TS values were calculated for 12 proteins decoy sets (each set has 31 conformations) generated by the Rosetta software package. Again, good correlations were achieved for all decoy sets: the mean, maximum, and minimum of R² were 0.73, 0.89, and 0.55, respectively. Finally, binding free energies were calculated for 6 protein systems (the numbers of inhibitors range from 4 to 18) using four scoring functions. Compared to the measured binding free energies, the mean R² of the six protein systems were 0.51, 0.47, 0.40, and 0.43 for MM-GBSA-WSAS, MM-GBSA-NMA, MM-PBSA-WSAS, and MM-PBSA-NMA, respectively. The mean rms errors of prediction were 1.19, 1.24, 1.41, 1.29 kcal/mol for the four scoring functions, correspondingly. Therefore, the two scoring functions employing WSAS achieved a comparable prediction performance to that of the scoring functions using NMA. It should be emphasized that no minimization was performed prior to the WSAS calculation in the last test. Although WSAS is not as rigorous as physical models such as quasi-harmonic analysis and thermodynamic integration (TI), it is computationally very efficient as only surface area calculation is involved and no structural minimization is required. Moreover, WSAS has achieved a comparable performance to normal-mode analysis. We expect that this model could find its applications in the fields like high throughput screening (HTS), molecular docking, and rational protein design. In those fields, efficiency is crucial since there are a large number of compounds, docking poses, or protein models to be evaluated. A list of acronyms and abbreviations used in this work is provided for quick reference.

Enthalpy-entropy compensation and cooperativity as thermodynamic epiphenomena of structural flexibility in ligand-receptor interactions

Ligand binding is a thermodynamically cooperative process in many biochemical systems characterized by the conformational flexibility of the reactants. However, the contribution of conformational entropy to cooperativity of ligation needs to be elucidated. Here, we perform kinetic and thermodynamic analyses on a panel of cycle-mutated peptides, derived from influenza H3 HA(306-319), interacting with wild type and a mutant HLA-DR. We observe that, within a certain range of peptide affinity, this system shows isothermal entropy-enthalpy compensation (iEEC). The incremental increases in conformational entropy measured as disruptive mutations are added in the ligand or receptor are more than sufficient in magnitude to account for the experimentally observed lack of free-energy decrease cooperativity. Beyond this affinity range, compensation is not observed, and therefore, the ability of the residual interactions to form a stable complex decreases in an exponential fashion. Taken together, our results indicate that cooperativity and iEEC constitute the thermodynamic epiphenomena of the structural fluctuation that accompanies ligand-receptor complex formation in flexible systems. Therefore, ligand binding affinity prediction needs to consider how each source of binding energy contributes synergistically to the folding and kinetic stability of the complex in a process based on the trade-off between structural tightening and restraint of conformational mobility.

Substrate-modulated thermal fluctuations affect long-range allosteric signaling in protein homodimers: exemplified in CAP

The role of conformational dynamics in allosteric signaling of proteins is increasingly recognized as an important and subtle aspect of this ubiquitous phenomenon. Cooperative binding is commonly observed in proteins with twofold symmetry that bind two identical ligands. We construct a coarse-grained model of an allosteric coupled dimer and show how the signal can be propagated between the distant binding sites via change in slow global vibrational modes alone. We demonstrate that modulation on substrate binding of as few as 5-10 slow modes can give rise to cooperativity observed in biological systems and that the type of cooperativity is given by change of interaction between the two monomers upon ligand binding. To illustrate the application of the model, we apply it to a challenging test case: the catabolite activator protein (CAP). CAP displays negative cooperativity upon association with two identical ligands. The conformation of CAP is not affected by the binding, but its vibrational spectrum undergoes a strong modification. Intriguingly, the first binding enhances thermal fluctuations, yet the second quenches them. We show that this counterintuitive behavior is, in fact, necessary for an optimal anticooperative system, and captured within a well-defined region of the model's parameter space. From analyzing the experimental results, we conclude that fast local modes take an active part in the allostery of CAP, coupled to the more-global slow modes. By including them into the model, we elucidate the role of the modes on different timescales. We conclude that such dynamic control of allostery in homodimers may be a general phenomenon and that our model framework can be used for extended interpretation of thermodynamic parameters in other systems.

Bridging timescales and length scales: from macroscopic flux to the molecular mechanism of antibiotic diffusion through porins

Our aim in this study was to provide an atomic description of ampicillin translocation through OmpF, the major outer membrane channel in Escherichia coli and main entry point for beta-lactam antibiotics. By applying metadynamics simulations, we also obtained the energy barriers along the diffusion pathway. We then studied the effect of mutations that affect the charge and size at the channel constriction zone, and found that in comparison to the wild-type, much lower energy barriers are required for translocation. The expected higher translocation rates were confirmed on the macroscopic scale by liposome-swelling assays. A microscopic view on the millisecond timescale was obtained by analysis of temperature-dependent ion current fluctuations in the presence of ampicillin and provide the enthalpic part of the energy barrier. By studying antibiotic translocation over various timescales and length scales, we were able to discern its molecular mechanism and rate-limiting interactions, and draw biologically relevant conclusions that may help in the design of drugs with enhanced permeation rates.

Analysis of cooperativity by isothermal titration calorimetry

Cooperative binding pervades Nature. This review discusses the use of isothermal titration calorimetry (ITC) in the identification and characterisation of cooperativity in biological interactions. ITC has broad scope in the analysis of cooperativity as it determines binding stiochiometries, affinities and thermodynamic parameters, including enthalpy and entropy in a single experiment. Examples from the literature are used to demonstrate the applicability of ITC in the characterisation of cooperative systems.

Structure of ristocetin A in complex with a bacterial cell-wall mimetic

Antimicrobial drug resistance is a serious public health problem and the development of new antibiotics has become an important priority. Ristocetin A is a class III glycopeptide antibiotic that is used in the diagnosis of von Willebrand disease and which has served as a lead compound for the development of new antimicrobial therapeutics. The 1.0 A resolution crystal structure of the complex between ristocetin A and a bacterial cell-wall peptide has been determined. As is observed for most other glycopeptide antibiotics, it is shown that ristocetin A forms a back-to-back dimer containing concave binding pockets that recognize the cell-wall peptide. A comparison of the structure of ristocetin A with those of class I glycopeptide antibiotics such as vancomycin and balhimycin identifies differences in the details of dimerization and ligand binding. The structure of the ligand-binding site reveals a likely explanation for ristocetin A's unique anticooperativity between dimerization and ligand binding.

Influence of the charge state on the structures and interactions of vancomycin antibiotics with cell-wall analogue peptides: experimental and theoretical studies

Charge matters! The charge state significantly influences the conformation and the binding energy between vancomycin antibiotic and bacterial cell-wall analogue peptides (see figure). Surface-induced dissociation (SID) studies provide a quantitative comparison between the stabilities of different charge states of the complex.In this study we examined the effect of the charge state on the energetics and dynamics of dissociation of the noncovalent complex between the vancomycin and the cell-wall peptide analogue N(alpha),N(epsilon)-diacetyl-L-Lys-D-Ala-D-Ala (V-Ac(2)LKdAdA). The binding energies between the vancomycin and the peptide were obtained from the RRKM (Rice, Ramsperger, Kassel, Marcus) modeling of the time- and energy-resolved surface-induced dissociation (SID) experiments. Our results demonstrate that the stability of the complex towards fragmentation increases in the order: doubly protonated<singly protonated<deprotonated. Dissociation of the singly protonated and singly deprotonated complex is characterized by very large entropy effects, which indicate a substantial increase in the conformational flexibility of the resulting products. The experimental threshold energies of (1.75+/-0.08) eV ((40.3+/-1.8) kcal mol(-1)) and (1.34+/-0.08) eV ((30.9+/-1.8) kcal mol(-1)) obtained for the deprotonated and singly protonated complexes, respectively, are in excellent agreement with the results of density functional theory calculations. The increased stability of the deprotonated complex observed experimentally is attributed to the presence of three charged sites in the deprotonated complex, as compared with only one charged site in the singly protonated complex. The low binding energy of (0.93+/-0.04) eV ((21.4+/-0.9) kcal mol(-1)) obtained for the doubly protonated complex suggests that this ion is destabilized by Coulomb repulsion between the singly protonated vancomycin and the singly protonated peptide comprising the complex.

Allostery: absence of a change in shape does not imply that allostery is not at play

Allostery is essential for controlled catalysis, signal transmission, receptor trafficking, turning genes on and off, and apoptosis. It governs the organism's response to environmental and metabolic cues, dictating transient partner interactions in the cellular network. Textbooks taught us that allostery is a change of shape at one site on the protein surface brought about by ligand binding to another. For several years, it has been broadly accepted that the change of shape is not induced; rather, it is observed simply because a larger protein population presents it. Current data indicate that while side chains can reorient and rewire, allostery may not even involve a change of (backbone) shape. Assuming that the enthalpy change does not reverse the free-energy change due to the change in entropy, entropy is mainly responsible for binding.

A computational investigation of allostery in the catabolite activator protein

The catabolite activator protein is a dimer that consists of two cAMP-binding subunits, each containing a C-terminus DNA-binding module and a N-terminus ligand binding domain. The system is well-known to exhibit negative cooperativity, whereby the binding of one cAMP molecule reduces the binding affinity of the other cAMP molecule by 2 orders of magnitude, despite the large separation between the cAMP binding pockets. Here we use extensive explicit-solvent molecular dynamics simulations (135 ns) to investigate the allosteric mechanism of CAP. Six trajectories were carried out for apo, singly liganded, and doubly liganded CAP, both in the presence and absence of DNA. Thorough analyses of the dynamics through the construction of dynamical cross-correlated maps, as well as essential dynamics analyses, indicated that the system experienced a switch in motion as a result of cAMP binding, in accordance with recent NMR experiments carried out on a truncated form of the protein. Analyses of conformer structures collected from the simulations revealed a remarkable event: the DNA-binding module was found to dissociate from the N-terminus ligand binding domain. An interesting aspect of this structural change is that it only occurred in unoccupied subunits, suggesting that the binding of cAMP provides additional stability to the system, consistent with the increase in entropy that was observed in our calculations and from isothermal titration calorimetry. Analysis of the distribution of intrinsic disorder propensities in CAP amino acid sequence using PONDR VLXT and VSL1 predictors revealed that the region connecting ligand-binding and DNA-binding domains of CAP have the potential to exhibit increased flexibility. We complemented these trajectories with free energy calculations following the MM-PBSA approach on more than 2000 snapshots that included 880 normal mode analysis. The resulting free energy differences between the singly liganded and doubly liganded states were in excellent agreement with isothermal titration calorimetry data. When the free energy calculations were carried out in the presence of DNA, we discovered that a switch in cooperativity occurred, so that the binding of the first cAMP promoted the binding of the other cAMP. The components of the free energy reveal that this effect is mainly entropic in nature, whereby the DNA reduces the degree of tightening that is observed in its absence, thereby promoting binding of the second cAMP. This finding prompted us to propose a new mechanism by which CAP triggers the transcription activation that is based on an order to disorder transition mediated by cAMP binding as well as DNA.

Nonequilibrium dynamic mechanism for allosteric effect

Allosteric regulation is often viewed as thermodynamic in nature. However, protein internal motions during an enzymatic reaction cycle can slow the hoping processes over numerous potential barriers. We propose that regulating molecules may function by modifying the nonequilibrium protein dynamics. The theory predicts that an enzyme under the new mechanism has a different temperature dependence, waiting time distribution of the turnover cycle, and dynamic fluctuation patterns with and without an effector. Experimental tests of the theory are proposed.

Molecular evolution of affinity and flexibility in the immune system

The immune system responds to the introduction of foreign antigens by rapidly evolving antibodies with increasing affinity for the antigen (i.e., maturation). To investigate the factors that control this process at the molecular level, we have assessed the changes in flexibility that accompany ligand binding at four stages of maturation in the 4-4-20 antibody. Our studies, based on molecular dynamics, indicate that increased affinity for the target ligand is associated with a decreased entropic cost to binding. The entropy of binding is unfavorable, opposing favorable enthalpic contributions that arise during complex formation. Computed binding free energies for the various antibody-ligand complexes qualitatively reproduce the trends observed in the experimentally derived values, although the absolute magnitude of free-energy differences is overestimated. Our results support the existence of a correlation between high-affinity interactions and decreased protein flexibility in this series of antibody molecules. This observation is likely to be a general feature of molecular association processes and key to the molecular evolution of antibody responses.

Guanylate kinase, induced fit, and the allosteric spring probe

Since the introduction of the induced-fit theory by D. E. Koshland Jr., it has been established that conformational motion invariably accompanies the execution of protein function. The catalytic activity of kinases, specifically, is associated with large conformational changes ( approximately 1 nm amplitude). In the case of guanylate kinase, upon substrate binding, the LID and nucleotide-monophosphate-binding domains are brought together and toward the CORE with large concerted movements about the alpha3 (helix 3) axis. However, whether the change in conformation mostly affects the catalytic rate or mostly increases binding affinities for one or the other substrate is unclear. We investigate this question using a nanotechnology approach based on mechanical stress. Using an "allosteric spring probe", we bias conformational states in favor of the "open" (substrate-free) conformation of the enzyme; the result is that the binding constant for the substrate guanosine monophosphate (GMP) is reduced by up to a factor of 10, whereas the binding constant for adenosine triphosphate (ATP) and the catalytic rate are essentially unaffected. The results show that the GMP-induced conformational change, which promotes catalysis, does not promote ATP binding, consistent with previous mutagenesis studies. Furthermore, they show that this conformational change is of the induced-fit type with respect to GMP binding (but not ATP binding). We elaborate on this point by proposing a quantitative criterion for the classification of conformational changes with respect to the induced-fit theory. More generally, these results show that the allosteric spring probe can be used to affect enzymatic activity in a continuously controlled manner, and also to affect specific steps of the reaction mechanism while leaving others unaffected. It is presumed that this will enable informative comparisons with the results of future molecular dynamics or statistical mechanics computations.

Chemical double-mutant cycles: dissecting non-covalent interactions

Thermodynamic double-mutant cycles and triple-mutant boxes are widely employed for the experimental quantification of non-covalent interactions and cooperative effects in proteins. This review describes the application of these powerful methodologies to the study of non-covalent interactions in synthetic systems.

Dynamic allostery of protein alpha helical coiled-coils

Alpha helical coiled-coils appear in many important allosteric proteins such as the dynein molecular motor and bacteria chemotaxis transmembrane receptors. As a mechanism for transmitting the information of ligand binding to a distant site across an allosteric protein, an alternative to conformational change in the mean static structure is an induced change in the pattern of the internal dynamics of the protein. We explore how ligand binding may change the intramolecular vibrational free energy of a coiled-coil, using parameterized coarse-grained models, treating the case of dynein in detail. The models predict that coupling of slide, bend and twist modes of the coiled-coil transmits an allosteric free energy of approximately 2kBT, consistent with experimental results. A further prediction is a quantitative increase in the effective stiffness of the coiled-coil without any change in inherent flexibility of the individual helices. The model provides a possible and experimentally testable mechanism for transmission of information through the alpha helical coiled-coil of dynein.

Coupling of global and local vibrational modes in dynamic allostery of proteins

It is now recognized that internal global protein dynamics play an important role in the allosteric function of many proteins. Alterations of protein flexibility on effector binding affect the entropic cost of binding at a distant site. We present a coarse-grained model for a potential amplification of such entropic allostery due to coupling of fast, localized modes to the slow, global modes. We show how such coupling can give rise to large compensating entropic and enthalpic terms. The model corresponds to the pattern of calorimetry and NMR data from experiments on the Met repressor.

Cooperativity in the self-assembly of porphyrin ladders

Cooperativity is a general feature of intermolecular interactions in biomolecular systems, but there are many different facets of the phenomenon that are not well understood. Positive cooperativity stabilizes a system as progressively more interactions are added, and the origin of the beneficial free energy may be entropic or enthalpic in origin. An "enthalpic chelate effect" has been proposed to operate through structural tightening that improves all of the functional group interactions in a complex, when it is more strongly bound. Here, we present direct calorimetric evidence that no such enthalpic effects exist in the cooperative assembly of supramolecular ladder complexes composed of metalloporphyrin oligomers coordinated to bipyridine ligands. The enthalpic contributions of the individual coordination interactions are 35 kJ.mol(-1) and constant over a range of free energies of self-assembly of -35 to -111 kJ.mol(-1). In rigid well defined systems of this type, the enthalpies of individual interactions are additive, and no enthalpic cooperative effects are apparent. The implication is that in more flexible, less well defined systems such as biomolecular assemblies, the enthalpy contributions available from specific functional group interactions are well defined and constant parameters.

Conformational and quantitative characterization of oritavancin-peptidoglycan complexes in whole cells of Staphylococcus aureus by in vivo 13C and 15N labeling

Solid-state NMR has been used to examine the cell walls of intact whole cells of Staphyloccus aureus grown on media containing D-[1-(13)C]alanine, [(15)N]glycine, and the alanine racemase inhibitor, alaphosphin. The results of in situ site-selective, four-frequency NMR experiments show directly for the first time that (i) 54% of the cell-wall peptidoglycan stems have D-alanine termini and 46%, D-alanine-D-alanine termini; (ii) the molar ratio of stems ending in D-alanine to esterified alditol repeats of cell-wall teichoic and lipoteichoic acids is 3:2; and (iii) 50% of the mature cell-wall binding sites for a fluorinated oritavancin analogue consist of two nearest-neighbor peptide stems of different glycan strands. The drug is bound to the D-Ala-D-Ala terminus of one stem and is proximate to the bridging pentaglycyl segment that cross-links the two stems. Structural details of the binding site are revealed in a model of the glycopeptide-peptidoglycan interaction produced by molecular dynamics simulations with internuclear distance restraints determined by NMR.

Evaluating the conformational entropy of macromolecules using an energy decomposition approach

We have developed a novel method to compute the conformational entropy of any molecular system via conventional simulation techniques. This method only requires that the total energy of the system is available and that the Hamiltonian is separable, with individual energy terms for the various degrees of freedom. Consequently the method, which we call the energy decomposition (Edcp) approach, is general and applicable to any large polymer in implicit solvent. Edcp is applied to estimate the entropy differences due to the peptide and ester groups in polyalanine and polyalanil ester. Ensembles over a wide range of temperatures were generated by replica exchange molecular dynamics, and densities of states were estimated using the weighted histogram analysis method. The results are compared with those obtained via evaluating the P ln P integral or employing the quasiharmonic approximation, other approaches widely employed to evaluate the entropy of molecular systems. Unlike the former method, Edcp can accommodate the correlations present between separate degrees of freedom. In addition, the Edcp model assumes no specific form for the underlying fluctuations present in the system, in contrast to the quasiharmonic approximation. For the molecules studied, the quasiharmonic approximation is observed to produce a good estimate of the vibrational entropy, but not of the conformational entropy. In contrast, our energy decomposition approach generates reasonable estimates for both of these entropy terms. We suggest that this approach embodies a simple yet effective solution to the problem of evaluating the conformational entropy of large macromolecules in implicit solvent.

Allosteric control through mechanical tension

Conformational changes of proteins modulate their function. In allosteric control, the conformational change is induced by the binding of a signaling molecule. Here we insert a "molecular spring" on the enzyme guanylate kinase, to control the conformation of this protein. The stiffness of the spring can be varied externally, which allows one to exert a controlled mechanical tension between the two points on the protein's surface where the spring is attached. We show that by applying and releasing the tension we can reversibly turn the enzyme off and on.

Affinity Capture Using Vancomycin-Bound Magnetic Nanoparticles for the MALDI-MS Analysis of Bacteria

Matrix-assisted laser desorption/ionization mass spectrometry (MALDI-MS) provides a straightforward means to differentiate microorganism species based on mass spectral fingerprinting. The pathogen cell concentration in an infected sample, however, is generally lower than that capable of being detected directly by MALDI-MS. Furthermore, the presence of proteins or metabolites in biological fluids always causes unavoidable interference for the identification of microorganism species. Vancomycin, which binds to D-Ala-D-Ala moieties on the cell walls of Gram-positive bacteria and, therefore, inhibits peptidoglycan synthesis, is one of the most potent antibiotics. Thus, we have employed vancomycin-modified magnetic nanoparticles as affinity probes to selectively trap Gram-positive pathogens from sample solutions; i.e., these bacteria can be isolated from sample solutions by applying a magnetic field. The isolated cells could then be characterized by MALDI-MS. This approach effectively reduces the interference of protein and metabolite signals in the mass spectra of Gram-positive bacteria because vancomycin has such high specificity for the D-Ala-D-Ala units of the cell walls. The lowest cell concentration we detected for both Staphylococcus saprophyticus and Staphylococcus aureus in a urine sample (3 mL) was approximately 7 x 10(4) cfu/mL.

Coarse-grained model of entropic allostery

Many signaling functions in molecular biology require proteins to bind to substrates such as DNA in response to environmental signals such as the simultaneous binding to a small molecule. Examples are repressor proteins which may transmit information via a conformational change in response to the ligand binding. An alternative entropic mechanism of "allostery" suggests that the inducer ligand changes the intramolecular vibrational entropy, not just the mean static structure. We present a quantitative, coarse-grained model of entropic allostery, which suggests design rules for internal cohesive potentials in proteins employing this effect. It also addresses the issue of how the signal information to bind or unbind is transmitted through the protein. The model may be applicable to a wide range of repressors and also to signaling in trans-membrane proteins.

The role of dynamics in allosteric regulation

The biomolecular conformational changes often associated with allostery are, by definition, dynamic processes. Recent publications have disclosed the role of pre-existing equilibria of conformational substates in this process. In addition, the role of dynamics as an entropic carrier of free energy of allostery has been investigated. Recent work thus shows that dynamics is pivotal to allostery, and that it constitutes much more than just the move from the 'T'-state to the 'R'-state. Emerging computational studies have described the actual pathways of allosteric change.

Free-Energy Calculations of Protein−Ligand Cation−π and Amino−π Interactions: From Vacuum to Proteinlike Environments

To probe the role of cation-pi and amino-pi interactions in the context of protein-ligand interactions, the stability of 55 X-ray cation/amino-pi motifs involving the Ade moieties of cofactor molecules and Arg, Lys, Asn, or Gln side chains of their host protein was evaluated using quantum chemistry calculations. The conjunction of vacuum interaction energies, vibrational entropy, and solvation contributions led to identify Arg-Ade as the most favorable cation/amino-pi complex in the solvents considered, followed by Asn/Gln-Ade and Lys-Ade: their minimum interaction free energies are approximately equal to -7, -4, and -2 kcal/mol, respectively, in the solvents of dielectric constant similar to that estimated for proteins (i.e., acetone, THF, and CCl(4)). Remarkably, these free-energy values of cation/amino-pi interactions correlate well with their frequency of occurrences in protein-ligand structures, which corroborates our approach in the absence of experimental data.