Next Article in Journal
Structural Analysis of the Michael-Michael Ring Closure (MIMIRC) Reaction Products
Next Article in Special Issue
The RHIM of the Immune Adaptor Protein TRIF Forms Hybrid Amyloids with Other Necroptosis-Associated Proteins
Previous Article in Journal
Identification of Key Aroma Compounds Responsible for the Floral Ascents of Green and Black Teas from Different Tea Cultivars
Previous Article in Special Issue
Molecular Dynamics Simulation Studies on the Aggregation of Amyloid-β Peptides and Their Disaggregation by Ultrasonic Wave and Infrared Laser Irradiation
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Reverse Engineering Analysis of the High-Temperature Reversible Oligomerization and Amyloidogenicity of PSD95-PDZ3

1
Department of Biotechnology and Life Science, Tokyo University of Agriculture and Technology, 2-24-16, Naka-cho, Koganei-shi 184-8588, Tokyo, Japan
2
Institute of Global Innovation Research, Tokyo University of Agriculture and Technology, 3-8-1, Harumi-cho, Fuchu-shi 183-8538, Tokyo, Japan
3
Department of Bioengineering, Nagaoka University of Technology, 1603-1, Kamitomioka-cho, Nagaoka-shi 940-2188, Niigata, Japan
4
Graduate School of Medical Life Science, Yokohama City University, 1-7-29 Suehiro, Yokohama 230-0045, Kanagawa, Japan
5
Department of Physical Chemistry, Institute of Biotechnology, Faculty of Sciences, University of Granada, 18071 Granada, Spain
*
Author to whom correspondence should be addressed.
These authors contributed equally to this work.
Molecules 2022, 27(9), 2813; https://doi.org/10.3390/molecules27092813
Submission received: 14 March 2022 / Revised: 9 April 2022 / Accepted: 18 April 2022 / Published: 28 April 2022

Abstract

:
PSD95-PDZ3, the third PDZ domain of the post-synaptic density-95 protein (MW 11 kDa), undergoes a peculiar three-state thermal denaturation (N In  D) and is amyloidogenic. PSD95-PDZ3 in the intermediate state (I) is reversibly oligomerized (RO: Reversible oligomerization). We previously reported a point mutation (F340A) that inhibits both ROs and amyloidogenesis and constructed the PDZ3-F340A variant. Here, we “reverse engineered” PDZ3-F340A for inducing high-temperature RO and amyloidogenesis. We produced three variants (R309L, E310L, and N326L), where we individually mutated hydrophilic residues exposed at the surface of the monomeric PDZ3-F340A but buried in the tetrameric crystal structure to a hydrophobic leucine. Differential scanning calorimetry indicated that two of the designed variants (PDZ3-F340A/R309L and E310L) denatured according to the two-state model. On the other hand, PDZ3-F340A/N326L denatured according to a three-state model and produced high-temperature ROs. The secondary structures of PDZ3-F340A/N326L and PDZ3-wt in the RO state were unfolded according to circular dichroism and differential scanning calorimetry. Furthermore, PDZ3-F340A/N326L was amyloidogenic as assessed by Thioflavin T fluorescence. Altogether, these results demonstrate that a single amino acid mutation can trigger the formation of high-temperature RO and concurrent amyloidogenesis.

1. Introduction

The two-state thermal denaturation process is a biophysical hallmark for a natively folded single-domain globular protein. A two-state thermal denaturation process exhibits a sharp endothermic peak as observed by micro-calorimetry [1,2,3,4], and the two-state unfolding can be formally confirmed by thermodynamically analyzing the heat capacity from differential scanning calorimetry (DSC) [5]. Exceptions to the two-state thermal unfolding are observed when a molten globule (MG) state forms upon thermal denaturation [6,7,8,9,10,11], usually under non-physiological conditions (acidic/high salt conditions) [8]. The equilibrium MG is a state where the secondary structure is retained, but the tertiary structure is loosely packed, similar to the kinetic intermediates observed during protein folding [7].
Post-synaptic density-95 protein (PSD-95) is a member of the membrane-associated guanylate kinase (MAGUK) family. Like other MAGUK proteins, PSD-95 consists of three PDZ domains, one SH3 domain, and one guanylate kinase [12,13,14,15]. PDZ3 is the third PDZ domain of PSD-95, containing three α-helices and six β-strands. It is small, globular, and it has a molecular weight of 11 kDa. PDZ3 undergoes a three-state denaturation process. The intermediate state is not a MG state but is an oligomer formed reversibly at temperatures as high as 60 to 70 °C, which we coined high-temperature reversible oligomer (RO) [16,17,18,19].
Recently, we found that the single mutations (F340A and L342A), which do not affect the physicochemical properties of PDZ3 at ambient temperatures, could inhibit the formation of RO and amyloids [16]. Both variants (PDZ3-F340A and PDZ3-L342A) undergoing a two-state unfolding process were designed by replacing hydrophobic residues at the interface of the tetrameric crystal structure of PDZ3 to alanine. Namely, the residues were identified by their large buried surface area (BSA) and accessible surface area (ASA). In particular, F340A mutation inhibited not only high-temperature RO but also amyloidogenesis; however, it would be of interest to understand how a single mutation could induce the formation of an RO state at a high-temperature.
This study applies a reverse engineering strategy to design point mutations to reintroduce high-temperature ROs in PDZ3-F340A. Reverse engineering test our hypothesis on the mechanisms underlying the formation of high-temperature RO in natural sequences. The PDZ3-F340A, which is RO-free, was used as a template protein. We hypothesized RO is produced by a hydrophobic residue on the surface of the monomeric protein and buried at the interface of the tetrameric structure. We thus selected three such hydrophilic residues and replaced them with Leucine (Leu) to enhance the hydrophobic interaction between the monomeric proteins. DSC analysis indicated that two variants unfolded according to the two-state model and did not form high-temperature ROs; however, RO was successfully reintroduced by the N326L mutation, and PDZ3-F340A/N326L was strongly amyloidogenic, confirming a correlation between the appearance of RO and amyloidogenesis.

2. Results and Discussion

2.1. Reverse Engineering of PDZ3

This study used PDZ3-F340A, which undergoes a two-state thermal denaturation as a template for examining whether one can find point mutations yielding RO at high temperatures. PDZ3-F340A itself was designed by mutating F340A in the wild-type of PSD95-PDZ3, which undergoes a three-state thermal denaturation with the concurrent production of high-temperature RO. The RO-producing mutations are assumed to be hydrophilic residues located on the surface of PDZ3 but at the interface of the tetrameric structural unit. The mutation site was determined by calculating the ASA and RSA using DSSP and by applying the following rules (i) high total monomeric accessible surface area (ASA) and tetrameric buried surface area (BSA), (ii) high total tetrameric relative solvent accessible area (RSA), (iii) non-hydrophobic residue, and, (iv) hydrophobic interaction between a side-chain without a steric clash (Table 1 and Table S1). We thus computed the ASA and BSA of PDZ3-F340A, modeled from the X-ray structure of PSD95-PDZ3-wt (PDB ID: 3I4W) using COOT. A large ASA indicates that the residue is located on the surface of the monomeric protein, while a large BSA indicates that the residue is in the interface of the tetrameric unit cell. We selected three hydrophilic residues based on their RSA and BSA (Figure 1a), namely R309, E310, N326, and visually confirmed their location using PyMol (Figure S1, see Supplementary Materials). Note that the calculation from DSSP was in line with those from PDBePISA [20], which we used in our previous reports [16,21]. Three variants, where the residues mentioned above were individually substituted to a Leucine (Leu) with the aim of inducing high-temperature RO by increasing hydrophobic interactions, were produced in E. coli, purified, and characterized as described in the following sections.

2.2. Biophysical Characterization of the PDZ3 Variants

The substitution of the candidate residues to Leu did not affect PDZ3 variants’ native structure and physiochemical properties at ambient temperatures. Sedimentation velocity analysis indicated that all PDZ3 variants were monomeric at 25 °C (Table S2, Figure S2). CD spectra of PDZ3 variants showed that the secondary structure contents of PDZ3 variants were identical to that of PDZ3-wt at temperatures up to 60 °C (Figure S3), and the denaturation was reversible as assessed by measuring the spectra after cooling the heated sample to 25 °C. The CD spectra of PDZ3 variants contained similar fractions of antiparallel and parallel β-strands more than α-helix structures in line with the secondary structure content calculated from the X-ray structure of PDZ3-wt (Table S3). The thermal denaturation curves measured by CD at a concentration of 0.5 mg/mL and pH 7.5 were sigmoidal (Figure 2), indicating an apparent two-state denaturation of the secondary structures (ND). For discussion, we estimated the apparent thermodynamic parameters, and we calculated the melting temperature (Tm) and van’t Hoff enthalpy (ΔHvan’t Hoff (Tm)) from the fitting of the CD denaturation curves (Table 2). The apparent Tm of PDZ3-F340A/N326L was the highest at 72.40 °C, whereas R309L and E310L slightly decreased the apparent Tm.

2.3. DSC Analysis and Thermodynamic Parameters

The DSC thermograms of reversely engineered PDZ3 variants were measured at a 0.5–1 mg/mL concentration in pH 7.5 with a +1 °C/min scan rate. A single endothermic peak was observed in DSC thermograms of PDZ3-F340A, as reported earlier [16], and for PDZ3-F340A/R309L and PDZ3-F340A/E310L. On the other hand, PDZ3-F340A/N326L exhibited two distinct endothermic peaks similar to those observed for PDZ3-wt (Figure 3).
A detailed analysis of the DSC curves was performed using DDCL3 with a two- and a three-state model (Table 3, Figure S4). The template PDZ3-F340A, PDZ3-F340A/E310L, and PDZ3-F340A/R309L were well fitted with a two-state model (ND). Global fitting of curves with protein concentrations of 0.5–1.0 mg/mL showed that PDZ3-F340A/N326L and PDZ3-wt formed tetrameric and pentameric ROs (N1/4(I4)D and N1/5(I5)D), respectively (Figure 4 and Figure S5). In addition, we observed a strong correlation between Tmid (N In+D) and the apparent Tm determined by CD (Figure 5 and Figure S6), but not between Tmid (N + InD) and Tm. This result suggests that the secondary structures of PDZ3-F340A/N326L and PDZ3-wt in the intermediate state is unfolded.

2.4. High-Temperature ROs and Amyloidogenesis of the Variants

To gain insight into the concurrent formation of high-temperature ROs with amyloidogenicity in PDZ3, we monitored the ThT and ANS fluorescence (Figure 6 and Figure S7). We used ThT to monitor amyloidogenesis because there is a strong relationship between beta-cross structures’ formation in fiber and oligomer forms [24,25,26]. ANS indicates molten globule-like properties and binds to partially exposed hydrophobic surfaces and cavities but also to aggregates with molten globule-like properties [24,27,28].
The ThT fluorescence intensity of PDZ3-F340A/N326L upon incubation at pH 7.5 and 1mg/mL at 60 °C and 70 °C for 3 h (Figure 6a and Figure S7a) increased within 5 min and became 5 times higher than that of PDZ3-wt. On the other hand, the ThT fluorescence of PDZ3-F340A/E310L, PDZ3-F340A/R309L, and PDZ3-F340A was essentially negligible. Thus, a single mutation, N326L, not only induced RO but strongly increased amyloidogenicity. For the purpose of discussion, let us note that the simultaneous increase of ANS and ThT fluorescence was also observed for Lysozyme [29].
ANS fluorescence of PDZ3-F340A/N326L slightly increased while PDZ3-wt increased within 10 min when incubated at 60 °C and 70 °C at pH 7.5 (Figure 6b and Figure S7b). After 3 h of incubation, the fluorescent intensity of PDZ3-F340A/N326L was 2.5-fold higher than the PDZ3-wt signal. The ANS fluorescence of the other variants was small. A similar phenomenon appeared at 70 °C incubation, but the fluorescence intensity of PDZ3-F340A/N326L and PDZ3-wt rapidly increased at the beginning and with higher intensity than at 60 °C incubation. This result suggests that PDZ3-F340A/N326L in the RO state has a molten globule-like property.
Finally, we measured the hydrodynamic radii (Rh) of the PDZ3 oligomers at temperatures from 25 °C to 90 °C (Figure S8). First, the Rh of PDZ3 variants were two to three-times larger at high than ambient temperatures (Table 4). In contrast, PDZ3-F340A/N326L showed a hydrodynamic radius (Rh) of 6.74 ± 0.13 nm at 70 °C, and the oligomerization was not fully reversible, as assessed by Rh measured after cooling the sample down to 25 °C. However, at a protein concentration of 0.5 mg/mL, the oligomerization was fully reversible (Figure S9 and Table S4), strongly suggesting the reversibility of the oligomer formed by N326L mutant, which we defined as a basic property of the RO state (R stands for reversible).

2.5. The N326L Mutation Induces RO and Amyloidogenesis

Despite being still rare, the number of high-temperature RO observations is gradually increasing, albeit under non-physiological conditions (Cytochrome c) [30], artificially redesigned protein for controlling solubility (tagged BPTI) [19], or globular domains with natural sequences (Dengue4 envelope domain 3 (DEN4ED3) [31]. In addition, we previously showed that the intermediate state (or RO) could be fully inhibited by replacing a single hydrophobic residue at the crystal interface with an alanine [16,21,31], which correlated with the inhibition of amyloidogenicity.
Here, we used a reversed engineering strategy to assess whether increasing the hydrophobicity of the crystal interface induces RO and concurrent amyloidogenicity. We constructed three variants and confirmed that their structure and physicochemical properties were conserved at low and ambient temperatures (Figures S1 and S6, and Table S5). Our experiment unambiguously indicated that one of the variants, PDZ3-F340A/N326L undergoes a three-state denaturation (N 1/4(I4) D), indicating that a single mutation can induce the formation of RO at high-temperature. In contrast, PDZ3-F340A/E310L and PDZ3-F340A/R309L do not induce RO as shown by a single endothermic peak in DSC thermograms and undergo a two-state denaturation (ND) like our template PDZ3-F340A.
It remains unclear why the N326L mutation-induced RO, but not R309L nor E310L, despite all three mutations being located at the crystallographic tetramer interface. Further inspection indicated that N326 is located on the β2 strand of PDZ3-F340A, whereas E310 and R309 are in a loop close to the β1-strand. In addition, the aggregation-prone region of PDZ3-F340A calculated by TANGO [32] indicated that N326 is placed in the 323–328 β2 strand, which has a high amyloidogenicity tendency (Figure 1b). On the other hand, R309 and E310 are in nonamyloidogenic regions. These regions were confirmed in PDZ3-wt by heteronuclear NMR experiments, being β2 engaged in the intermediate arrangement of the unfolding intermediate, whereas the β1 region does not [33]. The F340A and L342A mutations, which abolished RO and amyloidogenicity in the PDZ3-wt, were also located in the 335–343 region crossing the β3 strand amyloidogenic regions (data not shown). Thus, although E310 and R309 are located at the interface, which we assumed would increase the hydrophobic interaction and induce RO formation, they did not induce RO nor amyloidogenicity. In addition, the accessible surface area model created by PyMOL showed that the 323–328 region is buried when proteins arrange into the tetramer, but the 335–343 region is not (Figure 1c).

3. Materials and Methods

3.1. Protein Expression, Purification, and Identification

The proteins were prepared according to our previously reported protocol. In short, single mutations were introduced using a Quikchange protocol and a synthetic gene encoding PDZ3 cloned into a pBAT4 vector as the template. All variants were overexpressed in Escherichia coli strain BL21(DE3) with 1 L of LB medium. Protein expression was induced by adding 0.2 mM IPTG when the OD at 590 nm reached 0.6, and the culture was further incubated at 37 °C, 120 rpm for 4 h. The harvested cells were dissolved in 20 mL of 50 mM Tris-HCl (pH 8.7) and lysed by ultrasonication. The supernatant fraction of the cell lysate was then acidified to pH 3 by adding around 1 mL of 1 M HCl and ultracentrifuged. The recombinant proteins were purified from the supernatant by reverse-phase HPLC, lyophilized, and stored at −30 °C until use, as we reported in previous reports [34,35].
The molecular weight of the protein was confirmed by matrix-assisted laser desorption/ionization-time of flight (MALDI-TOF) MS measurements using the plate with Autoflex speed TOF/TOF (Bruker Daltonics, Fremont, CA, USA). The matrix solution was prepared by dissolving 10 mg of sinapic acid in 1 mL of a solution containing 300 μL of acetonitrile, 100 μL of 1% trifluoroacetic acid, and 600 μL of Milli Q water. Protein samples were prepared by mixing 1 μL of protein solution with 9 μL of the matrix solution. 1 μL of 10 µM, 1 µM, and 0.1 µM sample mixtures were spotted and air-dried on the MALDI-TOF MS plate. The molecular weights of all proteins were within 7 Da of the theoretical values (Table S6) computed with ExPASy’s ProtParam tool (https://web.expasy.org/protparam/ accessed on 30 May 2020).
Samples were prepared by dissolving lyophilized proteins in Milli Q water, and the protein concentrations of samples were adjusted to 0.2, 0.5, and 1.0 mg/mL in 50 mM potassium phosphate buffer (pH 7.5). The protein concentrations were determined by measuring absorbance at 280 nm (ε = 2980 M−1 cm−1) using a Nanodrop (Thermo Fisher Scientific, Waltham, MA, USA). The pH of the samples was confirmed just before performing the experiments. Freshly prepared samples were used for CD, DLS, and fluorescence spectroscopy measurements.

3.2. Differential Scanning Calorimetry (DSC) Measurements

Samples were prepared by dissolving lyophilized proteins in Milli Q water and dialyzed for 18 h at 4 °C in 50 mM potassium phosphate buffer (pH 7.5) using a Spectra/Por 3 membrane (MWCO of 3.5 kDa) with one buffer exchange. After dialysis, the protein concentration of samples was adjusted to 1 mg/mL and filtered with Disposable Ultrafiltration Unit, 200 K MWCO (ADVANTEC®, Tokyo, Japan), to remove aggregates. Protein concentrations and pHs of the samples were confirmed just before performing experiments.
DSC measurements were performed using a VP-DSC microcalorimeter (Malvern Panalytical Ltd, Malvern, UK) at a scan rate of +1.0 °C/min in the temperature range of 20 to 100 °C, essentially in line with our previous reports [36,37]. Baselines were recorded before measurements using a 50 mM potassium phosphate buffer (pH 7.5). The reversibility of the thermal unfolding was checked by repeating scans of the same sample. Thermodynamic parameters (Tmid and ΔH(Tmid)) were determined by analyzing the apparent heat capacity curves using a non-linear, least-squares fitting algorithm, DDCL3, and assuming a linear temperature dependence of the heat capacity of the native and denatured states [5,38].

3.3. Circular Dichroism (CD) Measurements

CD measurements was performed using a Jasco-J820 spectropolarimeter (Tokyo, Japan). A quartz cuvette with 2-mm optical path length was used. The secondary structure contents from CD spectra were calculated using BeStSel [39]. Thermal stability was measured at a protein concentration of 0.5 mg/mL in 50 mM potassium phosphate buffer (pH 7.5), at a +1.0 °C/min scan rate, and monitored between 25 °C and 90 °C using the CD value at 220 nm. Melting temperatures (Tm) were computed through least-squares fittings of experimental data to a two-state model using Origin 2020b (OriginLab Corp, Northampton, MA, USA) [40].

3.4. Dynamic Light Scattering (DLS) Measurements

DLS measurements were performed using a glass cuvette with a Zeta-nanosizer (Nano S, Malvern, UK). The sample was measured at 25–90 °C and reversed to 25 °C. The hydrodynamic radius (Rh) was calculated using the Stokes-Einstein equation from size-number plots [41].

3.5. Fluorescence Spectroscopy Measurement

ANS fluorescence was measured at an excitation wavelength of 380 nm at 60, 70 °C for 53 h. The emission spectra were monitored from 400 to 600 nm. ThT fluorescence was measured with an excitation wavelength of 444 nm, and the emission spectra were observed from 460 to 640 nm. The final concentration of Thioflavin T (ThT) and 8-Anilino-1-naphthalenesulfonate (ANS) was 12 and 20 μM, respectively. The dye was mixed with 300 μL of the samples, from which 60 μL were transferred to a 3.00 mm Hellma® micro cuvette.

3.6. Analytical Ultracentrifugation (AUC) Measurements

Samples were prepared by dissolving lyophilized proteins in Milli Q water and dialyed for 18 h at 4 °C against 50 mM potassium phosphate buffer pH 7.5. The protein concentrations of the samples were adjusted to 1 mg/mL by diluting with the dialyzed buffer. Samples were filtered with a 0.20 μm membrane filter (MilliporeSigma, Burlington, VT, USA) for removing aggregates. Protein concentrations and pH values of the samples were confirmed just before performing the experiments.
Sedimentation velocity experiments were carried out using an Optima XL-I analytical ultracentrifuge (Beckman-Coulter) with An-50 Ti analytical 8-place titanium rotor at 25 °C. Samples were transferred to a 12-mm double-sector epon charcoal-filled centerpiece and centrifuged at a rotor speed of 50,000 rpm, and the absorbance was monitored at 280 nm. Sedimentation velocity data were analyzed using the continuous distribution c(s) analysis module in the SEDFIT software [42]. The range of sedimentation coefficients, where the main peak was present, was integrated to obtain the weighted average sedimentation coefficient. The c(s) distribution was converted into c(M), a molar mass distribution. Solvent density, viscosity, and protein partial specific volumes were calculated using SEDTERP [43].

4. Conclusions

The reverse engineering strategy confirmed that hydrophobic interaction at the interface of the monomeric unit in the crystal tetramer induces the RO formation at a high-temperature. Thermodynamic analysis of the DSC denaturation curves indicated that the molar fraction of RO becomes maximal at 60 to 70 °C and that the secondary structures of PDZ3-F340A/N326L and PDZ3-wt in the intermediate state (RO state) are unfolded. Furthermore, the reversed engineering strategy confirmed the relationship between RO appearance and amyloidogenicity. Kinetic experiments need to establish whether ROs are on or off-pathway for amyloidogenesis, but the present results strongly favor the on-pathway hypothesis.

Supplementary Materials

The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/molecules27092813/s1, Figure S1: Accessible surface models of PDZ3 variants and PDZ3-wt by Pymol; Figure S2: Sedimentation velocity analysis of PDZ3 variants by AUC measurements; Figure S3: CD spectra of PDZ3 variants; Figure S4: Concentration dependence of DSC thermograms of PDZ3 variants fitting by DDCL3 analysis; Figure S5: Molar fraction of PDZ3 variants calculated by DDCL3 analysis of DSC thermograms; Figure S6: Fitting curve calculated by CD thermal denaturation superimposed to the molar fraction of In + D state calculated from DSC thermograms; Figure S7: Time-course fluorescence spectra of PDZ3 variants at 70°C by monitoring ThT and ANS fluorescence; Figure S8: Hydrodynamic radii (Rh) of PDZ3 variants using DLS measurements; Figure S9: Bar plot of Hydrodynamic radius (nm, Rh) by DLS measurements; Table S1: Accessible Surface Area values of artificial crystallographic PDZ3-F340A calculated by DSSP; Table S2: Molecular weight of PDZ3 variants determined by MALDI-TOF MS; Table S3: Secondary structure contents of PDZ3 variants calculated by BeStSel; Table S4: Sedimentation velocity analysis of PDZ3 variants at 25°C analyzed using SEDFIT and SEDNTERP; Table S5: Residues of PDZ3 variants between DSC raw data and fitting curves calculated with DDCL3; Table S6: Hydrodynamic radius (nm, Rh) by DLS.

Author Contributions

S.O., T.S. and Y.K. designed the study and wrote the manuscript. S.O. designed and purified the recombinant proteins and performed the spectroscopic measurements. T.S. and S.-i.K. performed DSC measurements and the data analysis. S.O., K.M. and J.R.H.T. carried out AUC analysis. J.C.M. provided the materials and co-wrote the manuscript. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by a MEXT scholarship to SO, a JSPS grant-in-aid for scientific research (KAKENHI: 18H02385, 21K05288, and 21K15049) and the TUAT’s Institute of Global Innovation Research.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Data is contained within the article or Supplementary Materials.

Acknowledgments

We are grateful to Brindha Subbaian and Kuroda lab members for suggestions with site-directed mutagenesis and technical assistance.

Conflicts of Interest

The authors declare no conflict of interest.

Sample Availability

The PDZ3-wt plasmid is available from JCM.

References

  1. Privalov, P.L. Stability of proteins: Small globular proteins. Adv. Protein Chem. 1979, 33, 167–241. [Google Scholar] [CrossRef] [PubMed]
  2. Pfeil, W. Thermodynamics of alpha-lactalbumin unfolding. Biophys. Chem. 1981, 13, 181–186. [Google Scholar] [CrossRef]
  3. Privalov, P.L.; Gill, S.J. Stability of protein structure and hydrophobic interaction. Adv. Protein Chem. 1988, 39, 191–234. [Google Scholar] [CrossRef]
  4. Makhatadze, G.I.; Privalov, P.L. Protein interactions with urea and guanidinium chloride. A calorimetric study. J. Mol. Biol. 1992, 226, 491–505. [Google Scholar] [CrossRef]
  5. Kidokoro, S.-I.; Wada, A. Determination of thermodynamic functions from scanning calorimetry data. Biopolymers 1987, 26, 213–229. [Google Scholar] [CrossRef]
  6. Ohgushi, M.; Wada, A. ‘Molten-globule state’: A compact form of globular proteins with mobile side-chains. FEBS Lett. 1983, 164, 21–24. [Google Scholar] [CrossRef] [Green Version]
  7. Kuwajima, K. The molten globule state as a clue for understanding the folding and cooperativity of globular-protein structure. Proteins Struct. Funct. Bioinform. 1989, 6, 87–103. [Google Scholar] [CrossRef]
  8. Kuroda, Y.; Kidokoro, S.-I.; Wada, A. Thermodynamic characterization of cytochrome c at low pH: Observation of the molten globule state and of the cold denaturation process. J. Mol. Biol. 1992, 223, 1139–1153. [Google Scholar] [CrossRef]
  9. Bychkova, V.E.; Semisotnov, G.V.; Balobanov, V.A.; Finkelstein, A.V. The molten globule concept: 45 years later. Biochemistry 2018, 83, S33–S47. [Google Scholar] [CrossRef]
  10. Kuroda, Y.; Endo, S.; Nakamura, H. How a novel scientific concept was coined the “molten globule state”. Biomolecules 2020, 10. [Google Scholar] [CrossRef] [Green Version]
  11. Vassilenko, K.S.; Uversky, V.N. Native-like secondary structure of molten globules. Biochim. Biophys. Acta 2002, 1594, 168–177. [Google Scholar] [CrossRef]
  12. Chen, X.; Levy, J.M.; Hou, A.; Winters, C.; Azzam, R.; Sousa, A.A.; Leapman, R.D.; Nicoll, R.A.; Reese, T.S. PSD-95 family MAGUKs are essential for anchoring AMPA and NMDA receptor complexes at the postsynaptic density. Proc. Natl. Acad. Sci. USA 2015, 112, E6983. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  13. Murciano-Calles, J.; Marin-Argany, M.; Cobos, E.S.; Villegas, S.; Martinez, J.C. The impact of extra-domain structures and post-translational modifications in the folding/misfolding behaviour of the third PDZ domain of MAGUK neuronal protein PSD-95. PLoS ONE 2014, 9, e98124. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  14. Savioz, A.; Leuba, G.; Vallet, P.G. A framework to understand the variations of PSD-95 expression in brain aging and in Alzheimer’s disease. Ageing Res. Rev. 2014, 18, 86–94. [Google Scholar] [CrossRef] [PubMed]
  15. Zhang, J.; Petit, C.M.; King, D.S.; Lee, A.L. Phosphorylation of a PDZ domain extension modulates binding affinity and interdomain interactions in postsynaptic density-95 (PSD-95) protein, a membrane-associated guanylate kinase (MAGUK). J. Biol. Chem. 2011, 286, 41776–41785. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  16. Saotome, T.; Mezaki, T.; Brindha, S.; Unzai, S.; Martinez, J.C.; Kidokoro, S.I.; Kuroda, Y. Thermodynamic analysis of point mutations inhibiting high-temperature reversible oligomerization of PDZ3. Biophys. J. 2020, 119, 1391–1401. [Google Scholar] [CrossRef]
  17. Murciano-Calles, J.; Cobos, E.S.; Mateo, P.L.; Camara-Artigas, A.; Martinez, J.C. An oligomeric equilibrium intermediate as the precursory nucleus of globular and fibrillar supramacromolecular assemblies in a PDZ domain. Biophys. J. 2010, 99, 263–272. [Google Scholar] [CrossRef] [Green Version]
  18. Murciano-Calles, J.; Cobos, E.S.; Mateo, P.L.; Camara-Artigas, A.; Martinez, J.C. A comparative analysis of the folding and misfolding pathways of the third PDZ domain of PSD95 investigated under different pH conditions. Biophys. Chem. 2011, 158, 104–110. [Google Scholar] [CrossRef] [Green Version]
  19. Nakamura, S.; Kibria, M.G.; Unzai, S.; Kuroda, Y.; Kidokoro, S.I. Reversible oligomerization and reverse hydrophobic effect induced by Isoleucine tags attached at the C-terminus of a simplified BPTI variant. Biochemistry 2020, 59, 3660–3668. [Google Scholar] [CrossRef]
  20. Krissinel, E.; Henrick, K. Inference of macromolecular assemblies from crystalline state. J. Mol. Biol 2007, 372, 774–797. [Google Scholar] [CrossRef]
  21. Saotome, T.; Onchaiya, S.; Brindha, S.; Mezaki, T.; Unzai, S.; Noguchi, K.; Martinez, J.C.; Kidokoro, S.I.; Kuroda, Y. Blocking PSD95-PDZ3’s amyloidogenesis through point mutations that inhibit high-temperature reversible oligomerization (RO). FEBS J. 2021. [Google Scholar] [CrossRef] [PubMed]
  22. Emsley, P.; Lohkamp, B.; Scott, W.G.; Cowtan, K. Features and development of Coot. Acta Cryst. D Biol. Cryst. 2010, 66, 486–501. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  23. Tien, M.Z.; Meyer, A.G.; Sydykova, D.K.; Spielman, S.J.; Wilke, C.O. Maximum allowed solvent accessibilities of residues in proteins. PLoS ONE 2013, 8, e80635. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  24. Hawe, A.; Sutter, M.; Jiskoot, W. Extrinsic fluorescent dyes as tools for protein characterization. Pharm. Res. 2008, 25, 1487–1499. [Google Scholar] [CrossRef] [Green Version]
  25. Dasgupta, M.; Kishore, N. Characterization and analysis of binding of Thioflavin T with partially folded and native states of alpha-lactalbumin protein by calorimetric and spectroscopic techniques. Int. J. Biol. Macromol. 2017, 95, 376–384. [Google Scholar] [CrossRef]
  26. Krebs, M.R.; Bromley, E.H.; Donald, A.M. The binding of thioflavin-T to amyloid fibrils: Localisation and implications. J. Struct. Biol. 2005, 149, 30–37. [Google Scholar] [CrossRef]
  27. Guliyeva, A.J.; Gasymov, O.K. ANS fluorescence: Potential to discriminate hydrophobic sites of proteins in solid states. Biochem. Biophys. Rep. 2020, 24, 100843. [Google Scholar] [CrossRef]
  28. Kibria, M.G.; Akazawa-Ogawa, Y.; Rahman, N.; Hagihara, Y.; Kuroda, Y. The immunogenicity of an anti-EGFR single domain antibody (VHH) is enhanced by misfolded amorphous aggregation but not by heat-induced aggregation. Eur. J. Pharm. Biopharm. 2020, 152, 164–174. [Google Scholar] [CrossRef]
  29. Sulatsky, M.I.; Sulatskaya, A.I.; Povarova, O.I.; Antifeeva, I.A.; Kuznetsova, I.M.; Turoverov, K.K. Effect of the fluorescent probes ThT and ANS on the mature amyloid fibrils. Prion 2020, 14, 67–75. [Google Scholar] [CrossRef]
  30. Nakamura, S.; Seki, Y.; Katoh, E.; Kidokoro, S.-I. Thermodynamic and structural properties of the acid molten globule state of horse cytochrome c. Biochemistry 2011, 50, 3116–3126. [Google Scholar] [CrossRef]
  31. Saotome, T.; Nakamura, S.; Islam, M.M.; Nakazawa, A.; Dellarole, M.; Arisaka, F.; Kidokoro, S.; Kuroda, Y. Unusual reversible oligomerization of unfolded dengue envelope protein domain 3 at high temperatures and its abolition by a point mutation. Biochemistry 2016, 55, 4469–4475. [Google Scholar] [CrossRef]
  32. Fernandez-Escamilla, A.M.; Rousseau, F.; Schymkowitz, J.; Serrano, L. Prediction of sequence-dependent and mutational effects on the aggregation of peptides and proteins. Nat. Biotechnol. 2004, 22, 1302–1306. [Google Scholar] [CrossRef] [PubMed]
  33. Marin-Argany, M.; Candel, A.M.; Murciano-Calles, J.; Martinez, J.C.; Villegas, S. The interconversion between a flexible beta-sheet and a fibril beta-arrangement constitutes the main conformational event during misfolding of PSD95-PDZ3 domain. Biophys. J. 2012, 103, 738–747. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  34. Khan, M.A.; Islam, M.M.; Kuroda, Y. Analysis of protein aggregation kinetics using short amino acid peptide tags. Biochim. Biophys. Acta (BBA)-Proteins Proteom. 2013, 1834, 2107–2115. [Google Scholar] [CrossRef] [PubMed]
  35. Rathnayaka, T.; Tawa, M.; Nakamura, T.; Sohya, S.; Kuwajima, K.; Yohda, M.; Kuroda, Y. Solubilization and folding of a fully active recombinant Gaussia luciferase with native disulfide bonds by using a SEP-Tag. Biochim. Biophys. Acta (BBA)-Proteins Proteom. 2011, 1814, 1775–1778. [Google Scholar] [CrossRef] [PubMed]
  36. Kato, A.; Yamada, M.; Nakamura, S.; Kidokoro, S.-I.; Kuroda, Y. Thermodynamic properties of BPTI variants with highly simplified amino acid sequences. J. Mol. Biol. 2007, 372, 737–746. [Google Scholar] [CrossRef]
  37. Islam, M.M.; Sohya, S.; Noguchi, K.; Kidokoro, S.-I.; Yohda, M.; Kuroda, Y. Thermodynamic and structural analysis of highly stabilized BPTIs by single and double mutations. Proteins Struct. Funct. Bioinform. 2009, 77, 962–970. [Google Scholar] [CrossRef]
  38. Kidokoro, S.-I.; Uedaira, H.; Wada, A. Determination of thermodynamic functions from scanning calorimetry data. II. For the system that includes self-dissociation / association process. Biopolymers 1988, 27, 271–297. [Google Scholar] [CrossRef]
  39. Micsonai, A.; Wien, F.; Kernya, L.; Lee, Y.H.; Goto, Y.; Refregiers, M.; Kardos, J. Accurate secondary structure prediction and fold recognition for circular dichroism spectroscopy. Proc. Natl Acad Sci USA 2015, 112, E3095–E3103. [Google Scholar] [CrossRef] [Green Version]
  40. Greenfield, N.J. Using circular dichroism collected as a function of temperature to determine the thermodynamics of protein unfolding and binding interactions. Nat. Protoc. 2006, 1, 2527–2535. [Google Scholar] [CrossRef]
  41. Stetefeld, J.; McKenna, S.A.; Patel, T.R. Dynamic light scattering: A practical guide and applications in biomedical sciences. Biophys Rev. 2016, 8, 409–427. [Google Scholar] [CrossRef] [PubMed]
  42. Schuck, P. Size-distribution analysis of macromolecules by sedimentation velocity ultracentrifugation and lamm equation modeling. Biophys J. 2000, 78, 1606–1619. [Google Scholar] [CrossRef] [Green Version]
  43. Laue, T.M. Computer-aided interpretation of analytical sedimentation data for proteins. Analytical Ultracentrifugation in Biochemistry and Polymer Science; Royal Society of Chemistry: London, UK, 1992; pp. 90–125. [Google Scholar]
Figure 1. Secondary structure and three-dimensional accessible surface structures of PDZ3-F340A (a) Total monomeric accessible surface area (ASA), total tetrameric buried surface area (BSA), and total relative solvent accessible surface area (RSA) values of R309, E310, and N326 residues are shown in square boxes; (b) TANGO analysis of β-aggregation (amyloid) prone regions in PDZ3-F340A. The 323–328 residue region crossing the β2 strand is represented in red. The 335–342 residue region crossing the β3 strand is represented in orange. The 384–393 residue region crossing the β6 strand and α-helix 3 at C-terminus is represented in yellow; (c) Three images of 90° counter-clockwise views of the PDZ3-F340A tetrameric structure represent the aggregation-prone region on the interface of a tetramer unit cell using Pymol. The color code is the same as in the amino acid sequence.
Figure 1. Secondary structure and three-dimensional accessible surface structures of PDZ3-F340A (a) Total monomeric accessible surface area (ASA), total tetrameric buried surface area (BSA), and total relative solvent accessible surface area (RSA) values of R309, E310, and N326 residues are shown in square boxes; (b) TANGO analysis of β-aggregation (amyloid) prone regions in PDZ3-F340A. The 323–328 residue region crossing the β2 strand is represented in red. The 335–342 residue region crossing the β3 strand is represented in orange. The 384–393 residue region crossing the β6 strand and α-helix 3 at C-terminus is represented in yellow; (c) Three images of 90° counter-clockwise views of the PDZ3-F340A tetrameric structure represent the aggregation-prone region on the interface of a tetramer unit cell using Pymol. The color code is the same as in the amino acid sequence.
Molecules 27 02813 g001
Figure 2. CD thermal denaturation curves of PDZ3 variants at 0.5 mg/mL and pH 7.5 at a scan rate of +1.0 °C/min. The CD values were monitored at 220 nm. Black dots are the experimental data, and red lines are the fitting curves.
Figure 2. CD thermal denaturation curves of PDZ3 variants at 0.5 mg/mL and pH 7.5 at a scan rate of +1.0 °C/min. The CD values were monitored at 220 nm. Black dots are the experimental data, and red lines are the fitting curves.
Molecules 27 02813 g002
Figure 3. Concentration dependence of DSC thermograms of PDZ3 variants at 0.5–1 mg/mL, pH 7.5, and 1 °C/min scan rate. Black and red dots show DSC thermograms at a protein concentration of 1 and 0.5 mg/mL, respectively.
Figure 3. Concentration dependence of DSC thermograms of PDZ3 variants at 0.5–1 mg/mL, pH 7.5, and 1 °C/min scan rate. Black and red dots show DSC thermograms at a protein concentration of 1 and 0.5 mg/mL, respectively.
Molecules 27 02813 g003
Figure 4. Molar fraction of PDZ3-F340A/E310L (N-D), F340A/N326L (N-I4-D), F340A/R309L (N-D), PDZ3-F340A (N-D model) and PDZ3-wt (N-I5-D model) by DDCL3 analysis of DSC thermograms at 1 mg/mL, pH 7.5, and 1 °C/min scan rate. The lines represent natively folded monomers (blue), intermediate oligomers (green), and unfolded monomers (red).
Figure 4. Molar fraction of PDZ3-F340A/E310L (N-D), F340A/N326L (N-I4-D), F340A/R309L (N-D), PDZ3-F340A (N-D model) and PDZ3-wt (N-I5-D model) by DDCL3 analysis of DSC thermograms at 1 mg/mL, pH 7.5, and 1 °C/min scan rate. The lines represent natively folded monomers (blue), intermediate oligomers (green), and unfolded monomers (red).
Molecules 27 02813 g004
Figure 5. Correlation plots between Tm calculated by CD denaturation curve and Tmid (N↔In + D), Tmid (N + In↔D) calculated from the molar fractions as determined by DSC analysis. The CD denaturation curve and DSC measurement were measured at 0.5 mg/mL and pH 7.5 at a scan rate of +1.0 °C/min. The red circles (●) represent Tmid (N↔In+D) and the blue open triangles (∆) present Tmid (N + In↔D).
Figure 5. Correlation plots between Tm calculated by CD denaturation curve and Tmid (N↔In + D), Tmid (N + In↔D) calculated from the molar fractions as determined by DSC analysis. The CD denaturation curve and DSC measurement were measured at 0.5 mg/mL and pH 7.5 at a scan rate of +1.0 °C/min. The red circles (●) represent Tmid (N↔In+D) and the blue open triangles (∆) present Tmid (N + In↔D).
Molecules 27 02813 g005
Figure 6. Time-course of RO formation of PDZ3 variants monitored by (a) ThT fluorescence at λEm 480 nm. The dye concentration was 12 μM; (b) ANS fluorescence at λEm 480 nm. The dye concentration was 20 μM. Protein concentration was 1 mg/mL in pH 7.5. Mixture samples were measured at 60 °C for 180 min. PDZ3-F340A/E310L (blue); PDZ3-F340A/N326L (orange); PDZ3-F340A/R309L (gray); PDZ3-F340A (yellow); and PDZ3-wt (green).
Figure 6. Time-course of RO formation of PDZ3 variants monitored by (a) ThT fluorescence at λEm 480 nm. The dye concentration was 12 μM; (b) ANS fluorescence at λEm 480 nm. The dye concentration was 20 μM. Protein concentration was 1 mg/mL in pH 7.5. Mixture samples were measured at 60 °C for 180 min. PDZ3-F340A/E310L (blue); PDZ3-F340A/N326L (orange); PDZ3-F340A/R309L (gray); PDZ3-F340A (yellow); and PDZ3-wt (green).
Molecules 27 02813 g006
Table 1. ASA, BSA, and total RSA of candidate residues of PDZ3-F340A.
Table 1. ASA, BSA, and total RSA of candidate residues of PDZ3-F340A.
ResidueMonomeric ASA (Å2)Tetrameric ASA (Å2)Tetrameric BSA (Å2)Total RSA (%)
ABCDABCDABCD
GLU 3101421451401432814514046114009795
ASN 326109107110107109232110708489089
ARG 30921621822021827218220202189001675
ASN 3691541601571601541715716001430073
GLU 3311751751731741751217317401630073
THR 32176777680762088005768073
ARG 3681771501741538413117474931907970
LEU 342858784858517178507067068
ALA 343888691888847488803943064
GLY 333961009698963496980660063
SER 3201311321321311311006613103266063
ARG 313169169164172901691648179009162
PRO 311676965672069651747005061
GLU 3951501461481493814614814911200050
LYS 3801381381381401381135014002588048
GLY 3222524222525202502222042
ILE 389504849481048491140003739
VAL 3283435333334223303331037
ASP 3661089110788756810777332301135
GLN 3914344434624443941003735
Residues with ASA > 0 Å2 of at least one of the chains are listed. The monomeric and tetrameric ASA values of four polypeptide chains (A, B, C, D) of PDZ3-F340A were calculated by DSSP. PDZ3-F340A and other variants were modeled using X-ray crystallographic data of PSD95-PDZ3 (PDB ID: 3I4W) consisting of four monomeric chains in each asymmetric unit cell using COOT (crystallographic object-oriented toolkit) [22] The modeled structures were used for calculating accessible surface area (ASA) by DSSP. Buried surface area (BSA) was calculated by subtracting the ASA in the tetrameric structure from the calculated ASA of the monomeric structure. Relative solvent accessibility (RSA) is the total tetrameric ASA of each residue divided by the maximum amino acid solvent accessibility from theoretical normalization values [23]. The selected residues are underlined.
Table 2. Apparent Tm and van’t Hoff enthalpy (ΔHvan’t Hoff (Tm)).
Table 2. Apparent Tm and van’t Hoff enthalpy (ΔHvan’t Hoff (Tm)).
NameTm (°C)ΔHvan’t Hoff (Tm) (kJ/mol)
PDZ3-F340A/E310L69.5 ± 0.0375.6 ± 3.3
PDZ3-F340A/N326L72.4 ± 0.0309.5 ± 2.6
PDZ3-F340A/R309L67.2 ± 0.0363.1 ± 2.7
PDZ3-F340A71.3 ± 0.0367.4 ± 3.6
PDZ3-wt71.1 ± 0.1255.3 ± 3.2
The parameters were calculated by fitting the CD denaturation curves using a two-state model in Origin 2020b software. CD denaturation curves were measured at 0.5 mg/mL, pH 7.5, 25–90 °C and +1.0 °C/min scan rate.
Table 3. Thermodynamics parameters were calculated by DSC measurements at 0.5–1 mg/mL of protein concentration.
Table 3. Thermodynamics parameters were calculated by DSC measurements at 0.5–1 mg/mL of protein concentration.
NameConcentration (mg/mL)TransitionTmid (°C) ΔHcal (Tmid) (kJ/mol)
PDZ3-F340A/E310L1N-D68.9 ± 0.1341.8 ± 2.7
0.5N-D69.3 ± 0.1343.1 ± 1.9
PDZ3-F340A/N326L1N-I468.1 ± 0.2233.4 ± 10.1
N-D73.4 ± 1.2250.8 ± 16.8
0.5N-I470.5 ± 0.1239.0 ± 7.3
N-D72.8 ± 0.7242.7 ± 4.4
PDZ3-F340A/R309L1N-D67.0 ± 0.1364.6 ± 2.6
0.5N-D66.8 ± 0.1359.9 ± 2.7
PDZ3-F340A1N-D70.0 ± 0.1369.0 ± 2.6
0.5N-D70.2 ± 0.2358.7 ± 2.7
PDZ3-wt1N-I564.9 ± 0.1251.4 ± 3.7
N-D67.9 ± 0.3296.0 ± 4.4
0.5N-I567.8 ± 0.1250.4 ± 8.2
N-D69.5 ± 0.5316.2 ± 8.5
Midpoint temperature (Tmid), calorimetric enthalpy (ΔHcal (Tmid)) of PDZ3 variants at pH 7.5 determined by fitting the DSC thermogram using DDCL3. PDZ3-F340A/E310L (N-D model); PDZ3-F340A/N326L (N-I4-D model); PDZ3-F340A/R309L (N-D model); PDZ3-F340A (N-D model); PDZ3-wt (N-I5-D model).
Table 4. Hydrodynamic radius (nm, Rh) by DLS.
Table 4. Hydrodynamic radius (nm, Rh) by DLS.
Name25 °C40 °C60 °C70 °C80 °C90 °C25 °C (Reverse)
PDZ3-F340A/E310L1.76 ± 0.041.46 ± 0.240.99 ± 0.110.94 ± 0.551.88 ± 0.021.63 ± 0.101.58 ± 0.07
PDZ3-F340A/N326L1.15 ± 0.621.60 ± 0.030.82 ± 0.066.74 ± 0.137.94 ± 0.227.36 ± 0.985.68 ± 0.24
PDZ3-F340A/R309L1.62 ± 0.031.62 ± 0.051.70 ± 0.012.00 ± 0.052.08 ± 0.012.08 ± 0.021.48 ± 0.02
PDZ3-F340A1.54 ± 0.041.56 ± 0.021.66 ± 0.001.81 ± 0.072.23 ± 0.022.41 ± 0.131.61 ± 0.02
PDZ3-wt1.65 ± 0.001.46 ± 0.021.65 ± 0.033.61 ± 0.073.06 ± 0.032.41 ± 0.051.66 ± 0.06
DLS was measured at 1 mg/mL, pH 7.5, and 25–90 °C, and after cooling the sample back to 25 °C after heating. Rh values were calculated from size-volume graphs. The errors are the standard deviation of three-times measurements with the same sample.
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Onchaiya, S.; Saotome, T.; Mizutani, K.; Martinez, J.C.; Tame, J.R.H.; Kidokoro, S.-i.; Kuroda, Y. Reverse Engineering Analysis of the High-Temperature Reversible Oligomerization and Amyloidogenicity of PSD95-PDZ3. Molecules 2022, 27, 2813. https://doi.org/10.3390/molecules27092813

AMA Style

Onchaiya S, Saotome T, Mizutani K, Martinez JC, Tame JRH, Kidokoro S-i, Kuroda Y. Reverse Engineering Analysis of the High-Temperature Reversible Oligomerization and Amyloidogenicity of PSD95-PDZ3. Molecules. 2022; 27(9):2813. https://doi.org/10.3390/molecules27092813

Chicago/Turabian Style

Onchaiya, Sawaros, Tomonori Saotome, Kenji Mizutani, Jose C. Martinez, Jeremy R. H. Tame, Shun-ichi Kidokoro, and Yutaka Kuroda. 2022. "Reverse Engineering Analysis of the High-Temperature Reversible Oligomerization and Amyloidogenicity of PSD95-PDZ3" Molecules 27, no. 9: 2813. https://doi.org/10.3390/molecules27092813

APA Style

Onchaiya, S., Saotome, T., Mizutani, K., Martinez, J. C., Tame, J. R. H., Kidokoro, S. -i., & Kuroda, Y. (2022). Reverse Engineering Analysis of the High-Temperature Reversible Oligomerization and Amyloidogenicity of PSD95-PDZ3. Molecules, 27(9), 2813. https://doi.org/10.3390/molecules27092813

Article Metrics

Back to TopTop