Subscriber access provided by UNIVERSITY OF ADELAIDE LIBRARIES
Article
Pathogenic Mutations Induce Partial Structural Changes in Native #-Sheet Structure of Transthyretin and Accelerate Aggregation Kwang Hun Lim, Anvesh K. R. Dasari, Renze Ma, Ivan Hung, Zhehong Gan, Jeffery W. Kelly, and Michael C Fitzgerald Biochemistry, Just Accepted Manuscript • DOI: 10.1021/acs.biochem.7b00658 • Publication Date (Web): 18 Aug 2017 Downloaded from http://pubs.acs.org on August 19, 2017
Just Accepted “Just Accepted” manuscripts have been peer-reviewed and accepted for publication. They are posted online prior to technical editing, formatting for publication and author proofing. The American Chemical Society provides “Just Accepted” as a free service to the research community to expedite the dissemination of scientific material as soon as possible after acceptance. “Just Accepted” manuscripts appear in full in PDF format accompanied by an HTML abstract. “Just Accepted” manuscripts have been fully peer reviewed, but should not be considered the official version of record. They are accessible to all readers and citable by the Digital Object Identifier (DOI®). “Just Accepted” is an optional service offered to authors. Therefore, the “Just Accepted” Web site may not include all articles that will be published in the journal. After a manuscript is technically edited and formatted, it will be removed from the “Just Accepted” Web site and published as an ASAP article. Note that technical editing may introduce minor changes to the manuscript text and/or graphics which could affect content, and all legal disclaimers and ethical guidelines that apply to the journal pertain. ACS cannot be held responsible for errors or consequences arising from the use of information contained in these “Just Accepted” manuscripts.
Biochemistry is published by the American Chemical Society. 1155 Sixteenth Street N.W., Washington, DC 20036 Published by American Chemical Society. Copyright © American Chemical Society. However, no copyright claim is made to original U.S. Government works, or works produced by employees of any Commonwealth realm Crown government in the course of their duties.
Page 1 of 36
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60
Biochemistry
Pathogenic Mutations Induce Partial Structural Changes
in
Native
β-Sheet
Structure
of
Transthyretin and Accelerate Aggregation Kwang Hun Lim*,†, Anvesh K. R. Dasari†, Renze Ma , Ivan Hung‡, Zhehong Gan‡, Jeffery W. Kelly§, Michael C. Fitzgerald †
Department of Chemistry, East Carolina University, Greenville, NC 27858, USA. Department
of Chemistry, Duke University, 124 Science Drive, Durham, NC 27708-0346, USA. ‡ Center of Interdisciplinary Magnetic Resonance (CIMAR), National High Magnetic Field Laboratory (NHMFL), 1800 East, Paul Dirac Dr., Tallahassee, FL 32310, USA. § Department of Molecular and Experimental Medicine and the Skaggs Institute for Chemical Biology, The Scripps Research Institute, La Jolla, CA 92037, USA.
Corresponding author: Kwang Hun Lim, Department of Chemistry, East Carolina University, Greenville, NC 27858, USA. E-mail:
[email protected]; (T) 252-328-9805 KEYWORDS TTR, Amyloid, Solid-state NMR, PDSD, amyloidosis, misfolding, V30M, L55P, SPROX, mass spectrometry
ACS Paragon Plus Environment
1
Biochemistry
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60
Page 2 of 36
ABBREVIATIONS NMR, nuclear magnetic resonance; PDSD, proton-driven spin-diffusion; TTR, transthyretin; CP, cross-polarization; WT, wild-type; CD, circular dichroism; SPROX, stability of proteins from rates of oxidation; Leu, leucine; Val, Valine; Gly, glycine; Tyr, tyrosine; MAS, magic-angle spinning; HNSB, dimethyl(2-hydroxy-5-nitrobenzyl)sulfonium bromide; GdmCl, guanidinium chloride; H/D exchange, Hydrogen/deuterium (H/D) exchange.
ACS Paragon Plus Environment
2
Page 3 of 36
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60
Biochemistry
Abstract Amyloid formation of natively folded proteins involves global and/or local unfolding of the native state to form aggregation-prone intermediates. Here we report solid-state NMR structural studies of amyloid derived from wild-type (WT) and more aggressive mutant forms of transthyretin (TTR) to investigate the structural changes associated with effective TTR aggregation. We employed selective 13C-labeling schemes to investigate structural features of βstructured core regions in amyloid states of WT and two mutant forms (V30M and L55P) of TTR. Analyses of the
13
C-13C correlation solid-state NMR spectra revealed that WT TTR
aggregates contain an amyloid core consisting of native-like CBEF and DAGH β-sheet structures and the mutant TTR amyloids adopt a similar amyloid core structure with native-like CBEF and AGH β-structures. However, the V30M mutant amyloid was shown to have a different DA βstructure. In addition, strand D is more disordered even in the native state of L55P TTR, indicating that the pathogenic mutations affect the DA β-structure, leading to more effective amyloid formation. The NMR results are consistent with our mass spectrometry-based thermodynamic analyses that showed the amyloidogenic precursor states of WT and mutant TTRs adopt folded structures, but the mutant precursor states are less stable than that of WT TTR. Analyses of the oxidation rate of methionine sidechain also revealed that the sidechain of residue Met-30 pointing between strands D and A is not protected from the oxidation in V30M mutant, while protected in the native state, supporting that the DA β-structure might be disrupted in V30M mutant amyloid.
ACS Paragon Plus Environment
3
Biochemistry
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60
Page 4 of 36
Introduction Transthyretin (TTR) is one of the few globular proteins that undergo aberrant conformational transitions from natively folded states into insoluble β-structured aggregates in vivo.1-3 Extracelluar deposition of TTR aggregates in various organs including heart, lung, and peripheral nerves is associated with senile systemic amyloidosis, familial amyloidotic polyneuropathy, familial amyloid cardiomyopathy, and rarely central nervous system selective amyloidosis. 4-6 The TTR amyloidoses exhibit extreme variations of the disease phenotype.4,7-9 For example, aggregation of wild-type (WT) TTR affects primarily the heart and lung, causing senile systemic amyloidosis that affects nearly 25% of the population over age 80.10 Various single-point mutations including V30M and L55P cause exclusively neurological disorders (polyneuropathy). On the other hand, many of the mutations (T60A, I68L, L111M and V122I) are associated with mainly cardiac phenotype, while some of the mutations develop mixed symptoms.11 Such heterogeneity suggests that mutant forms of TTR may have different misfolding pathways and adopt distinct amyloid conformations.12 TTR is a 55 kDa homo-tetrameric protein that binds and transports holoretinol binding protein in the bloodstream. TTR is also the major carrier of thyroxine (T4) in the plasma and cerebrospinal fluid.13 The native state of the 127-residue TTR monomer adopts a primarily βsheet tertiary structure with two four-β-stranded anti-parallel sheets (CBEF and DAGH) that are arranged into a β-sandwich (Fig. 1).14-16 The TTR monomers form a dimer through anti-parallel β-sheet interactions between the H and H’ strands, and between the strands F and F’. Association of two dimers mainly by hydrophobic interactions between the residues in AB and GH loops leads to the formation of tetramer.
ACS Paragon Plus Environment
4
Page 5 of 36
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60
Biochemistry
Figure 1. A schematic diagram of the misfolding and amyloid formation mechanism of TTR.
Previous biophysical studies showed that dissociation of the tetrameric TTR to monomers is a rate-determining step in TTR amyloidogenesis.17-19 The dissociated monomers may then undergo a conformational transition to a partially unfolded intermediate, which self-assembles into amyloid (Figure 1).17-19 It was also shown that most pathogenic single-point mutations destabilize the native TTR tetramer, subsequently facilitating the dissociation of the tetramer and amyloid formation.8,9 However, structural studies using X-ray crystallography revealed that the pathogenic mutations do not induce major changes in the β-sheet tertiary structure that might be responsible for the enhanced aggregation propensity, 20,21 and thereby the precise effect of the mutations on the TTR aggregation remains elusive. The single-point mutations may accelerate only the tetramer dissociation to monomers, and the mutant TTR monomer may have the same misfolding pathway as the WT TTR, which will lead to identical structures of the end product, amyloid. Local structural perturbations by the mutations may also initiate a distinct misfolding transition, resulting in different amyloid structures relative to WT TTR.12 Structural investigation of the soluble precursor states and the final product amyloid derived from WT and mutant forms of TTR would be essential to gain insights into the misfolding pathways. Previous structural studies of TTR aggregates using various biophysical techniques revealed that TTR amyloids contain extensive native-like β-structures. However, the extent of the native-
ACS Paragon Plus Environment
5
Biochemistry
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60
Page 6 of 36
like β-structures has been controversial.22-28 Structural investigations of TTR amyloid using electron paramagnetic resonance (EPR) spectroscopy with site-directed spin labeling suggested that strands C and D dislocate from the main β-sheet structures, exposing strands B and A, respectively, for intermolecular association to amyloid.25,29,30 Hydrogen/deuterium (H/D) exchange NMR experiments on Y114C mutant amyloid also indicated that strands C and D fold away from the β-sheet structure on the basis of relatively lower solvent protection factors of the strands C and D regions.22 However, recent protease digestion experiments suggested a considerably different structure for WT amyloid. It was proposed that strand D in the DAGH βsheet makes contacts with strand E in the amyloid core, implying that the native β-sheet is completely disrupted and rearranged into non-native β-sheet structures in amyloid. 31 Very recently, we reported solid-state NMR structural studies of wild-type (WT) TTR amyloid formed at pH 4.4, revealing that the two native β-sheets (CBEF and DAGH) remain largely unaltered after amyloid formation. 32 We also demonstrated that 13C-13C correlation solidstate NMR experiments with selective
13
CO and
13
Cα labeling schemes are effective for
exploring native-like β-sheet structures in amyloid. 32,33 Here we extended our solid-state NMR structural studies to mutant forms of TTR (V30M and L55P) to investigate the effect of the pathogenic mutations on TTR amyloid structure. Our results reveal that the mutant forms of TTR amyloid have similar native-like structural features (CBEF and AGH) as WT amyloid, which is consistent with our circular dichroism (CD) and mass spectrometry-based thermodynamic analyses of the amyloidogenic precursor states of WT and TTR variants. However, the DA βstructure appears more disordered in the two mutant TTR amyloids, which may be linked to more effective amyloid formation of the mutant TTRs.
ACS Paragon Plus Environment
6
Page 7 of 36
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60
Biochemistry
Materials and Methods Protein Expression and Purification Recombinant wild type and mutant forms of TTR were expressed and purified from an Escherichia coli cells transformed from the pMMHa plasmid as described previously.34 Briefly, overexpressed proteins in the soluble fraction of the sonicated lysates were extracted using ammonium sulfate precipitation methods. The precipitates were resuspended with deionized water and purified using anion exchange Q column, followed by size exclusion gel-filtration chromatography using Superdex HR 200 column (GE Health care). The expression and purification methods yielded a large quantity of TTR proteins (~ 80 mg from a 1L M9 medium culture). The protein concentration was calculated using an extinction coefficient of 7.76 × 104 M-1 cm-1 at 280 nm. Amyloid samples were obtained by incubating the protein (0.2 mg/mL) in 200 mM acetate buffer (100 mM KCl, 1 mM EDTA, pH 4.4) for a period of 30 days at 37 °C. The insoluble amyloid was spun down and washed twice with deionized water to remove remaining tetramers and soluble aggregates. The TTR amyloid was examined by transmission electron microscopy (TEM) and thioflavin T (ThT) binding assay in our previous studies.33
CD spectroscopy The CD spectra were recorded by scanning from 250 nm to 190 nm on a Jasco J-815 spectropolarimeter (Easton, MD) using 1 and 3 mm path length Suprasil quartz cells depending on the protein concentration. Very low protein concentrations of 0.09 and 0.02 mg/ml in 15 mM sodium acetate buffer (pH 4.4) were used for the CD experiments.
ACS Paragon Plus Environment
7
Biochemistry
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60
Page 8 of 36
Solid-state NMR Solid-state NMR spectra were recorded using Bruker 600- and 830-MHz spectrometers equipped with a 3.2 mm MAS probe. Two-dimensional
13
C-13C correlation NMR spectra were
recorded using a proton-driven spin diffusion (PDSD) mixing scheme at spinning frequencies of 17 kHz (600 MHz) and 11 – 12 kHz (830 MHz). The spinning speed was set close to the ∆ωiso = ωr (600 MHz) and ∆ωiso = 2ωr rotational resonance (RR) conditions 35,36 for efficient polarization transfer. Since the second spinning sideband from 13CO carbons is still observable under the RR conditions, the spinning speed was set slightly below for each 13CO/13Cα spin pairs such that the sideband of the carbonyl carbon appears ~ 5 ppm away from the NMR resonance of the carbon. The 90° pulse-lengths for 1H and
13
13
Cα
C were 3.0 and 2.5 µs, respectively, and the two-
pulse phase-modulated (TPPM) decoupling scheme was employed with a radio-frequency field strength of 85 kHz. For the 2D PDSD spectra, complex data points of 1024 × 362 and 1024 × 145 were collected for the native and amyloid states, respectively, with an acquisition delay of 2 sec, and 48 – 64 FIDs were accumulated for each t1 data point. The cross-peaks observed in all of the WT PDSD spectra were assigned in our previous solid-state NMR study.32 For the selectively
13
CO- and
13
Cα-labeled samples, the M9 media supplemented with
unlabeled amino acids (100 mg per liter culture) and
13
CO- and
13
Cα-labeled amino acids (50
mg/L) were used for protein expression. In our bacterial expression system, serine was synthesized from glycine, and thus 13Cα-Gly was used to label serine. The selective labeling of the specific amino acid was confirmed by
13
C solid-state NMR spectra of the native state of
TTR. For more extensive amino acid selective labeling, 15N-NH4Cl and 13C-glucose were used in M9 media supplemented with the other unlabeled amino acids (100 mg/L) by utilizing biosynthetic pathways of the amino acids in bacteria. 37 For example, in order to selectively label
ACS Paragon Plus Environment
8
Page 9 of 36
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60
Biochemistry
eleven amino acids (P, L, M, V, K, D, A, G, S, H and T), M9 medium was supplemented with 15
N-NH4Cl, 13C-glucose and the other nine unlabeled amino acids (100 mg/L).
Chemical Modification and Mass Spectrometry-based Analyses The chemical denaturation curves generated here were obtained using a previously established tryptophan
modification
protocol
involving
dimethyl(2-hydroxy-5-nitrobenzyl)sulfonium
bromide (HNSB). 38 Briefly, the HNSB experiments were performed by reacting the protein (1 µM) with HNSB (3.68mM) in a series of different buffers (pH 4.4) that contained increasing concentration of guanidinium chloride (GdmCl). The HNSB modification reactions in each buffer were allowed to proceed for 10 min at 37 °C. The reactions were quenched during a desalting step using Millipore C18 Zip-Tips. The desalted protein samples were each analyzed by MALDI-TOF mass spectrometry using sinapinic acid as the matrix. The MALDI mass spectra were acquired on an UltraFlex II TOF/TOF (Bruker Daltonics) mass spectrometer in the reflective and positive ion mode using a smartbeam Nd:YAG laser (355 nm). Spectra were collected using the following instrument parameters: 25kV ion source 1 voltage, 21.9kV ion source 2 voltage, 9.9kV lens voltage, 26.3kV reflector voltage, 13.7kV reflector 2 voltage, 100ns pulsed ion extraction, and matrix gating to 4000 Da. The protein ion signal in the resulting mass spectra were used to extract ∆Masswt,av values as we have previously described. 38 The time course studies of Met-30 oxidation in V30M TTR were performed by reacting the protein (1 µM) with H2O2 (100 mM) at pH 4.4 (200 mM acetate buffer) or pH 7.4 (200 mM phosphate buffer) for a specific of time at 37 °C. The reaction was quenched by adding a 6-fold molar excess of L-methionine over H2O2. A buffer exchange was performed using Amicon Ultra 10K centrifugal filters to obtain protein solutions in 100 mM triethyammonium bicarbonate
ACS Paragon Plus Environment
9
Biochemistry
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60
Page 10 of 36
buffer (pH 8.5). The protein solutions were digested with Lys-C protease, and the resulting peptide solutions were analyzed by MALDI-TOF mass spectrometry using α-cyano-4hydroxycinnamic acid (HCCA) as the matrix for MALDI-MS analysis. For determination of the oxidation status of Met-30 in the V30M mutant, the MALDI analysis was focused on the ion signals from the Met-30 containing peptide VLDAVRGSPAINVAMHVFRK in both its wildtype and oxidized forms. The MALDI ion signals of the wild-type and oxidized forms of the VLDAVRGSPAINVAMHVFRK peptide were identified by their masses. MALDI mass spectra for the time-course experiment were acquired on a Voyager-DE PRO (Applied Biosystems) mass spectrometer in the linear and positive ion mode using a nitrogen laser (337 nm). Spectra were collected using the following instrument parameters: manual control mode, 20kV accelerating voltage 93% grid, 0.05% guide wire, 100 nsec delay time, 60 shots/spectrum, and low mass gate 500 Da. The time-course data were analyzed as previously described. 38 Briefly, the increases in mass as a function of reaction times were obtained by calculating ∆Masswt,av values, using a weighted average of the intensities of the ion signals from the wild-type and oxidized peptide, and by fitting the data to a single exponential equation.
Results β-structured amyloidogenic precursor states The native state of TTR monomer adopts a primarily β-sheet tertiary structure, in which eight β-strands denoted A-H are arranged into a β-sandwich with two β-sheets (CBEF and DAGH, Figure 1). Thus, a local unfolding transition of the outer β-strands may be able to trigger aggregations to cross-β structured amyloid containing substantial native-like structures. Circular
ACS Paragon Plus Environment
10
Page 11 of 36
dichroism (CD) spectroscopy was used to examine structural features of amyloidogenic precursor states of WT and mutant forms of TTR (V30M and L55P) at the most amyloidogenic pH of 4.4. Our previous CD analysis of WT TTR showed that the amyloidogenic precursor state of WT TTR has almost identical secondary structures to those of native state of TTR. 33 The structural feature of the WT amyloidogenic state was compared to those of the TTR variants using CD spectroscopy (Figure 2). Notably, the amyloidogenic states of the more aggressive TTR variants exhibit similar CD spectra with the same minimum signal intensity at 212 nm, indicative of β-sheet conformations. The CD spectra suggest that the WT and TTR variants share similar native-like secondary structures.
15
ellipticity (mdeg)
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60
Biochemistry
Figure 2. CD spectra of WT and two mutant WT V30M L55P
10 5
forms of TTR (0.09 mg/ml) in 10 mM sodium acetate buffer (pH 4.4) recorded at 20
0
o
C.
Identical spectra were obtained at a lower protein
-5
concentration (0.02 mg/ml), suggesting that the
-10 190
200
210
220
230
240
wavelength (nm)
250 CD signals originate mainly from monomeric
precursor states of TTR. Four scans were collected and averaged for the CD spectra.
It was shown that TTR tetramers are dissociated to monomers at the amyloidogenic condition of pH 4.4. Previous ultracentrifugation experiments showed that tetrameric TTR is in a dynamic equilibrium with monomeric forms of TTR at pH 4.4 (approximately 25 % of the tetramers are present at a protein concentration of 0.2 mg/ml).24 Thus TTR would exist predominantly as a
ACS Paragon Plus Environment
11
Biochemistry
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60
Page 12 of 36
monomeric form under our experimental conditions (0.09 – 0.02 mg/ml), which indicates that the CD signals recorded at pH 4.4 mostly originate from amyloidogenic TTR monomers. Thermodynamic stabilities of the amyloidogenic precursor states were examined using a chemical modification and mass spectrometry-based approach, which has been previously described. 38 In this mass spectrometry-based analysis, we monitored the unfolding transition of the precursor states by evaluating the chemical denaturant dependence of the reaction of dimethyl(2-hydroxy-5-nitrobenzyl)sulfonium bromide (HNSB) with tryptophan (Trp) residues.38 At increasing concentrations of chemical denaturant, the TTR protein is unfolded and the Trp residues are more exposed to solvent and more readily react with HNSB, resulting in an increased mass of the protein that can be detected by using mass spectrometry. In the natively folded state of TTR, both Trp residues (W41 and W79) appear to be largely protected from the reaction with HNSB (Supplemental Figure S1), as evidenced by the relatively low ∆Masswt,av values of the pre-transition baselines of the chemical denaturation curves in Figure 3. In the presence of increasing concentrations of chemical denaturant (GdmCl in this work), the protected Trp residues in the native state become available for the chemical reaction with HNSB, leading to an increase in the protein mass, as demonstrated in Figure 3. Thus, the chemical denaturant dependence of the HNSB modification reaction can report on the thermodynamic stability of the WT and mutant forms of TTR. The cooperative unfolding transitions observed in Figure 3 suggest that the amyloidogenic precursor states adopt folded conformations. A qualitative analysis of the transition midpoints of the chemical denaturation curves generated in Figure 3 also indicates that the precursor states of the TTR variants are less stable than that of WT TTR. The less stable amyloidogenic precursor states of the TTR variants might be linked to the stronger aggregation propensity of the TTR variants.
ACS Paragon Plus Environment
12
Page 13 of 36
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60
Biochemistry
Figure 3. Chemical denaturation curves generated for WT, V30M, and L55P TTR using the tryptophan modification protocol. The vertical dotted lines indicate the concentrations of guanidinium chloride (GdmCl) at the transitions midpoints, which were 3.8, 1.6 and 1.4 M for WT, V30M, and L55P, respectively. The data points and error bars represent the average and standard deviation of the ∆Masswt,av values determined from the 8 mass spectra collected on each sample.
Structural features of the amyloid states derived from WT and mutant forms of TTR (V30M and L55P) The CD and mass spectrometry data described above reveal that the amyloidogenic precursor states of WT and mutant forms of TTR adopt overall very similar folded native-like conformations. TTR effectively forms amyloid at pH 4.4 via a downhill mechanism where major conformational changes are not necessary.39 Thus, the native-like structural features observed in the CD and mass spectrometry-based analyses may be retained in TTR aggregates. Solid-state NMR was used to probe the native-like β-structures in amyloid state of WT TTR. Twodimensional (2D) cross-polarization (CP) based
13
C-13C and
15
N-13C correlation NMR spectra
ACS Paragon Plus Environment
13
Biochemistry
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60
Page 14 of 36
were acquired for amyloid state of WT TTR and compared to that of the native state (Supplemental Figure S2). For the 13C-13C correlation experiments, proton-driven spin diffusion (PDSD) mixing scheme40 with a short mixing time (25 ms) was used to record intra-residue correlation NMR spectra (Supplemental Figure S2a and S2b). The
15
N-13C correlation NMR
spectra shown in Supplemental Figure S2c were acquired for the amide 15N and the 13Cα (NCA). The 2D 13C-13C correlation and NCA spectra show that the NMR resonances from the amyloid state (red) are overlapped well with those of the native state (black), suggesting the presence of native-like structures in amyloid state. It is notable that some of the NMR resonances in the native state are not observed in the amyloid state spectrum. Previous solid-state NMR studies of amyloids derived from various polypeptides revealed that strong signals in the dipolar-based NMR spectrum of amyloid mostly originate from highly structured amyloid core regions, while NMR signals are not observable from other unstructured regions due to extensive linebroadening.41-44 Our previous structural analyses of WT TTR amyloid also showed that loop regions such as AB and EF loop are more disordered, 33 suggesting that NMR signals from the disordered regions are not observable in the 2D solid-state NMR spectrum for the amyloid state of WT TTR in Supplemental Figure S2. The solid-state NMR studies were extended to investigate structural features of mutant forms of TTR amyloid. Figure 4a shows
13
C cross-polarization (CP) magic angle spinning (MAS)
NMR spectra of amyloid states of WT and TTR variants (V30M and L55P). The three amyloid states exhibit a similar pattern in the chemical shifts, while NMR signals at aliphatic regions (20 – 70 ppm) are weaker in intensity for the two mutant forms of TTR amyloid. These NMR results indicate that amyloid states of the two TTR variants might be slightly more disordered than that of WT TTR. Structural features of the three amyloids were also examined with 2D
13
C-13C
ACS Paragon Plus Environment
14
Page 15 of 36
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60
Biochemistry
correlation NMR experiments (Figure 4b). Overlay of the aliphatic regions of the 2D spectra for the three amyloid state suggest that the amyloid state of WT and mutant forms of TTR share a similar structural feature in amyloid core regions that produce strong NMR signals in the dipolarbased 2D correlation experiments.
Figure 4. (a)
13
C CPMAS spectra of the amyloid states of WT (black), V30M (red), L55P
(green) TTR where eleven amino acids (P, L, M, V, K, D, A, G, S, H and T) are uniformly 13
C/15N labeled. (b) Overlay of the 2D
13
C-13C correlation NMR spectra of the three amyloid
states with a PDSD mixing time of 25 ms. Each PDSD spectrum was compared with the other two PDSD spectra separately [WT (black), V30M (red), L55P (green)].
ACS Paragon Plus Environment
15
Biochemistry
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60
Page 16 of 36
DAGH β -sheet structure in amyloid states The combined structural analyses described above indicate that the amyloid sates of TTR contain extensive native-like β-sheet structures. Solid-state NMR has proven to be a powerful tool for detailed structural characterization of non-crystalline amyloid at atomic resolution. 43-55 However, structural analyses of the tertiary β-structure within this 127-residue TTR amyloid are a daunting task. We recently demonstrated that specific labeling schemes, which generate isolated 13CO-13Cα dipolar-coupled spin pairs in the native β-sheets (Figure 5a), are of great use for probing the native β-sheet structure in amyloid state.32,33 For example, selective
13
CO-Leu
and 13Cα-Tyr labeling will produce three dipolar-coupled 13CO-13Cα spin pairs at distances of 4– 6 Å in the AGH strands. Solid-state NMR experiments that detect spin-pairs with a separation of up to 6 Å can then be used to investigate any structural changes in the β-structure. Figure 5b shows 13C cross-polarization (CP) magic-angle spinning (MAS) NMR spectra of the native and amyloid states of TTR. The NMR peaks at ~175 and ~60 ppm in the 13C CPMAS spectrum of the L/Y labeled TTR confirms the selective enrichment of
13
CO and
13
Cα carbons, respectively
(Figure 5b and Supplemental Figure S3). Although the NMR resonances in the amyloid states are substantially broader than those of the native state, our previous solid-state NMR experiments showed that TTR amyloid used in our studies is a uniform protein assembly consisting of a single monomeric core conformer on the basis of a single set of cross-peaks in extensive
13
C-13C correlation solid-state NMR spectra.
32,33
It is also notable that the NMR
resonance for residue Y78 in EF loop is more significantly broadened in the amyloid states, indicating that the loop regions might be disordered in the amyloid states.
ACS Paragon Plus Environment
16
Page 17 of 36
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60
Biochemistry
a
C
B
K48 G47
A29
E
F
V30 5.0 I73 4.9 A91 5.0Å S46 H31 E72 E92 V32 A45 V71 V93 K70 4.9 V94 F44 F33 P43 R34 4.4 Y69 F95 E42 K35 I68 T96 W41 A36 G67 A97
A D
G 5.7 R104
H T123
L12 Y105 5.0 V122 4.7Å M13 T106 V121 L55 V14 4.3 I107 A120 E54
K15 V16 L17 D18
A108 A109
T119 4.3 T118 S117 L110
L111
Y116
S112 5.2 S115
b
Figure 5. (a) Selective labeling schemes for the solid-state NMR experiments. The (green) and
13
13
CO
Cα (red) carbons of the amino acids are labeled if the internuclear distance in the
native TTR tetramer is within 6 Å. (b)
13
C CPMAS spectra of the native and amyloid states of
TTR prepared from 13CO-Leu and 13Cα-Tyr labeled protein. The NMR resonance of the residue Y78 in the EF helix is slightly shifted upfield and becomes much broader in amyloid state, suggesting that the EF helix becomes disordered in the amyloid state.
The 2D 13C-13C correlation experiments using a long PDSD mixing time of 500 ms were used to probe the weakly dipolar-coupled spin pairs at distances of 5 – 6 Å in AGH β-structure in the amyloid states of WT and TTR variants (Figure 5a and Figure 6). Figure 6a (the left spectrum) shows overlay of the PDSD spectra for the L/Y (black) and L/M (red) spin pairs in WT amyloid state. The two cross-peaks for the L110-Y116 and L111-Y116, and one cross-peak for L12-Y105 (black in Figure 6a) were observed in the WT TTR amyloid. The 2D NMR experiments were
ACS Paragon Plus Environment
17
Biochemistry
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60
Page 18 of 36
also conducted to probe the L/M spin pair in the DA β-structure for WT amyloid (Figure 5a and red in the left spectrum in Figure 6a). The two cross-peaks for the L55-M13 and L12-M13 were detected in the WT amyloid state. These NMR results indicate the presence of native-like DAGH β-structure in WT amyloid state. Additional evidence for the DAGH β-structure was obtained with the 2D PDSD experiments on a
13
CO-Val/13Cα-Tyr labeled TTR amyloid (Figure 5a and
right spectrum in Figure 6a) that supports a close proximity between strands G and H. The solid-state NMR experiments were conducted on amyloid states of the two TTR variants (Figure 6b and 6c). The NMR cross-peaks for the L/Y and V/Y spin pairs in the AGH β-structure (black in Figure 6) were all observed in the mutant forms of TTR amyloid, suggesting that the mutant amyloids also contain the native-like AGH β-structure. Similar structural features in AGH structure is supported by additional 2D experiments on I/V spin pairs (AG strands in Figure 5a and Supplemental Figure S4). In contrast to the spin-pairs in the AGH β-structure, the cross-peak from the L55-M13 spin pair in DA β-structure is not observed in V30M amyloid (Figure 6b). The different 2D spectrum for V30M mutant indicates that the V30M TTR may have a different DA substructure in amyloid state. In the stronger amyloidogenic mutant L55P, the distance between P55 and M13 is higher (6.5 Å) even in the native state, in comparison to the native tetrameric TTR (5 Å), as shown in Supplemental Figure S5. The strand D is also more disordered in the native state of the more amyloidogenic mutant L55P TTR. Theses results indicate that destabilization of the DA β-structure may be critical to the effective amyloid formation, as suggested by previous computational studies. 56,57
ACS Paragon Plus Environment
18
Page 19 of 36
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60
Biochemistry
Figure 6. 2D PDSD spectra of amyloid states of the 13CO-Leu/13Cα-Tyr (black, left; labeled on AGH strands),
13
CO-Leu/13Cα-Met (red; labeled on DA), and
13
CO-Val/13Cα-Tyr (black, right;
labeled on GH) TTR with a contour level of 3.0 % with respect to the diagonal peak obtained using a mixing time of 500 ms at a 1H frequency of 830 MHz. L/Y and L/M spin pairs were unambiguously assigned in our previous solid-state NMR studies.32 In order to confirm that the cross-peaks originate from the spin pairs within the TTR monomer, the 2D PDSD experiments were conducted on the mixture of singly 13CO- and 13Cα- labeled TTRs (1:1 ratio, for example a mixture of
13
CO-Leu-TTR and
13
Cα-Tyr-TTR) and no cross-peak was observed from the
mixture.
ACS Paragon Plus Environment
19
Biochemistry
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60
Page 20 of 36
Native-like CBEF β -sheet structure The solid-state NMR experiments were extended to probe CBEF β-structure. Our previous solid-state NMR experiments showed that cross-peaks for I73-M30 and I73-A91 spin pairs in the BEF β-structure (Figure 5a) were observed in the 2D PDSD spectra obtained for V30M amyloid,32 which is in good agreement with the NMR spectra for the I/V labeling scheme (Supplemental Figure S4). Additional F/Y spin pairs in the BEF β-structure were investigated with the
13
C-13C correlation experiments (Figure 7). Our recent solid-state NMR experiments
combined with single-point mutagenesis successfully assigned the two NMR cross-peaks for F33-Y69 and F95-Y69 that are overlapped in the 2D spectrum of WT TTR amyloid sample.32 The two cross-peaks for the F33-Y69 and F95-Y69 observed in the WT amyloid (Figure 7) were also detected with the same line width and intensity in those of the two mutant forms, suggesting the two single-point pathogenic mutations do not affect the BEF β-structure in mutant amyloids.
Figure 7. 2D PDSD spectra of amyloid states of the
13
CO-Phe and
13
Cα-Tyr (labeled on BEF
strands) TTR obtained using a mixing time of 500 ms at a 1H frequency of 830 MHz. * denotes spinning sidebands.
The 2D PDSD experiments were also conducted to investigate CB β-structure in the amyloid states of the TTR variants (Figure 5a and Figure 8). For the structural study of the CB
ACS Paragon Plus Environment
20
Page 21 of 36
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60
Biochemistry
substructure, three amino acids (13CO-His,
13
Cα-Gly and
13
Cα-Ser) were selectively labeled,
generating the H31-S46 spin pair in strands B and C, respectively (Figure 8). The 2D
13
C-13C
correlation NMR spectra were collected for the amyloid states of the two TTR variants and compared to that of WT TTR amyloid. The two cross-peaks in the WT amyloid PDSD spectrum were assigned in our previous solid-state NMR studies.32 The presence of the similar cross-peak for H31-S46 in both V30M and L55P amyloid suggest that the two mutant amyloids have a similar CB substructure. Additional evidence for the native-like CB β-structure in the amyloid states was obtained with the 2D PDSD experiments of
13
CO-Gly and
13
Cα-Val TTR samples
(Figure 5a and 9). The cross-peak from the G47-V30 spin pair that is absent in the V30M spectrum supports the intact CB β-structure in the WT and L55P mutant amyloids.
Figure 8. 2D PDSD spectra of amyloid states of the 13CO-His, 13Cα-Gly, and 13Cα-Ser (labeled on CB strands) TTR obtained using a mixing time of 500 ms at a 1H frequency of 830 MHz. * denotes spinning sidebands, and the same contour level was used for the three spectra.
ACS Paragon Plus Environment
21
Biochemistry
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60
Page 22 of 36
Figure 9. 2D PDSD spectra of amyloid states of the 13COGly and
13
Cα-Val (labeled on CB strands) TTR using a
mixing time of 500 ms at a 1H frequency of 830 MHz [WT (black), L55P (green), and V30M (red)]. The cross-peak marked by ? could not be assigned. The cross-peak indicated by the dashed line is absent in the V30M amyloid (red), and thus it was assigned as G47-V30.
Time-course analysis of Met-30 oxidation in V30M TTR The solid-state NMR experiments of V30M amyloid described above suggest that strand D might be unfolded in amyloid state, which may fully expose strand A for intermolecular interactions. The sidechain of residue Met-30 points toward the DA substructure (Supplemental Figure S6) in the native V30M TTR, and the longer sidechain of methionine than that of valine may destabilize the DA structure, inducing a structural change under the amyloidogenic condition. In order to probe the local environment of the methionine sidechain, oxidation rates of the methionine sidechain were analyzed using mass spectrometry. 58 In this mass spectrometrybased approach we measured the time-course of the H2O2 oxidation of the Met-30 sidechain. In the native tetrameric V30M TTR at pH 7.4, the sidechain of residue Met-30 is protected from the oxidation (Figure 10). On the contrary, the Met-30 sidechain is readily oxidized at pH 4.4, suggesting that the M30 sidechain protected in the native state is exposed to solvent in the amyloidogenic precursor state at pH 4.4. We note that the H2O2 oxidation of an unprotected Met
ACS Paragon Plus Environment
22
Page 23 of 36
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60
Biochemistry
residue is not significantly different at pHs 4.5 and 7.4 (see Supplemental Table S1). Thus, the different reaction rates observed in Figure 10 appear to be due to structural differences in the protein. We also note that the second-order reaction rate constant calculated from the time-course data on V30M TTR at pH 4.4 and 37 oC (63 M-1h-1) was only slightly smaller than that measured for an unprotected methionine in the model peptide under the identical condition (82 M-1h-1)59 (see Supplemental Table S1).
Figure 10. Time course analysis of Met-30 sidechain oxidation with H2O2 in V30M TTR. Data obtained on the Lys-C generated TTR peptide VLDAVRGSPAINVAMHVFRK, which contains residue Met-30, are shown. The data points and error bars represent the average and standard deviation of the ∆Masswt,av values determined from the five mass spectra collected on each sample. The dotted blue line represents the best fit of the pH 7.4 data to a line. The solid red line represents the best fit of the data to a single exponential equation.
ACS Paragon Plus Environment
23
Biochemistry
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60
Page 24 of 36
Discussion Despite the extensive mechanistic studies of TTR misfolding and amyloid formation, detailed molecular mechanism has remained largely unknown. Structural characterizations of amyloidogenic precursor states and the final product amyloid are essential to the understanding of the amyloid formation mechanism. However, structural studies of amyloidogenic precursor states are of great challenge due to their dynamic, aggregation-prone properties. In this study, we employed CD spectroscopy and mass spectrometry to investigate structural features of amyloidogenic precursor states of TTR using a very low protein concentration (1 µM) to minimize aggregation. Our structural analyses using CD spectroscopy showed that the amyloidogenic precursor states of WT TTR adopt a native-like β-sheet structure at the most amyloidogenic pH of 4.4. The CD results are in good agreements with our mass spectrometrybased thermodynamic analyses of WT TTR at pH 4.4 that showed the precursor state undergoes a cooperative unfolding transition upon additions of a denaturant, characteristic of folded proteins. The CD spectra, in conjunction with mass spectrometry-based unfolding studies, clearly demonstrated that amyloidogenic precursor states of WT TTR adopt largely folded native-like βsheet conformations. TTR amyloidogenesis is characterized by more than 100 single-point mutations that cause early onset familial amyloidoses. Most of the pathogenic mutations tested in vitro were shown to kinetically or thermodynamically destabilize the native tetrameric state of TTR, accelerating the dissociation of the tetramers to monomers, which is the rate-determining step in TTR amyloid formation process. 8,9,16,24 However, the effect of the single-point pathogenic mutations on the molecular structure of the monomeric precursor states and amyloid has not been reported. Previous structural studies using X-ray crystallography revealed that the pathogenic TTR
ACS Paragon Plus Environment
24
Page 25 of 36
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60
Biochemistry
variants have almost identical crystal structures to that of the tetrameric state of WT TTR. Only local structural perturbations at the vicinity of the mutation were observed. Recent relaxation NMR experiments identified flexible regions undergoing conformational fluctuations that might be linked to the dissociation of the TTR tetramer to aggregation-prone monomers.60 The pathogenic mutations (V30M and L55P) were shown to induce more extensive conformational fluctuations of the same flexible regions as those observed for WT TTR. However, a notable difference was observed in strand D. The short strand with little flexibility in WT TTR undergoes substantial conformational fluctuations in the pathogenic mutants.60 These NMR results are largely consistent with our biophysical studies that showed the amyloidogenic states of the TTR variants are less stable than that of WT TTR. In addition, our solid-state NMR results reveal that amyloid states of WT and two mutant forms of TTR adopt a similar native-like CBEF and AGH β-sheet structure, but different DA substructures were observed in the mutant amyloids. The solid-state NMR results described here suggest that the native-like DA substructure is maintained in the WT amyloid state on the basis of the cross-peak from the L55/M13 spin pair. However, the NMR signal for the spin pair in the DA β-strands was not observed in V30M amyloid, indicating that strand D may undergo a structural change. Although the single-point mutation (V30M) is observed in strand C, the side chain of the methionine points toward the strands D and A in the other DAGH β-sheet (Supplemental Figure S6). The longer side chain of methionine than that of valine may crash with other side chains in strands D and A, inducing a structural change in the DA substructure, as was suggested by the previous relaxation dispersion NMR experiments. 60 In the stronger amyloidogenic mutant L55P, the distance between residues P55 and M13 is higher (6.5 Å) even in the native state, in comparison to the native WT
ACS Paragon Plus Environment
25
Biochemistry
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60
Page 26 of 36
tetrameric TTR (5 Å), as shown in Supplementary Figure S5. The strand D is also more disordered in the more amyloidogenic mutant L55P TTR, which is consistent with our mass spectra-based analysis that showed the amyloidogenic state of the mutant is less stable than that of WT TTR. In addition, the sidechain of residue Met-30 that points toward DA substructure is protected from oxidations in the native V30M, while the sidechain is readily available for the oxidation in the amyloidogenic state of V30M TTR. The combined results indicate that the DA substructure becomes disrupted in the amyloidogenic states of the TTR variants. Previous H/D exchange experiments also reported the structural change of strand D in another TTR variant (Y114C) amyloid. 22 In addition, triple mutations on strand D (G53S/E54D/L55S) were shown to slightly unfold the strand D, accelerating the aggregation.61 Taken all together, the destabilization of the DA β-structure under amyloidogenic conditions may trigger more effective amyloid formation of the TTR variants, as suggested by other biophysical studies. 22,34,56,62 The structural studies of the amyloidogenic precursor states and the final product, amyloid, provide valuable insights into TTR aggregation pathway. Our previous solid-state studies revealed that the native-like CBEF and DAGH β-structure is retained in WT TTR amyloid state, while AB loop regions that are hydrogen-bonded with strand A undergo a structural transition. 33 The partly exposed strand A might be involved in intermolecular interactions essential to TTR amyloid formation. In our present solid-state NMR studies reported here, the amyloid states of the two TTR variants are also shown to adopt the similar native-like CBEF and AGH β-sheet structure. However, strand D appears undergo a structural change in the amyloidogenic states of the TTR variants, which may fully expose strand A for intermolecular interactions. The structural change of the strand D may accelerate amyloid formation of the pathogenic TTR variants, which is in good agreements with previous computational studies. 56,62
ACS Paragon Plus Environment
26
Page 27 of 36
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60
Biochemistry
The structural information of the monomeric conformer in amyloid state will also help to develop a quaternary structural model of TTR amyloid. The combined solid-state NMR studies with CD and mass spectrometry-based analyses show that the monomeric conformer in WT TTR amyloid adopts a native-like β-sheet structure, suggesting that only outer strands (C, F, D, H, and part of the strand A due to the structural change of the AB loop region) are available for the intermolecular interactions. In addition, previous relaxation dispersion NMR experiments of the monomeric and tetrameric TTR showed that CBEF β-sheet remains stable while the other βsheet (DAGH) undergoing extensive conformational fluctuations on millisecond time scales associated with misfolding and amyloid formation. 60,63 These NMR results suggest that DAGH β-sheet may be involved in intermolecular associations. Indeed, recent mechanistic studies of TTR aggregation showed that intermolecular interactions between strands H play a key role in TTR amyloid formation, and strand A has a strong aggregation propensity. 64 In addition, TTR dimers cross-linked between the strands H and H’ were shown to form amyloid. 65 On the basis of the previous biophysical studies and our solid-state NMR results, we propose that the nativelike monomeric TTR forms a dimer through H-H’ interactions, and the dimers may oligomerize via A-A’ interactions (Figure 11), as was suggested by previous EPR studies.25 The more exposed strand A due to the structural change of strand D in the TTR variants may not only facilitate the A-A’ interactions, accelerating TTR aggregation, but also lead to different intermolecular association patterns and thus distinct quaternary structures. The mutant amyloids with distinct quaternary structures may interfere with tissue-specific cellular components such as receptors and signaling proteins, resulting in phenotype diversities observed in TTR amyloidogenesis.
ACS Paragon Plus Environment
27
Biochemistry
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60
Page 28 of 36
Figure 11. Proposed structural models for the amyloid states of WT and mutant forms of TTR. On the basis of our solution and solid-state NMR results, we propose that the CBEF β-sheet remains unchanged in the amyloid states. However, we cannot completely rule out that the TTR amyloid used in this work may contain a minor conformer where the strand C is unfolded and involved in intermolecular interactions, as was proposed by the EPR studies based on the EPR signal from strands B and B’.25 It is also plausible that the EPR signal may come from interactions between protofibrils that contain native-like CBEF β-sheet structures.
Finally, there is growing evidence that amyloid formation of natively folded proteins can proceed via native-like intermediates populated by conformational fluctuations of the natively folded state. 66-69 Although structural characterization of the amyloidogenic intermediate state is essential to understanding molecular mechanism of the amyloid formation process, the structural analysis is of great challenge due to transient and aggregation-prone nature of the intermediate states. Our CD and mass spectrometry-based thermodynamic analyses used relatively low protein concentrations (~ 1 µM, monomer concentration) to prevent protein aggregation, which allowed the structural characterization of the amyloidogenic precursor states. Solid-state NMR was then used to investigate the native-like structures in the final product, amyloid, using the selectively labeling schemes. The combined structural studies of the initial and final states of the protein in amyloid formation process would provide valuable insights into molecular
ACS Paragon Plus Environment
28
Page 29 of 36
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60
Biochemistry
mechanisms of misfolding and amyloid formation of the other natively folded proteins such as superoxide dismutase 1 (SOD1) and β-microglobulin.
ASSOCIATED CONTENT Supporting Information. 2D PDSD NMR spectra for the mutant TTR amyloids. Schematic diagrams of monomeric TTR displaying sidechains of Trp and Met. Oxidation rates for the Metcontaining model peptides. The following files are available free of charge.
AUTHOR INFORMATION Corresponding Author Email:
[email protected] Funding Sources This work was supported by NIH Grants NS084138 (KHL), NS097490 (KHL) and DK46335 (JWK), and the Skaggs Institute of Chemical Biology (JWK). The solid-state NMR spectra were acquired at the National High Magnetic Field Laboratory, which is supported by NSF Cooperative Agreement No. DMR-1157490 and the State of Florida. Notes The authors declare no competing financial interests.
ACS Paragon Plus Environment
29
Biochemistry
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60
Page 30 of 36
ACKNOWLEDGMENT We thank Profs. P. E. Wright (Scripps) and D. E. Wemmer (University of California at Berkeley) for helpful discussion. We also acknowledge Prof. Kenney (ECU) for assistance in the CD experiments.
REFERENCES (1) Dobson, C. M. (2003) Protein folding and misfolding. Nature 426, 884-890. (2) Kelly, J. W. (1998) The alternative conformations of amyloidogenic proteins and their multi-step assembly pathways. Curr. Opin. Struct. Biol. 8, 101-106. (3) Jahn, T. R., and Radford, S. E. (2008) Folding versus aggregation: polypeptide conformations on competing pathways. Arch Biochem Biophys. 469, 100-117. (4) Jacobson, D. R., Pastore, R. D., Yaghoubian, R., Kane, I., Gallo, G., Buck, F. S., and Buxbaum, J. N. (1997) Variant-sequence transthyretin (isoleucine 122) in late-onset cardiac amyloidosis in black Americans. N. Engl. J. Med. 336, 466-473. (5) Connors, L. H., Lim, A., Prokaeva, T., Roskens, V. A., and Costello, C. E. (2003) Tabulation of human transthyretin (TTR) variants, 2003. Amyloid 10, 160-184. (6) Saraiva, M. J. (1995) Transthyretin mutations in health and disease. Hum. Mutat. 5, 191-196. (7) Rapezzi, C., Quarta, C. C., Riva, L., Longhi, S., Gallelli, I., Lorenzini, M., Ciliberti, P., Biagini, E., Salvi, F., and Branzi, A. (2010) Transthyretin-related amyloidoses and the heart: a clinical overview. Nat Rev Cardiol. 7, 398-408. (8) Hammarstrom, P., Jiang, X., Hurshman, A. R., Powers, E. T., and Kelly, J. W. (2002) Sequencedependent denaturation energetics: A major determinant in amyloid disease diversity. Proc Natl Acad Sci U S A. 99 Suppl 4, 16427-16432. (9) Sekijima, Y., Wiseman, R. L., Matteson, J., Hammarstrom, P., Miller, S. R., Sawkar, A. R., Balch, W. E., and Kelly, J. W. (2005) The biological and chemical basis for tissue-selective amyloid disease. Cell 121, 73-85. (10) Westermark, P., Sletten, K., Johansson, B., and Cornwell, G. G.,3rd. (1990) Fibril in senile systemic amyloidosis is derived from normal transthyretin. Proc. Natl. Acad. Sci. USA. 87, 2843-2845. (11) Sousa, M. M., and Saraiva, M. J. (2003) Neurodegeneration in familial amyloid polyneuropathy: from pathology to molecular signaling. Prog. Neurobiol. 71, 385-400.
ACS Paragon Plus Environment
30
Page 31 of 36
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60
Biochemistry
(12) Ihse, E., Ybo, A., Suhr, O., Lindqvist, P., Backman, C., and Westermark, P. (2008) Amyloid fibril composition is related to the phenotype of hereditary transthyretin V30M amyloidosis. J. Pathol. 216, 253-261. (13) Buxbaum, J. N., and Reixach, N. (2009) Transthyretin: the servant of many masters. Cell Mol. Life Sci. 66, 3095-3101. (14) Blake, C. C., Geisow, M. J., Oatley, S. J., Rerat, B., and Rerat, C. (1978) Structure of prealbumin: secondary, tertiary and quaternary interactions determined by Fourier refinement at 1.8 A. J. Mol. Biol. 121, 339-356. (15) Hamilton, J. A., Steinrauf, L. K., Braden, B. C., Liepnieks, J., Benson, M. D., Holmgren, G., Sandgren, O., and Steen, L. (1993) The x-ray crystal structure refinements of normal human transthyretin and the amyloidogenic Val-30-->Met variant to 1.7-A resolution. J. Biol. Chem. 268, 2416-2424. (16) Connelly, S., Choi, S., Johnson, S. M., Kelly, J. W., and Wilson, I. A. (2010) Structure-based design of kinetic stabilizers that ameliorate the transthyretin amyloidoses. Curr. Opin. Struct. Biol. 20, 54-62. (17) Lai, Z. H., Colon, W., and Kelly, J. W. (1996) The acid-mediated denaturation pathway of transthyretin yields a conformational intermediate that can self-assemble into amyloid. Biochemistry 35, 6470-6482. (18) Lashuel, H. A., Lai, Z. H., and Kelly, J. W. (1998) Characterization of the transthyretin acid denaturation pathways by analytical ultracentrifugation: Implications for wild-type, V30M, and L55P amyloid fibril formation. Biochemistry 37, 17851-17864. (19) Quintas, A., Vaz, D. C., Cardoso, I., Saraiva, M. J., and Brito, R. M. (2001) Tetramer dissociation and monomer partial unfolding precedes protofibril formation in amyloidogenic transthyretin variants. J. Biol. Chem. 276, 27207-27213. (20) Palaninathan, S. K. (2012) Nearly 200 X-ray crystal structures of transthyretin: what do they tell us about this protein and the design of drugs for TTR amyloidoses? Curr. Med. Chem. 19, 2324-2342. (21) Hornberg, A., Eneqvist, T., Olofsson, A., Lundgren, E., and Sauer-Eriksson, A. E. (2000) A comparative analysis of 23 structures of the amyloidogenic protein transthyretin. J. Mol. Biol. 302, 649669. (22) Olofsson, A., Ippel, J. H., Wijmenga, S. S., Lundgren, E., and Ohman, A. (2004) Probing solvent accessibility of transthyretin amyloid by solution NMR spectroscopy. J. Biol. Chem. 279, 5699-5707. (23) Liu, K., Cho, H. S., Lashuel, H. A., Kelly, J. W., and Wemmer, D. E. (2000) A glimpse of a possible amyloidogenic intermediate of transthyretin. Nat. Struct. Biol. 7, 754-757. (24) Lai, Z. H., Colon, W., and Kelly, J. W. (1996) The acid-mediated denaturation pathway of transthyretin yields a conformational intermediate that can self-assemble into amyloid. Biochemistry 35, 6470-6482. (25) Serag, A. A., Altenbach, C., Gingery, M., Hubbell, W. L., and Yeates, T. O. (2002) Arrangement of subunits and ordering of beta-strands in an amyloid sheet. Nat. Struct. Biol. 9, 734-739.
ACS Paragon Plus Environment
31
Biochemistry
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60
Page 32 of 36
(26) Bateman, D. A., Tycko, R., and Wickner, R. B. (2011) Experimentally derived structural constraints for amyloid fibrils of wild-type transthyretin. Biophys. J. 101, 2485-2492. (27) Liu, K., Kelly, J. W., and Wemmer, D. E. (2002) Native state hydrogen exchange study of suppressor and pathogenic variants of transthyretin. J. Mol. Biol. 320, 821-832. (28) Liu, K., Cho, H. S., Hoyt, D. W., Nguyen, T. N., Olds, P., Kelly, J. W., and Wemmer, D. E. (2000) Deuterium-proton exchange on the native wild-type transthyretin tetramer identifies the stable core of the individual subunits and indicates mobility at the subunit interface. J. Mol. Biol. 303, 555-565. (29) Elam, J. S., Taylor, A. B., Strange, R., Antonyuk, S., Doucette, P. A., Rodriguez, J. A., Hasnain, S. S., Hayward, L. J., Valentine, J. S., Yeates, T. O., and Hart, P. J. (2003) Amyloid-like filaments and water-filled nanotubes formed by SOD1 mutant proteins linked to familial ALS. Nat. Struct. Biol. 10, 461-467. (30) Laidman, J., Forse, G. J., and Yeates, T. O. (2006) Conformational change and assembly through edge beta strands in transthyretin and other amyloid proteins. Acc. Chem. Res. 39, 576-583. (31) Bateman, D. A., Tycko, R., and Wickner, R. B. (2011) Experimentally derived structural constraints for amyloid fibrils of wild-type transthyretin. Biophys. J. 101, 2485-2492. (32) Lim, K. H., Dasari, A. K., Hung, I., Gan, Z., Kelly, J. W., Wright, P. E., and Wemmer, D. E. (2016) Solid-State NMR Studies Reveal Native-like beta-Sheet Structures in Transthyretin Amyloid. Biochemistry 55, 5272-5278. (33) Lim, K. H., Dasari, A. K., Hung, I., Gan, Z., Kelly, J. W., and Wemmer, D. E. (2016) Structural Changes Associated with Transthyretin Misfolding and Amyloid Formation Revealed by Solution and Solid-State NMR. Biochemistry 55, 1941-1944. (34) Colon, W., and Kelly, J. W. (1992) Partial Denaturation of Transthyretin Is Sufficient for Amyloid Fibril Formation Invitro. Biochemistry 31, 8654-8660. (35) Raleigh, D. P., Levitt, M. H., and Griffin, R. G. (1988) Rotational resonance in solid state NMR. Chem. Phys. Lett. 147, 71-76. (36) Seidel, K., Lange, A., Becker, S., Hughes, C. E., Heise, H., and Baldus, M. (2004) Protein solid-state NMR resonance assignments from (13C,13C) correlation spectroscopy. Phys. Chem. Chem. Phys. 6, 5090-5093. (37) Hong, M., and Jakes, K. (1999) Selective and extensive 13C labeling of a membrane protein for solid-state NMR investigations. J. Biomol. NMR. 14, 71-74. (38) Xu, Y., Strickland, E. C., and Fitzgerald, M. C. (2014) Thermodynamic analysis of protein folding and stability using a tryptophan modification protocol. Anal. Chem. 86, 7041-7048. (39) Hurshman, A. R., White, J. T., Powers, E. T., and Kelly, J. W. (2004) Transthyretin aggregation under partially denaturing conditions is a downhill polymerization. Biochemistry 43, 7365-7381.
ACS Paragon Plus Environment
32
Page 33 of 36
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60
Biochemistry
(40) Szeverenyi, N. M., Sullivan, M. J., and Maciel, G. E. (1982) Observation of spin exchange by twodimensional fourier transform 13C cross polarization-magic-angle spinning. J. Magn. Reson. 47, 462-475. (41) Helmus, J. J., Surewicz, K., Surewicz, W. K., and Jaroniec, C. P. (2010) Conformational flexibility of Y145Stop human prion protein amyloid fibrils probed by solid-state nuclear magnetic resonance spectroscopy. J. Am. Chem. Soc. 132, 2393-2403. (42) Debelouchina, G. T., Platt, G. W., Bayro, M. J., Radford, S. E., and Griffin, R. G. (2010) Magic angle spinning NMR analysis of beta2-microglobulin amyloid fibrils in two distinct morphologies. J. Am. Chem. Soc. 132, 10414-10423. (43) Daebel, V., Chinnathambi, S., Biernat, J., Schwalbe, M., Habenstein, B., Loquet, A., Akoury, E., Tepper, K., Muller, H., Baldus, M., Griesinger, C., Zweckstetter, M., Mandelkow, E., Vijayan, V., and Lange, A. (2012) beta-Sheet core of tau paired helical filaments revealed by solid-state NMR. J. Am. Chem. Soc. 134, 13982-13989. (44) Barbet-Massin, E., Ricagno, S., Lewandowski, J. R., Giorgetti, S., Bellotti, V., Bolognesi, M., Emsley, L., and Pintacuda, G. (2010) Fibrillar vs crystalline full-length beta-2-microglobulin studied by high-resolution solid-state NMR spectroscopy. J. Am. Chem. Soc. 132, 5556-5557. (45) Tang, M., Comellas, G., and Rienstra, C. M. (2013) Advanced solid-state NMR approaches for structure determination of membrane proteins and amyloid fibrils. Acc. Chem. Res. 46, 2080-2088. (46) Lu, J. X., Qiang, W., Yau, W. M., Schwieters, C. D., Meredith, S. C., and Tycko, R. (2013) Molecular structure of beta-amyloid fibrils in Alzheimer's disease brain tissue. Cell 154, 1257-1268. (47) Debelouchina, G. T., Bayro, M. J., Fitzpatrick, A. W., Ladizhansky, V., Colvin, M. T., Caporini, M. A., Jaroniec, C. P., Bajaj, V. S., Rosay, M., Macphee, C. E., Vendruscolo, M., Maas, W. E., Dobson, C. M., and Griffin, R. G. (2013) Higher order amyloid fibril structure by MAS NMR and DNP spectroscopy. J. Am. Chem. Soc. 135, 19237-19247. (48) Wasmer, C., Lange, A., Van Melckebeke, H., Siemer, A. B., Riek, R., and Meier, B. H. (2008) Amyloid fibrils of the HET-s(218-289) prion form a beta solenoid with a triangular hydrophobic core. Science 319, 1523-1526. (49) Helmus, J. J., Surewicz, K., Apostol, M. I., Surewicz, W. K., and Jaroniec, C. P. (2011) Intermolecular alignment in Y145Stop human prion protein amyloid fibrils probed by solid-state NMR spectroscopy. J. Am. Chem. Soc. 133, 13934-13937. (50) Linser, R., Dasari, M., Hiller, M., Higman, V., Fink, U., Lopez del Amo, J. M., Markovic, S., Handel, L., Kessler, B., Schmieder, P., Oesterhelt, D., Oschkinat, H., and Reif, B. (2011) Proton-detected solid-state NMR spectroscopy of fibrillar and membrane proteins. Angew. Chem. Int. Ed Engl. 50, 45084512. (51) Hong, M., and Schmidt-Rohr, K. (2013) Magic-angle-spinning NMR techniques for measuring longrange distances in biological macromolecules. Acc. Chem. Res. 46, 2154-2163. (52) Yan, S., Suiter, C. L., Hou, G., Zhang, H., and Polenova, T. (2013) Probing structure and dynamics of protein assemblies by magic angle spinning NMR spectroscopy. Acc. Chem. Res. 46, 2047-2058.
ACS Paragon Plus Environment
33
Biochemistry
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60
Page 34 of 36
(53) Parthasarathy, S., Inoue, M., Xiao, Y., Matsumura, Y., Nabeshima, Y., Hoshi, M., and Ishii, Y. (2015) Structural Insight into an Alzheimer's Brain-Derived Spherical Assembly of Amyloid beta by Solid-State NMR. J. Am. Chem. Soc. 137, 6480-6483. (54) Weingarth, M., and Baldus, M. (2013) Solid-state NMR-based approaches for supramolecular structure elucidation. Acc. Chem. Res. 46, 2037-2046. (55) Kotler, S. A., Brender, J. R., Vivekanandan, S., Suzuki, Y., Yamamoto, K., Monette, M., Krishnamoorthy, J., Walsh, P., Cauble, M., Holl, M. M., Marsh, E. N., and Ramamoorthy, A. (2015) High-resolution NMR characterization of low abundance oligomers of amyloid-beta without purification. Sci. Rep. 5, 11811. (56) Rodrigues, J. R., Simoes, C. J., Silva, C. G., and Brito, R. M. (2010) Potentially amyloidogenic conformational intermediates populate the unfolding landscape of transthyretin: insights from molecular dynamics simulations. Protein Sci. 19, 202-219. (57) Lei, M., Yang, M., and Huo, S. (2004) Intrinsic versus mutation dependent instability/flexibility: a comparative analysis of the structure and dynamics of wild-type transthyretin and its pathogenic variants. J. Struct. Biol. 148, 153-168. (58) West, G. M., Tang, L., and Fitzgerald, M. C. (2008) Thermodynamic analysis of protein stability and ligand binding using a chemical modification- and mass spectrometry-based strategy. Anal. Chem. 80, 4175-4185. (59) Yin, J., Chu, J. W., Ricci, M. S., Brems, D. N., Wang, D. I., and Trout, B. L. (2004) Effects of antioxidants on the hydrogen peroxide-mediated oxidation of methionine residues in granulocyte colonystimulating factor and human parathyroid hormone fragment 13-34. Pharm. Res. 21, 2377-2383. (60) Das, J. K., Mall, S. S., Bej, A., and Mukherjee, S. (2014) Conformational flexibility tunes the propensity of transthyretin to form fibrils through non-native intermediate states. Angew. Chem. Int. Ed Engl. 53, 12781-12784. (61) Eneqvist, T., Andersson, K., Olofsson, A., Lundgren, E., and Sauer-Eriksson, A. E. (2000) The betaslip: a novel concept in transthyretin amyloidosis. Mol. Cell. 6, 1207-1218. (62) Yang, M., Lei, M., and Huo, S. (2003) Why is Leu55-->Pro55 transthyretin variant the most amyloidogenic: insights from molecular dynamics simulations of transthyretin monomers. Protein. Sci. 12, 1222-1231. (63) Lim, K. H., Dyson, H. J., Kelly, J. W., and Wright, P. E. (2013) Localized structural fluctuations promote amyloidogenic conformations in transthyretin. J. Mol. Biol. 425, 977-988. (64) Saelices, L., Johnson, L. M., Liang, W. Y., Sawaya, M. R., Cascio, D., Ruchala, P., Whitelegge, J., Jiang, L., Riek, R., and Eisenberg, D. S. (2015) Uncovering the Mechanism of Aggregation of Human Transthyretin. J. Biol. Chem. 290, 28932-28943. (65) Serag, A. A., Altenbach, C., Gingery, M., Hubbell, W. L., and Yeates, T. O. (2001) Identification of a subunit interface in transthyretin amyloid fibrils: evidence for self-assembly from oligomeric building blocks. Biochemistry 40, 9089-9096.
ACS Paragon Plus Environment
34
Page 35 of 36
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60
Biochemistry
(66) Chiti, F., and Dobson, C. M. (2009) Amyloid formation by globular proteins under native conditions. Nat. Chem. Biol. 5, 15-22. (67) Jahn, T. R., Parker, M. J., Homans, S. W., and Radford, S. E. (2006) Amyloid formation under physiological conditions proceeds via a native-like folding intermediate. Nat. Struct. Mol. Biol. 13, 195201. (68) Bemporad, F., and Chiti, F. (2009) "Native-like aggregation" of the acylphosphatase from Sulfolobus solfataricus and its biological implications. FEBS Lett. 583, 2630-2638. (69) Nordlund, A., and Oliveberg, M. (2006) Folding of Cu/Zn superoxide dismutase suggests structural hotspots for gain of neurotoxic function in ALS: parallels to precursors in amyloid disease. Proc. Natl. Acad. Sci. U. S. A. 103, 10218-10223.
ACS Paragon Plus Environment
35
Biochemistry
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60
Page 36 of 36
Pathogenic Mutations Induce Partial Structural Changes
in
Native
β-Sheet
Structure
of
Transthyretin and Accelerate Aggregation Kwang Hun Lim*,†, Anvesh K. R. Dasari†, Renze Ma , Ivan Hung‡, Zhehong Gan‡, Jeffery W. Kelly§, Michael C. Fitzgerald
For Table of Contents Use Only
ACS Paragon Plus Environment
36