Dynamics and Thermodynamics of Transthyretin Association from Molecular Dynamics Simulations

Hint approach on Transthyretin folding/unfolding mechanism comprehension

Non-covalent intramolecular interactions play a key role in the protein folding process. Aminoacidic mutations or changes in physiological conditions such as pH and/or temperature variations can compromise intramolecular stability generating misfolding or unfolding proteins with consequent impairment of functionality and the triggering of pathological states. The intramolecular HINT scoring function recently implemented and validated, is proposed as a rapid and sensitive method for the evaluation of different conformational states characterizing destabilization processes. In this work, the stability of Transthyretin, whose denaturation is related to amyloid fibril formation, is evaluated by generating multiple structural mutated models under different pH conditions in comparison with experimental data. These results suggest that the HINT scoring function can be used for an accurate and rapid evaluation and computational prediction of the effects of structural changes on any protein system.

Solvent quality and solvent polarity in polypeptides

Using molecular dynamics and thermodynamic integration, we report on the solvation process of seven polypeptides (GLY, ALA, ILE, ASN, LYS, ARG, GLU) in water and in cyclohexane. The polypeptides are selected to cover the full hydrophobic scale while varying their chain length from tri- to undeca-homopeptides, providing indications on possible non-additivity effects as well as the role of the peptide backbone in the overall stability of the polypeptides. The use of different solvents and different polypeptides allows us to investigate the relation between solvent quality - the capacity of a given solvent to fold/unfold a given biopolymer often described on a scale ranging from "good" to "poor"; and solvent polarity - related to the specific interactions of any solvent with respect to a reference solvent. Undeca-glycine is found to be the only polypeptide to have a stable collapse in water (polar solvent), with the other hydrophobic polypeptides displaying repeated folding and unfolding events in water, with polar polypeptides presenting even more complex behavior. By contrast, all polypeptides are found to keep an extended conformation in cyclohexane, irrespective of their polarity. All considered polypeptides are also found to have favorable solvation free energy independent of the solvent polarity and their intrinsic hydrophobicity, clearly highlighting the prominent stabilizing role of the peptide backbone - with the solvation process largely enthalpically dominated in polar polypeptides and partially entropically driven for hydrophobic polypeptides. Our study thus reveals the complexity of the solvation process of polypeptides defying the common view "like dissolves like", with the solute polarity playing the most prominent role. The absence of mirror symmetry upon the inversion of polarities of both the solvent and the polypeptides is confirmed.

Atomic Mechanisms of Transthyretin Tetramer Dissociation Studied by Molecular Dynamics Simulations

The dissociation of the transthyretin (TTR) tetramer into a monomer is closely related to various TTR amyloidoses in humans. While the tetramer dissociation has been reported to be the rate-limiting step for TTR aggregation, few details are known about the mechanism. Here, molecular dynamics (MD) simulations were performed by combining conventional MD and biased metadynamics to investigate the mechanism for the wild-type (WT) and mutant (T119M) structures. Both were found to have a great deal in common. Conventional MD simulations reveal that interfacial hydrophobic interactions contribute significantly to stabilize the tetramer. Interfacial residues including L17, V20, L110, and V121 with close contacts form a hydrophobic channel. Metadynamics simulations indicate that the mouth opening of the hydrophobic channel is the first and the most difficult step for dissociation. Interactions of V20 between opposing dimers lock four monomers into the tetramer, and disruption of the interactions is found to be involved in the final step. During the dissociation, an increasing extent of solvation was observed by calculating the radial distribution functions of water around interfacial hydrophobic residues, suggesting that water plays a role in driving the tetramer dissociation. Moreover, compared to T119, residue M119 has a longer side chain that extends into the hydrophobic channel, making solvation more difficult, consistent with a higher energy barrier for dissociation of the T119M tetramer. This result provides a good explanation for the protective role of the T119M mutation. Overall, this study can provide atomic-level insights to better understand the pathogenesis of TTR amyloidosis and guide rational drug design in the future.

Transthyretin: From Structural Stability to Osteoarticular and Cardiovascular Diseases

Transthyretin (TTR) is a tetrameric protein transporting hormones in the plasma and brain, which has many other activities that have not been fully acknowledged. TTR is a positive indicator of nutrition status and is negatively correlated with inflammation. TTR is a neuroprotective and oxidative-stress-suppressing factor. The TTR structure is destabilized by mutations, oxidative modifications, aging, proteolysis, and metal cations, including Ca2+. Destabilized TTR molecules form amyloid deposits, resulting in senile and familial amyloidopathies. This review links structural stability of TTR with the environmental factors, particularly oxidative stress and Ca2+, and the processes involved in the pathogenesis of TTR-related diseases. The roles of TTR in biomineralization, calcification, and osteoarticular and cardiovascular diseases are broadly discussed. The association of TTR-related diseases and vascular and ligament tissue calcification with TTR levels and TTR structure is presented. It is indicated that unaggregated TTR and TTR amyloid are bound by vicious cycles, and that TTR may have an as yet undetermined role(s) at the crossroads of calcification, blood coagulation, and immune response.

Transthyretin Misfolding, A Fatal Structural Pathogenesis Mechanism

Transthyretin (TTR) is an essential transporter of a thyroid hormone and a holo-retinol binding protein, found abundantly in human plasma and cerebrospinal fluid. In addition, this protein is infamous for its amyloidogenic propensity, causing various amyloidoses in humans, such as senile systemic amyloidosis, familial amyloid polyneuropathy, and familial amyloid cardiomyopathy. It has been known for over two decades that decreased stability of the native tetrameric conformation of TTR is the main cause of these diseases. Yet, mechanistic details on the amyloidogenic transformation of TTR were not clear until recent multidisciplinary investigations on various structural states of TTR. In this review, we discuss recent advancements in the structural understanding of TTR misfolding and amyloidosis processes. Special emphasis has been laid on the observations of novel structural features in various amyloidogenic species of TTR. In addition, proteolysis-induced fragmentation of TTR, a recently proposed mechanism facilitating TTR amyloidosis, has been discussed in light of its structural consequences and relevance to acknowledge the amyloidogenicity of TTR.

Divergence Entropy-Based Evaluation of Hydrophobic Core in Aggressive and Resistant Forms of Transthyretin

The two forms of transthyretin differing slightly in the tertiary structure, despite the presence of five mutations, show radically different properties in terms of susceptibility to the amyloid transformation process. These two forms of transthyretin are the object of analysis. The search for the sources of these differences was carried out by means of a comparative analysis of the structure of these molecules in their native and early intermediate stage forms in the folding process. The criterion for assessing the degree of similarity and differences is the status of the hydrophobic core. The comparison of the level of arrangement of the hydrophobic core and its initial stages is possible thanks to the application of divergence entropy for the early intermediate stage and for the final forms. It was shown that the minimal differences observed in the structure of the hydrophobic core of the forms available in PDB, turned out to be significantly different in the early stage (ES) structure in folding process. The determined values of divergence entropy for both ES forms indicate the presence of the seed of hydrophobic core only in the form resistant to amyloid transformation. In the form of aggressively undergoing amyloid transformation, the structure lacking such a seed is revealed, being a stretched one with a high content of β-type structure. In the discussed case, the active presence of water in the structural transformation of proteins expressed in the fuzzy oil drop model (FOD) is of decisive importance for the generation of the final protein structure. It has been shown that the resistant form tends to generate a centric hydrophobic core with the possibility of creating a globular structure, i.e., a spherical micelle-like form. The aggressively transforming form reveals in the structure of its early intermediate, a tendency to form the ribbon-like micelle as observed in amyloid.

Missense NR2F1 variant in monozygotic twins affected with the Bosch-Boonstra-Schaaf optic atrophy syndrome

BACKGROUND

The Bosch-Boonstra-Schaaf optic atrophy syndrome (BBSOAS) is an autosomal-dominant disorder (OMIM615722) mostly characterized by optic atrophy and/or hypoplasia, mild intellectual disability, hypotonia, seizures/infantile epilepsy. This disorder is caused by loss-of-function alterations of NR2F1 (i.e., either whole gene deletions or single nucleotide variants) and, to date, 40 patients have been identified with deletions or mutations in this gene. Here we describe two monozygotic twins harboring a de novo missense variant in the DNA-binding domain of NR2F1 (c.313G>A, p.Gly105Ser), with well-characterized features associated to BBSOAS.

METHODS

Patients' DNA was analyzed by exome sequencing identifying the missense variant c.313G>A in NR2F1 (NM_005654.4). Furthermore, molecular modeling was performed to evaluate putative differences in DNA binding between wild-type and mutated NR2F1.

RESULTS

The missense variant is predicted to be likely pathogenetic following the ACMG (American College of Medical Genetics and Genomics)/AMP (Association for Molecular Pathology) guidelines. Indeed, dynamic simulation experiments highlighted that the Gly105Ser substitution let the formation of a hydrogen bond between the S105 side chain and R142 and a base (G5) of the DNA sequence, allowing us to hypothesize that the G105 residue might be evolutionary conserved due to the absence of a side chain, besides glycine conformational features. Therefore, the G105S variation seems to cause a stiffening and a possible deformation in the protein-DNA complex due to the interaction of residues R142-S105 and G5 on the DNA, compared to the wild-type.

CONCLUSION

In summary, we described two monozygotic twins harboring a novel Gly105Ser mutation in NR2F1 DNA binding domain, displaying the classical phenotype of BBSOAS-affected patients. Our computational data suggest a dominant negative effect of this newly characterized missense variant. To date, this is the first genetic report analyzing in silico structural consequences of NR2F1 Gly105Ser substitution.