BioMed Research International

BioMed Research International / 2018 / Article

Research Article | Open Access

Volume 2018 |Article ID 7480749 |

Cedrix J. Dongmo Foumthuim, Alessandra Corazza, Rodolfo Berni, Gennaro Esposito, Federico Fogolari, "Dynamics and Thermodynamics of Transthyretin Association from Molecular Dynamics Simulations", BioMed Research International, vol. 2018, Article ID 7480749, 14 pages, 2018.

Dynamics and Thermodynamics of Transthyretin Association from Molecular Dynamics Simulations

Academic Editor: Carmen Domene
Received16 Feb 2018
Accepted06 May 2018
Published05 Jun 2018


Molecular dynamics simulations are used in this work to probe the structural stability and the dynamics of engineered mutants of transthyretin (TTR), i.e., the double mutant F87M/L110M (MT-TTR) and the triple mutant F87M/L110M/S117E (3M-TTR), in relation to wild-type. Free energy analysis from end-point simulations and statistical effective energy functions are used to analyze trajectories, revealing that mutations do not have major impact on protein structure but rather on protein association, shifting the equilibria towards dissociated species. The result is confirmed by the analysis of 3M-TTR which shows dissociation within the first 10 ns of the simulation, indicating that contacts are lost at the dimer-dimer interface, whereas dimers (formed by monomers which pair to form two extended -sheets) appear fairly stable. Overall the simulations provide a detailed view of the dynamics and thermodynamics of wild-type and mutant transthyretins and a rationale of the observed effects.

1. Introduction

1.1. TTR: Origin and Function

Formerly known as prealbumin, transthyretin (TTR) is a 55 kDa globular oligomeric protein made up by four identical monomeric units (I-IV, see Figure 1(a)) each composed of 127 amino acid residues. Before being released into plasma, TTR is principally produced in liver, choroids plexus, and retina [14]. It is mainly a carrier protein that binds to retinol binding protein (RBP) to transport vitamin A in plasma and represents the leading transporter of thyroxine (T4) in cerebrospinal fluid (CSF) [57].

1.2. TTR: Structure

The 3D-conformation of all TTR-tetramers displays high degree of symmetry as shown at high resolution by Blake and coworkers [8] and Hörnberg et al. [9]. TTR is an overall -sheet protein with a small -helix domain between strands E and F (see Figure 1(a)). Its structure is highly ordered with a flexible N-terminal region that could not be resolved in many crystallographic studies. In each TTR monomer, about 45 % of residues are arranged in a sandwich immunoglobulin-like topology made up by two four antiparallel-strand -sheets; the inner sheet DAGH is opposed to the outer sheet CBEF. Its secondary structure also exhibits a short -strand portion A (see Figure 1) which is antiparallel (folded back) to strand A through a -turn (i+5) and which is involved in dimer-dimer contact [10]. The dimer of monomers I and II (Figures 1(b) and 1(c)), also known as primary dimer, is the crystallographic asymmetric unit structurally symmetric to the dimer of monomers III-IV. It is stabilized by a network of six main chain hydrogen bond interactions established by three pairs of residues (A120(NH)(OC)Y114, T118(NH)(OC)Y116, Y116(NH)(OC)T118) at the HH’ interface (Figure 2). Monomer I inner sheet DAGH forms with its homologue H’G’A’D’ in monomer II a kind of symmetric antiparallel pseudo-continuous eight-strand -sheet centered on HH’ (DAGH-H’G’A’D’), Figure 1(b) [10, 11]. At the opposite side, the pseudo-continuity in outer sheets CBEF-F’E’B’C’(Figure 1(c)) is rather loose and only 4 backbone hydrogen bonds (F87(CO)(HN)T96 and E89(NH)(OC)V94) are established (Figure 2). While dimers are mainly kept together by hydrogen bonds at monomer-monomer interface, tetramers become the dimers of dimers essentially stabilized by hydrophobic contacts established within neighbouring subunits at the dimer-dimer interface via AB-GH loop interactions involving mainly residues L17, A19, V20, L110, P113, T119, and V121.

1.3. TTR-Related Amyloidosis

TTR-related amyloidosis can be inherited in the case of genetic mutations or can be nonhereditary when it is due to wild-type (WT-TTR) [12]. WT-TTR amyloidosis also termed senile systemic amyloidosis (SSA) is the prevalent form of TTR amyloidosis which principally affects the heart. It develops with ageing and requires TTR tetramer dissociation and partial unfolding. Mutation-induced TTR amyloidosis is instead associated with familial amyloid polyneuropathy (FAP) which impacts essentially on peripheral nervous system and familial amyloid cardiomyopathy (FAC) with broad implication on heart [7, 12]. The interconnection between protein structure and stability and its capability to form amyloid fibrils have motivated several structural studies. To date, about 80 point mutations have been correlated to human inherited amyloidosis [13].

1.4. Designed Mutants

WT-TTR displays a very stable tetrameric molecular assembly. It is established that its structure is resistant to dissociation at physiological concentrations within the pH range 5-7 [14]. Several amyloid fibril formation models have been proposed and reviewed [1517]. However, there is almost a consensus that the common mechanism involves the dissociation of the native tetrameric TTR into unstable but folded monomers, followed by local unfolding of the latter into multiple nonnative amyloidogenic intermediate states that self-assemble in solution [2, 6, 1518]. Furthermore, evidence has been provided pointing to the dissociation of the native tetramer as the rate-limiting step towards aggregation [5, 1922].

Therefore knowledge of the dynamics and thermodynamics of dissociation of TTR native structure is an important issue as tetramer dissociation, monomer misfolding, and self-assembly of amyloidogenic monomers into amyloid and other aggregate morphologies are known to be linked to several human degenerative diseases. Specifically, TTR represents one of the few examples whereby it was possible to devise a drug for TTR amyloidogenesis [23] that acts through the stabilization of a native structure, thus indicating that thermodynamics studies and modeling of tetramer dissociation are a very important target.

To mimic the pathological situation, previous site-directed mutagenesis, namely, mutation F87M/L110M (MT-TTR), was carried out to promote the dissociation of TTR tetramer into monomers [5]. To further shift the equilibrium towards monomers to obtain a more homogeneous dissociated species, an additional amino acid replacement has been introduced in MT-TTR molecule (S117E) [24]. Indeed, the latter mutant containing three mutations (3M-TTR) proved to be markedly more prone to in vitro monomerization in comparison with the double mutant MT-TTR, whose tetrameric states could be recovered by Tafamidis binding, at variance with 3M-TTR [24].

The three mutations are not reported in the ClinVar database [25], although pathological mutations do occur in neighbouring residues.

1.5. TTR Amyloidosis Therapeutics

As mentioned previously, the dissociation of TTR tetramer is believed to be the initial step into its fibrillation pathway [26]. Therefore, many strategies to TTR amyloidosis prevention exploit its ability to bind small molecules in the T4 binding channel (see Figure 1(a)), mimicking its hormone binding capability, thereby producing kinetic stabilization of the tetramer [27, 28]. Other clinical remedies to FAC and FAP amyloidosis often employ organs transplantation (e.g., liver), even though not all the affected organs can be transplanted (e.g., choroids plexus, where TTR is produced as well [29]). Furthermore, additional treatments against TTR amyloidosis have been experimented, among which there are gene therapy [30] and the resorption of amyloid deposits [31].

The aim of this work is to provide a rationale of the observed dynamics and thermodynamics of the wild-type and mutant transthyretin based on atomistic molecular dynamics simulations. The methods applied here show the relative stability of monomers, dimers, and tetramers and provide a description of both enthalpic and entropic contributions taking advantage of recently developed methods.

2. Materials and Methods

2.1. Free Energy from End-Point MD Simulation

The free energy of a system, with respect to some reference state, can be in principle deduced from conformational samples obtained, e.g., in a molecular dynamics simulation.

For a well equilibrated system, the enthalpy can be estimated as an ensemble average of the potential energy over a set of conformations from the MD trajectory. The bottleneck in obtaining the Gibbs free energy is the proper modeling of the entropic part and in particular of the solvation entropic effects.

In order to obtain free energy estimates all molecular dynamics snapshots have been analyzed using the GBSA implicit solvent model as implemented in a home-written version of the software Bluues [32, 33].

The Gibbs free energy in the context of implicit solvent MD simulation has been discussed in reviews [3437]. We follow here our recent perspective on the issue [37]. The free energy of a system may be written aswhere is the solute potential, is the potential of mean forces (with the solvation energy), is the Boltzmann constant, and is the temperature.

The above expression can be rearranged and written as the sum of an entropic term involving only solute coordinates and an ensemble average including also solvent enthalpy and entropy:where is the configurational entropy of the solute.

The above formula constitutes the basis of free energy estimation from end-point simulations detailed in the following sections. The solute vacuum potential energy was obtained by selecting the solute from the trajectories and recalculating the energy involving only solute terms.

The solute potential energy and the solvation energy were calculated from molecular dynamics snapshots of the solute. First the molecular surface was generated using the program MSMS by Sanner [38]; then the vertices and normals to the surface were read by the program Bluues [32, 33] and used to compute Generalized Born radii and to compute solvation energy according to the GB model [39].

The solvent accessible surface outputted by MSMS was used to compute the apolar contribution to the solvation energy using a surface tension constant of 5 cal/( mol) [40, 41].

For entropic contributions is rewritten in terms of solute conformational probability density:where is the density in configurational space. The latter is estimated using the nearest neighbour method [4250]. The basic idea of the method is to provide a description of the configuration in a space where a metric is defined (e.g., the torsional angle space) and then use the distance to the nearest neighbour of each sample to estimate the density of probability around that sample. When this is done with the caveats discussed in depth by Demchuk and coworkers [42], the average of the estimated density at each sample provides the estimate of the configurational entropy by the equation above.

In order to cope with the high dimensionality of the space, with largely decoupled degrees of freedom, the mutual information expansion (MIE) method [44, 45] and the maximum information spanning tree method (MIST) [51, 52] were proposed.

The reader is referred to the cited literature for details on the methods.

Entropic contributions due to changes in the internal degrees of freedom were calculated using the nearest neighbour and the MIST methods as implemented in our program PDB2ENTROPY (URL: which will be described elsewhere.

Rotational and orientational entropy was computed using the nearest neighbour method in the Euclidean space approximation, which is an excellent approximation for the number of samples at hand, as described previously [50]. The calculation is implemented in our program PDB2TRENT (URL: which will be described elsewhere. The same approach has been used by us before [53].

2.2. Molecular Models

The X-ray structures used as starting configuration are the wild-type transthyretin (WT-TTR) taken from the RCSB Protein Data Bank (PDB) [54], PDB id: 1F41 [9], solved at the resolution of 1.3 Å; the double point mutant (MT-TTR) F87M/L110M (PDB id: 1GKO [5]) resolved at 2.10 Å; and the in silico engineered triple mutant F87M/L110M/S117E (3M-TTR). The structure of triple mutant was obtained from that of MT-TTR mutating the S117 by E (S117E) using the protein modeling software Swiss-PdbViewer [55]. The N-terminal residues 1-9 and C-terminals 126-127 are not present in the structure and were not modeled. Tetramers were built from deposited asymmetric units (dimers) by applying crystal symmetry operations. All the crystallization water molecules were removed prior to run the simulations.

2.3. Molecular Dynamics Simulations

A set of independent simulation runs of different length ranging from 50 to 250 ns were performed on monomer, dimer, and tetramer of wild-type and mutant structures of TTR employing either CHARMM27 all atoms force field with CMAP correction [56] or amber99sb-ildn molecular mechanics force field [57]. The protein atoms were placed at the center of a cubic box; the system was solvated using the 3-site rigid water model TIP3P [5860]. In each system an equivalent number of solvent molecules were replaced by counterions to obtain a neutral system.

The systems were first minimized using the steepest descent minimization algorithm with a minimization step size of 0.1 nm and a maximum convergence force of 1000.0 kJmol−1nm−1. The equilibration phase was done in 2 steps: 100 ps in NVT ensemble followed by 100 ps in NPT ensemble. During the first equilibration stage, the leap-frog integrator with integration time-step of 0.002 ps was used to update the changes in the system. Particle Mesh Ewald summation [61, 62] accounted for long-range electrostatics interactions. The temperature was equilibrated to a reference value of 300 K using the velocity rescaling (modified Berendsen thermostat) [63] with a coupling constant of 0.1 ps. Short-range electrostatics and van der Waals interactions were truncated with a 10 Å cutoff. All bonds were constrained with the LINCS algorithm [64]. In NPT equilibration stage, the previous parameters were still used and the pressure was stabilized to 1.0 bar using the Parrinello-Rahman pressure coupling [65, 66] with a coupling constant of 2.0 ps.

2.4. Molecular Dynamics Simulation Analysis

MD trajectories were analyzed with available structural-based tools in Gromacs-5.0.4 [41, 67]. The thermodynamic stability of the tetrameric systems was further processed using the Academic License version of the bioinformatics tool FoldX [68], using the commands Stability and AnalyseComplex, respectively, to gain information on protein stability and interacting interface free energies. In all the cases, the values were averaged over the whole simulation trajectory. AnalyseComplex command outputs the Gibbs interacting free energy of binding for a complex formation (folding), say AB (), computed as , where is the free energy of folding. Similarly the statistical effective energy function proposed by Berrera et al. [69], referred hereafter as BMF, has been used with home-written routines.

The computation of conformational entropy employs the nearest neighbour formalism whose rationale is the estimation of the local probability density around each sample by counting its number of neighbours within a hypersphere of radius equal the distance from that sample to its nearest neighbours. Thus, the values discussed in the results section are averaged for the nearest neighbour, set as default for the routines previously cited. Conformational entropies were calculated over 5000 frames.

Pictures were either collected with PyMOL [70] or VMD [71], secondary structures were assigned using DSSP program [72], and H-bonds occupancy was computed using the tool reporting only H-bonds with occupancy greater than 10 %.

3. Results and Discussion

3.1. Monomer to Dimer and Dimer to Tetramer Associations

The enthalpic contributions to the association of monomers and dimers were assessed considering both single species simulations and detaching from the tetramer simulations dimers first and then monomers. In both cases favorable protein-protein interactions appear largely overestimated by the force-field used, as recently observed also for other systems [73, 74]. The “enthalpic” contributions encompass the solvation energy which includes also the entropy of the solvent. For this reason the term enthalpy is used rather freely only with the purpose of separating the terms related to the potential of mean force and those computed from the analysis of conformational samples.

The enthalpy of two wild-type (mutant) monomers associating in a dimer computed by the ensemble averages of the potential is -49.7 (-49.8) kcal/mol and that of two dimers associating in a tetramer is -81.0 (-77.2) kcal/mol leading to a very favorable -180.4 (-176.8) kcal/mol association enthalpy for tetramerization.

The reduction in conformational entropy (i.e., changes in entropy due to internal degrees of freedom) from two monomers to dimer and from two dimers to tetramer, computed considering only correlations within each residue, is -6.7 (-12.8) e.u. and -27.2 (-24.8) e.u., respectively, thus opposing tetramerization, at 300K, by 24.3 (30.1) kcal/mol.

The rotational-translational entropy () was computed for monomer-monomer and dimer-dimer association, using the Euclidean approximation for the rotational-translational space, resulting in 19.0 (17.3) e.u. for the dimer-dimer association and 17.5 (17.5) e.u. for the monomer-monomer association.

Rotational-translational entropy thus opposes tetramerization by 32.2 (31.2) kcal/mol.

Overall, therefore, the large favorable enthalpy (likely to be overestimated here) is opposed by conformational and rotational-translational entropy which we can estimate to result in total 56.5 (61.3) kcal/mol. These results are summarized in Table 1.

(kcal mol−1)-T (kcal mol−1)-T (kcal mol−1)

WT-TTR2MD-49.7 3.74.0 0.810.4
2DT-81.0 2.916.2 0.511.3

MT-TTR2MD-49.8 2.67.6 0.810.4
2DT-77.2 6.514.8 1.110.3

The detailed analysis of the terms contributing to the enthalpy shows that covalent energy terms contribute very little to the energy difference and Lennard-Jones terms make the largest favorable contribution to association, whereas favorable electrostatic and unfavorable desolvation contributions provide an overall unfavorable but smaller (about one-third in absolute value) contribution.

The analysis provides thus a lower stability of the tetramer assembly for the mutant relative to the wild-type protein. The uncertainties of the methodology are however comparable to this difference.

3.2. Statistical Effective Energy Functions Analysis

The structural and thermodynamic stability of both mutant structures in relation to wild-type was further probed thanks to two independent bioinformatic tools, i.e., FoldX [68] (Figure 3(a)) and the statistical effective energy function BMF [69] (Figure 3(b)). This analysis provides a complementary view of the thermodynamic picture presented above.

We computed the folding free energy differences between mutants and wild-type (G=-) throughout the tetramer-unfolded monomers equilibrium, Figure 3. Some words of caution are due when presenting the following data because the energy functions used take into account implicitly the entropy loss from internal degrees of freedom but not that arising from external degrees of freedom as implied by the tetramer to dimers and the dimers to monomers transitions.

Finally, notwithstanding the fluctuations and differences in underlying principles and datasets, the two approaches used here are consistent with each other as detailed below.

The panels (A) in Figures 3(a) and 3(b) point out that both mutants have positive difference in free energy as compared to WT-TTR and thus are less stable. Their average free energy computed using BMF (FoldX) is +0.3 and +7.8 kcal/mol (+2.5 and +15.6 kcal/mol), respectively, for MT-TTR and 3M-TTR, confirming the decreased tetramer thermodynamic stability upon mutations, significant in the case of 3M-TTR. This result is consistent with previously determined experimental data, which indicated that monomers of 3M-TTR, at variance with MT-TTR, could not reform tetramers upon the addition of Tafamidis, a known stabilizer of the tetrameric form [24].

In order to track down all contributions, G’s values were computed and reported in subsets (B, C, and D) of Figures 3(a) and 3(b) for all the steps in the equilibrium pathway from tetramers to unfolded monomers (required before monomers self-assemble into amyloid fibrils) through dimers. We called , , and the free energy for the dissociation of tetramers into dimers (TD), for the formation of folded monomers from dimers (DM), and for the unfolding of individual monomers (M→), respectively. Based on simulations we guessed that tetramers first break down into dimers I/II and III/IV, i.e., along the C2 crystallographic axis, instead of I/III and II/IV. Thus, (dissociation free energy at I/II-III/IV dimer-dimer contact), (dissociation free energy at I-II and III-IV monomer-monomer contacts), and .

using BMF (FoldX) average is -1.1 and -8.7 kcal/mol (-2.5 and -26.3 kcal/mol) for MT-TTR and 3M-TTR, respectively, relative to WT-TTR. This means both equilibriums are shifted towards the right (subsets (B) in Figures 3(a) and 3(b)), i.e., the formation of dimers. The quite high (absolute) value displayed by 3M-TTR indicates the propensity of tetrameric assembly of the latter structure to dissociate into dimers I/II and III/IV as compared to the wild-type protein. Considering the dimer-monomer equilibrium (subsets (C) in Figures 3(a) and 3(b)), ’s average using BMF (FoldX) is -1.2 and -2.5 kcal/mol (-3.3 and -5.2 kcal/mol) over the trajectory for MT-TTR and 3M-TTR, respectively, relative to WT-TTR. The values are however rather limited compared to the overall stability computed using BMF for the process of four WT-TTR monomers associating into 2 dimers (-14.2 kcal/mol).

Finally the change with respect to WT-TTR in folding free energy of monomers has been computed. Both BMF and FoldX predict a shift towards the formation of more stable monomers for both mutants 0; i.e., the equilibrium is shifted towards the folded form (panels (D) in Figures 3(a) and 3(b)). The unfavorable formation of unfolded monomers in both mutants strongly correlates with the secondary structure analysis displayed in Figure 5 showing that notwithstanding the enhanced dissociation of tetramers into monomers through dimers monomers remain stable (folded). Besides, this result would strongly suggest that both monomers of our mutant variants are nonamyloidogenic, as is known as far as MT-TTR is concerned [5].

3.3. Localised Structural Transitions

Local conformational transitions were assessed by computing the dihedral phi () and psi () angles of individual residues involved in the mutation and represented as Ramachandran plots, Figure 4.

The X-ray () dihedral angles for residue 87 are (-84,-55), (-84,-40), and (-85,-41) in WT-TTR, MT-TTR, and 3M-TTR, respectively. The same values reported from simulation are as follows: (-7512, -469), (-8016, -4413), and (-9115, -3815), Figure 4(a). These latter lie in the allowed regions for -helical secondary structure, and [75], confirming that no secondary structural transition occurred at this position.

The / dihedral angles of L110 and M110 residues Figure 4(b) look even more similar than the previously observed ones. No major deviation was observed in individual monomers and the dihedral angles for residue 110 in each variant keep close to the equilibrium value. The average reported simulation values are (-12211, 1228), (-11612, 12210), and (-11711, 11910), respectively, in WT-TTR, MT-TTR, and 3M-TTR. Indeed, these fall into the region corresponding to -sheet secondary structure, and [75], and imply that no structural transition is seen in the process of point mutation L110M.

In Figure 4(c), while both TTR variants are displaying nearly the same / averages, (-13710, 1409), (-13616, 13412), and (-13117, 13511), it is seen that a distinct region is being accessed by mutant variants. The region is defined by and corresponding to left handed helix conformation. The transition is observed only once for both MT-TTR and 3M-TTR within the first 10 ns of the simulation possibly as a transient adaptation to in silico mutation of residue 110.

3.4. Secondary Structural Changes

The illustration of dynamical processes in our systems during the simulations was done by computing the secondary structure, Figure 5. The latter analysis was performed thanks to the DSSP program which determines the existence of hydrogen bonds as a criterion for the presence of secondary structure [72]. While dihedral () angles (Figure 4) were useful for assessment of local structural changes, secondary structure is able to highlight the global structural changes.

In Figure 5 we can see that the helical region exhibits large structural fluctuations and even gets disrupted from time to time, like in monomer II for all the TTR variants. F-strand undergoes some structural fluctuations at its beginning, particularly in monomers I, III, and IV. In G-H loop, some part of the structure is being converted from bend to -turn and conversely. D-strand displays quite high structural fluctuations in monomers II and IV of WT-TTR, in monomers II and III of MT-TTR, and in all monomers for 3M-TTR. A- and C-strands make the most stable parts of the molecules. E-strand in the vicinity of -helix is showing some fluctuations in monomers I, III, and IV of all the variants. In spite of structural fluctuations exhibited by major parts of the tetrameric molecular assemblies, Figure 5 clearly points out that monomers are nevertheless remaining folded; i.e., they preserve essentially their secondary structures. This observation is supported by the fact that no evident disruption and secondary structure conversion are seen in the core region (made by the -strands) of the different systems studied.

3.5. Intra- and Interchains Hydrogen Bonds Occupancy

In order to understand the mechanism of 3M-TTR dissociation we computed the H-bonds occupancy along the MD trajectory considering both main chain-main chain, main chain-side chain, and side chain-side chain hydrogen bonds, at monomer-monomer and dimer-dimer interfaces. Figure 6 summarizes the occupancies of the interactions that were identified to significantly perturb the tetrameric assembly of TTR upon mutations. It should be noted that for 3M-TTR the occupancy reflects an average of both (starting) associated and dissociated configurations.

In Figure 6, seven major H-bonds (with 80 % occupancy) can be seen almost symmetrically distributed in both dimers, three at I-II interface and four at III-IV interface. These are being lost in 3M-TTR with their average occupancy dropping below 40 %. These include the interactions Y114(CO)-A120(HN) and A120(HN)-Y114(O) (main chain-main chain) and T119()-S115() and S115()-T119() (side chain-side chain). Most of them are located in the H-strand and some in the G-H loop.

At the interface of symmetric units (I-III and II-IV) there is a significant loss in H-bond occupancy in 3M-TTR, in particular at II-IV interface. The occupancies of 22GLY(O)-122VAL(HN) and 122VAL(HN)-22GLY(O) (main chain-main chain) located in H-strand and A-B loops lie under the 10 % (Figure 6). Interestingly, it is worth mentioning that, in 3M-TTR, several side chain-side chain interactions are being formed following dissociation and subsequent temporary association. These involve many inner sheets residues like K15, R104, and E117, most of which show significant deviation from the starting WT-TTR X-ray structure.

Analysis of hydrogen bonds connectivity confirms the pivotal role of such networks in preserving the protein integrity, typically in some dedicated portions, in this case H-strand, G-H, and A-B loops. These susceptible portions of TTR that build up the monomer-monomer and dimer-dimer interfaces are definitely playing a key role in destabilization of tetrameric subunits of TTR, especially in the case of 3M-TTR. Indeed, the disappearance of H-bonds in the previously mentioned domains along the timescale of the simulations brings a possible explanation on the integrity loss of tetramers in 3M-TTR.

3.6. Mechanism of Tetramer Dissociation

The path of dissociation of dimers in 3M-TTR may be followed during the simulation (Figure 7). The first step in conformational transition is the dissociation at the interface I/IV involving residues 17-24 and 110-123, including position 117 mutated to glutamic acid in 3M-TTR at variance with MT-TTR where the corresponding residue is a serine. The transition appears driven by the electrostatic repulsion of the pairs of acidic residues E117 close in each dimer. First the interface is weakened (up to 15 ns) and disrupted (20 ns); then both I/II and III/IV dimers remain rigid. Dimer III/IV rotates for most of the simulation about a hinge centered on salt bridges E51 (I)-R104 (III) and R104 (I)-E51 (III) and hydrogen bonds E51 (III)-T123 (I) and at the end only about the latter two interactions. The final (possibly transient) conformation is stabilized by salt bridges E117 (I)-R21 (III) and K15-E54 at I/III interface. Other interactions at II/III interface are mostly hydrophobic. No sign of dissociation or conformational rearrangement at dimers is observed during the simulation. For both dimers the proximity of E117 acidic groups should result in repulsion which is however reduced by ionic interactions.

4. Conclusions

Molecular dynamics simulations were used in this work to probe the structural stability and the dynamics of engineered mutants of transthyretin (TTR), i.e., the double mutant F87M/L110M (MT-TTR) and the triple mutant F87M/L110M/S117E (3M-TTR), in relation to wild-type. The analysis of trajectories reveals that mutations do not have major impact on protein structure, and the thermodynamic analysis confirms this picture.

Estimation of the free energy from end-point simulations shows that the two mutations F87M/L110M on monomers/dimer and dimers/tetramer equilibria favor the dissociation. The results are however within the uncertainties of the methodology.

Consistent with the latter analysis evaluation of stability with the statistical effective energy functions FoldX and BMF shows that the mutations and the triple mutations F87M/L110M/S117E affect differently the same equilibria and the stability of the monomers. The latter is almost not perturbed by mutations, whereas the equilibria are shifted towards the dissociated species relative to wild-type TTR. This is confirmed by the analysis of 3M-TTR which shows partial dissociation within the first 20 ns of the simulation, implying that contacts are lost at the dimer-dimer interface, whereas dimers appear fairly stable.

Overall the simulations provide a detailed view of the dynamics and thermodynamics of wild-type and mutant transthyretin and a rationale of the observed effects.

Data Availability

The data used to support the findings of this study are available from the corresponding author upon request.

Conflicts of Interest

The authors declare that they have no conflicts of interest.


Cedrix J. Dongmo Foumthuim acknowledges the CINECA Projects HP10CCSKB9 and HP10CIWPJ4 for the availability of high performance computing resources and the Ph.D. fellowship provided by the “Ministero dell’Istruzione, dell’Universita e della Ricerca”.


  1. J. A. Hamilton and M. D. Benson, “Transthyretin: A review from a structural perspective,” Cellular and Molecular Life Sciences, vol. 58, no. 10, pp. 1491–1521, 2001. View at: Publisher Site | Google Scholar
  2. J. R. Rodrigues, C. J. V. Simões, C. G. Silva, and R. M. M. Brito, “Potentially amyloidogenic conformational intermediates populate the unfolding landscape of transthyretin: Insights from molecular dynamics simulations,” Protein Science, vol. 19, no. 2, pp. 202–219, 2010. View at: Publisher Site | Google Scholar
  3. P. Ferreira, O. Sant'Anna, N. Varejão et al., “Correction: Structure-Based Analysis of A19D, a Variant of Transthyretin Involved in Familial Amyloid Cardiomyopathy,” PLoS ONE, vol. 8, no. 12, 2013. View at: Publisher Site | Google Scholar
  4. J. B. Reinés, T. R. Vera, M. U. Martín et al., “Epidemiology of transthyretin-associated familial amyloid polyneuropathy in the Majorcan area: Son Llàtzer Hospital descriptive study,” Orphanet Journal of Rare Diseases, vol. 9, no. 1, article no. 29, 2014. View at: Publisher Site | Google Scholar
  5. X. Jiang, C. S. Smith, H. M. Petrassi et al., “An engineered transthyretin monomer that is nonamyloidogenic, unless it is partially denatured,” Biochemistry, vol. 40, no. 38, pp. 11442–11452, 2001. View at: Publisher Site | Google Scholar
  6. K. H. Lim, H. J. Dyson, J. W. Kelly, and P. E. Wright, “Localized structural fluctuations promote amyloidogenic conformations in transthyretin,” Journal of Molecular Biology, vol. 425, no. 6, pp. 977–988, 2013. View at: Publisher Site | Google Scholar
  7. A. Banerjee, H. R. Bairagya, B. P. Mukhopadhyay, T. K. Nandi, and D. K. Mishra, “Conserved water mediated H-bonding dynamics of Ser117 and Thr119 residues in human transthyretinthyroxin complexation: Inhibitor modeling study through docking and molecular dynamics simulation,” Journal of Molecular Graphics and Modelling, vol. 44, pp. 70–80, 2013. View at: Google Scholar
  8. C. C. F. Blake, M. J. Geisow, S. J. Oatley, B. Rérat, and C. Rérat, “Structure of prealbumin: Secondary, tertiary and quaternary interactions determined by Fourier refinement at 1.8 Å,” Journal of Molecular Biology, vol. 121, no. 3, pp. 339–356, 1978. View at: Publisher Site | Google Scholar
  9. A. Hörnberg, T. Eneqvist, A. Olofsson, E. Lundgren, and A. E. Sauer-Eriksson, “A comparative analysis of 23 structures of the amyloidogenic protein transthyretin,” Journal of Molecular Biology, vol. 302, no. 3, pp. 649–669, 2000. View at: Publisher Site | Google Scholar
  10. M. Haupt, M. P. Blakeley, S. J. Fisher et al., “Binding site asymmetry in human transthyretin: Insights from a joint neutron and X-ray crystallographic analysis using perdeuterated protein,” IUCrJ, vol. 1, pp. 429–438, 2014. View at: Publisher Site | Google Scholar
  11. M. Cianci, C. Folli, F. Zonta, P. Florio, R. Berni, and G. Zanotti, “Structural evidence for asymmetric ligand binding to transthyretin,” Acta Crystallographica Section D: Biological Crystallography, vol. 71, pp. 1582–1592, 2015. View at: Publisher Site | Google Scholar
  12. S. Dubrey, E. Ackermann, and J. Gillmore, “The transthyretin amyloidoses: Advances in therapy,” Postgraduate Medical Journal, vol. 91, no. 1078, pp. 439–448, 2015. View at: Publisher Site | Google Scholar
  13. N. Pasquato, R. Berni, C. Folli, B. Alfieri, L. Cendron, and G. Zanotti, “Acidic pH-induced Conformational Changes in Amyloidogenic Mutant Transthyretin,” Journal of Molecular Biology, vol. 366, no. 3, pp. 711–719, 2007. View at: Publisher Site | Google Scholar
  14. Z. Lai, W. Colón, and J. W. Kelly, “The acid-mediated denaturation pathway of transthyretin yields a conformational intermediate that can self-assemble into amyloid,” Biochemistry, vol. 35, no. 20, pp. 6470–6482, 1996. View at: Publisher Site | Google Scholar
  15. F. Chiti and C. M. Dobson, “Protein misfolding, functional amyloid, and human disease,” Annual Review of Biochemistry, vol. 75, pp. 333–366, 2006. View at: Publisher Site | Google Scholar
  16. G. Invernizzi, E. Papaleo, R. Sabate, and S. Ventura, “Protein aggregation: Mechanisms and functional consequences,” The International Journal of Biochemistry & Cell Biology, vol. 44, no. 9, pp. 1541–1554, 2012. View at: Publisher Site | Google Scholar
  17. F. Chiti and C. M. Dobson, “Protein misfolding, amyloid formation, and human disease: A summary of progress over the last decade,” Annual Review of Biochemistry, vol. 86, pp. 27–68, 2017. View at: Publisher Site | Google Scholar
  18. M. J. Greene, E. S. Klimtchuk, D. C. Seldin, J. L. Berk, and L. H. Connors, “Cooperative stabilization of transthyretin by clusterin and diflunisal,” Biochemistry, vol. 54, no. 2, pp. 268–278, 2015. View at: Publisher Site | Google Scholar
  19. H. A. Lashuel, Z. Lai, and J. W. Kelly, “Characterization of the transthyretin acid denaturation pathways by analytical ultracentrifugation: Implications for wild-type, V30M, and L55P amyloid fibril formation,” Biochemistry, vol. 37, no. 51, pp. 17851–17864, 1998. View at: Publisher Site | Google Scholar
  20. A. Quintas, M. J. M. Saraiva, and R. M. M. Britot, “The tetrameric protein transthyretin dissociates to a non-native monomer in solution. A novel model for amyloidogenesis,” The Journal of Biological Chemistry, vol. 274, no. 46, pp. 32943–32949, 1999. View at: Publisher Site | Google Scholar
  21. A. Quintas, D. C. Vaz, I. Cardoso, M. J. M. Saraiva, and R. M. M. Brito, “Tetramer Dissociation and Monomer Partial Unfolding Precedes Protofibril Formation in Amyloidogenic Transthyretin Variants,” The Journal of Biological Chemistry, vol. 276, no. 29, pp. 27207–27213, 2001. View at: Publisher Site | Google Scholar
  22. X. Jiang, J. N. Buxbaum, and J. W. Kelly, “The V122I cardiomyopathy variant of transthyretin increases the velocity of rate-limiting tetramer dissociation, resulting in accelerated amyloidosis,” Proceedings of the National Acadamy of Sciences of the United States of America, vol. 98, no. 26, pp. 14943–14948, 2001. View at: Publisher Site | Google Scholar
  23. S. M. Johnson, S. Connelly, C. Fearns, E. T. Powers, and J. W. Kelly, “The transthyretin amyloidoses: From delineating the molecular mechanism of aggregation linked to pathology to a regulatory-agency-approved drug,” Journal of Molecular Biology, vol. 421, no. 2-3, pp. 185–203, 2012. View at: Publisher Site | Google Scholar
  24. G. Zanotti, F. Vallese, A. Ferrari et al., “Structural and dynamics evidence for scaffold asymmetric flexibility of the human transthyretin tetramer,” PLoS ONE, vol. 12, no. 12, pp. 1–14, 2017. View at: Google Scholar
  25. M. J. Landrum, J. M. Lee, M. Benson et al., “ClinVar: improving access to variant interpretations and supporting evidence,” Nucleic Acids Research, vol. 46, no. D1, pp. D1062–D1067, 2018. View at: Publisher Site | Google Scholar
  26. J. W. Kelly and P. T. Lansbury, “A chemical approach to elucidate tin mechanism of transthyretin and βprotein amyloid fibril formation,” Amyloid, vol. 1, no. 3, pp. 186–205, 1994. View at: Publisher Site | Google Scholar
  27. S. A. Peterson, T. Klabunde, H. A. Lashuel, H. Purkey, J. C. Sacchettini, and J. W. Kelly, “Inhibiting transthyretin conformational changes that lead to amyloid fibril formation,” Proceedings of the National Acadamy of Sciences of the United States of America, vol. 95, no. 22, pp. 12956–12960, 1998. View at: Publisher Site | Google Scholar
  28. M. R. Almeida, L. Gales, A. M. Damas, I. Cardoso, and M. J. Saraiva, “Small transthyretin (TTR) ligands as possible therapeutic agents in TTR amyloidoses,” CNS & Neurological Disorders - Drug Targets , vol. 4, no. 5, pp. 587–596, 2005. View at: Publisher Site | Google Scholar
  29. Y. Ando and M. Ueda, “Diagnosis and therapeutic approaches to transthyretin amyloidosis,” Current Medicinal Chemistry, vol. 19, no. 15, pp. 2312–2323, 2012. View at: Publisher Site | Google Scholar
  30. M. Nakamura and Y. Ando, “Applications of gene therapy for familial amyloidotic polyneuropathy,” Expert Opinion on Biological Therapy, vol. 4, no. 10, pp. 1621–1627, 2004. View at: Publisher Site | Google Scholar
  31. M. P. Sebastião, G. Merlini, M. J. Saraiva, and A. M. Damas, “The molecular interaction of 4’-iodo-4’-deoxydoxorubicin with Leu-55Pro transthyretin ‘amyloid-like’ oligomer leading to disaggregation,” Biochemical Journal, vol. 351, no. 1, pp. 273–279, 2000. View at: Publisher Site | Google Scholar
  32. F. Fogolari, A. Corazza, V. Yarra, A. Jalaru, P. Viglino, and G. Esposito, “Bluues: A program for the analysis of the electrostatic properties of proteins based on generalized Born radii,” BMC Bioinformatics, vol. 13, no. 4, 2012. View at: Publisher Site | Google Scholar
  33. I. Walsh, G. Minervini, A. Corazza, G. Esposito, S. C. E. Tosatto, and F. Fogolari, “Bluues server: Electrostatic properties of wild-type and mutated protein structures,” Bioinformatics, vol. 28, no. 16, Article ID bts343, pp. 2189-2190, 2012. View at: Publisher Site | Google Scholar
  34. M. K. Gilson, J. A. Given, B. L. Bush, and J. A. McCammon, “The statistical-thermodynamic basis for computation of binding affinities: A critical review,” Biophysical Journal, vol. 72, no. 3, pp. 1047–1069, 1997. View at: Publisher Site | Google Scholar
  35. B. Roux and T. Simonson, “Implicit solvent models,” Biophysical Chemistry, vol. 78, no. 1-2, pp. 1–20, 1999. View at: Publisher Site | Google Scholar
  36. J. Wereszczynski and J. A. McCammon, “Statistical mechanics and molecular dynamics in evaluating thermodynamic properties of biomolecular recognition,” Quarterly Reviews of Biophysics, vol. 45, no. 1, pp. 1–25, 2012. View at: Publisher Site | Google Scholar
  37. F. Fogolari, A. Corazza, and G. Esposito, “Free Energy, Enthalpy and Entropy from Implicit Solvent End-Point Simulations,” Frontiers in Molecular Biosciences, vol. 5, 2018. View at: Publisher Site | Google Scholar
  38. M. F. Sanner, A. J. Olson, and J.-C. Spehner, “Reduced surface: An efficient way to compute molecular surfaces,” Biopolymers, vol. 38, no. 3, pp. 305–320, 1996. View at: Publisher Site | Google Scholar
  39. D. Bashford and D. A. Case, “Generalized born models of macromolecular solvation effects,” Annual Review of Physical Chemistry, vol. 51, no. 1, pp. 129–152, 2000. View at: Publisher Site | Google Scholar
  40. D. Sitkoff, K. A. Sharp, and B. Honig, “Accurate calculation of hydration free energies using macroscopic solvent models,” The Journal of Physical Chemistry C, vol. 98, no. 7, pp. 1978–1988, 1994. View at: Publisher Site | Google Scholar
  41. M. J. Abraham, D. van der Spoel, E. Lindahl, and B. Hess, the GROMACS development team: GROMACS User Manual version 5.0.4, 2014,
  42. H. Singh, N. Misra, V. Hnizdo, A. Fedorowicz, and E. Demchuk, “Nearest neighbor estimates of entropy,” American Journal of Mathematical and Management Sciences, vol. 23, no. 3-4, pp. 301–321, 2003. View at: Publisher Site | Google Scholar | MathSciNet
  43. V. Hnizdo, A. Fedorowicz, H. Singh, and E. Demchuk, “Statistical thermodynamics of internal rotation in a hindering potential of mean force obtained from computer simulations,” Journal of Computational Chemistry, vol. 24, no. 10, pp. 1172–1183, 2003. View at: Publisher Site | Google Scholar
  44. B. J. Killian, J. Yundenfreund Kravitz, and M. K. Gilson, “Extraction of configurational entropy from molecular simulations via an expansion approximation,” The Journal of Chemical Physics, vol. 127, no. 2, Article ID 024107, 2007. View at: Publisher Site | Google Scholar
  45. V. Hnizdo, T. Jun, B. J. Killian, and M. K. Gilson, “Efficient calculation of configurational entropy from molecular simulations by combining the mutual-information expansion and nearest-neighbor methods,” Journal of Computational Chemistry, vol. 29, no. 10, pp. 1605–1614, 2008. View at: Publisher Site | Google Scholar
  46. A. T. Fenley, B. J. Killian, V. Hnizdo, A. Fedorowicz, D. S. Sharp, and M. K. Gilson, “Correlation as a determinant of configurational entropy in supramolecular and protein systems,” The Journal of Physical Chemistry B, vol. 118, no. 24, pp. 6447–6455, 2014. View at: Publisher Site | Google Scholar
  47. D. J. Huggins, “Comparing distance metrics for rotation using the k-nearest neighbors algorithm for entropy estimation,” Journal of Computational Chemistry, vol. 35, no. 5, pp. 377–385, 2014. View at: Publisher Site | Google Scholar
  48. D. J. Huggins, “Quantifying the entropy of binding for water molecules in protein cavities by computing correlations,” Biophysical Journal, vol. 108, no. 4, pp. 928–936, 2015. View at: Publisher Site | Google Scholar
  49. F. Fogolari, A. Corazza, S. Fortuna et al., “Distance-based configurational entropy of proteins from molecular dynamics simulations,” PLoS ONE, vol. 10, no. 7, Article ID e0132356, pp. 1–26, 2015. View at: Publisher Site | Google Scholar
  50. F. Fogolari, C. J. Dongmo Foumthuim, S. Fortuna, M. A. Soler, A. Corazza, and G. Esposito, “Accurate Estimation of the Entropy of Rotation-Translation Probability Distributions,” Journal of Chemical Theory and Computation, vol. 12, no. 1, pp. 1–8, 2016. View at: Publisher Site | Google Scholar
  51. B. M. King and B. Tidor, “MIST: Maximum Information Spanning Trees for dimension reduction of biological data sets,” Bioinformatics, vol. 25, no. 9, pp. 1165–1172, 2009. View at: Publisher Site | Google Scholar
  52. B. M. King, N. W. Silver, and B. Tidor, “Efficient calculation of molecular configurational entropies using an information theoretic approximation,” The Journal of Physical Chemistry B, vol. 116, no. 9, pp. 2891–2904, 2012. View at: Publisher Site | Google Scholar
  53. C. J. Dongmo Foumthuim, A. Corazza, G. Esposito, and F. Fogolari, “Molecular dynamics simulations of β2-microglobulin interaction with hydrophobic surfaces,” Molecular BioSystems, vol. 13, no. 12, pp. 2625–2637, 2017. View at: Publisher Site | Google Scholar
  54. H. M. Berman, J. Westbrook, Z. Feng et al., “The protein data bank,” Nucleic Acids Research, vol. 28, no. 1, pp. 235–242, 2000. View at: Publisher Site | Google Scholar
  55. M. U. Johansson, V. Zoete, O. Michielin, and N. Guex, “Defining and searching for structural motifs using DeepView/Swiss-PdbViewer,” BMC Bioinformatics, vol. 13, no. 1, article 173, 2012. View at: Publisher Site | Google Scholar
  56. E. Lindahl, P. Bjelkmar, P. Larsson, M. A. Cuendet, and B. Hess, “Implementation of the charmm force field in GROMACS: Analysis of protein stability effects from correction maps, virtual interaction sites, and water models,” Journal of Chemical Theory and Computation, vol. 6, no. 2, pp. 459–466, 2010. View at: Publisher Site | Google Scholar
  57. K. Lindorff-Larsen, S. Piana, K. Palmo et al., “Improved side-chain torsion potentials for the Amber ff99SB protein force field,” Proteins: Structure, Function, and Bioinformatics, vol. 78, no. 8, pp. 1950–1958, 2010. View at: Publisher Site | Google Scholar
  58. W. L. Jorgensen, J. Chandrasekhar, J. D. Madura, R. W. Impey, and M. L. Klein, “Comparison of simple potential functions for simulating liquid water,” The Journal of Chemical Physics, vol. 79, no. 2, pp. 926–935, 1983. View at: Publisher Site | Google Scholar
  59. W. L. Jorgensen and J. D. Madura, “Temperature and size dependence for monte carlo simulations of TIP4P water,” Molecular Physics, vol. 56, no. 6, pp. 1381–1392, 1985. View at: Publisher Site | Google Scholar
  60. M. W. Mahoney and W. L. Jorgensen, “A five-site model for liquid water and the reproduction of the density anomaly by rigid, nonpolarizable potential functions,” The Journal of Chemical Physics, vol. 112, no. 20, pp. 8910–8922, 2000. View at: Publisher Site | Google Scholar
  61. T. Darden, D. York, and L. Pedersen, “Particle mesh Ewald: an N·log(N) method for Ewald sums in large systems,” The Journal of Chemical Physics, vol. 98, no. 12, pp. 10089–10092, 1993. View at: Publisher Site | Google Scholar
  62. U. Essmann, L. Perera, M. L. Berkowitz, T. Darden, H. Lee, and L. G. Pedersen, “A smooth particle mesh Ewald method,” The Journal of Chemical Physics, vol. 103, no. 19, pp. 8577–8593, 1995. View at: Publisher Site | Google Scholar
  63. G. Bussi, D. Donadio, and M. Parrinello, “Canonical sampling through velocity rescaling,” The Journal of Chemical Physics, vol. 126, no. 1, Article ID 014101, 2007. View at: Publisher Site | Google Scholar
  64. B. Hess, H. Bekker, H. J. C. Berendsen, and J. G. E. M. Fraaije, “LINCS: a linear Constraint Solver for molecular simulations,” Journal of Computational Chemistry, vol. 18, no. 12, pp. 1463–1472, 1997. View at: Publisher Site | Google Scholar
  65. M. Parrinello and A. Rahman, “Polymorphic transitions in single crystals: a new molecular dynamics method,” Journal of Applied Physics, vol. 52, no. 12, pp. 7182–7190, 1981. View at: Publisher Site | Google Scholar
  66. S. Nosé and M. L. Klein, “Constant pressure molecular dynamics for molecular systems,” Molecular Physics, vol. 50, no. 5, pp. 1055–1076, 1983. View at: Publisher Site | Google Scholar
  67. B. Hess, C. Kutzner, D. van der Spoel, and E. Lindahl, “GRGMACS 4: algorithms for highly efficient, load-balanced, and scalable molecular simulation,” Journal of Chemical Theory and Computation, vol. 4, no. 3, pp. 435–447, 2008. View at: Publisher Site | Google Scholar
  68. J. Schymkowitz, J. Borg, F. Stricher, R. Nys, F. Rousseau, and L. Serrano, “The FoldX web server: an online force field,” Nucleic Acids Research, vol. 33, no. 2, pp. W382–W388, 2005. View at: Publisher Site | Google Scholar
  69. M. Berrera, H. Molinari, and F. Fogolari, “Amino acid empirical contact energy definitions for fold recognition in the space of contact maps,” BMC Bioinformatics, vol. 4, article no. 8, 2003. View at: Publisher Site | Google Scholar
  70. W. L. Delano, “The PyMOL Molecular Graphics System Version Schrödinger, LLC,” 2002, View at: Google Scholar
  71. W. Humphrey, A. Dalke, and K. Schulten, “VMD: visual molecular dynamics,” Journal of Molecular Graphics, vol. 14, no. 1, pp. 33–38, 1996. View at: Publisher Site | Google Scholar
  72. W. Kabsch and C. Sander, “Dictionary of protein secondary structure: pattern recognition of hydrogen-bonded and geometrical features,” Biopolymers, vol. 22, no. 12, pp. 2577–2637, 1983. View at: Publisher Site | Google Scholar
  73. D. Petrov and B. Zagrovic, “Are Current Atomistic Force Fields Accurate Enough to Study Proteins in Crowded Environments?” PLoS Computational Biology, vol. 10, no. 5, Article ID e1003638, 2014. View at: Publisher Site | Google Scholar
  74. L. A. Abriata and M. Dal Peraro, “Assessing the potential of atomistic molecular dynamics simulations to probe reversible protein-protein recognition and binding,” Scientific Reports, vol. 5, Article ID 10549, 2015. View at: Publisher Site | Google Scholar
  75. S. Hovmöller, T. Zhou, and T. Ohlson, “Conformations of amino acids in proteins,” Acta Crystallographica Section D: Biological Crystallography, vol. 58, no. 5, pp. 768–776, 2002. View at: Publisher Site | Google Scholar

Copyright © 2018 Cedrix J. Dongmo Foumthuim et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

More related articles

 PDF Download Citation Citation
 Download other formatsMore
 Order printed copiesOrder

Related articles

We are committed to sharing findings related to COVID-19 as quickly as possible. We will be providing unlimited waivers of publication charges for accepted research articles as well as case reports and case series related to COVID-19. Review articles are excluded from this waiver policy. Sign up here as a reviewer to help fast-track new submissions.