Theoretical-Studies on the Modes of Binding of Some of the Beta-Glycosidically Linked Glucobioses to Concanavalin-A

The probable modes of binding for methyl-α-d-sophoroside, methyl-β-d-sophoroside, laminariboise and cellobiose to concanavalin A have been determined using theoretical methods. Methyl-d-sophorosides can bind to concanavalin A in two modes, i.e. by placing their reducing as well as non-reducing sugar units in the carbohydrate specific binding site, whereas laminaribiose and cellobiose can reach the binding site only with their non-reducing glucose units. However, the probability for methyl-α-d-sophoroside to bind to concanavalin A with its reducing sugar residue as the occupant of the binding site is much higher than it is with its non-reducing sugar residue as the occupant of the sugar binding site. A few of the probable conformers of methyl-β-d-sophoroside can bind to concanavalin A with either the reducing or non-reducing glucose unit. Higher energy conformers of cellobiose or laminaribiose can reach the binding site with their non-reducing residues alone. The relative differences in the binding affinities of these disaccharides are mainly due to the differences in the availability of proper conformers which can reach the binding site and to non-covalent interactions between the sugar and the protein. This study also suggests that though the sugar binding site of concanavalin A accommodates a single sugar residue, the residue outwards from the binding site also interacts with concanavalin A, indicating the existence of extended concanavalin A carbohydrate interactions.

Effect of Crowding Agents, Signal Peptide, and Chaperone SecB on the Folding and Aggregation of E. coli Maltose Binding Protein

SecB, a soluble cytosolic chaperone component of the Secexport pathway, binds to newly synthesized precursor proteins and prevents their premature aggregation and folding and subsequently targets them to the translocation machinery on the membrane. PreMBP, the precursor form of maltose binding protein, has a 26-residue signal sequence attached to the N-terminus of MBP and is a physiological substrate of SecB. We examine the effect of macromolecular crowding and SecB on the stability and refolding of denatured preMBP and MBP. PreMBP was less stable than MBP (ΔTm =7( 0.5 K) in both crowded and uncrowded solutions. Crowding did not cause any substantial changes in the thermal stability ofMBP(ΔTm=1(0.4 K) or preMBP (ΔTm=0(0.6 K), as observed in spectroscopically monitored thermal unfolding experiments. However, both MBP and preMBP were prone to aggregation while refolding under crowded conditions. In contrast to MBP aggregates, which were amorphous, preMBP aggregates form amyloid fibrils.Under uncrowded conditions, a molar excess of SecB was able to completely prevent aggregation and promote disaggregation of preformed aggregates of MBP. When a complex of the denatured protein and SecB was preformed, SecB could completely prevent aggregation and promote folding of MBP and preMBP even in crowded solution. Thus, in addition to maintaining substrates in an unfolded, export-competent conformation, SecB also suppresses the aggregation of its substrates in the crowded intracellular environment. SecB is also able to promote passive disaggregation of macroscopic aggregates of MBP in the absence of an energy source such as ATP or additional cofactors. These experiments also demonstrate that signal peptide can reatly influence protein stability and aggregation propensity.

Effect of configuration of the inhibitors on the mode of binding to the enzyme, thermolysin

Preferred conformations of the competitive inhibitors glycyl-L-phenylalanine and glycyl-D-phenylalanine and their mode of binding to thermolysin have been studied. The difference in configuration is shown to affect significantly the mode of binding to thermolysin. Gly-D-Phe prefers to enter the active site in the global minimum conformation whereas Gly-L-Phe may enter in a higher energy conformation. Moreover, D-enantiomer is shown to have a better fit than the L-counterpart in the active site.

Complex carbohydrates: 2. The modes of binding of complex carbohydrates to concanavalin A - a computer modelling approach

The probable modes of binding of some complex carbohydrates, which have the trimannosidic core structure (Man3GlcNAc2), to concanavalin A (Con A) have been determined using a computer modelling technique. These studies show that Con a can bind to the terminal mannose residues of the trimannosidic core structure and to the internal mannosyl as well as to the terminal N-acetylglucosamine residues of the N-acetylglucosamine substituted trimannosidic core structure. The oligosaccharide with terminal mannose residues can bind in its minimum energy conformers, whereas the oligosaccharide with internal mannosyl and terminal N-acetylglucosamine residues can bind only in higher energy conformers. In addition the former oligosaccharide forms more hydrogen bonds with Con A than the latter. These results suggest that, for these oligosaccharides, the terminal mannose residue has a much higher probability of reaching the binding site than either the internal mannosyl or the terminal N-acetylglucosamine residues. The substitution of a bisecting N-acetylglucosamine residue on these oligosaccharides, affects significantly the accessibility of the residues which bind to Con A and thereby reduces their binding affinity. It thus seems that the binding affinity of an oligosaccharide to Con A depends not only on the number of sugar residues which possess free 3-, 4- and 6-hydroxyl groups but also on the accessibility of these sugar residues to Con A. This study also reveals that the sugar binding site of Con A is small and that the interactions between Con A and carbohydrates are extended slightly beyond the single sugar residue that is placed in the binding site.

Theoretical studies on penicillins, penicillin sulphones and their 1-oxa-1-dethia and 1-carba-1-dethia analogues: binding specificities of transeptidases and penicillinases

Empirical potential energy calculations have been carried out to determine the preferred conformations of penicillins and penicillin sulphones and their 1-oxa-1-dethia and 1-carba-1-dethia analogues. With the exception of 1-oxa-1-dethia penicillins, all the other compounds favour C2 and the C3 puckered conformations of their five-membered rings. Replacement of C2 methyl groups by hydrogen atoms as in bisnorpenicillin V or oxidation of sulphur in position 1 as in sulphones, makes the C3 puckered form much less favourable. Addition of an amino-acyl group at the C6 atom, however, makes the C3 puckered form more favoured in penicillin G or V and in 1-carba-1-dethia penicillins. Through the replacement of the sulphur atom at position 1 by an oxygen atom or by a -CH2 group increases the non-planarity of the lactam peptide bond, it significantly affects the relative disposition of the C3 carboxyl group with respect to the β-lactam ring. These conformational differences have been correlated with the biological activities of these compounds. The present study suggests that the conformation of the bicyclic ring system may be more important for initial binding with the crosslinking enzyme(s) involved in the biosynthesis of bacterial cell-wall peptidoglycan and that the mode of binding is influenced by the nature of the side-group at the C6 atom. These studies predict, in agreement with experimental results, that the 1-oxa-1-dethia penicillin nulceus is an inhibitor of penicillianses. The study also suggests that the stereospecificities of the crosslinking enzyme(s) and penicillinases are very similar with regard to the nature of the side-group at the 6 atom and the confirmation of the bicyclic ring system. However, the confirmational requirement for the bicyclic ring system appears to be more specific in the former enzyme than in the latter.

New ternary copper(II) complexes of L-alanine and heterocyclic bases: DNA binding and oxidative DNA cleavage activity

Four new ternary copper(II) complexes of alpha-amino acid having polypyridyl bases of general formulation [Cu(L-ala)(B)(H2O)](X)(1-4), where L-ala is L-alanine, B is an N,N-donor heterocyclic base, viz. 2,2'-bipyridine (bpy, 1), 1,10-phenanthroline (phen, 2) and 5,6-phenanthroline dione (dione, 3), dipyrido[3,2:2',3'-f] quinoxaline (dpq, 4), and X = ClO4-/NO3- are synthesized, characterized by various spectroscopic and X-ray crystallographic methods. The complexes show a distorted square-pyramidal (4 + 1) CuN3O2 coordination geometry. The one-electron paramagnetic complexes (1-4) display a low energy d-d band near 600 nm in aqueous medium and show a quasi-reversible cyclic voltammetric response due to one-electron Cu(II)/Cu(I) reduction near - 100 mV (versus SCE) in DMF-0.1 M TBAP. Binding interactions of the complexes with calf thymus DNA (CT-DNA) were investigated by UV-Vis absorption titration, ethidium bromide displacement assay, viscometric titration experiment and DNA melting studies. All the complexes barring the complexes 1 and 3 are avid binder to the CT-DNA in the DNA minor groove giving an order: 4 > 2 >>>1, 3. The complexes 2 and 4 show appreciable chemical nuclease activity in the presence of 3-mercaptopropionic acid (MPA) as a reducing agent. Hydroxyl radical was investigated to be the DNA cleavage active species. Control experiments in the presence of distamycin-A show primarily minor groove-binding propensity for the complexes 2 and 4 to the DNA.

Free energy calculations of glycosaminoglycan - Protein interactions

Glycosaminoglycans (GAGs) are complex highly charged linear polysaccharides that have a variety of roles in biological processes. We report the first use of molecular dynamics (MD) free energy calculations using the MM/PBSA method to investigate the binding of GAGs to protein molecules, namely the platelet endothelial cell adhesion molecule 1 (PECAM-1) and annexin A2. Calculations of the free energy of the binding of heparin fragments of different sizes reveal the existence of a region of low GAG-binding affinity in domains 5-6 of PECAM-1 and a region of high affinity in domains 2-3, consistent with experimental data and ligand-protein docking studies. A conformational hinge movement between domains 2 and 3 was observed, which allows the binding of heparin fragments of increasing size (pentasaccharides to octasaccharides) with an increasingly higher binding affinity. Similar simulations of the binding of a heparin fragment to annexin A2 reveal the optimization of electrostatic and hydrogen bonding interactions with the protein and protein-bound calcium ions. In general, these free energy calculations reveal that the binding of heparin to protein surfaces is dominated by strong electrostatic interactions for longer fragments, with equally important contributions from van der Waals interactions and vibrational entropy changes, against a large unfavorable desolvation penalty due to the high charge density of these molecules.

Free energy calculations of the interactions of c-Jun-based synthetic peptides with the c-Fos protein

The c-Fos–c-Jun complex forms the activator protein 1 transcription factor, a therapeutic target in the treatment of cancer. Various synthetic peptides have been designed to try to selectively disrupt the interaction between c-Fos and c-Jun at its leucine zipper domain. To evaluate the binding affinity between these synthetic peptides and c-Fos, polarizable and nonpolarizable molecular dynamics (MD) simulations were conducted, and the resulting conformations were analyzed using the molecular mechanics generalized Born surface area (MM/GBSA) method to compute free energies of binding. In contrast to empirical and semiempirical approaches, the estimation of free energies of binding using a combination of MD simulations and the MM/GBSA approach takes into account dynamical properties such as conformational changes, as well as solvation effects and hydrophobic and hydrophilic interactions. The predicted binding affinities of the series of c-Jun-based peptides targeting the c-Fos peptide show good correlation with experimental melting temperatures. This provides the basis for the rational design of peptides based on internal, van der Waals, and electrostatic interactions.

Calculations of the free energy of interaction of the c-Fos−c-Jun coiled coil: Effects of the solvation model and the inclusion of polarization effects

The leucine zipper region of activator protein-1 (AP-1) comprises the c-Jun and c-Fos proteins and constitutes a well-known coiled coil protein−protein interaction motif. We have used molecular dynamics (MD) simulations in conjunction with the molecular mechanics/Poisson−Boltzmann generalized-Born surface area [MM/PB(GB)SA] methods to predict the free energy of interaction of these proteins. In particular, the influence of the choice of solvation model, protein force field, and water potential on the stability and dynamic properties of the c-Fos−c-Jun complex were investigated. Use of the AMBER polarizable force field ff02 in combination with the polarizable POL3 water potential was found to result in increased stability of the c-Fos−c-Jun complex. MM/PB(GB)SA calculations revealed that MD simulations using the POL3 water potential give the lowest predicted free energies of interaction compared to other nonpolarizable water potentials. In addition, the calculated absolute free energy of binding was predicted to be closest to the experimental value using the MM/GBSA method with independent MD simulation trajectories using the POL3 water potential and the polarizable ff02 force field, while all other binding affinities were overestimated.

The GAF Domain of the cGMP-Binding, cGMP-Specific Phosphodiesterase (PDE5) Is a Sensor and a Sink for cGMP

We describe here a novel sensor for cGMP based on the GAF domain of the cGMP-binding, cGMP-specific phosphodiesterase 5 (PDE5) using bioluminescence resonance energy transfer (BRET). The wild type GAFa domain, capable of binding cGMP with high affinity, and a mutant (GAFaF163A) unable to bind cGMP were cloned as fusions between GFP and Rluc for BRET2 assays. BRET2 ratios of the wild type GAFa fusion protein, but not GAFaF163A, increased in the presence of cGMP but not cAMP. Higher basal BRET2 ratios were observed in cells expressing the wild type GAFa domain than in cells expressing GAFaF163A. This was correlated with elevated basal intracellular levels of cGMP, indicating that the GAF domain could act as a sink for cGMP. The tandem GAF domains in full length PDE5 could also sequester cGMP when the catalytic activity of PDE5 was inhibited. Therefore, these results describe a cGMP sensor utilizing BRET2 technology and experimentally demonstrate the reservoir of cGMP that can be present in cells that express cGMP-binding GAF domain-containing proteins. PDE5 is the target for the anti-impotence drug sildenafil citrate; therefore, this GAF-BRET2 sensor could be used for the identification of novel compounds that inhibit cGMP binding to the GAF domain, thereby regulating PDE5 catalytic activity.

Excitation energy transfer from a fluorophore to single-walled carbon nanotubes

We study the process of electronic excitation energy transfer from a fluorophore to the electronic energy levels of a single-walled carbon nanotube. The matrix element for the energy transfer involves the Coulombic interaction between the transition densities on the donor and the acceptor. In the Foumlrster approach, this is approximated as the interaction between the corresponding transition dipoles. For energy transfer from a dye to a nanotube, one can use the dipole approximation for the dye, but not for the nanotube. We have therefore calculated the rate using an approach that avoids the dipole approximation for the nanotube. We find that for the metallic nanotubes, the rate has an exponential dependence if the energy that is to be transferred, h is less than a threshold and a d(-5) dependence otherwise. The threshold is the minimum energy required for a transition other than the k(i,perpendicular to)=0 and l=0 transition. Our numerical evaluation of the rate of energy transfer from the dye pyrene to a (5,5) carbon nanotube, which is metallic leads to a distance of similar to 165 A degrees up to which energy transfer is appreciable. For the case of transfer to semiconducting carbon nanotubes, apart from the process of transfer to the electronic energy levels within the one electron picture, we also consider the possibility of energy transfer to the lowest possible excitonic state. Transfer to semiconducting carbon nanotubes is possible only if>=epsilon(g)-epsilon(b). The long range behavior of the rate of transfer has been found to have a d(-5) dependence if h >=epsilon(g). But, when the emission energy of the fluorophore is in the range epsilon(g)>h >=epsilon(g)-epsilon(b), the rate has an exponential dependence on the distance. For the case of transfer from pyrene to the semiconducting (6,4) carbon nanotube, energy transfer is found to be appreciable up to a distance of similar to 175 A degrees.

Thermodynamic and kinetic studies on the mechanism of binding of methylumbelliferyl galactosides to the basic agglutinin from winged bean (Psophocarpus tetragonolobus)

The binding of winged bean basic agglutinin (WBA I) to 4-methylumbelliferyl (MeUmb) galactosides was examined by extrinsic fluorescence titration and stopped-flow spectrofluorimetry. Upon binding to WBA I, MeUmb alpha-galactosides show quenching in fluorescence intensity, decrease in UV absorbance with a concomitant blue shift, and decrease in fluorescence excited-state lifetimes. However, their beta-analogues show enhancement in fluorescence intensity, increase in UV absorbance with a red shift, and an increase in fluorescence excited-state lifetimes. This implies that the umbelliferyl groups of alpha- and beta-galactosides experience non-polar and polar microenvironments, respectively, upon binding to WBA I. Replacement of the anomeric hydroxyl group of galactose by 4-methylumbelliferyl moiety increases the affinity of resulting saccharides. Substitution of C-2 hydroxyl of galactose by an acetamido group leads to increased affinity due to a favorable entropy change. This suggests that acetamido group of MeUmb-alpha/beta-GalNAc binds to a relatively non-polar subsite of WBA I. Most interestingly, this substitution also reduces the association rate constants dramatically. Inspection of the activation parameters reveals that the enthalpy of activation is the limiting factor for the differences in the forward rate constants for these saccharides and the entropic contribution to the activation energy is small

Computer modelling approach to study the modes of binding of alpha- and beta-anomers of D-galactose, D-fucose and D-glucose to L-arabinose-binding protein

The modes of binding of alpha- and beta-anomers of D-galactose, D-fucose and D-glucose to L-arabinose-binding protein (ABP) have been studied by energy minimization using the low resolution (2.4 A) X-ray data of the protein. These studies suggest that these sugars preferentially bind in the alpha-form to ABP, unlike L-arabinose where both alpha- and beta-anomers bind almost equally. The best modes of binding of alpha- and beta-anomers of D-galactose and D-fucose differ slightly in the nature of the possible hydrogen bonds with the protein. The residues Arg 151 and Asn 232 of ABP from bidentate hydrogen bonds with both L-arabinose and D-galactose, but not with D-fucose or D-glucose. However in the case of L-arabinose, Arg 151 forms hydrogen bonds with the hydroxyl group at the C-4 atom and the ring oxygen, whereas in case of D-galactose it forms bonds with the hydroxyl groups at the C-4 and C-6 atoms of the pyranose ring. The calculated conformational energies also predict that D-galactose is a better inhibitor than D-fucose and D-glucose, in agreement with kinetic studies. The weak inhibitor D-glucose binds preferentially to one domain of ABP leading to the formation of a weaker complex. Thus these studies provide information about the most probable binding modes of these sugars and also provide a theoretical explanation for the observed differences in their binding affinities.

Computer simulation of protein—carbohydrate complexes: application to arabinose-binding protein and pea lectin

The CCEM method (Contact Criteria and Energy Minimisation) has been developed and applied to study protein-carbohydrate interactions. The method uses available X-ray data even on the native protein at low resolution (above 2.4 Å) to generate realistic models of a variety of proteins with various ligands.The two examples discussed in this paper are arabinose-binding protein (ABP) and pea lectin. The X-ray crystal structure data reported on ABP-β-l-arabinose complex at 2.8, 2.4 and 1.7 Å resolution differ drastically in predicting the nature of the interactions between the protein and ligand. It is shown that, using the data at 2.4 Å resolution, the CCEM method generates complexes which are as good as the higher (1.7 Å) resolution data. The CCEM method predicts some of the important hydrogen bonds between the ligand and the protein which are missing in the interpretation of the X-ray data at 2.4 Å resolution. The theoretically predicted hydrogen bonds are in good agreement with those reported at 1.7 Å resolution. Pea lectin has been solved only in the native form at 3 Å resolution. Application of the CCEM method also enables us to generate complexes of pea lectin with methyl-α-d-glucopyranoside and methyl-2,3-dimethyl-α-d-glucopyranoside which explain well the available experimental data in solution.

Elucidation of the conformational free energy landscape in H.pylori LuxS and its implications to catalysis

Background: One of the major challenges in understanding enzyme catalysis is to identify the different conformations and their populations at detailed molecular level in response to ligand binding/environment. A detail description of the ligand induced conformational changes provides meaningful insights into the mechanism of action of enzymes and thus its function. Results: In this study, we have explored the ligand induced conformational changes in H. pylori LuxS and the associated mechanistic features. LuxS, a dimeric protein, produces the precursor (4,5-dihydroxy-2,3-pentanedione) for autoinducer-2 production which is a signalling molecule for bacterial quorum sensing. We have performed molecular dynamics simulations on H. pylori LuxS in its various ligand bound forms and analyzed the simulation trajectories using various techniques including the structure network analysis, free energy evaluation and water dynamics at the active site. The results bring out the mechanistic details such as co operativity and asymmetry between the two subunits, subtle changes in the conformation as a response to the binding of active and inactive forms of ligands and the population distribution of different conformations in equilibrium. These investigations have enabled us to probe the free energy landscape and identify the corresponding conformations in terms of network parameters. In addition, we have also elucidated the variations in the dynamics of water co-ordination to the Zn2+ ion in LuxS and its relation to the rigidity at the active sites. Conclusions: In this article, we provide details of a novel method for the identification of conformational changes in the different ligand bound states of the protein, evaluation of ligand-induced free energy changes and the biological relevance of our results in the context of LuxS structure-function. The methodology outlined here is highly generalized to illuminate the linkage between structure and function in any protein of known structure.