Monosialogangliosides and Their Interaction with Cholera Toxin—Investigation by Molecular Modeling and Molecular Mechanics

Abstract

Molecular mechanics and molecular dynamics studies are performed to investigate the conformational preference of cell surface monosialogangliosides (GM3, GM2 and GM1) in aqueous environment. Water mediated hydrogen bonding network plays a significant role in the structural stabilization of GM3, GM2 and GM1. The spatial flexibility of NeuNAc of gangliosides at the binding site of cholera toxin reveals a limited allowed eulerian space of 2.4% with a much less allowed eulerian space (1.4%) for external galactose of GM1. The molecular mechanics of monosialoganglioside-cholera toxin complex reveals that cholera toxin can accommodate the monosialogangliosides in three different modes. Direct and water mediated hydrogen bonding interactions stabilize these binding modes and play an essential role in defining the order of specificity for different monosialogangliosides towards cholera toxin. This study identifies the NeuNAc binding site as a site for design of inhibitors that can restrict the pathogenic activity of cholera toxin.     

Disialogangliosides and their interaction with cholera toxin - investigation by molecular modeling, molecular mechanics and molecular dynamics

Molecular mechanics and molecular dynamics studies are performed to investigate the conformational preference of cell surface disialogangliosides (GD1A, GD1B and GD3) in aqueous environment. The molecular mechanics calculation reveals that water mediated hydrogen bonding network plays a significant role in the structural stabilization of GD1A, GD1B and GD3. These water mediated hydrogen bonds not only exist between neighboring residues but also exist between residues that are separated by 2 to 3 residues in between. The conformational energy difference between different conformational states of gangliosides correlates very well with the number of water mediated and direct hydrogen bonds. The spatial flexibility of NeuNAc of gangliosides at the binding site of cholera toxin is worked out. The NeuNAc has a limited allowed eulerian space at the binding site of Cholera Toxin (2.4%). The molecular modeling, molecular mechanics and molecular dynamics of disialoganglioside-cholera toxin complex reveal that cholera toxin can accommodate the disialoganglioside GD1A in three different modes. A single mode of binding is permissible for GD1B and GD3. Direct and water mediated hydrogen bonding interactions stabilizes these binding modes and play an essential role in defining the order of specificity for different disialogangliosides towards cholera toxin. This study not only provides models for the disialoganglioside-cholera toxin complexes but also identifies the NeuNAc binding site as a site for design of inhibitors that can restrict the pathogenic activity of cholera toxin.     

Disialogangliosides and Their Interaction with Cholera Toxin—Investigation by Molecular Modeling,Molecular Mechanics and Molecular Dynamics

Abstract

Molecular mechanics and molecular dynamics studies are performed to investigate the conformational preference of cell surface disialogangliosides (GD1A, GD1B and GD3) in aqueous environment. The molecular mechanics calculation reveals that water mediated hydrogen bonding network plays a significant role in the structural stabilization of GD1A, GD1B and GD3. These water mediated hydrogen bonds not only exist between neighboring residues but also exist between residues that are separated by 2 to 3 residues in between. The conformational energy difference between different conformational states of gangliosides correlates very well with the number of water mediated and direct hydrogen bonds. The spatial flexibility of NeuNAc of gangliosides at the binding site of cholera toxin is worked out. The NeuNAc has a limited allowed eulerian space at the binding site of Cholera Toxin (2.4%). The molecular modeling, molecular mechanics and molecular dynamics of disialo- ganglioside-cholera toxin complex reveal that cholera toxin can accommodate the disialo- ganglioside GD1A in three different modes. A single mode of binding is permissible for GD1B and GD3. Direct and water mediated hydrogen bonding interactions stabilizes these binding modes and play an essential role in defining the order of specificity for different disialogangliosides towards cholera toxin. This study not only provides models for the disialoganglioside-cholera toxin complexes but also identifies the NeuNAc binding site as a site for design of inhibitors that can restrict the pathogenic activity of cholera toxin.     

Conformations of higher gangliosides and their binding with cholera toxin - investigation by molecular modeling, molecular mechanics, and molecular dynamics

Molecular mechanics and molecular dynamics studies are performed to investigate the conformational preference of cell surface higher gangliosides (GT1A and GT1B) and their interaction with Cholera Toxin. The water mediated hydrogen bonding network exists between sugar residues in gangliosides. An integrated molecular modeling, molecular mechanics, and molecular dynamics calculation of cholera toxin complexed with GT1A and GT1B reveal that, the active site of cholera toxin can accommodate these higher gangliosides. Direct and water mediated hydrogen bonding interactions stabilize these binding modes and play an essential role in defining the order of specificity for different higher ganglioside towards cholera toxin. This study identifies that the binding site of cholera toxin is shallow and can accommodate a maximum of two NeuNAc residues. The NeuNAc binding site of cholera toxin may be crucial for the design of inhibitors that can prevent the infection of cholera.     

Conformations of Higher Gangliosides and Their Binding with Cholera Toxin—Investigation by Molecular Modeling,Molecular Mechanics,and Molecular Dynamics

Abstract

Molecular mechanics and molecular dynamics studies are performed to investigate the conformational preference of cell surface higher gangliosides (GT1A and GT1B) and their interaction with Cholera Toxin. The water mediated hydrogen bonding network exists between sugar residues in gangliosides. An integrated molecular modeling, molecular mechanics, and molecular dynamics calculation of cholera toxin complexed with GT1A and GT1B reveal that, the active site of cholera toxin can accommodate these higher gangliosides. Direct and water mediated hydrogen bonding interactions stabilize these binding modes and play an essential role in defining the order of specificity for different higher ganglioside towards cholera toxin. This study identifies that the binding site of cholera toxin is shallow and can accommodate a maximum of two NeuNAc residues. The NeuNAc binding site of cholera toxin may be crucial for the design of inhibitors that can prevent the infection of cholera.     

Molecular simulation of N-acetylneuraminic acid analogs and molecular dynamics studies of cholera toxin-Neu5Gc complex

Cholera toxin (CT) is an AB₅ protein complex secreted by the pathogen Vibrio cholera, which is responsible for cholera infection. N-acetylneuraminic acid (NeuNAc) is a derivative of neuraminic acid with nine-carbon backbone. NeuNAc is distributed on the cell surface mainly located in the terminal components of glycoconjugates, and also plays an important role in cell–cell interaction. In our current study, molecular docking and molecular dynamic (MD) simulations were implemented to identify the potent NeuNAc analogs with high-inhibitory activity against CT protein. Thirty-four NeuNAc analogs, modified in different positions C-1/C-2/C-4/C-5/C-7/C-8/C-9, were modeled and docked against the active site of CT protein. Among the 34 NeuNAc analogs, the analog Neu5Gc shows the least extra precision glide score of ?9.52 and glide energy of ?44.71?kcal/mol. NeuNAc analogs block the CT active site residues HIS:13, ASN:90, LYS:91, GLN:56, GLN:61, and TRP:88 through intermolecular hydrogen bonding. The MD simulation for CT-Neu5Gc docking complex was performed using Desmond. MD simulation of CT-Neu5Gc complex reveals the stable nature of docking interaction.     

Theoretical studies on the conformation of monosialogangliosides and disialogangliosides

The preferred conformation of gangliosides GM3, GM2, GM1, GD1a and GD1b have been studied by computing their potential energies. The conformation of NeuNAc in GM3 differs from that expected for the same residue in GM2 and GM1. The NeuNAc residues in GM2 and GM1 exhibit identical conformations. Theory predicts that the terminal NeuNAc of GD1a is conformationally similar to that of GM3 and that the internal one is similar in conformation to those present in GM2 and GM1 in agreement with NMR studies. The differences in chemical shifts of the C2 and C3 carbons of the internal and terminal NeuNAc of GD1a have been attributed to differences in orientation. The present studies suggest that the binding site of cholera toxin is much smaller than that of tetanus toxin. The preferred shape of these gangliosides correlate well with their biological properties.     

Molecular dynamics of sialic acid analogues complex with cholera toxin and DFT optimization of ethylene glycol-mediated zinc nanocluster conjugation

Cholera is an infectious disease caused by cholera toxin (CT) protein of bacterium Vibrio cholerae. A sequence of sialic acid (N-acetylneuraminic acid, NeuNAc or Neu5Ac) analogues modified in its C-5 position is modelled using molecular modelling techniques and docked against the CT followed by molecular dynamics simulations. Docking results suggest better binding affinity of NeuNAc analogue towards the binding site of CT. The NeuNAc analogues interact with the active site residues GLU:11, TYR:12, HIS:13, GLY:33, LYS:34, GLU:51, GLN:56, HIE:57, ILE:58, GLN:61, TRP:88, ASN:90 and LYS:91 through intermolecular hydrogen bonding. Analogues N-glycolyl-NeuNAc, N-Pentanoyl-NeuNAc and N-Propanoyl-NeuNAc show the least XPGscore (docking score) of ?9.90, ?9.16, and ?8.91, respectively, and glide energy of ?45.99, ?42.14 and ?41.66 kcal/mol, respectively. Stable nature of CT-N-glycolyl-NeuNAc, CT-N-Pentanoyl-NeuNAc and CT-N-Propanoyl-NeuNAc complexes was verified through molecular dynamics simulations, each for 40 ns using the software Desmond. All the nine NeuNAc analogues show better score for drug-like properties, so could be considered as suitable candidates for drug development for cholera infection. To improve the enhanced binding mode of NeuNAc analogues towards CT, the nine NeuNAc analogues are conjugated with Zn nanoclusters through ethylene glycol (EG) as carriers. The NeuNAc analogues conjugated with EG-Zn nanoclusters show better binding energy towards CT than the unconjugated nine NeuNAc analogues. The electronic structural optimization of EG-Zn nanoclusters was carried out for optimizing their performance as better delivery vehicles for NeuNAc analogues through density functional theory calculations. These sialic acid analogues may be considered as novel leads for the design of drug against cholera and the EG-Zn nanocluster may be a suitable carrier for sialic acid analogues.     

Molecular dynamics simulation of GM1 in phospholipid bilayer

We report molecular dynamics simulation of fully hydrated lipid bilayer of dimyristoyl phosphatidyl choline (DMPC) at room temperature with ganglioside GM1 attached to it in the upper layer under periodic boundary conditions. The simulation results indicate that the presence of a single GM1 molecule has local effects on the bilayer. Three sugar residues (GalNAc-Gal-Glc) of the pentasaccharide head group of GM1 remain on the lipid surface where as the NeuNAc residue extends out in the aqueous layer. The radial distribution functions suggest ordering of water molecules near the glycerol and carboxyl group of the sialic acid in the upper layer. One of the ceramide chains of GM1, the sphingosine chain, folds up and is stacked under the sugar residues lying on the surface. The other ceramide chain is inserted into the lipid bilayer. The arrangement of the polar head group as well as the acyl chains of the lipids which are immediate neighbours of the GM1 are modified compared to the non-neighbour ones and others at the lower layer. The time average conformation of GM1-pentasaccharide is stabilized by a number of inter residue hydrogen bonds that were observed experimentally. The trajectory average conformation of GM1-pentasaccharide was docked on to the cholera toxin molecule and the minimized complex reveals alternative binding modes between the toxin and the GM1-pentasaccharide moiety. The results of these simulation studies might help to understand the structure and nature of the effects of GM1 on the membrane at atomic resolution.     

Conformational analysis of GT1B ganglioside and its interaction with botulinum neurotoxin type B: a study by molecular modeling and molecular dynamics

The conformational property of oligosaccharide GT1B in aqueous environment was studied by molecular dynamics (MD) simulation using all-atom model. Based on the trajectory analysis, three prominent conformational models were proposed for GT1B. Direct and water-mediated hydrogen bonding interactions stabilize these structures. The molecular modeling and 15 ns MD simulation of the Botulinum Neuro Toxin/B (BoNT/B) - GT1B complex revealed that BoNT/B can accommodate the GT1B in the single binding mode. Least mobility was seen for oligo-GT1B in the binding pocket. The bound conformation of GT1B obtained from the MD simulation of the BoNT/B-GT1B complex bear a close conformational similarity with the crystal structure of BoNT/A-GT1B complex. The mobility noticed for Arg 1268 in the dynamics was accounted for its favorable interaction with terminal NeuNAc. The internal NeuNAc1 tends to form 10 hydrogen bonds with BoNT/B, hence specifying this particular site as a crucial space for the therapeutic design that can restrict the pathogenic activity of BoNT/B.     

Crystal structure of cholera toxin B-pentamer bound to receptor GM1 pentasaccharide.

Cholera toxin (CT) is an AB5 hexameric protein responsible for the symptoms produced by Vibrio cholerae infection. In the first step of cell intoxication, the B-pentamer of the toxin binds specifically to the branched pentasaccharide moiety of ganglioside GM1 on the surface of target human intestinal epithelial cells. We present here the crystal structure of the cholera toxin B-pentamer complexed with the GM1 pentasaccharide. Each receptor binding site on the toxin is found to lie primarily within a single B-subunit, with a single solvent-mediated hydrogen bond from residue Gly 33 of an adjacent subunit. The large majority of interactions between the receptor and the toxin involve the 2 terminal sugars of GM1, galactose and sialic acid, with a smaller contribution from the N-acetyl galactosamine residue. The binding of GM1 to cholera toxin thus resembles a 2-fingered grip: the Gal(beta 1-3)GalNAc moiety representing the "forefinger" and the sialic acid representing the "thumb." The residues forming the binding site are conserved between cholera toxin and the homologous heat-labile enterotoxin from Escherichia coli, with the sole exception of His 13. Some reported differences in the binding affinity of the 2 toxins for gangliosides other than GM1 may be rationalized by sequence differences at this residue. The CTB5:GM1 pentasaccharide complex described here provides a detailed view of a protein:ganglioside specific binding interaction, and as such is of interest not only for understanding cholera pathogenesis and for the design of drugs and development of vaccines but also for modeling other protein:ganglioside interactions such as those involved in GM1-mediated signal transduction.     

Conformational dynamics of sialyl Lewisx in aqueous solution and its interaction with selectinE. A study by molecular dynamics

Three dimensional structures of sialyl Lewis(x) (SLe(x)) in aqueous solution and bound to selectinE are described based on an exhaustive conformational analysis and several long molecular dynamics simulations using different glycosidic regions as starting conformations. It appears from this study that when the oligosaccharide is free in solution the NeuNAcalpha(2-3)Gal segment favors glycosidic conformation in three different regions in the (Phi,Psi) plane with propensity of populations in the ratio 1:8:1. Each one of these structures is characteristically stabilized by specific hydrogen bonding interaction between NeuNAc and Gal. On the other hand, the Gal-GlcNAc-Fuc segment can exist in four different conformational states. Based on the topology of SLe(x) we are able to predict that out of all the allowed conformations in solution only two of these structures possess a geometry that would fit without steric clashes into the binding location of selectinE. In both of these binding modes, segment Gal-GlcNAc-Fuc adopts a unique conformation. The only difference between the two SLe(x) conformers that can successfully bind to selectinE is given by two possible regions in glycosidic space in the fragment NeuNAcalpha(2-3)Gal. A large conformational departure from the crystallographic data is observed for two lysine residues at the binding site of selectinE. These two residues play an important role when SLe(x) binds selectinE in aqueous solution. These findings help reconcile the X-ray data, in which these residues appear to be 1 nm away from SLe(x), with recent liquid NMR data reporting couplings between these protein residues and the sugar.     

Comparison of the glycolipid-binding specificities of cholera toxin and porcineEscherichia coli heat-labile enterotoxin: identification of a receptor-active non-ganglioside glycolipid for the heat-labile toxin in infant rabbit small intestine

The binding specificities of cholera toxin andEscherichia coli heat-labile enterotoxin were investigated by binding of¹²⁵I-labelled toxins to reference glycosphingolipids separated on thin-layer chromatograms and coated in microtitre wells. The binding of cholera toxin was restricted to the GM1 ganglioside. The heat-labile toxin showed the highest affinity for GM1 but also bound, though less strongly, to the GM2, GD2 and GD1b gangliosides and to the non-acid glycosphingolipids gangliotetraosylceramide and lactoneotetraosylceramide. The infant rabbit small intestine, a model system for diarrhoea induced by the toxins, was shown to contain two receptor-active glycosphingolipids for the heat-labile toxin, GM1 ganglioside and lactoneotetraosylceramide, whereas only the GM1 ganglioside was receptor-active for cholera toxin. Preliminary evidence was obtained, indicating that epithelial cells of human small intestine also contain lactoneotetraosylceramide and similar sequences. By computer-based molecular modelling, lactoneotetraosylceramide was docked into the active site of the heat-labile toxin, using the known crystal structure of the toxin in complex with lactose. Interactions which may explain the relatively high toxin affinity for this receptor were found.Abbreviations CT cholera toxin - CT-B B-subunits of cholera toxin - LT Escherichia coli heat-labile enterotoxin - hLT humanEscherichia coli heat-labile enterotoxin - pLT porcineEscherichia coli heat-labile enterotoxin - EI electron ionization     

Insights into the binding specificity of wild type and mutated wheat germ agglutinin towards Neu5Acα(2‐3)Gal: a study by in silico mutations and molecular dynamics simulations

Wheat germ agglutinin (WGA) is a plant lectin, which specifically recognizes the sugars NeuNAc and GlcNAc. Mutated WGA with enhanced binding specificity can be used as biomarkers for cancer. In silico mutations are performed at the active site of WGA to enhance the binding specificity towards sialylglycans, and molecular dynamics simulations of 20 ns are carried out for wild type and mutated WGAs (WGA1, WGA2, and WGA3) in complex with sialylgalactose to examine the change in binding specificity. MD simulations reveal the change in binding specificity of wild type and mutated WGAs towards sialylgalactose and bound conformational flexibility of sialylgalactose. The mutated polar amino acid residues Asn114 (S114N), Lys118 (G118K), and Arg118 (G118R) make direct and water mediated hydrogen bonds and hydrophobic interactions with sialylgalactose. An analysis of possible hydrogen bonds, hydrophobic interactions, total pair wise interaction energy between active site residues and sialylgalactose and MM‐PBSA free energy calculation reveals the plausible binding modes and the role of water in stabilizing different binding modes. An interesting observation is that the binding specificity of mutated WGAs (cyborg lectin) towards sialylgalactose is found to be higher in double point mutation (WGA3). One of the substituted residues Arg118 plays a crucial role in sugar binding. Based on the interactions and energy calculations, it is concluded that the order of binding specificity of WGAs towards sialylgalactose is WGA3 > WGA1 > WGA2 > WGA. On comparing with the wild type, double point mutated WGA (WGA3) exhibits increased specificity towards sialylgalactose, and thus, it can be effectively used in targeted drug delivery and as biological cell marker in cancer therapeutics. Copyright © 2014 John Wiley & Sons, Ltd.     

A comparative study of trypsin specificity based on QM/MM molecular dynamics simulation and QM/MM GBSA calculation

Hydrogen bonding and polar interactions play a key role in identification of protein-inhibitor binding specificity. Quantum mechanics/molecular mechanics molecular dynamics (QM/MM MD) simulations combined with DFT and semi-empirical Hamiltonian (AM1d, RM1, PM3, and PM6) methods were performed to study the hydrogen bonding and polar interactions of two inhibitors BEN and BEN1 with trypsin. The results show that the accuracy of treating the hydrogen bonding and polar interactions using QM/MM MD simulation of PM6 can reach the one obtained by the DFT QM/MM MD simulation. Quantum mechanics/molecular mechanics generalized Born surface area (QM/MM-GBSA) method was applied to calculate binding affinities of inhibitors to trypsin and the results suggest that the accuracy of binding affinity prediction can be significantly affected by the accurate treatment of the hydrogen bonding and polar interactions. In addition, the calculated results also reveal the binding specificity of trypsin: (1) the amidinium groups of two inhibitors generate favorable salt bridge interaction with Asp189 and form hydrogen bonding interactions with Ser190 and Gly214, (2) the phenyl of inhibitors can produce favorable van der Waals interactions with the residues His58, Cys191, Gln192, Trp211, Gly212, and Cys215. This systematic and comparative study can provide guidance for the choice of QM/MM MD methods and the designs of new potent inhibitors targeting trypsin.     

Theoretical investigation on the glycan‐binding specificity of Agrocybe cylindracea galectin using molecular modeling and molecular dynamics simulation studies

Galectins are β‐galactoside binding proteins which have the ability to serve as potent antitumor, cancer biomarker, and induce tumor cell apoptosis. Agrocybe cylindracea galectin (ACG) is a fungal galectin which specifically recognizes α(2,3)‐linked sialyllactose at the cell surface that plays extensive roles in the biological recognition processes. To investigate the change in glycan‐binding specificity upon mutations, single point and double point site‐directed in silico mutations are performed at the binding pocket of ACG. Molecular dynamics (MD) simulation studies are carried out for the wild‐type (ACG) and single point (ACG1) and double point (ACG2) mutated ACGs to investigate the dynamics of substituted mutants and their interactions with the receptor sialyllactose. Plausible binding modes are proposed for galectin–sialylglycan complexes based on the analysis of hydrogen bonding interactions, total pair‐wise interaction energy between the interacting binding site residues and sialyllactose and binding free energy of the complexes using molecular mechanics–Poisson–Boltzmann surface area. Our result shows that high contribution to the binding in different modes is due to the direct and water‐mediated hydrogen bonds. The binding specificity of double point mutant Y59R/N140Q of ACG2 is found to be high, and it has 26 direct and water‐mediated hydrogen bonds with a relatively low‐binding free energy of −47.52 ± 5.2 kcal/mol. We also observe that the substituted mutant Arg59 is crucial for glycan‐binding and for the preference of α(2,3)‐linked sialyllactose at the binding pocket of ACG2 galectin. When compared with the wild‐type and single point mutant, the double point mutant exhibits enhanced affinity towards α(2,3)‐linked sialyllactose, which can be effectively used as a model for biological cell marker in cancer therapeutics. Copyright © 2015 John Wiley & Sons, Ltd.     

Solution dynamics of the oligosaccharide moiety of ganglioside GM1: Comparison of solution conformations with the bound state conformation in association with cholera toxin B-pentamer

The solution dynamics of the oligosaccharide moiety of ganglioside G_M1 have been determined by use of a combination of ¹H rotating frame Overhauser effect measurements and restrained molecular dynamics simulations, It is found that the Galβ1-3 and NeuNAc moieties which are primarily recognized by cholera toxin both exhibit considerable torsional flexibility about their respective glycosidic linkages. A comparison with the bound state conformation of the ganglioside in association with cholera toxin B-pentamer, shows that a low energy conformation of the oligosaccharide, which closely approximates the globel minimum, is selected upon binding.     

Protein reagent modification of cholera toxin: characterization of effects on antigenic, receptor-binding and toxic properties.

The effects of protein modification procedures on the biologically most important properties of cholera toxin, i.e. the toxic activity, the GM1 receptor-binding capacity and the antigenic (antibody-fixing) properties, have been studied quantitatively using microgram amounts or less of toxin protein. Most of the 24 group-specific reagents used had either no inhibitory effect on the toxic or the combination of GM1-binding and antibody-fixing properties of cholera toxin, or they had a concomitant inhibitory effect on these activities. Separate testing of GM1- and antibody-binding revealed a close, but not absolute, structural association between these properties, Amino group reactive substances were particularly effective in decreasing the GM1-binding activity, while leucine aminopeptidase had no effect. This suggests that lysine residues may be involved in binding toxin to the acidic GM1 receptor. Sodium dodecylsulphate and mercaptoethanol, which caused dissociation of the subunits of cholera toxin as indicated by polyacrylamide gel electrophoresis, abolished toxicity without inhibiting the concomitant GM1- and antibody-binding properties of the toxin. Similar differential effects were also obtained with three reagents which did not seem to change the aggregation state of the toxin. These substances all had specificity for arginine, suggesting that arginyl residues of the toxin molecule may be involved in a 'toxic site' distinct from the receptor-binding site(s). A selective effect on the toxic site was also found by treating the toxin with carboxypeptidase or trypsin in the presence of urea; in the absence of urea no enzymic effect on any toxin property was noted.     

Membrane traffic and the cellular uptake of cholera toxin.

In nature, cholera toxin (CT) and the structurally related E. coli heat labile toxin type I (LTI) must breech the epithelial barrier of the intestine to cause the massive diarrhea seen in cholera. This requires endocytosis of toxin-receptor complexes into the apical endosome, retrograde transport into Golgi cisternae or endoplasmic reticulum (ER), and finally transport of toxin across the cell to its site of action on the basolateral membrane. Targeting into this pathway depends on toxin binding ganglioside GM1 and association with caveolae-like membrane domains. Thus to cause disease, both CT and LTI co-opt the molecular machinery used by the host cell to sort, move, and organize their cellular membranes and substituent components.     

Analysis of cholera toxin-ganglioside interactions by flow cytometry

Cholera toxin entry into mammalian cells is mediated by binding of the pentameric B subunit (CTB) to ganglioside GM(1) in the cell membrane. We used flow cytometry to quantitatively measure in real time the interactions of fluorescently labeled pentameric cholera toxin B-subunit (FITC-CTB) with its ganglioside receptor on microsphere-supported phospholipid membranes. A model that describes the multiple steps of this mode of recognition was developed to guide our flow cytometric experiments and extract relevant equilibrium and kinetic rate constants. In contrast to previous studies, our approach takes into account receptor cross-linking, an important feature for multivalent interactions. From equilibrium measurements, we determined an equilibrium binding constant for a single subunit of FITC-CTB binding monovalently to GM(1) presented in bilayers of approximately 8 x 10(7) M(-1) while that for binding to soluble GM(1)-pentasaccharide was found to be approximately 4 x 10(6) M(-1). From kinetic measurements, we determined the rate constant for dissociation of a single site of FITC-CTB from microsphere-supported bilayers to be (3.21 +/- 0.03) x 10(-3) s(-1), and the rate of association of a site on FITC-CTB in solution to a GM(1) in the bilayer to be (2.8 +/- 0.4) x 10(4) M(-1) s(-1). These values yield a lower estimate for the equilibrium binding constant of approximately 1 x 10(7) M(-1). We determined the equilibrium surface cross-linking constant [(1.1 +/- 0.1) x 10(-12) cm(2)] and from this value and the value for the rate constant for dissociation derived a value of approximately 3.5 x 10(-15) cm(2) s(-1) for the forward rate constant for cross-linking. We also compared the interaction of the receptor binding B-subunit with that of the whole toxin (A- and B-subunits). Our results show that the whole toxin binds with approximately 100-fold higher avidity than the pentameric B-subunit alone which is most likely due to the additional interaction of the A(2)-subunit with the membrane surface. Interaction of cholera toxin B-subunit and whole cholera toxin with gangliosides other than GM(1) revealed specific binding only to GD1(b) and asialo-GM(1). These interactions, however, are marked by low avidity and require high receptor concentrations to be observed.