Stabilization of the structure of horse plasma vitamin D binding protein by disulfide bonds.

Vitamin D binding protein (DBP) was isolated from horse plasma in a four-step procedure that involved Affi-Gel Blue affinity chromatography, gel filtration, hydroxylapatite chromatography, and anion exchange high-pressure liquid chromatography. The yield of DBP from 80 mL of plasma was 6-7 mg. Horse plasma DBP closely resembles other plasma DBPs, being a tryptophan-free protein of Mr 53,000. It is able to bind to and block the polymerization of monomeric actin. The secondary structure of DBP was calculated from circular dichroism measurements to be 39% alpha-helix, 42% beta-sheet, and 19% random coil. Circular dichroism and fluorescence studies revealed that the disulfide bonds of DBP contribute substantial structural stabilization to the molecule with respect to thermal denaturation. The thermal stability of DBP can be used to advantage. Incorporation of a brief treatment at 70 degrees C into the preparative scheme enables omission of one chromatographic step, without detectable alteration of the purified product.     

Paramyxovirus membrane protein enhances antibody production to new antigenic determinants in the actin molecule: a model for virus-induced autoimmunity.

Paramyxovirus membrane (M) protein specifically binds to cellular actin but not to bovine serum albumin or myoglobin, as determined by affinity chromatography and enzyme-linked immunosorbent assay. The binding site for M protein on actin is different from the binding sites for antiactin antibodies. The interaction of M protein with actin resulted in production of antibodies to several new antigenic sites on the actin molecule. Five rabbits immunized with actin alone produced antibodies against the N-terminal sequence (residues 1 to 39). Another five rabbits immunized with a mixture of M protein and actin produced antibodies against a C-terminal fragment and a central region as well as the N-terminal fragment. By immunoblotting with proteolytic fragments of actin, the new antigenic sites were located between amino acid residues 40 to 113, 114 to 226, and 227 to 375. Antisera taken from some patients with recent measles virus infections demonstrated antiactin antibodies to sites other than the N-terminal fragment of actin (amino acids 1 to 39). The interaction of paramyxovirus M protein with actin and the subsequent production of antibodies to new antigenic sites may serve as a model for one of the mechanisms of virus-induced autoimmunity.     

DNA-induced partial unfolding of prion protein leads to its polymerisation to amyloid

The full-length mouse recombinant prion protein (23-231 amino acid residues) contains all of its structural elements viz. three alpha-helices and a short two-stranded antiparallel beta-sheet in its C-terminal fragment comprising 121-231 amino acid residues. The incubated mixture of this prion protein fragment and nucleic acid results in the formation of amyloid fibres evidenced from electron microscopy, birefringence and fluorescence of the fibre bound Congo Red and Thioflavin T dyes, respectively. The secondary structure of the amyloid formed in nucleic acid solution is similar to the in vivo isolated prion protein 27-30 amyloid but unlike in it, a hydrophobic milieu is absent in the 121-231 amyloid. Thermal denaturation study demonstrates a partial unfolding of the protein fragment in nucleic acid solution. We propose that nucleic acid catalyses unfolding of prion protein helix 1 followed by a nucleation-dependent polymerisation of the protein to amyloid.     

The carboxy-terminal third of dystrophin enhances actin binding activity

Dystrophin is an actin binding protein that is thought to stabilize the cardiac and skeletal muscle cell membranes during contraction. Here, we investigated the contributions of each dystrophin domain to actin binding function. Cosedimentation assays and pyrene-actin fluorescence experiments confirmed that a fragment spanning two-thirds of the dystrophin molecule [from N-terminal actin binding domain (ABD) 1 through ABD2] bound actin filaments with high affinity and protected filaments from forced depolymerization, but was less effective in both assays than full-length dystrophin. While a construct encoding the C-terminal third of dystrophin displayed no specific actin binding activity or competition with full-length dystrophin, our data show that it confers an unexpected regulation of actin binding by the N-terminal two-thirds of dystrophin when present in cis. Time-resolved phosphorescence anisotropy experiments demonstrated that the presence of the C-terminal third of dystrophin in cis also influences actin interaction by restricting actin rotational amplitude. We propose that the C-terminal region of dystrophin allosterically stabilizes an optimal actin binding conformation of dystrophin.     

Recombinant 70-kDa protein from the amino-terminal region of rat fibronectin inhibits binding of fibronectin to cells and bacteria

Binding of fibronectin to substrate-attached cells and to Staphylococcus aureus is mediated by the amino-terminal 70-kDa portion of fibronectin. The 70-kDa amino-terminus is composed of nine type I and two type II internal homology units, each containing two intrachain disulfide bonds. The exact structural features of the 70-kDa amino-terminus that are necessary for binding to cells and bacteria are not known. We characterized a recombinant 70-kDa protein from the amino-terminus of rat fibronectin using a baculovirus expression system. Recombinant 70-kDa (r70kDa) protein was easily purified in high amounts from the conditioned medium by affinity chromatography on gelatin-agarose. Secretion was much less when N-linked glycosylation was blocked by tunicamycin. Like the native fragment, the r70kDa protein contains intrachain disulfide bonds. In addition, the r70kDa protein was indistinguishable from the nonrecombinant 70-kDa fragment in its ability to compete for binding sites on fibroblasts and S. aureus. Thus, the r70kDa protein retains the important functional characteristics of the native fragment. This expression system is well adapted to studying the structural features important for the interaction of 70-kDa protein with cells.     

Solution structure and interdomain interactions of the Saccharomyces cerevisiae "TATA binding protein" (TBP) probed by radiolytic protein footprinting

Although atomic-resolution crystal structures of the conserved C-terminal domain of several species of TBP and their complexes with DNA have been determined, little information is available concerning the structure in solution of full-length TBP containing both the conserved C-terminal and nonconserved N-terminal domains. Quantitation of the amino acid side chain oxidation products generated by synchrotron X-ray radiolysis by mass spectrometry has been used to determine the solvent accessibility of individual residues in monomeric Saccharomyces cerevisiae TATA binding protein (TBP) free in solution and in the TBP-DNA complex. Amino acid side chains within the C-terminal domain of unliganded full-length TBP that are predicted to be accessible from crystal structures of the isolated domain are protected from oxidation. Residues within the N-terminal domain are also protected from oxidation in both the absence and presence of DNA. Some residues within the DNA-binding "saddle" of the C-terminal domain are protected upon formation of a TBP-DNA complex as expected, while others are protected in both the absence and presence of bound DNA. In addition, residues on the upper side of the beta-sheets undergo reactivity changes as a function of DNA binding. These data suggest that the DNA-binding saddle of monomeric unliganded yeast TBP is only partially accessible to solvent, the N-terminal domain is partially structured, and the N- and C-terminal domains form a different set of contacts in the free and DNA-bound protein. The functional implications of these results are discussed.     

The immunoglobulin D-binding part of the outer membrane protein MID from Moraxella catarrhalis comprises 238 amino acids and a tetrameric structure

Moraxella catarrhalis IgD-binding protein (MID), a 200-kDa outer membrane protein comprising 2,139 amino acids, has recently been isolated and shown to display a unique and specific affinity for human IgD. To identify the IgD-binding region, MID was digested with proteases. In addition, a series of truncated fragments of MID were manufactured and expressed in Escherichia coli followed by analysis for IgD binding in Western and dot blots. The smallest fragment with essentially preserved IgD binding was comprised of 238 amino acid residues (MID(962-1200)). Shorter recombinant proteins gradually lost IgD-binding capacity, and the shortest IgD-binding fragment comprising 157 amino acids (MID(985-1142)) displayed a 1,000-fold reduced IgD binding compared with the full-length molecule. The truncated MID(962-1200) was efficiently attracted to a standard IgD serum and to purified myeloma IgD(kappa) and IgD(lambda) sera but not to IgG, IgM, or IgA myeloma sera. Furthermore, the fragment specifically bound to peripheral blood B lymphocytes, and the binding was inhibited by preincubation with anti-IgD-Fab polyclonal antibodies. Results obtained by introducing five amino acids randomly into MID(962-1200) using transposons suggested that alpha-helix structures were important for IgD binding. Ultracentrifugation experiments and gel electrophoresis revealed that native MID(962-1200) was a tetramer. Interestingly, tetrameric MID(962-1200) attracted IgD more than 20-fold more efficiently than the monomeric form. Thus, a tetrameric structure of MID(962-1200) is crucial for optimal IgD-binding capacity.     

Biologic significance of disulfide bonds in human IgE molecules.

E myeloma protein, PS, was reduced in different concentrations of dithiothreitol (DTT) for 1 hr followed by alkylation with 14C-iodoacetamide. The affinity of the reduced-alkylated molecules for target cells was evaluated by their ability 1) to sensitize primate skin in a reversed P-K reaction, 2) to sensitize human basophils in a reversed-type histamine release and 3) to block passive sensitization with reaginic antibody. Antibody-epsilon0 antibody was employed for reversed type reactions to avoid participation of cell-bound normal IgE in the reactions. The sensitizing activity of IgE did not change following reduction in 1 mM DTT, which split inter-heavy-light chain disulfide bond. The activity of IgE significantly diminished after reduction in 2 mM DTT followed by alkylation. This treatment resulted in the cleavage of two intra-epsilon-chain disulfide bonds, which are present between the hinge and the Fd portion of the molecules. The reduced-alkylated protein was capable of sensitizing primate skin and human basophils, however, a much higher concentration of the reduced-alkylated protein than the native protein was required for passive sensitization. The optimal sensitization period for the reversed P-K reaction was 3 hr with the reduced-alkylated protein. The protein had the ability to block passive sensitization with reaginic antibody. The reduced-alkylated protein and the native protein were labeled with 125I, and binding of these proteins with human basophils was examined by autoradiography. The results showed that affinity of the reduced-alkylated protein for basophils was less than that of native protein. Since the disulfide bonds split by 2 mM DTT were not included in the Fc portion of the molecules, the Fc fragment was obtained from the reduced-alkylated protein and was tested for affinity for basophils. It was found that the Fc fragment had higher affinity than the reduced-alkylated protein. Recovery of the affinity by papain digestion strongly suggested that cleavage of disulfide bonds in the Fab portion of the molecules induced conformational changes in the Fc portion which is involved in binding to the target cells. Reduction of IgE with 10 mM DTT followed by alkylation resulted in cleavage of 5 disulfide bonds, which is accompanied by a loss of both sensitizing and blocking activities. The fifth disulfide bond which was cleaved by 10 mM DTT, but not by 2 mM DTT, appears to be an inter-heavy chain disulfide bond in the Fc portion of the epsilon-chains. Neither epsilon1 nor epsilon2 determinants in the Fc portion of epsilon-chains were degraded by this treatment.     

Interaction of myosin subfragment 1 with forms of monomeric actin

The ability of myosin subfragment 1 to interact with monomeric actin complexed to sequestering proteins was tested by a number of different techniques such as affinity absorption, chemical cross-linking, fluorescence titration, and competition procedures. For affinity absorption, actin was attached to agarose immobilized DNase I. Both chymotryptic subfragment 1 isoforms (S1A1 and S1A2) were retained by this affinity matrix. Fluorescence titration employing pyrenyl-actin in complex with deoxyribonuclease I (DNase I) or thymosin beta4 demonstrated S1 binding to these actin complexes. A K(D) of 5 x 10(-8) M for S1A1 binding to the actin-DNase I complex was determined. Fluorescence titration did not indicate binding of S1 to actin in complex with gelsolin segment 1 (G1) or vitamin D-binding protein (DBP). However, fluorescence competition experiments and analysis of tryptic cleavage patterns of S1 indicated its interaction with actin in complex with DBP or G1. Formation of the ternary DNase I-acto-S1 complex was directly demonstrated by sucrose density sedimentation. S1 binding to G-actin was found to be sensitive to ATP and an increase in ionic strength. Actin fixed in its monomeric state by DNase I was unable to significantly stimulate the Mg2+-dependent S1-ATPase activity. Both wild-type and a mutant of Dictyostelium discoideum myosin II subfragment 1 containing 12 additional lysine residues within an insertion of 20 residues into loop 2 (K12/20-Q532E) were found to also interact with actin-DNase I complex. Binding of the K12/20-Q532E mutant to the actin-DNase I complex occurred with higher affinity than wild-type S1 and was less sensitive to mono- and divalent cations.     

Tropomyosin ends determine the stability and functionality of overlap and troponin T complexes

Tropomyosin binds end to end along the actin filament. Tropomyosin ends, and the complex they form, are required for actin binding, cooperative regulation of actin filaments by myosin, and binding to the regulatory protein, troponin T. The aim of the work was to understand the isoform and structural specificity of the end-to-end association of tropomyosin. The ability of N-terminal and C-terminal model peptides with sequences of alternate alpha-tropomyosin isoforms, and a troponin T fragment that binds to the tropomyosin overlap, to form complexes was analyzed using circular dichroism spectroscopy. Analysis of N-terminal extensions (N-acetylation, Gly, AlaSer) showed that to form an overlap complex between the N-terminus and the C-terminus requires that the N-terminus be able to form a coiled coil. Formation of a ternary complex with the troponin T fragment, however, effectively takes place only when the overlap complex sequences are those found in striated muscle tropomyosins. Striated muscle tropomyosins with N-terminal modifications formed ternary complexes with troponin T that varied in affinity in the order: N-acetylated > Gly > AlaSer > unacetylated. The circular dichroism results were corroborated by native gel electrophoresis, and the ability of the troponin T fragment to promote binding of full-length tropomyosins to filamentous actin.     

cDNA cloning and 1.75 A crystal structure determination of PPL2, an endochitinase and N-acetylglucosamine-binding hemagglutinin from Parkia platycephala seeds

Parkia platycephala lectin 2 was purified from Parkia platycephala (Leguminosae, Mimosoideae) seeds by affinity chromatography and RP-HPLC. Equilibrium sedimentation and MS showed that Parkia platycephala lectin 2 is a nonglycosylated monomeric protein of molecular mass 29 407+/-15 Da, which contains six cysteine residues engaged in the formation of three intramolecular disulfide bonds. Parkia platycephala lectin 2 agglutinated rabbit erythrocytes, and this activity was specifically inhibited by N-acetylglucosamine. In addition, Parkia platycephala lectin 2 hydrolyzed beta(1-4) glycosidic bonds linking 2-acetoamido-2-deoxy-beta-D-glucopyranose units in chitin. The full-length amino acid sequence of Parkia platycephala lectin 2, determined by N-terminal sequencing and cDNA cloning, and its three-dimensional structure, established by X-ray crystallography at 1.75 A resolution, showed that Parkia platycephala lectin 2 is homologous to endochitinases of the glycosyl hydrolase family 18, which share the (betaalpha)8 barrel topology harboring the catalytic residues Asp125, Glu127, and Tyr182.     

Unfolding intermediates of the extracellular domain of the interferon gamma receptor.

Reduction of proteins which require disulfide bonds to be stable in the folded state is accompanied by step-wise unfolding. A soluble human interferon gamma receptor produced in Escherichia coli was used to investigate the kinetics of formation of unfolding intermediates. The protein includes 8 cysteine residues forming four disulfide bonds. It was reduced by using either dithiothreitol or the thioredoxin reduction system. Reduction with dithiothreitol resulted in formation of mainly four monomeric unfolding species as visualized by sodium dodecyl sulfate-polyacrylamide gels. The enzymatically catalyzed reaction produced only small amounts of two monomeric products and mostly delivered oligomeric and polymeric forms. In both cases, the ligand binding capacity of the receptor was significantly reduced immediately after appearance of the first intermediate. The intermediates involved interchange of disulfide bonds and did not show ligand binding capacity. Some of them were recognized by specific antibodies which detect conformational epitopes on the native interferon gamma receptor. On the basis of the antibody binding, a preliminary characterization of the formed intermediates was attempted. When the soluble receptor was reduced in the presence of denaturing agents, the reduction products were different from the unfolding intermediates generated in the absence of denaturants.     

Disulfide structure of the leucine-rich repeat C-terminal cap and C-terminal stalk region of Nogo-66 receptor

Nogo-66 receptor (NgR1) is a leucine-rich repeat (LRR) protein that forms part of a signaling complex modulating axon regeneration. Previous studies have shown that the entire LRR region of NgR1, including the C-terminal cap of the LRR, LRRCT, is needed for ligand binding, and that the adjacent C-terminal region (CT stalk) of the NgR1 contributes to interaction with its coreceptors. To provide structure-based information for these interactions, we analyzed the disulfide structure of full-length NgR1. Our analysis revealed a novel disulfide structure in the C-terminal region of the NgR1, wherein the two Cys residues, Cys-335 and Cys-336, in the CT stalk are disulfide-linked to Cys-266 and Cys-309 in the LRRCT region: Cys-266 is linked to Cys-335, and Cys-309 to Cys-336. The other two Cys residues, Cys-264 and Cys-287, in the LRRCT region are disulfide-linked to each other. The analysis also showed that Cys-419 and Cys-429, in the CT stalk region, are linked to each other by a disulfide bond. Although published crystal structures of a recombinant fragment of NgR1 had revealed a disulfide linkage between Cys-266 and Cys-309 in the LRRCT region and we verified its presence in the corresponding fragment, this is artificially caused by the truncation of the protein, since this linkage was not detected in intact NgR1 or a slightly larger fragment containing Cys-335 and Cys-336. A structural model of the LRRCT with extended residues 311-344 from the CT stalk region is proposed, and its function in coreceptor binding is discussed.     

Assistance of maltose binding protein to the in vivo folding of the disulfide-rich C-terminal fragment from Plasmodium falciparum merozoite surface protein 1 expressed in Escherichia coli

The C-terminal fragment of Plasmodium falciparum merozoite surface protein 1 (F19) is a leading candidate for the development of a malaria vaccine. Successful vaccination trials on primates, immunochemistry, and structural studies have shown the importance of its native conformation for its protective role against infection. F19 is a disulfide-rich protein, and the correct pairing of its 12 half-cystines is required for the native state of the protein. F19 has been produced in the Escherichia coli periplasm, which has an oxidative environment favorable for the formation of disulfide bonds. F19 was either expressed as a fusion with the maltose binding protein (MBP) or directly addressed to the periplasm by fusing it with the MBP signal peptide. Direct expression of F19 in the periplasm led to a misfolded protein with a heterogeneous distribution of disulfide bridges. On the contrary, when produced as a fusion protein with E. coli MBP, the F19 moiety was natively folded. Indeed, after proteolysis of the fusion protein, the resulting F19 possesses the structural characteristics and the immunochemical reactivity of the analogous fragment produced either in baculovirus-infected insect cells or in yeast. These results demonstrate that the positive effect of MBP in assisting the folding of passenger proteins extends to the correct formation of disulfide bridges in vivo. Although proteins or protein fragments fused to MBP have been frequently expressed with success, our comparative study evidences for the first time the helping property of MBP in the oxidative folding of a disulfide-rich protein.     

Involvement of disulfide bonds and histidine 172 in a unique beta-sheet to alpha-helix transition of alpha 1-acid glycoprotein at the biomembrane interface

Human alpha(1)-acid glycoprotein (AGP), which is comprised of 183 amino acid residues and 5 carbohydrate chains, is a major plasma protein that binds to basic and neutral drugs as well as to steroid hormones. It has a beta-sheet-rich structure in aqueous solution. Our previous findings suggest that AGP forms an alpha-helix structure through an interaction with biomembranes. We report herein on a study of the mechanism of alpha-helix formation in AGP using various modified AGPs. The disulfide reduced AGP (R-AGP) was extensively unfolded, whereas asialylated AGP (A-AGP) maintained the native structure. Intriguingly, reduced and asialylated AGP (RA-AGP) increased the alpha-helix content as observed in the presence of biomembrane models, and showed a significant decrease in ligand binding capacity. This suggests that AGP has an innate tendency to form an alpha-helix structure, and disulfide bonds are a key factor in the conformational transition between the beta-sheet and alpha-helix structures. However, RA-AGP with all histidine residues chemically modified (HRA-AGP) was found to lose the intrinsic ability to form an alpha-helix structure. Furthermore, disulfide reduction of the H172A mutant expressed in Pichia pastoris also caused a similar loss of folding ability. The present results indicate that disulfide bonds and the C-terminal region, including H172 of AGP, play important roles in alpha-helix formation in the interaction of the protein with biomembranes.     

Functional modulation of the isolated glycoprotein Ib binding domain of von Willebrand factor expressed in Escherichia coli

We have expressed in Escherichia coli the domain of von Willebrand factor (vWF) containing the binding site for platelet glycoprotein (GP) Ib and used it to study the regulation of vWF-platelet interaction. The recombinant fragment, comprising residues 445-733 of the mature vWF subunit and designated rvWF445-733, did not have the native conformation of the corresponding domain in the intact molecule because, in order to prevent formation of random aggregates, the seven cysteine residues in the sequence were reduced and alkylated. Unlike native vWF, rvWF445-733 bound to GP Ib in the absence of any modulator, suggesting that the lack of disulfide bonds and/or carbohydrate side chains within this domain may expose platelet interaction sites. In the presence of two modulators, the glycopeptide ristocetin and the snake protein botrocetin, rvWF445-733 inhibited native vWF binding to GP Ib as well as platelet aggregation mediated by vWF, suggesting that both the fragment and the native molecule interact with the same site on platelets. This conclusion was also supported by the observation that the recombinant fragment competed with the binding to platelets of an anti-GP Ib monoclonal antibody known to inhibit vWF binding. Botrocetin formed a complex with rvWF445-733, but the affinity of this interaction was approximately 25-fold lower than with native vWF. However, the complexes of botrocetin with either rvWF445-733 or multimeric native vWF bound to GP Ib with similar dissociation constant. Therefore, conformational attributes of vWF regulate its affinity for botrocetin, but once the complex is formed, interaction with GP Ib is independent of native vWF conformation. These findings provide insights into the regulation of vWF-platelet interaction.     

Mapping the cofilin binding site on yeast G-actin by chemical cross-linking

Cofilin is a major cytoskeletal protein that binds to both monomeric actin (G-actin) and polymeric actin (F-actin) and is involved in microfilament dynamics. Although an atomic structure of the G-actin-cofilin complex does not exist, models of the complex have been built using molecular dynamics simulations, structural homology considerations, and synchrotron radiolytic footprinting data. The hydrophobic cleft between actin subdomains 1 and 3 and, alternatively, the cleft between actin subdomains 1 and 2 have been proposed as possible high-affinity cofilin binding sites. In this study, the proposed binding of cofilin to the subdomain 1/subdomain 3 region on G-actin has been probed using site-directed mutagenesis, fluorescence labeling, and chemical cross-linking, with yeast actin mutants containing single reactive cysteines in the actin hydrophobic cleft and with cofilin mutants carrying reactive cysteines in the regions predicted to bind to G-actin. Mass spectrometry analysis of the cross-linked complex revealed that cysteine 345 in subdomain 1 of mutant G-actin was cross-linked to native cysteine 62 on cofilin. A cofilin mutant that carried a cysteine substitution in the α3-helix (residue 95) formed a cross-link with residue 144 in actin subdomain 3. Distance constraints imposed by these cross-links provide experimental evidence for cofilin binding between actin subdomains 1 and 3 and fit a corresponding docking-based structure of the complex. The cross-linking of the N-terminal region of recombinant yeast cofilin to actin residues 346 and 374 with dithio-bis-maleimidoethane (12.4 Å) and via disulfide bond formation was also documented. This set of cross-linking data confirms the important role of the N-terminal segment of cofilin in interactions with G-actin.     

Isolation and characterisation of a sialic acid-binding lectin (carcinoscorpin) from Indian horseshoe crab Carcinoscorpius rotunda cauda

A sialic acid-binding lectin, carcinoscorpin, has been purified to apparent homogeneity in 40% yield from the Indian horseshoe carb, Carcinoscorpius rotunda cauda. This glycoprotein lectin of molecular weight 420,000 was composed of two non-identical subunits of molecular weights 27,000 and 28,000 as determined by gel electrophoresis in the presence of sodium dodecyl sulfate. The hemagglutination activity of the lectin was susceptible to guanidine-HCl; modification of tyrosyl and tryptophanyl residues also inhibited the activity although alkylation of the -SH group, reduction of disulfide bonds or modification of amino and carboxyl groups were without any effect. The monomeric form of the lectin produced by succinylation of native protein was inactive in binding to sialoglycoconjugates.     

Computer simulations of glycolytic enzyme interactions with F-actin

Muscle actin and fructose-1,6-bisphosphate aldolase (aldolase) were chemically crosslinked to produce an 80 kDa product representing one subunit of aldolase linked to one subunit of actin. Hydroxylamine digestion of the crosslinked product resulted in two 40.5 kDa fragments, one that was aldolase linked to the 12 N-terminal residues of actin. Brownian dynamics simulations of muscle aldolase and GAPDH with F-actin (muscle, yeast, and various mutants) estimated the association free energy. Mutations of residues 1-4 of muscle actin to Ala individually or two in combination of the first four residues reduced the estimated binding free energy. Simulations showed that muscle aldolase binds with the same affinity to the yeast actin as to the double mutated muscle actin; these mutations make the N-terminal of muscle actin identical to yeast, supporting the conclusion that the actin N-terminus participates in binding. Because the depth of free energy wells for yeast and the double mutants is less than for native rabbit actin, the simulations support experimental findings that muscle aldolase and GAPDH have a higher affinity for muscle actin than for yeast actin. Furthermore, Brownian dynamics revealed that the lower affinity of yeast actin for aldolase and GAPDH compared to muscle actin, was directly related to the acidic residues at the N-terminus of actin.