The nonlinear structure of the desmoplakin plakin domain and the effects of cardiomyopathy-linked mutations

Desmoplakin is a cytoplasmic desmosomal protein that plays a vital role in normal intercellular adhesion. Mutations in desmoplakin can result in devastating skin blistering diseases and arrhythmogenic right ventricular cardiomyopathy, a heart muscle disorder associated with ventricular arrhythmias, heart failure, and sudden death. The desmoplakin N-terminal region is a 1056-amino-acid sequence of unknown structure. It mediates interactions with other desmosomal proteins, is found in a variety of plakin proteins, and spans what has been termed the “plakin domain,” which includes residues 180-1022 and consists of six spectrin repeats (SRs) and an Src homology 3 domain. Herein we elucidate the architecture of desmoplakin's plakin domain, as well as its constituent tandem SRs. Small-angle X-ray scattering analysis shows that the entire plakin domain has an “L” shape, with a long arm and a short arm held at a perpendicular angle. The long arm is 24.0 nm long and accommodates four stably folded SRs arranged in tandem. In contrast, the short arm is 17.9 nm in length and accommodates two independently folded repeats and an extended C-terminus. We show that mutations linked to arrhythmogenic right ventricular cardiomyopathy (K470E and R808C) cause local conformational alterations, while the overall folded structure is maintained. This provides the first structural and mechanistic insights into an entire plakin domain and provides a basis for understanding the critical role of desmoplakin in desmosome function.     

The structure of a tandem pair of spectrin repeats of plectin reveals a modular organization of the plakin domain

Plectin is a large and versatile cytoskeletal linker and member of the plakin protein family. Plakins share a conserved region called the plakin domain located near their N terminus. We have determined the crystal structure of an N-terminal fragment of the plakin domain of plectin to 2.05 A resolution. This region is adjacent to the actin-binding domain and is required for efficient binding to the integrin alpha6beta4 in hemidesmosomes. The structure is formed by two spectrin repeats connected by an alpha-helix that spans these two repeats. While the first repeat is very similar to other known structures, the second repeat is structurally different with a hydrophobic core, narrower than that in canonical spectrin repeats. Sequence analysis of the plakin domain revealed the presence of up to nine consecutive spectrin repeats organized in an array of tandem modules, and a Src-homology 3 domain inserted in the central spectrin repeat. The structure of the plakin domain is reminiscent of the modular organization of members of the spectrin family. The architecture of the plakin domain suggests that it forms an elongated and flexible structure, and provides a novel molecular explanation for the contribution of plectin and other plakins to the elasticity and stability of tissues subjected to mechanical stress, such as the skin and striated muscle.     

Structures of two intermediate filament-binding fragments of desmoplakin reveal a unique repeat motif structure

Desmosomes are intercellular junctions in which cadherin cell adhesion molecules are linked to the intermediate filament (IF) system. Desmoplakin is a member of the plakin family of IF-binding proteins. The C-terminal domain of desmoplakin (DPCT) mediates binding to IFs in desmosomes. The DPCT sequence contains three regions, termed A, B and C, consisting of 4.5 copies of a 38-amino acid repeat motif. We demonstrate that these regions form discrete subdomains that bind to IFs and report the crystal structures of domains B and C. In contrast to the elongated structures formed by other kinds of repeat motifs, the plakin repeats form a globular structure with a unique fold. A conserved basic groove found on the domain may represent an IF-binding site.     

Structural basis for myopathic defects engendered by alterations in the myosin rod

While mutations in the myosin subfragment 1 motor domain can directly disrupt the generation and transmission of force along myofibrils and lead to myopathy, the mechanism whereby mutations in the myosin rod influences mechanical function is less clear. Here, we used a combination of various imaging techniques and molecular dynamics simulations to test the hypothesis that perturbations in the myosin rod can disturb normal sarcomeric uniformity and, like motor domain lesions, would influence force production and propagation. We show that disrupting the rod can alter its nanomechanical properties and, in vivo, can drive asymmetric myofilament and sarcomere formation. Our imaging results indicate that myosin rod mutations likely disturb production and/or propagation of contractile force. This provides a unifying theory where common pathological cascades accompany both myosin motor and specific rod domain mutations. Finally, we suggest that sarcomeric inhomogeneity, caused by asymmetric thick filaments, could be a useful index of myopathic dysfunction.     

Biophysical simulations and structure-based modeling of residue interaction networks in the tumor suppressor proteins reveal functional role of cancer mutation hotspots in molecular communication

In the current study, we have combined molecular simulations and energetic analysis with dynamics-based network modeling and perturbation response scanning to determine molecular signatures of mutational hotspot residues in the p53, PTEN, and SMAD4 tumor suppressor proteins. By examining structure, energetics and dynamics of these proteins, we have shown that inactivating mutations preferentially target a group of structurally stable residues that play a fundamental role in global propagation of dynamic fluctuations and mediating allosteric interaction networks. Through integration of long-range perturbation dynamics and network-based approaches, we have quantified allosteric potential of residues in the studied proteins. The results have revealed that mutational hotspot sites often correspond to high centrality mediating centers of the residue interaction networks that are responsible for coordination of global dynamic changes and allosteric signaling. Our findings have also suggested that structurally stable mutational hotpots can act as major effectors of allosteric interactions and mutations in these positions are typically associated with severe phenotype. Modeling of shortest inter-residue pathways has shown that mutational hotspot sites can also serve as key mediating bridges of allosteric communication in the p53 and PTEN protein structures. Multiple regression models have indicated that functional significance of mutational hotspots can be strongly associated with the network signatures serving as robust predictors of critical regulatory positions responsible for loss-of-function phenotype. The results of this computational investigation are compared with the experimental studies and reveal molecular signatures of mutational hotspots, providing a plausible rationale for explaining and localizing disease-causing mutations in tumor suppressor genes.     

Role of E542 and E545 missense mutations of PIK3CA in breast cancer: a comparative computational approach

Recent statistics describe breast cancer as the leading cause of death among women across the world with varied causes and reasons. Lifestyle, diet, genetic and environmental factors introduce their generous contributions towards breast cancer, among which genetic factors have lately become one of the most important aspects in understanding the mechanism. Although various genes have already been reported in causing breast cancer, PIK3CA stands second on the list. Mutations observed in this gene have the ability to trigger the different activities of the cell, thereby bypassing the regular cellular cycle. Among the mutations in PIK3CA, three hotspot mutations were commonly reported, one in the catalytic domain (position HIS1047) and other two in the helical domain (position GLU542 and GLU545). In the helical domain of PIK3CA, the lysine substitution at 542–545 positions was significantly studied in causing breast cancer. To compare the deleterious effect of these mutations, in silico prediction tools along with molecular dynamics simulations and molecular docking approach was initiated to analyse the change in binding landscape upon mutation. In this comparative analysis, we report that the mere existence of mutant E545K can trigger the function of the protein but may not be as harmful as H1047R. Among the two mutations E542K and E545K, the latter shows the most deleterious effect that correlates with the previous reported experimental studies. We assume the results observed in this combinatorial computational study might further pave a better way for providing better treatment procedures.     

Structural Basis of Type 2A von Willebrand Disease Investigated by Molecular Dynamics Simulations and Experiments

The hemostatic function of von Willebrand factor is downregulated by the metalloprotease ADAMTS13, which cleaves at a unique site normally buried in the A2 domain. Exposure of the proteolytic site is induced in the wild-type by shear stress as von Willebrand factor circulates in blood. Mutations in the A2 domain, which increase its susceptibility to cleavage, cause type 2A von Willebrand disease. In this study, molecular dynamics simulations suggest that the A2 domain unfolds under tensile force progressively through a series of steps. The simulation results also indicated that three type 2A mutations in the C-terminal half of the A2 domain, L1657I, I1628T and E1638K, destabilize the native state fold of the protein. Furthermore, all three type 2A mutations lowered in silico the tensile force necessary to undock the C-terminal helix 6 from the rest of the A2 domain, the first event in the unfolding pathway. The mutations F1520A, I1651A and A1661G were also predicted by simulations to destabilize the A2 domain and facilitate exposure of the cleavage site. Recombinant A2 domain proteins were expressed and cleavage assays were performed with the wild-type and single-point mutants. All three type 2A and two of the three predicted mutations exhibited increased rate of cleavage by ADAMTS13. These results confirm that destabilization of the helix 6 in the A2 domain facilitates exposure of the cleavage site and increases the rate of cleavage by ADAMTS13.     

Hinged Plakin Domains Provide Specialized Degrees of Articulation in Envoplakin,Periplakin and Desmoplakin

Envoplakin, periplakin and desmoplakin are cytoskeletal proteins that provide structural integrity within the skin and heart by resisting shear forces. Here we reveal the nature of unique hinges within their plakin domains that provides divergent degrees of flexibility between rigid long and short arms composed of spectrin repeats. The range of mobility of the two arms about the hinge is revealed by applying the ensemble optimization method to small-angle X-ray scattering data. Envoplakin and periplakin adopt ‘L’ shaped conformations exhibiting a ‘helicopter propeller’-like mobility about the hinge. By contrast desmoplakin exhibits essentially unrestricted mobility by ‘jack-knifing’ about the hinge. Thus the diversity of molecular jointing that can occur about plakin hinges includes ‘L’ shaped bends, ‘U’ turns and fully extended ‘I’ orientations between rigid blocks of spectrin repeats. This establishes specialised hinges in plakin domains as a key source of flexibility that may allow sweeping of cellular spaces during assembly of cellular structures and could impart adaptability, so preventing irreversible damage to desmosomes and the cell cytoskeleton upon exposure to mechanical stress.     

Integrating NMR,SAXS, and Atomistic Simulations: Structure and Dynamics of a Two-Domain Protein

Multidomain proteins with two or more independently folded functional domains are prevalent in nature. Whereas most multidomain proteins are linked linearly in sequence, roughly one-tenth possess domain insertions where a guest domain is implanted into a loop of a host domain, such that the two domains are connected by a pair of interdomain linkers. Here, we characterized the influence of the interdomain linkers on the structure and dynamics of a domain-insertion protein in which the guest LysM domain is inserted into a central loop of the host CVNH domain. Expanding upon our previous crystallographic and NMR studies, we applied SAXS in combination with NMR paramagnetic relaxation enhancement to construct a structural model of the overall two-domain system. Although the two domains have no fixed relative orientation, certain orientations were found to be preferred over others. We also assessed the accuracies of molecular mechanics force fields in modeling the structure and dynamics of tethered multidomain proteins by integrating our experimental results with microsecond-scale atomistic molecular dynamics simulations. In particular, our evaluation of two different combinations of the latest force fields and water models revealed that both combinations accurately reproduce certain structural and dynamical properties, but are inaccurate for others. Overall, our study illustrates the value of integrating experimental NMR and SAXS studies with long timescale atomistic simulations for characterizing structural ensembles of flexibly linked multidomain systems.     

Molecular mechanism of intermediate filament recognition by plakin proteins

The plakin family of cytolinkers interacts with intermediate filaments (IFs) through plakin repeat domain (PRD) and linker modules. Recent structure/function studies have established the molecular basis of envoplakin-PRD and periplakin-linker interactions with vimentin. Both plakin modules share a broad basic groove which recognizes acidic rod elements on IFs, a mechanism that is applicable to other plakin family members. This review postulates a universal IF engagement mechanism that illuminates the specific effects of pathogenic mutations associated with diseases including arrhythmogenic right ventricular cardiomyopathy, and reveals how diverse plakin proteins offer tailored IF tethering to ensure stable, dynamic and regulated cellular structures.     

Molecular dynamics as a tool to detect protein foldability. A mutant of domain B1 of protein G with non-native secondary structure propensities

The usefulness of molecular dynamics to assess the structural integrity of mutants containing several mutations has been investigated. Our goal was to determine whether molecular dynamics would be able to discriminate mutants of a protein having a close-to-wild-type fold, from those that are not folded under the same conditions. We used as a model the B1 domain of protein G in which we replaced the unique central alpha-helix by the sequence of the second beta-hairpin, which has a strong intrinsic propensity to form this secondary structure in solution. In the resulting protein, one-third of the secondary structure has been replaced by a non-native one. Models of the mutants were built based on the three-dimensional structure of the wild-type GB1 domain. During 2 ns of molecular dynamics simulations on these models, mutants containing up to 10 mutations in the helix retained the native fold, while another mutant with an additional mutation unfolded. This result is in agreement with our circular dichroism and NMR experiments, which indicated that the former mutants fold into a structure similar to the wild-type, as opposed to the latter mutant which is partly unfolded. Additionally, a mutant containing six mutations scattered through the surface of the domain, and which is unfolded, was also detected by the simulation. This study suggests that molecular dynamics calculations could be performed on molecular models of mutants of a protein to evaluate their foldability, prior to a mutagenesis experiment.     

The impact of dyskeratosis congenita mutations on the structure and dynamics of the human telomerase RNA pseudoknot domain

The pseudoknot domain is a functionally crucial part of telomerase RNA and influences the activity and stability of the ribonucleoprotein complex. Autosomal dominant dyskeratosis congenita (DKC) is an inherited disease that is linked to mutations in telomerase RNA and impairs telomerase function. In this paper, we present a computational prediction of the influence of two base DKC mutations on the structure, dynamics, and stability of the pseudoknot domain. We use molecular dynamics simulations, MM-GBSA free energy calculations, static analysis, and melting simulations analysis. Our results show that the DKC mutations stabilize the hairpin form and destabilize the pseudoknot form of telomerase RNA. Moreover, the P3 region of the predicted DKC-mutated pseudoknot structure is unstable and fails to form as a defined helical stem. We directly compare our predictions with experimental observations by calculating the enthalpy of folding and melting profiles for each structure. The enthalpy values are in very good agreement with values determined by thermal denaturation experiments. The melting simulations and simulations at elevated temperatures show the existence of an intermediate structure, which involves the formation of two UU base pairs observed in the hairpin form of the pseudoknot domain.     

Molecular dynamics simulation of the SH3 domain aggregation suggests a generic amyloidogenesis mechanism

We use molecular dynamics simulation to study the aggregation of Src SH3 domain proteins. For the case of two proteins, we observe two possible aggregation conformations: the closed form dimer and the open aggregation state. The closed dimer is formed by "domain swapping"-the two proteins exchange their RT-loops. All the hydrophobic residues are buried inside the dimer so proteins cannot further aggregate into elongated amyloid fibrils. We find that the open structure-stabilized by backbone hydrogen bond interactions-packs the RT-loops together by swapping the two strands of the RT-loop. The packed RT-loops form a beta-sheet structure and expose the backbone to promote further aggregation. We also simulate more than two proteins, and find that the aggregate adopts a fibrillar double beta-sheet structure, which is formed by packing the RT-loops from different proteins. Our simulations are consistent with a possible generic amyloidogenesis scenario.     

The Impact of Dyskeratosis Congenita Mutations on the Structure and Dynamics of the Human Telomerase RNA Pseudoknot Domain

Abstract

The pseudoknot domain is a functionally crucial part of telomerase RNA and influences the activity and stability of the ribonucleoprotein complex. Autosomal dominant dyskeratosis congenita (DKC) is an inherited disease that is linked to mutations in telomerase RNA and impairs telomerase function. In this paper, we present a computational prediction of the influence of two base DKC mutations on the structure, dynamics, and stability of the pseudoknot domain. We use molecular dynamics simulations, MM-GBSA free energy calculations, static analysis, and melting simulations analysis. Our results show that the DKC mutations stabilize the hairpin form and destabilize the pseudoknot form of telomerase RNA. Moreover, the P3 region of the predicted DKC-mutated pseudoknot structure is unstable and fails to form as a defined helical stem. We directly compare our predictions with experimental observations by calculating the enthalpy of folding and melting profiles for each structure. The enthalpy values are in very good agreement with values determined by thermal denaturation experiments. The melting simulations and simulations at elevated temperatures show the existence of an intermediate structure, which involves the formation of two UU base pairs observed in the hairpin form of the pseudoknot domain.     

Molecular dynamics simulations of Alzheimer's beta-amyloid protofilaments

Filamentous amyloid aggregates are central to the pathology of Alzheimer's disease. We use all-atom molecular dynamics (MD) simulations with explicit solvent and multiple force fields to probe the structural stability and the conformational dynamics of several models of Alzheimer's beta-amyloid fibril structures, for both wild-type and mutated amino acid sequences. The structural models are based on recent solid state NMR data. In these models, the peptides form in-register parallel beta-sheets along the fibril axis, with dimers of two U-shaped peptides located in layers normal to the fibril axis. Four different topologies are explored for stacking the beta-strand regions against each other to form a hydrophobic core. Our MD results suggest that all four NMR-based models are structurally stable, and we find good agreement with dihedral angles estimated from solid-state NMR experiments. Asp23 and Lys28 form buried salt-bridges, resulting in an alternating arrangement of the negatively and positively charged residues along the fibril axis that is reminiscent of a one-dimensional ionic crystal. Interior water molecules are solvating the buried salt-bridges. Based on data from NMR measurements and MD simulations of short amyloid fibrils, we constructed structural models of long fibrils. Calculated X-ray fiber diffraction patterns show the characteristics of packed beta-sheets seen in experiments, and suggest new experiments that could discriminate between various fibril topologies.     

The allosteric role of the Ca2+ switch in adhesion and elasticity of C-cadherin

Modular proteins such as titin, fibronectin, and cadherin are ubiquitous components of living cells. Often involved in signaling and mechanical processes, their architecture is characterized by domains containing a variable number of heterogeneous “repeats” arranged in series, with either flexible or rigid linker regions that determine their elasticity. Cadherin repeats arranged in series are unique in that linker regions also feature calcium-binding motifs. While it is well known that the extracellular repeats of cadherin proteins mediate cell-cell adhesion in a calcium-dependent manner, the molecular mechanisms behind the influence of calcium in adhesion dynamics and cadherin's mechanical response are not well understood. Here we show, using molecular dynamics simulations, how calcium ions control the structural integrity of cadherin's linker regions, thereby affecting cadherin's equilibrium dynamics, the availability of key residues involved in cell-cell adhesion, and cadherin's mechanical response. The all-atom, multi-nanosecond molecular dynamics simulations involved the entire C-cadherin extracellular domain solvated in water (a 345,000 atom system). Equilibrium simulations show that the extracellular domain maintains its crystal conformation (elongated and slightly curved) when calcium ions are present. In the absence of calcium ions, however, it assumes a disordered conformation. The conserved residue Trp², which is thought to insert itself into a hydrophobic pocket of another cadherin molecule (thereby providing the basis for cell-cell adhesion), switches conformation from exposed to intermittently buried upon removal of calcium ions. Furthermore, the overall mechanical response of C-cadherin's extracellular domain is characterized at low force by changes in shape (tertiary structure elasticity), and at high force by unraveling of secondary structure elements (secondary structure elasticity). This mechanical response is modulated by calcium ions at both low and high force, switching from a stiff, rod-like to a soft, spring-like behavior upon removal of ions. The simulations provide an unprecedented molecular view of calcium-mediated allostery in cadherins, also illustrating the general principles of linker-mediated elasticity of modular proteins relevant not only for cell-cell adhesion and sound transduction, but also muscle elasticity.     

Structure and dynamics of the von Willebrand Factor C6 domain

Von Willebrand disease (VWD) is a bleeding disorder with different levels of severity. VWD-associated mutations are located in the von Willebrand factor (VWF) gene, coding for the large multidomain plasma protein VWF with essential roles in hemostasis and thrombosis. On the one hand, a variety of mutations in the C-domains of VWF are associated with increased bleeding upon vascular injury. On the other hand, VWF gain-of-function (GOF) mutations in the C4 domain have recently been identified, which induce an increased risk of myocardial infarction. Mechanistic insights into how these mutations affect the molecular behavior of VWF are scarce and holistic approaches are challenging due to the multidomain and multimeric character of this large protein. Here, we determine the structure and dynamics of the C6 domain and the single nucleotide polymorphism (SNP) variant G2705R in C6 by combining nuclear magnetic resonance spectroscopy, molecular dynamics simulations and aggregometry. Our findings indicate that this mutation mostly destabilizes VWF by leading to a more pronounced hinging between both subdomains of C6. Hemostatic parameters of variant G2705R are close to normal under static conditions, but the missense mutation results in a gain-of-function under flow conditions, due to decreased VWF stem stability. Together with the fact that two C4 variants also exhibit GOF characteristics, our data underline the importance of the VWF stem region in VWF’s hemostatic activity and the risk of mutation-associated prothrombotic properties in VWF C-domain variants due to altered stem dynamics.     

Structural Characterization of N-WASP Domain V Using MD Simulations with NMR and SAXS Data

Because of their large conformational heterogeneity, structural characterization of intrinsically disordered proteins (IDPs) is very challenging using classical experimental methods alone. In this study, we use NMR and small-angle x-ray scattering (SAXS) data with multiple molecular dynamics (MD) simulations to describe the conformational ensemble of the fully disordered verprolin homology domain of the neural Aldrich syndrome protein involved in the regulation of actin polymerization. First, we studied several back-calculation software of SAXS scattering intensity and optimized the adjustable parameters to accurately calculate the SAXS intensity from an atomic structure. We also identified the most appropriate force fields for MD simulations of this IDP. Then, we analyzed four conformational ensembles of neural Aldrich syndrome protein verprolin homology domain, two generated with the program flexible-meccano with or without NMR-derived information as input and two others generated by MD simulations with two different force fields. These four conformational ensembles were compared to available NMR and SAXS data for validation. We found that MD simulations with the AMBER-03w force field and the TIP4P/2005s water model are able to correctly describe the conformational ensemble of this 67-residue IDP at both local and global level.     

Force Field Bias in Protein Folding Simulations

Long timescale (>1 μs) molecular dynamics simulations of protein folding offer a powerful tool for understanding the atomic-scale interactions that determine a protein's folding pathway and stabilize its native state. Unfortunately, when the simulated protein fails to fold, it is often unclear whether the failure is due to a deficiency in the underlying force fields or simply a lack of sufficient simulation time. We examine one such case, the human Pin1 WW domain, using the recently developed deactivated morphing method to calculate free energy differences between misfolded and folded states. We find that the force field we used favors the misfolded states, explaining the failure of the folding simulations. Possible further applications of deactivated morphing and implications for force field development are discussed.     

Abnormal connexin43 in arrhythmogenic right ventricular cardiomyopathy caused by plakophilin-2 mutations

Arrhythmogenic right ventricular cardiomyopathy (ARVC) is a disorder of cardiomyocyte intercalated disk proteins causing sudden death. Heterozygous mutations of the desmosomal protein plakophilin-2 (PKP-2) are the commonest genetic cause of ARVC. Abnormal gap junction connexin43 expression has been reported in autosomal dominant forms of ARVC (Naxos and Carvajal disease) caused by homozygous mutations of desmosomal plakoglobin and desmoplakin. In tissue culture, suppression of PKP-2 results in decreased expression of connexin43. We sought to characterize the expression and localization of connexin43 in patients with ARVC secondary to heterozygous PKP-2 mutations. Complete PKP-2 gene sequencing of 27 ARVC patients was utilized to identify mutant genotypes. Endomyocardial biopsies of identified carriers were then assessed by immunofluorescence to visualize intercalated disk proteins. N-cadherin was targeted to highlight intercalated disks, followed by counterstaining for PKP-2 or connexin43 using confocal double immunofluorescence microscopy. Immunofluorescence was quantified using an AdobeA Photoshop protocol, and colocalization coefficients were determined. PKP-2 siRNA experiments were performed in mouse cardiomyocyte (HL1) cell culture with Western blot analysis to assess connexin43 expression following PKP-2 suppression. Missense and frameshift mutations of the PKP-2 gene were found in four patients with biopsy material available for analysis. Immunofluorescent studies showed PKP-2 localization to the intercalated disk despite mutations, but associated with decreased connexin43 expression and abnormal colocalization. PKP-2 siRNA in HL1 culture confirmed decreased connexin43 expression. Reduced connexin43 expression and localization to the intercalated disk occurs in heterozygous human PKP-2 mutations, potentially explaining the delayed conduction and propensity to develop arrhythmias seen in this disease.