Lessons from functional and structural analyses of disease-associated genetic variants in the complement alternative pathway

Complement is an essential component of innate immunity and a major trigger of inflammatory responses. A critical step in complement activation is the formation of the C3 convertase of the alternative pathway (AP), a labile bimolecular complex formed by activated fragments of the C3 and factor B components that is fundamental to provide exponential amplification of the initial complement trigger. Regulation of the AP C3 convertase is essential to maintain complement homeostasis in plasma and to protect host cells and tissues from damage by complement. During the last decade, several studies have associated genetic variations in components and regulators of the AP C3 convertase with a number of chronic inflammatory diseases and susceptibility to infection. The functional characterization of these protein variants has helped to decipher the critical pathogenic mechanisms involved in some of these complement related disorders. In addition, these functional data together with recent 3D structures of the AP C3 convertase have provided fundamental insights into the assembly, activation and regulation of the AP C3 convertase.     

Anion shielding of electrostatic repulsions in transthyretin modulates stability and amyloidosis: insight into the chaotrope unfolding dichotomy

The balance between stabilizing forces and the localized electrostatic repulsions destabilizing the transthyretin (TTR) tetramer is tunable via anion shielding. The two symmetrical anion interaction sites in TTR are comprised of residues Lys15 and Lys15' from opposing subunits on the periphery of the two thyroxine binding sites. These epsilon-ammonium groups repel one another and destabilize the tetramer, unless an appropriate anion is present, which stabilizes the tetramer. Chaotrope denaturation of TTR exhibits unusual behavior in that urea appears to be a stronger denaturant than GdmCl (guanidinium chloride), even though GdmCl is typically twice as powerful as a denaturant. The shift in the midpoint of the urea denaturation curve to higher concentrations as well as the increase in the mole fraction of tetramer that is highly resistant to denaturation with increasing KCl concentration provides strong evidence that anion shielding stabilizes the TTR tetramer. A consequence of tetramer stabilization is folding hysteresis, because the high GdmCl concentrations required to denature the anion-stabilized tetramer do not allow refolding of the unfolded monomers. The formation of amyloid fibrils by TTR requires that its normal tetrameric structure dissociate to alternatively folded monomers, a process mediated by acidification (pH 5-4). This process is inhibited by Cl(-) ions in a concentration-dependent fashion. Chloride ion may not be the relevant physiological TTR stability modulator, but it is the main focus of these studies explaining the hysteresis observed in the denaturation and refolding studies with GdmCl.     

Tertiary and quaternary structural differences between two genetic variants of bovine casein by small-angle X-ray scattering

The casein complexes of bovine milk consist of four major protein fractions, alpha s1, alpha s2, beta, and kappa. Colloidal particles of casein (termed micelles) contain inorganic calcium and phosphate; they are very roughly spherical with an average radius of 650 A. Removal of Ca2+ leads to the formation of smaller protein aggregates (submicelles) with an average radius of 94 A. Two genetic variants, A and B, of the predominant fraction, alpha s1-casein, result in milks with markedly different physical properties, such as solubility and heat stability. To investigate the molecular basis for these differences, small-angle X-ray scattering was performed on the respective colloidal micelles and submicelles. Scattering curves for submicelles of both variants showed multiple Gaussian character; data for the B variant were previously interpreted in terms of two concentric regions of different electron density, i.e., a "compact" core and a relatively "loose" shell. For the submicelle of A, there was a third Gaussian, reflecting a negative contribution due to interparticle interference. Molecular parameters for submicelles of both A and B are in agreement with hydrodynamic data in the literature. Data for the micelles, for which scattering yields cross-sectional information, were fitted by a sum of three Gaussians for both variants; for these, the corresponding two lower radii of gyration represent the two concentric regions of the submicelles, while the third reflects the average packing of submicelles within the micellar cross section. Most of the molecular parameters obtained showed small but consistent differences between A and B, but for submicelles within the micelle several differences were particularly notable: A has a greater molecular weight for the "compact" region of the constituent submicelle (82,000 vs 60,000) and a much greater submicellar packing number (6:1 vs 3:1). Reasons for these and other differences are to be sought in sequence differences and in differences in calcium-binding sites and charge distribution.     

Mass spectrometric immunoassay for quantitative determination of transthyretin and its variants

Transthyretin (TTR, or prealbumin) is a tetrameric protein found in plasma and cerebrospinal fluid. Its major role is to transport thyroid hormones (thyroxin-T4) and retinol (through association with retinol-binding protein). TTR has been studied extensively due to the great number of point mutations that result in sequence heterogeneity. Many of these variants are associated with pathological conditions that result in extracellular deposition of amyloid fibers in tissues. In this work, we have developed a rapid mass spectrometric immunoassay for determination and quantification of TTR and its variants from human serum and plasma samples. The assay was fully characterized in terms of its precision, linearity and recovery characteristics. The new assay was also compared with a conventional TTR ELISA. Furthermore, we have applied the optimized method to analyze TTR and its modifications in 44 human plasma samples, and in the process optimized a method for TTR proteolytic digestion and identification of point mutations.     

Role of the tertiary and quaternary structures in the stability of dimeric copper, zinc superoxide dismutases

The equilibrium unfolding process of human Cu,Zn superoxide dismutase has been quantitatively monitored through circular dichroism and fluorescence spectroscopy as a function of increasing guanidinium hydrochloride concentration. The process occurs through the formation of a monomeric intermediate species following a three-state transition equilibrium. Comparison with the stability of the prokaryotic Cu,Zn SOD from P. leiognathi shows that the eukaryotic enzyme is more stable than the prokaryotic enzyme by approximately 3 kcal/mol. This difference is due to the monomer-to-unfolded equilibrium, while the dimer-to-monomer equilibrium is comparable for the two enzymes despite their different intersubunit interactions. These results are confirmed by the unfolding of the copper-depleted derivatives. The Cu,Zn superoxide dismutase represents a good example of how evolution has found two independent quaternary assemblies maintaining the same dimer stability.     

Quantification of tertiary structural conservation despite primary sequence drift in the globin fold.

The globin family of protein structures was the first for which it was recognized that tertiary structure can be highly conserved even when primary sequences have diverged to a virtually undetectable level of similarity. This principle of structural inertia in molecular evolution is now evident for many other protein families. We have performed a systematic comparison of the sequences and structures of 6 representative hemoglobin subunits as diverse in origin as plants, clams, and humans. Our analysis is based on a 97-residue helical core in common to all 6 structures. Amino acid sequence identities range from 12.4% to 42.3% in pairwise comparisons, and, despite these variations, the maximal RMS deviation in alpha-carbon positions is 3.02 A. Overall, sequence similarity and structural deviation are significantly anticorrelated, with a correlation coefficient of -0.71, but for a set of structures having under 20% pairwise identity, this anticorrelation falls to -0.38, which emphasizes the weak connection between a specific sequence and the tertiary fold. There is substantial variability in structure outside the helical core, and functional characteristics of these globins also differ appreciably. Nevertheless, despite variations in detail that the sequence dissimilarities and functional differences imply, the core structures of these globins remain remarkably preserved.     

The ionization of a buried glutamic acid is thermodynamically linked to the stability of Leishmania mexicana triose phosphate isomerase.

The amino acid sequence of Leishmania mexicana triose phosphate isomerase is unique in having at position 65 a glutamic acid instead of a glutamine. The stability properties of LmTIM and the E65Q mutant were investigated by pH and guanidinium chloride-induced unfolding. The crystal structure of E65Q was determined. Three important observations were made: (a) there are no structural rearrangements as the result of the substitution; (b) the mutant is more stable than the wild-type; and (c) the stability of the wild-type enzyme shows strong pH dependence, which can be attributed to the ionization of Glu65. Burying of the Glu65 side chain in the uncharged environment of the dimer interface results in a shift in pKa of more than 3 units. The pH-dependent decrease in overall stability is due to weakening of the monomer-monomer interactions (in the dimer). The E65Q substitution causes an increase in stability as the result of the formation of an additional hydrogen bond in each subunit (DeltaDeltaG degrees of 2 kcal.mol-1 per monomer) and the elimination of a charged group in the dimer interface (DeltaDeltaG degrees of at least 9 kcal.mol-1 per dimer). The computated shift in pKa and the stability of the dimer calculated from the charge distribution in the protein structure agree closely with the experimental results. The guanidinium chloride dependence of the unfolding constant was smaller than expected from studies involving monomeric model proteins. No intermediates could be identified in the unfolding equilibrium by combining fluorescence and CD measurements. Study of a stable monomeric triose phosphate isomerase variant confirmed that the phenomenon persists in the monomer.     

Functional and structural stability are linked in phytoplankton metacommunities of different connectivity

The diversity–stability debate is a long‐standing issue in ecology, asking whether more diverse communities show higher stability over time and more rapid recovery from disturbances. Connection to undisturbed habitats is thought to affect compositional and functional stability after disturbances. Therefore, we established marine phytoplankton metacommunities consisting of three microcosms (local patches), which were connected by tubes opened for different time intervals to create 5 levels of connectivity. We performed two experiments differing by homogeneous (HOM) or heterogeneous (HET) supply of irradiance across patches. As disturbance we either removed 75% of the algal biomass locally from one randomly chosen patch, or 25 or 75% regionally from each local patch. By comparing these treatments to an undisturbed control, we analyzed resilience (rate of recovery) and final recovery (recovery ratio) with regard to biomass (functional stability) and species composition (structural stability). In both experiments (HET, HOM), functional and structural aspects of stability responded significantly to connectivity and disturbance treatments. Functional resilience was enhanced by increasing connectivity (HET and HOM), which partially also increased functional recovery (HET) and structural resilience (HOM). By contrast, the treatment‐induced gradients in diversity (species richness and evenness) had no clear effect on functional resilience or recovery. Instead, structural and functional resilience were strongly correlated regardless of patch quality, indicating that only a full recovery in community composition ensured functional stability. Our findings suggest that connectivity plays a pivotal role in maintaining ecosystem stability under pulse disturbance such that a more complete understanding of stability requires spatially explicit approaches.     

The relationship of the quaternary structure of allophycocyanin to its spectrum

Allophycocyanin II in its trimer form (α₃β₃) at pH 7.0 has an absorption maximum at 652 nm. This band is selectively reduced in intensity at pH 7.0 when various salts are added. The loss of 652 nm absorption follows the order: NaClO₄ ? NaNO₃ > NaBr > NaCl. When the NaClO₄ concentration is in the range 0.6-1.0 m the 652-nm band is entirely lost, and sedimentation equilibrium and velocity studies suggest that the trimer is completely dissociated to monomers (αβ). Hydrophobic interactions appear to be important in maintaining the trimer. The monomer absorption maximum is at 616 nm. A series of experiments using these salts demonstrated at intermediate 652-nm intensities and the two extrema that an isobestic point at 626 nm is present which indicates an equilibrium between two species. Corresponding to the loss of 652 nm absorption is the disappearance of 661 nm fluorescence emission and the appearance of a new band at 642 nm. Removal of the NaClO₄ by dialysis essentially restores the 652-nm absorption and 661-nm emission and the trimeric protein structure. The near ultraviolet region is only slightly perturbed during the loss of 652 nm absorption. In the absence of any additional salts these spectral changes also occur in pH 7.0 buffer at very low protein concentrations.     

Binding affects the tertiary and quaternary structures of the Shigella translocator protein IpaB and its chaperone IpgC

Shigella flexneri uses its type III secretion system (T3SS) to promote invasion of human intestinal epithelial cells as the first step in causing shigellosis, a life-threatening form of dysentery. The Shigella type III secretion apparatus (T3SA) consists of a basal body that spans the bacterial envelope and an exposed needle that injects effector proteins into target cells. The nascent Shigella T3SA needle is topped with a pentamer of the needle tip protein invasion plasmid antigen D (IpaD). Bile salts trigger recruitment of the first hydrophobic translocator protein, IpaB, to the tip complex where it senses contact with a host membrane. In the bacterial cytoplasm, IpaB exists in a complex with its chaperone IpgC. Several structures of IpgC have been determined, and we recently reported the 2.1 ? crystal structure of the N-terminal domain (IpaB(74.224)) of IpaB. Like IpgC, the IpaB N-terminal domain exists as a homodimer in solution. We now report that when the two are mixed, these homodimers dissociate and form heterodimers having a nanomolar dissociation constant. This is consistent with the equivalent complexes copurified after they had been co-expressed in Escherichia coli. Fluorescence data presented here also indicate that the N-terminal domain of IpaB possesses two regions that appear to contribute additively to chaperone binding. It is also likely that the N-terminus of IpaB adopts an alternative conformation as a result of chaperone binding. The importance of these findings within the functional context of these proteins is discussed.     

Legume lectin family, the 'natural mutants of the quaternary state', provide insights into the relationship between protein stability and oligomerization

Legume lectins family of proteins, despite having the same 'jelly roll' tertiary structural fold at monomeric level, exhibit considerable variation in their quaternary structure arising out of small changes in their sequence. Nevertheless, their folding behavior and stability correlates very well with their patterns of assembly into dimers and tetramers. A conservation of their fold during evolution, its wide distribution in many protein families together with the availability of structural information on them make them interesting as proteins to explore the effect of inter- versus intra-subunit interactions in the stability of multimeric proteins. Additionally, as 'natural mutants' of quaternary association, proteins of legume lectin family provide interesting paradigms for studies addressing the effect of subunit oligomerization on the stability, folding and function as well as the evolution of multimeric structures.     

Computation of conformational entropy from protein sequences using the machine-learning method--application to the study of the relationship between structural conservation and local structural stability

A complete protein sequence can usually determine a unique conformation; however, the situation is different for shorter subsequences--some of them are able to adopt unique conformations, independent of context; while others assume diverse conformations in different contexts. The conformations of subsequences are determined by the interplay between local and nonlocal interactions. A quantitative measure of such structural conservation or variability will be useful in the understanding of the sequence-structure relationship. In this report, we developed an approach using the support vector machine method to compute the conformational variability directly from sequences, which is referred to as the sequence structural entropy. As a practical application, we studied the relationship between sequence structural entropy and the hydrogen exchange for a set of well-studied proteins. We found that the slowest exchange cores usually comprise amino acids of the lowest sequence structural entropy. Our results indicate that structural conservation is closely related to the local structural stability. This relationship may have interesting implications in the protein folding processes, and may be useful in the study of the sequence-structure relationship.     

The relationship between molecular weight and quaternary structure of globular proteins

The relationship between the molecular weight and the number of subunits in oligomeric globular proteins consisting of identical subunits has been analyzed. It has been shown that the molecular weights of the subunits are distributed about a mean value of 48,000 and consequently that the molecular weights of the native oligomeric proteins are distributed in clearly distinguishable molecular weight regions. This observation allows the probability of a particular oligomeric structure to be predicted from a measurement of the oligomer molecular weight alone, which is useful in a number of types of study of protein structure, particularly comparative studies. Calculations have been performed which suggest that there is no thermodynamic limitation, in terms of the subunit interactions themselves, to the size of an oligomeric protein with a given number of subunits. Rather, an individual polypeptide chain itself has inherent size limitations, which consequently limits the molecular weight of the corresponding oligomer.     

Use of analytical gel chromatography to analyze tertiary and quaternary structural changes in E. coli aspartate transcarbamylase

E. coli aspartate transcarbamylase (ATCase) is a large (310 kDa) protein that undergoes major changes in quaternary structure when substrates and regulatory nucleotides bind. We have used analytical gel chromatography to detect quaternary structure changes in both the holoenzyme and its catalytic subunit (c3), to characterize the quaternary structure of single site mutant proteins and to monitor urea-induced dissociation and unfolding of c3. Binding of the bisubstrate analog PALA (N-(phosphonacetyl)-L-aspartate) to ATCase and c3 has been shown to alter s20.w by -3.3% and + 1.4%, respectively [Howlett, G.J. and Schachman, H.K. (1977), Biochemistry 23, 5077-5083]. The corresponding changes in the chromatographic partition coefficient (sigma) are -2.6 +/- 0.3% and 5.5 +/- 1.9% on Sephacryl S400HR and S200, respectively. Partition coefficients of mutant ATCases with single site mutations in the c chain differ from those of the wild-type protein by +/- 0.5% in small zone experiments; for example, mutations Arg 269----Gly and Glu 239----Gln alter the partition coefficient by 0.4% and -0.5%, respectively. The partition coefficient of mutant Glu 50----Gln is identical to the wild type enzyme. In the presence of saturating PALA, partition coefficients of Glu 50----Gln and Arg 269----Gly, but not Glu 239----Gln are identical to those of the wild type. Results for Glu 239----Gln are consistent with measurements of activity, small angle X-ray scattering and sedimentation coefficient that indicate that mutations at this site shift the quaternary structure towards the R state [Ladjimi and Kantrowitz (1988), Biochemistry 27, 276-83; Vachette and Hervé, cited by Kantrowitz and Lipscomb (1988), Science 241, 669-674; Newell and Schachman (1988), FASEB J. 2, A551]. Results for Glu 50----Gln are also consistent with measurements of activity (Ladjimi et al. (1988), Biochemistry 27, 268-276). The changes in tertiary and quaternary structure that result from urea-induced denaturation of c3 result in larger changes in the partition coefficient. Dissociation into folded monomers in 1-1.75 M urea is accompanied by a 4.6% increase in partition coefficient, while denaturation at greater than 5 M urea gives rise to a 43% decrease on S-300 Sephacryl. The bisubstrate analog PALA suppresses dissociation and increases the cooperativity of the unfolding reaction.     

Correlation between the structural stability and aggregation propensity of proteins

Protein aggregation, being an outcome of improper protein folding, is largely dependent on the folding kinetics of a protein. Previous studies have reported a positive correlation between the stability of the secondary structural elements of a protein and their rate of folding/unfolding. In this in silico study, the secondary and tertiary structures of proteins a) that form inclusion bodies on overexpression in Escherichia coli, b) that form amyloid fibrils and c) that are soluble on overexpression in E. coli are analyzed for certain features that are known to be associated with structural stability. The study revealed that the soluble proteins seem to have a higher rate of folding (based on contact order) and a lower percentage of exposed hydrophobic residues as compared to the inclusion body forming or amyloidogenic proteins. The soluble proteins also seem to have a more favored helix and strand composition (based on the known secondary structural propensities of amino acids). The secondary structure analyses also reveal that the evolutionary pressure is directed against protein aggregation. This understanding of the positive correlation between structural stability and solubility, along with the other parameters known to influence aggregation, could be exploited in the design of mutations aimed at reducing the aggregation propensity of the proteins.     

Intrinsic versus mutation dependent instability/flexibility: a comparative analysis of the structure and dynamics of wild-type transthyretin and its pathogenic variants

Transthyretin (TTR) is one of the about 20 known human proteins associated with amyloidosis which is characterized by the accumulation of amyloid fibrils in tissues or extracellular matrix surrounding vital organs. Unlike Alzheimer's fibrils that comprise a fragment of a large precursor protein, TTR amyloid fibrils are composed of both full-length protein and fragments of the molecule. The native state of TTR is a homotetramer with eight beta-strands organized into a beta-sandwich in each monomer. To elucidate the structural reorganization mechanisms preceding amyloid formation, it is important to characterize the dynamic features of the wild-type native state as well as to reveal the influence of disease-associated mutations on the structure and dynamics. Molecular dynamics (MD) simulations complement X-ray crystallography and D-H exchange to capture the intrinsically unstable/flexible sites of the wild-type as well as the mutation dependent unstable sites of the pathogenic variants. Our results of MD simulations have shown that the Leu55-->Pro (L55P) mutation occurs in an intrinsically unstable site, leading to substantial local and global structural changes. This observation supports the early speculation that the C-strand-loop-D-strand rearrangement leads to the formation of amyloidogenic intermediates. In addition to the D strand, the alpha-helical region and the strands at the monomer-monomer interface are also intrinsically unstable. The central channel of L55P-TTR undergoes opening and closing fluctuations, which may provide an explanation for the fact that while the mutation is far from the channel, the mutant shows a substantial low binding affinity of thyroxine.     

Involvement of single residue tryptophan 548 in the quaternary structural stability of pigeon cytosolic malic enzyme

Pigeon cytosolic malic enzyme has a double dimer quaternary structure with three tryptophanyl residues in each monomer distributed in different structural domains. The enzyme showed a three-state unfolding phenomenon upon increasing the urea concentration (Chang, H. C., Chou, W. Y., and Chang, G. G. (2002) J. Biol. Chem. 277, 4663-4671). At urea concentration of 4-4.5 m, where the intermediate form was detected, the enzyme existed as partially unfolded dimers, which were easily polymerized. Mn2+ provided full protection against the polymerization. To further characterize this phenomenon, three mutants of the enzyme (W129, W321, and W548), each with only one tryptophanyl residue left, were constructed. All these mutants were successfully overexpressed in Escherichia coli cells and purified to homogeneity. Changes in the circular dichroism spectra of all mutants revealed a three-state urea-unfolding process in the absence of Mn2+. In the presence of 4 mm Mn2+, W548 and wild type (WT) enzymes shifted to monophasic, while W129 and W321 were still biphasic. Similar results were obtained from the fluorescence spectral changes, except for W321, which showed monophasic denaturation curve with or without Mn2+. Analytical ultracentrifugation analysis indicated that the mutant enzymes were polymerized at 4.5 m urea, and Mn2+ provided protective effect on W548 and WT enzymes only. Other mutants with mutated Trp-548 polymerized at 4.5 m urea in the absence or presence of 4 mm Mn2+. The above results indicate that a single residue, Trp-548, in the subunit interface region, is responsible for the integrity of the quaternary structure of the pigeon cytosolic malic enzyme.