Molecular dynamics simulation of amyloid beta dimer formation

Recent experiments with amyloid beta (Abeta) peptide indicate that formation of toxic oligomers may be an important contribution to the onset of Alzheimer's disease. The toxicity of Abeta oligomers depends on their structure, which is governed by assembly dynamics. Due to limitations of current experimental techniques, a detailed knowledge of oligomer structure at the atomic level is missing. We introduce a molecular dynamics approach to study Abeta dimer formation. 1), We use discrete molecular dynamics simulations of a coarse-grained model to identify a variety of dimer conformations; and 2), we employ all-atom molecular mechanics simulations to estimate thermodynamic stability of all dimer conformations. Our simulations of a coarse-grained Abeta peptide model predicts 10 different planar beta-strand dimer conformations. We then estimate the free energies of all dimer conformations in all-atom molecular mechanics simulations with explicit water. We compare the free energies of Abeta(1-42) and Abeta(1-40) dimers. We find that 1), dimer conformations have higher free energies compared to their corresponding monomeric states; and 2), the free-energy difference between the Abeta(1-42) and the corresponding Abeta(1-40) dimer conformation is not significant. Our results suggest that Abeta oligomerization is not accompanied by the formation of thermodynamically stable planar beta-strand dimers.     

Interaction between alpha 1 and beta 1 subunits of human hemoglobin

We prepared normal and modified alpha and beta globulin chains in which C-terminal residues were enzymatically removed. The CD spectra of the deoxy form of these chains and the reconstituted modified Hb's were measured in the Soret region. The CD spectra of the modified Hb's were markedly different from the arithmetic means of respective spectra of their constituent chains. This difference was ascribed to the interaction between alpha 1 and beta 1 subunits to make the alpha 1 beta 1 dimer. The peak wavelength of the difference CD spectra could be classified into two groups, one was 433 +/- 1 nm and the other 437 +/- 1 nm. A comparison of this classification with the previously identified quaternary structures revealed that the R and T structures showed a maximum of the difference CD spectra at 437 +/- 1 nm and 433 +/- 1 nm, respectively. These results indicated that the R and T structures differed in the interaction between alpha 1 and beta 1 subunits.     

Differences in the activation mechanism between the alpha and beta subunits of human meprin

Meprins are zinc-endopeptidases of the astacin family, which are expressed as membrane-bound or secreted forms in renal and intestinal brush-border membranes of mouse, rat and man. There are two types of meprin subunits, alpha and beta, which form disulfide-bonded homo- and heterodimers; further oligomerization is mediated by non-covalent interactions. Both subunits are translated as proenzymes that have to be activated by removal of an N-terminal propeptide. In the gut, the most probable activator is trypsin. In addition, plasmin has been shown to activate the human alpha subunit in colorectal cancer tissue. In the present study we have overexpressed the human meprin alpha subunit and a His-tagged soluble tail-switch-mutant of meprin beta in Baculovirus-infected insect cells. The recombinant homo-oligomeric proteins were purified by gel filtration and affinity chromatography with yields of up to 10 mg/l cell culture medium and analyzed with regard to their activation mechanism. While both alpha and beta homo-oligomers are activated by trypsin, only meprin alpha homo-oligomers are processed to their mature form by plasmin. These results indicate a different accessibility of the propeptide in meprin homo-oligomers and suggest an explanation for the appearance of meprin hetero-oligomers consisting of active alpha, but latent beta subunits.     

C subunits binding to the protein kinase A RI alpha dimer induce a large conformational change

We present structural data on the RI alpha isoform of the cAMP-dependent protein kinase A that reveal, for the first time, a large scale conformational change within the RI alpha homodimer upon catalytic subunit binding. This result infers that the inhibition of catalytic subunit activity is not the result of a simple docking process but rather is a multi-step process involving local conformational changes both in the cAMP-binding domains as well as in the linker region of the regulatory subunit that impact the global structure of the regulatory homodimer. The results were obtained using small-angle neutron scattering with contrast variation and deuterium labeling. From these experiments we derived information on the shapes and dispositions of the catalytic subunits and regulatory homodimer within a holoenzyme reconstituted with a deuterated regulatory subunit. The scattering data also show that, despite extensive sequence homology between the isoforms, the overall structure of the type I alpha holoenzyme is significantly more compact than the type II alpha isoform. We present a model of the type I alpha holoenzyme, built using available high-resolution structures of the component subunits and domains, which best fits the neutron-scattering data. In this model, the type I alpha holoenzyme forms a flattened V shape with the RI alpha dimerization domain at the point of the V and the cAMP-binding domains of the RI alpha subunits with their bound catalytic subunits at the ends.     

Ligand recombination to the alpha and beta subunits of human hemoglobin

The rebinding of CO, O2, NO, methyl, ethyl, n-propyl, and n-butyl isocyanide to isolated alpha and beta chains and intact hemoglobin at pH 7, 20 degrees C was examined both during and after a 30-ns dye laser pulse. The resultant absorbance changes were analyzed in terms of a linear three-step reaction scheme: Hb + X in equilibrium with C in equilibrium with B in equilibrium with A or HbX, where A is the final bound state, and C and B are geminate states. Rate constants were assigned for each of the transitions in this mechanism using fitting procedures described previously for analyzing ligand rebinding to sperm whale myoglobin at room temperature (Gibson, Q. H., Olson, J. S., McKinnie, R. E., and Rohlfs, R. J. (1986) J. Biol. Chem. 261, 10228-10239). Five major conclusions were obtained. First, initial geminate recombination phases for the NO and O2 complexes of hemoglobin and its isolated subunits exhibit half-times equal to approximately 12 and approximately 440 ps, respectively. These values are in excellent agreement with more direct, picosecond measurements of the geminate recombination of HbNO (Cornelius, P. A., Hochstrasser, R. M., and Steele, A. W. (1983) J. Mol. Biol. 163, 119-128) and HbO2 (Friedman, J. M., Scott, T. W., Fisanick, G. J., Simon, S. R., Findsen, E. W., Ondrias, M. R., and MacDonald, V. W. (1985) Science 229, 187-229) following extremely short laser pulses. Second, the correspondence between our nanosecond measurements and the published picosecond data suggests strongly that the intrinsic photochemical yield of all ferrous, hexacoordinate heme complexes approaches one. Third, the major differences between the isolated alpha and beta chains involve the rate of ligand migration to the solvent, kC----X and the extent of recombination from the second geminate state, C, as measured by the ratio kC----B/kC----X. Fourth, for both isolated chains and intact hemoglobin, the rate and equilibrium constants for the formation of the initial O2 geminate state starting from ligand in the solvent (i.e. kX----B and KX----B) are 5-10 times greater than the corresponding parameters for the formation of the first CO geminate state. Fifth, the rate-limiting step for NO, O2, and isonitrile binding to hemoglobin and its isolated subunits is ligand migration up to the initial geminate state (i.e. kX----B). In the case of CO binding, both migration to state B and iron-ligand bond formation (kB----A) affect the overall, bimolecular association rate constant.     

Properties of the interspecies hybrids between haptoglobin alpha and beta subunits

1. Subunits alpha isolated from human haptoglobin were recombined with beta subunits of equine haptoglobin, and vice versa. Both hybrid proteins were separated on electrophoresis in polyacrylamide gel into four bands with mobilities corresponding to tetramers 2alpha.2beta, trimers 2alpha.beta, and dimers alpha.beta, in addition to free subunits beta. 2. The binding ability of haemoglobin and the antigenic specificity of tetramers depended on the origin of beta subunit. 3. Reduction of native and hybrid proteins with 2-mercaptoethanol led to gradual formation of alpha.beta, alpha, and beta; the components 2alpha.beta and 2alpha appeared in trace amounts.     

Site-specific mutagenesis of the alpha subunit of tryptophan synthase from Salmonella typhimurium. Changing arginine 179 to leucine alters the reciprocal transmission of substrate-induced conformational changes between the alpha and beta 2 subunits

Arginine 179 of the alpha subunit of tryptophan synthase of Salmonella typhimurium was changed to leucine by site-directed mutagenesis. The mutant alpha subunit was expressed in S. typhimurium, purified and crystallized as the alpha 2 beta 2 complex, and characterized by kinetic studies under steady-state reaction conditions. The rate of cleavage of indole 3-glycerol phosphate (alpha reaction) is reduced by 60% in the mutant alpha 2 beta 2 complex, whereas the rate of L-tryptophan synthesis from indole and L-serine (beta reaction) is unchanged. Thus, arginine 179 is not obligatory for catalysis, for binding of indole 3-glycerol phosphate, or for interaction of the alpha and beta 2 subunits. However, changing arginine 179 to leucine does have striking effects on ligand-dependent properties of this multienzyme complex. Ligands of the alpha subunit (DL-alpha-glycerophosphate and indole 3-propanol phosphate) which strongly inhibit the beta reaction of the native alpha 2 beta 2 complex have a slight stimulatory effect on the beta reaction of the mutant alpha 2 beta 2 complex. Likewise, L-serine, a ligand of the beta subunit which produces a 5-fold reduction in the Km for the alpha ligand indole 3-glycerol phosphate in the native alpha 2 beta 2 complex, has no effect on the mutant alpha 2 beta 2 complex. These results suggest that arginine 179 of the alpha subunit plays a role in the reciprocal transmission of substrate-induced conformational changes which occur between native alpha and beta 2 subunits in the alpha 2 beta 2 complex.     

Nanosecond-timescale conformational dynamics of the human alpha7 nicotinic acetylcholine receptor

We explore the conformational dynamics of a homology model of the human alpha7 nicotinic acetylcholine receptor using molecular dynamics simulation and analyses of root mean-square fluctuations, block partitioning of segmental motion, and principal component analysis. The results reveal flexible regions and concerted global motions of the subunits encompassing extracellular and transmembrane domains of the subunits. The most relevant motions comprise a bending, hinged at the beta10-M1 region, accompanied by concerted tilting of the M2 helices that widens the intracellular end of the channel. Despite the nanosecond timescale, the observations suggest that tilting of the M2 helices may initiate opening of the pore. The results also reveal direct coupling between a twisting motion of the extracellular domain and dynamic changes of M2. Covariance analysis of interresidue motions shows that this coupling arises through a network of residues within the Cys and M2-M3 loops where Phe135 is stabilized within a hydrophobic pocket formed by Leu270 and Ile271. The resulting concerted motion causes a downward shift of the M2 helices that disrupts a hydrophobic girdle formed by 9' and 13' residues.     

Role of disulfide bond formation in the folding of human chorionic gonadotropin beta subunit into an alpha beta dimer assembly-competent form.

The malignant trophoblastic cell line JAR was used as a model system to study protein folding in intact cells. We have used this model previously to identify conformational intermediates in the production of an assembly-competent form of the human chorionic gonadotropin beta subunit (Ruddon, R. W., Krzesicki, R. F., Norton, S. E., Beebe, J. S., Peters, B. P., and Perini, F. (1987) J. Biol. Chem. 262, 12533-12540). The earliest biosynthetic precursor of the human chorionic gonadotropin beta subunit detectable in JAR cells pulse labeled for 2 min is p beta 1, a form that lacks half of the six intrachain disulfide bonds observed in the fully processed dimer form of beta and that does not combine with the alpha subunit. p beta 1 is rapidly (t1/2 approximately 4 min) converted into p beta 2, which has a full complement of intrachain disulfide bonds and does combine with the alpha subunit. In this study, we have identified the three late forming disulfide bonds involved in the transition of p beta 1 into the assembly-compete form, p beta 2. The last three disulfide bonds to form are those between cysteines 9 and 90, 23 and 72, and 93 and 100. These were identified in JAR cell lysates that had been pulse labeled with [35S]cysteine for 2 or 5 min followed by trapping of the cysteine thiols with iodoacetic acid before immunopurification of the beta subunit forms. Immunopurified p beta 1 was treated with trypsin under nonreducing conditions to liberate [35S]cysteine-containing peptides from the disulfide-linked beta core polypeptide. These tryptic peptides were then separated by high performance liquid chromatography and sequenced to determine the location of the carboxymethyl-[35S]cysteine residues. The three late forming disulfide bonds are most likely the ones involved in stabilizing the conformation of the beta subunit that is required for combination with alpha to form the biologically functional alpha beta heterodimer.     

Molecular dynamics simulations and free energy analyses on the dimer formation of an amyloidogenic heptapeptide from human beta2-microglobulin: implication for the protofibril structure

Amyloid formation is associated with many neurodegenerative diseases. Recent findings suggest that early oligomeric aggregates could be major sources of toxicity. We present a computational investigation of the first step of amyloid initiation-dimer formation of a seven residue peptide (NHVTLSQ) from human beta2-microglobulin at pH 2.0, which renders +2.0 units charges to each peptide. A total of over 1.2 micros of simulations with explicit solvent and 1.0 micros of simulations with implicit solvent were conducted. Main-chain conformational restraint was applied to facilitate the formation of ordered dimers. An antiparallel beta-sheet with six main-chain hydrogen bonds was dominant in the implicit solvent simulations. In contrast, no stable dimers were observed in the two negative controls, the mouse heptapeptide (KHDSMAE, +3.0 units charges) and the scrambled human heptapeptide (QVLHTSN). Explicit solvent simulations presented a more complex scenario. The wild-type human heptapeptide formed predominantly antiparallel beta-sheets ( approximately 38%) although parallel ones ( approximately 12%) were also observed. Hydrophobic contacts preceded hydrogen bond saturation in the majority of the association events in the explicit solvent simulations, highlighting the important role of hydrophobic interaction in amyloid initiation. The fact that the mouse dimer dissociated immediately after the removal of conformational restraint suggests that the higher conformational entropy barrier, along with the stronger charge repulsion and weaker hydrophobic interaction, contributed to its inability to form amyloid fibril. The closeness of positive charge pairs in the dimers of the scrambled human heptapeptide may prohibit further beta-sheet extension and fibril growth. Combining the results from simulations and free energy analyses, we propose that the building block for this amyloid fibril is an antiparallel dimer with a two-residue register shift and six main-chain hydrogen bonds. A double-layer protofibril structure is also proposed in which two antiparallel beta-sheets face each other and are held together by hydrophobic staples and hydrogen bonds of the polar side-chains.     

Allosteric communication between alpha and beta subunits of tryptophan synthase: modelling the open-closed transition of the alpha subunit

Ligand binding to the alpha-subunit of the alpha2beta2 complex of tryptophan synthase induces the alphaloop6 closure over the alpha-active site. This conformational change is associated with the formation of a hydrogen bond between alphaGly181 NH group and betaSer178 carbonyl oxygen, a key event for the triggering of intersubunit allosteric signals. Mutation of betaSer178 to Pro and alphaGly181 to Pro, Ala, Phe and Val abolishes the ligand-induced intersubunit communication. Molecular dynamics methods were applied to simulate the conformation of the highly flexible and crystallographically undetectable open state of alphaloop6 in the wild type and in the alpha181 mutants. The open conformation of alphaloop6 is favoured in the wild type enzyme in the absence of alpha-ligands, and in the alpha181 mutants both in the presence and absence of bound ligands. A very good correlation was found between the extent of limited tryptic proteolysis and both the hydrogen bond distance between alphaX181 and betaSer178, obtained from the molecular dynamics simulation, and the hydrogen bond strength, evaluated by HINT, an empirical force field that takes into account both enthalpic and entropic contributions. Comparison of the open and closed conformations of alphaloop6 suggests a pathway for substrate entrance into the alpha-active site and provides an explanation for the limited catalytic efficiency of the open state.     

Packing-induced conformational and functional changes in the subunits of alpha -crystallin

The heteroaggregate alpha-crystallin and homoaggregates of its subunits, alphaA- and alphaB-crystallins, function like molecular chaperones and prevent the aggregation of several proteins. Although modulation of the chaperone-like activity of alpha-crystallin by both temperature and chaotropic agents has been demonstrated in vitro, the mechanism(s) of its regulation in vivo have not been elucidated. The subunits of alpha-crystallin exchange freely, resulting in its dynamic and variable quaternary structure. Mixed aggregates of the alpha-crystallins and other mammalian small heat shock proteins (sHSPs) have also been observed in vivo. We have investigated the time-dependent structural and functional changes during the course of heteroaggregate formation by the exchange of subunits between homoaggregates of alphaA- and alphaB-crystallins. Native isoelectric focusing was used to follow the time course of subunit exchange. Circular dichroism revealed large tertiary structural alterations in the subunits upon subunit exchange and packing into heteroaggregates, indicating specific homologous and heterologous interactions between the subunits. Subunit exchange also resulted in quaternary structural changes as demonstrated by gel filtration chromatography. Interestingly, we found time-dependent changes in chaperone-like activity against the dithiothreitol-induced aggregation of insulin, which correlated with subunit exchange and the resulting tertiary and quaternary structural changes. Heteroaggregates of varying subunit composition, as observed during eye lens epithelial cell differentiation, generated by subunit exchange displayed differential chaperone-like activity. It was possible to alter chaperone-like activity of preexisting oligomeric sHSPs by alteration of subunit composition by subunit exchange. Our results demonstrate that subunit exchange and the resulting structural and functional changes observed could constitute a mechanism of regulation of chaperone-like activity of alpha-crystallin (and possibly other mammalian sHSPs) in vivo.     

The interaction between subunits in the tubulin dimer.

Limited proteolysis and chemical cross-linking techniques have been used to study the interaction between alpha- and beta-tubulin subunits. Trypsin digestion of tubulin dimer resulted in the cleavage of the alpha-subunit into two fragments, whereas chymotrypsin cleaved the beta-subunit into two distinct fragments. All of these fragments have been mapped on the tubulin subunits by further proteolysis with formic acid. Cross-linking of trypsin- and chymotrypsin-cleaved subunits has been performed with two different cross-linker agents of different cross-linking distance. The addition of formaldehyde resulted in the cross-linking of the alpha-tubulin N-terminal fragment with beta-tubulin C-terminal domain. The same result was obtained when methyl 4-mercaptobutyrimidate was used.     

Antisera specific for the alpha 1, alpha 2, alpha 3, and beta subunits of the Na,K-ATPase: differential expression of alpha and beta subunits in rat tissue membranes

We have developed a panel of antibodies specific for the alpha 1, alpha 2, alpha 3, and beta subunits of the rat Na,K-ATPase. TrpE-alpha subunit isoform fusion proteins were used to generate three antisera, each of which reacted specifically with a distinct alpha subunit isotype. Western blot analysis of rat tissue microsomes revealed that alpha 1 subunits were expressed in all tissues while alpha 2 subunits were expressed in brain, heart, and lung. The alpha 3 subunit, a protein whose existence had been inferred from cDNA cloning, was expressed primarily in brain and copurified with ouabain-inhibitable Na,K-ATPase activity. An antiserum specific for the rat Na,K-ATPase beta subunit was generated from a TrpE-beta subunit fusion protein. Western blot analysis showed that beta subunits were present in kidney, brain, and heart. However, no beta subunits were detected in liver, lung, spleen, thymus, or lactating mammary gland. The distinct tissue distributions of alpha and beta subunits suggest that different members of the Na,K-ATPase family may have specialized functions.     

Reciprocal modulation between the alpha and beta 4 subunits of hSlo calcium-dependent potassium channels

Large conductance Ca(2+)-dependent potassium (K(Ca) or maxi K) channels are composed of a pore-forming alpha subunit and an auxiliary beta subunit. We have shown that the brain-specific beta4 subunit modulates the voltage dependence, activation kinetics, and toxin sensitivity of the hSlo channel (Weiger, T. M., Holmqvist, M. H., Levitan, I. B., Clark, F. T., Sprague, S., Huang, W. J., Ge, P., Wang, C., Lawson, D., Jurman, M. E., Glucksmann, M. A., Silos-Santiago, I., DiStefano, P. S., and Curtis, R. (2000) J. Neurosci. 20, 3563-3570). We investigated here the N-linked glycosylation of the beta4 subunit and its effect on the modulation of the hSlo alpha subunit. When expressed alone in HEK293 cells, the beta4 subunit runs as a single molecular weight band on an SDS gel. However, when coexpressed with the hSlo alpha subunit, the beta4 subunit appears as two different molecular weight bands. Enzymatic deglycosylation or mutation of the N-linked glycosylation residues in beta4 converts it to a single lower molecular weight band, even in the presence of the hSlo alpha subunit, suggesting that the beta4 subunit can be present as an immature, core glycosylated form and a mature, highly glycosylated form. Blockage of protein transport from the endoplasmic reticulum to the Golgi compartment with brefeldin A abolishes the mature, highly glycosylated beta4 band. Glycosylation of the beta4 subunit is not required for its binding to the hSlo channel alpha subunit. It also is not necessary for cell membrane targeting of the beta4 subunit, as demonstrated by surface biotinylation experiments. However, the double glycosylation site mutant beta4 (beta4 N53A/N90A) protects the channel less against toxin blockade, as compared with the hSlo channel coexpressed with wild type beta4 subunit. Taken together, these data show that the pore-forming alpha subunit of the hSlo channel promotes N-linked glycosylation of its auxiliary beta4 subunit, and this in turn influences the modulation of the channel by the beta4 subunit.     

Pressure dissociation and conformational drift of the beta dimer of tryptophan synthase

Micromolar solutions of tryptophan synthase beta 2 dimer dissociate into monomers in the pressure range of 800-1600 bars as shown by studies of the spectral shift of the intrinsic fluorescence and of the fluorescence polarization of dansyl conjugates. At 25 degrees C the standard change in volume on dissociation (dV0) of the holoprotein was -162 mL mol-1, and the dissociation constant at 1 bar was K0 = 3.7 10(-10) M. Pyridoxal-reduced holoprotein and apoprotein had, within 10%, the same dV0, but K0 was decreased in the reduced protein (6 X 10(-11) M) and increased in the apoprotein (3.6 X 10(-9) M). At 4 degrees C the free energy of association of the holoprotein was reduced by 1.4 kcal mol-1, but dV0 was unchanged. In all the protein forms the decompression curves differed from the respective compression curves, indicating the loss of some free energy of association following separation of the monomers. This hysteretic behavior was largest in the apoprotein and amounted to a loss of 2.6 kcal mol-1 in the free energy of association. When the pressure was rapidly raised to 2.2 kbars, half-dissociation of the reduced pyridoxal beta 2 dimer took approximately 12 min. Upon return to atmospheric pressure reassociation was complete in 2-3 min and half of the enzyme activity was regained in 10 min; pyridoxal fluorescence recovered more slowly with a biphasic course. The independent return of these properties and the hysteretic behavior indicate that subunit separation is followed by a conformational drift like that observed in lactate dehydrogenase dissociated by either pressure or temperature or in enolase dissociated by dilution.     

Chromosomal assignment of genes encoding the alpha, beta, and gamma subunits of human complement protein C8: identification of a close physical linkage between the alpha and the beta loci

The eighth component of human complement (C8) is a serum protein containing three nonidentical subunits (alpha, beta, gamma) that are arranged as a disulfide-linked alpha-gamma dimer and a noncovalently associated beta chain. In earlier genetic studies, electrophoretic analysis of C8 protein polymorphisms revealed several allelic variants of alpha-gamma and beta. These were governed by separate loci designated C8A and C8B for alpha-gamma and beta, respectively. Genetic linkage analyses indicated that these loci were linked to each other and to chromosome 1 marker loci PGM1 and Rh, but it was unclear at the time if C8A was a single locus coding for a single-chain precursor form of alpha-gamma or if separate loci existed for alpha and gamma. Since evidence now indicates that alpha, beta, and gamma are encoded by separate genes, cDNA probes corresponding to each subunit were used to make direct assignments of the individual loci. Analysis of somatic cell hybrids revealed that only the alpha and beta loci are located on chromosome 1. Parallel analysis of genomic DNA digests using 5' and 3'-specific cDNA probes showed they are physically linked (less than 2.5 kb) and oriented 5' alpha-beta 3'. Further probing of the hybrid panel revealed that gamma is located on chromosome 9q. Thus, the observed genetic linkage of alpha-gamma to beta must be determined solely by alpha. In accordance with these findings, the C8 loci should now be designated C8A, C8B, and C8G for alpha, beta and gamma, respectively.     

Biophysical characterization of interactions between the core binding factor alpha and beta subunits and DNA

Core binding factors (CBFs) play key roles in several developmental pathways and in human disease. CBFs consist of a DNA binding CBFalpha subunit and a non-DNA binding CBFbeta subunit that increases the affinity of CBFalpha for DNA. We performed sedimentation equilibrium analyses to unequivocally establish the stoichiometry of the CBFalpha:beta:DNA complex. Dissociation constants for all four equilibria involving the CBFalpha Runt domain, CBFbeta, and DNA were defined. Conformational changes associated with interactions between CBFalpha, CBFbeta, and DNA were monitored by nuclear magnetic resonance and circular dichroism spectroscopy. The data suggest that CBFbeta 'locks in' a high affinity DNA binding conformation of the CBFalpha Runt domain.