Critical Structural Defects Explain Filamin A Mutations Causing Mitral Valve Dysplasia.

Mitral valve diseases affect ∼3% of the population and are the most common reasons for valvular surgery because no drug-based treatments exist. Inheritable genetic mutations have now been established as the cause of mitral valve insufficiency, and four different missense mutations in the filamin A gene (FLNA) have been found in patients suffering from nonsyndromic mitral valve dysplasia (MVD). The filamin A (FLNA) protein is expressed, in particular, in endocardial endothelia during fetal valve morphogenesis and is key in cardiac development. The FLNA-MVD-causing mutations are clustered in the N-terminal region of FLNA. How the mutations in FLNA modify its structure and function has mostly remained elusive. In this study, using NMR spectroscopy and interaction assays, we investigated FLNA-MVD-causing V711D and H743P mutations. Our results clearly indicated that both mutations almost completely destroyed the folding of the FLNA5 domain, where the mutation is located, and also affect the folding of the neighboring FLNA4 domain. The structure of the neighboring FLNA6 domain was not affected by the mutations. These mutations also completely abolish FLNA's interactions with protein tyrosine phosphatase nonreceptor type 12, which has been suggested to contribute to the pathogenesis of FLNA-MVD. Taken together, our results provide an essential structural and molecular framework for understanding the molecular bases of FLNA-MVD, which is crucial for the development of new therapies to replace surgery.

A novel structural unit in the N-terminal region of filamins.

Immunoglobulin-like (Ig) domains are a widely expanded superfamily that act as interaction motifs or as structural spacers in multidomain proteins. Vertebrate filamins (FLNs), which are multifunctional actin-binding proteins, consist of 24 Ig domains. We have recently discovered that in the C-terminal rod 2 region of FLN, Ig domains interact with each other forming functional domain pairs, where the interaction with signaling and transmembrane proteins is mechanically regulated by weak actomyosin contraction forces. Here, we investigated if there are similar inter-domain interactions around domain 4 in the N-terminal rod 1 region of FLN. Protein crystal structures revealed a new type of domain organization between domains 3, 4, and 5. In this module, domains 4 and 5 interact rather tightly, whereas domain 3 has a partially flexible interface with domain 4. NMR peptide titration experiments showed that within the three-domain module, domain 4 is capable for interaction with a peptide derived from platelet glycoprotein Ib. Crystal structures of FLN domains 4 and 5 in complex with the peptide revealed a typical ß sheet augmentation interaction observed for many FLN ligands. Domain 5 was found to stabilize domain 4, and this could provide a mechanism for the regulation of domain 4 interactions.

Molecular mechanism of T-cell protein tyrosine phosphatase (TCPTP) activation by mitoxantrone.

T-cell protein tyrosine phosphatase (TCPTP) is a ubiquitously expressed non-receptor protein tyrosine phosphatase. It is involved in the negative regulation of many cellular signaling pathways. Thus, activation of TCPTP could have important therapeutic applications in diseases such as cancer and inflammation. We have previously shown that the α-cytoplasmic tail of integrin α1ß1 directly binds and activates TCPTP. In addition, we have identified in a large-scale high-throughput screen six small molecules that activate TCPTP. These small molecule activators include mitoxantrone and spermidine. In this study, we have investigated the molecular mechanism behind agonist-induced TCPTP activation. By combining several molecular modeling and biochemical techniques, we demonstrate that α1-peptide and mitoxantrone activate TCPTP via direct binding to the catalytic domain, whereas spermidine does not interact with the catalytic domain of TCPTP in vitro. Furthermore, we have identified a hydrophobic groove surrounded by negatively charged residues on the surface of TCPTP as a putative binding site for the α1-peptide and mitoxantrone. Importantly, these data have allowed us to identify a new molecule that binds to TCPTP, but interestingly cannot activate its phosphatase activity. Accordingly, we describe here mechanism of TCPTP activation by mitoxantrone, the cytoplasmic tail of α1-integrin, and a mitoxantrone-like molecule at the atomic level. These data provide invaluable insight into the development of novel TCPTP activators, and may facilitate the rational discovery of small-molecule cancer therapeutics.

The C-terminal rod 2 fragment of filamin A forms a compact structure that can be extended.

Filamins are large proteins that cross-link actin filaments and connect to other cellular components. The C-terminal rod 2 region of FLNa (filamin A) mediates dimerization and interacts with several transmembrane receptors and intracellular signalling adaptors. SAXS (small-angle X-ray scattering) experiments were used to make a model of a six immunoglobulin-like domain fragment of the FLNa rod 2 (domains 16-21). This fragment had a surprising three-branched structural arrangement, where each branch was made of a tightly packed two-domain pair. Peptides derived from transmembrane receptors and intracellular signalling proteins induced a more open structure of the six domain fragment. Mutagenesis studies suggested that these changes are caused by peptides binding to the CD faces on domains 19 and 21 which displace the preceding domain A-strands (18 and 20 respectively), thus opening the individual domain pairs. A single particle cryo-EM map of a nine domain rod 2 fragment (domains 16-24), showed a relatively compact dimeric particle and confirmed the three-branched arrangement as well as the peptide-induced conformation changes. These findings reveal features of filamin structure that are important for its interactions and mechanical properties.

Assembly of a filamin four-domain fragment and the influence of splicing variant-1 on the structure.

Filamins are scaffold proteins that bind to various proteins, including the actin cytoskeleton, integrin adhesion receptors, and adaptor proteins such as migfilin. Alternative splicing of filamin, largely constructed from 24 Ig-like domains, is thought to have a role in regulating its interactions with other proteins. The filamin A splice variant-1 (FLNa var-1) lacks 41 amino acids, including the last ß-strand of domain 19, FLNa(19), and the first ß-strand of FLNa(20) that was previously shown to mask a key binding site on FLNa(21). Here, we present a structural characterization of domains 18-21, FLNa(18-21), in the FLNa var-1 as well as its nonspliced counterpart. A model of nonspliced FLNa(18-21), obtained from small angle x-ray scattering data, shows that these four domains form an L-shaped structure, with one arm composed of a pair of domains. NMR spectroscopy reveals that in the splice variant, FLNa(19) is unstructured whereas the other domains retain the same fold as in their canonical counterparts. The maximum dimensions predicted by small angle x-ray scattering data are increased upon migfilin binding in the FLNa(18-21) but not in the splice variant, suggesting that migfilin binding is able to displace the masking ß-strand and cause a rearrangement of the structure. Possible function roles for the spliced variants are discussed.

Interaction mechanism of endogenous PP2A inhibitor protein ENSA with PP2A.

The vast diversity of protein phosphatase 2A (PP2A) holoenzyme composition ensures its multifaceted role in the regulation of cellular growth and signal transduction. In several pathological conditions, such as cancer, PP2A is inhibited by endogenous inhibitor proteins. Several PP2A inhibitor proteins have been identified, one of which is α-endosulfine (ENSA). ENSA inhibits PP2A activity when it is phosphorylated at Ser67 by Greatwall (Gwl) kinase. The role of ENSA in PP2A inhibition is rather well characterized, but knowledge of the mechanism of inhibition is scarce. In this study, we have performed comprehensive structural characterization of ENSA, and its interaction with PP2A A- and various B56-subunit isoforms by combining NMR spectroscopy, small-angle X-ray scattering (SAXS) and interaction assays. The results clearly indicate that ENSA is an intrinsically disordered protein containing three transient α-helical structures. ENSA was observed to interact PP2A mainly via A-subunit, as the affinity with the A-subunit is significantly stronger than with any of the B56 subunits. Based on our results, it seems that ENSA follows the dock-and-coalesce mechanism in associating with PP2A A-subunit. Taken together, our results provide an essential structural and molecular framework to understanding molecular bases of ENSA-mediated PP2A inhibition, which is crucial for the development of new therapies for diseases linked to PP2A inhibition.

The Interaction Mechanism of Intrinsically Disordered PP2A Inhibitor Proteins ARPP-16 and ARPP-19 With PP2A.

Protein phosphatase 2A (PP2A) activity is critical for maintaining normal physiological cellular functions. PP2A is inhibited by endogenous inhibitor proteins in several pathological conditions including cancer. A PP2A inhibitor protein, ARPP-19, has recently been connected to several human cancer types. Accordingly, the knowledge about ARPP-19-PP2A inhibition mechanism is crucial for the understanding the disease development and the therapeutic targeting of ARPP-19-PP2A. Here, we show the first structural characterization of ARPP-19, and its splice variant ARPP-16 using NMR spectroscopy, and SAXS. The results reveal that both ARPP proteins are intrinsically disordered but contain transient secondary structure elements. The interaction mechanism of ARPP-16/19 with PP2A was investigated using microscale thermophoresis and NMR spectroscopy. Our results suggest that ARPP-PP2A A-subunit interaction is mediated by linear motif and has modest affinity whereas, the interaction of ARPPs with B56-subunit is weak and transient. Like many IDPs, ARPPs are promiscuous binders that transiently interact with PP2A A- and B56 subunits using multiple interaction motifs. In summary, our results provide a good starting point for future studies and development of therapeutics that block ARPP-PP2A interactions.

1H, 13C and 15N NMR chemical shift assignments of cAMP-regulated phosphoprotein-19 and -16 (ARPP-19 and ARPP-16).

Protein Phosphatase 2A, PP2A, the principal Serine/threonine phosphatase, has major roles in broad range of signaling pathways that include regulation of cell cycle, cell proliferation and neuronal signaling. The loss of function of PP2A is linked with many human diseases, like cancer and neurodegenerative disorders. Protein phosphatase 2A (PP2A) functions as tumor suppressor and its tumor suppressor activity is inhibited by the overexpression of PP2A inhibitor proteins in most of the cancers. ARPP-19/ARPP-16 has been identified as one of the potential PP2A inhibitor proteins. Here, we report the resonance assignment of backbone 1H, 13C and 15N atoms of human ARPP-19 and ARPP-16 proteins. These chemical shift values can provide valuable information for the further study of the dynamics and interaction of ARPP-proteins to PP2A using NMR spectroscopy.

Non-syndromic Mitral Valve Dysplasia Mutation Changes the Force Resilience and Interaction of Human Filamin A.

Filamin A (FLNa), expressed in endocardial endothelia during fetal valve morphogenesis, is key in cardiac development. Missense mutations in FLNa cause non-syndromic mitral valve dysplasia (FLNA-MVD). Here, we aimed to reveal the currently unknown underlying molecular mechanism behind FLNA-MVD caused by the FLNa P637Q mutation. The solved crystal structure of the FLNa3-5 P637Q revealed that this mutation causes only minor structural changes close to mutation site. These changes were observed to significantly affect FLNa's ability to transmit cellular force and to interact with its binding partner. The performed steered molecular dynamics simulations showed that significantly lower forces are needed to split domains 4 and 5 in FLNA-MVD than with wild-type FLNa. The P637Q mutation was also observed to interfere with FLNa's interactions with the protein tyrosine phosphatase PTPN12. Our results provide a crucial step toward understanding the molecular bases behind FLNA-MVD, which is critical for the development of drug-based therapeutics.

Cooperative symmetric to asymmetric conformational transition of the apo-form of scavenger decapping enzyme revealed by simulations.

Decapping is a central step in eukaryotic mRNA turnover and in gene expression regulation. The human scavenger decapping enzyme, DcpS, catalyses cap hydrolysis following mRNA degradation. DcpS is a dimeric enzyme, with two active sites. Crystal structures suggest that DcpS must undergo significant conformational changes upon ligand binding, but the mechanism of this transition is unknown. Here, we report two long timescale (20 ns) molecular dynamics simulations of the apo-form of DcpS. The dimer is observed to undergo a strikingly cooperative motion, with one active site closing while the other opens. The amplitude of the conformational change is 6-21 A and the apparent timescale is 4-13 ns. These findings indicate that the crystallographically observed symmetric conformation of apo-form of DcpS is only a minor conformation in solution. The simulations also show that active sites are structurally connected via the domain-swapped dimer structure of the N-terminal domain, even in the absence of a bound ligand. These findings suggest a functional reason for the enzyme existing as a dimer, and may be widely relevant, also for other dimeric proteins.

Structural and Functional Insights Into Lysostaphin-Substrate Interaction.

Lysostaphin from Staphylococcus simulans and its family enzymes rapidly acquire prominence as the next generation agents in treatment of S. aureus infections. The specificity of lysostaphin is promoted by its C-terminal cell wall targeting domain selectivity toward pentaglycine bridges in S. aureus cell wall. Scission of these cross-links is carried out by its N-terminal catalytic domain, a zinc-dependent endopeptidase. Understanding the determinants affecting the efficiency of catalysis and strength and specificity of interactions lies at the heart of all lysostaphin family enzyme applications. To this end, we have used NMR, SAXS and molecular dynamics simulations to characterize lysostaphin structure and dynamics, to address the inter-domain interaction, the enzyme-substrate interaction as well as the catalytic properties of pentaglycine cleavage in solution. Our NMR structure confirms the recent crystal structure, yet, together with the molecular dynamics simulations, emphasizes the dynamic nature of the loops embracing the catalytic site. We found no evidence for inter-domain interaction, but, interestingly, the SAXS data delineate two preferred conformation subpopulations. Catalytic H329 and H360 were observed to bind a second zinc ion, which reduces lysostaphin pentaglycine cleaving activity. Binding of pentaglycine or its lysine derivatives to the targeting domain was found to be of very low affinity. The pentaglycine interaction site was located to the N-terminal groove of the domain. Notably, the targeting domain binds the peptidoglycan stem peptide Ala-d-γ-Glu-Lys-d-Ala-d-Ala with a much higher, micromolar affinity. Binding site mapping reveals two interaction sites of different affinities on the surface of the domain for this peptide.

Skeletal Dysplasia Mutations Effect on Human Filamins' Structure and Mechanosensing.

Cells' ability to sense mechanical cues in their environment is crucial for fundamental cellular processes, leading defects in mechanosensing to be linked to many diseases. The actin cross-linking protein Filamin has an important role in the conversion of mechanical forces into biochemical signals. Here, we reveal how mutations in Filamin genes known to cause Larsen syndrome and Frontometaphyseal dysplasia can affect the structure and therefore function of Filamin domains 16 and 17. Employing X-ray crystallography, the structure of these domains was first solved for the human Filamin B. The interaction seen between domains 16 and 17 is broken by shear force as revealed by steered molecular dynamics simulations. The effects of skeletal dysplasia associated mutations of the structure and mechanosensing properties of Filamin were studied by combining various experimental and theoretical techniques. The results showed that Larsen syndrome associated mutations destabilize or even unfold domain 17. Interestingly, those Filamin functions that are mediated via domain 17 interactions with other proteins are not necessarily affected as strongly interacting peptide binding to mutated domain 17 induces at least partial domain folding. Mutation associated to Frontometaphyseal dysplasia, in turn, transforms 16-17 fragment from compact to an elongated form destroying the force-regulated domain pair.

Flexible Structure of Peptide-Bound Filamin A Mechanosensor Domain Pair 20-21.

Filamins (FLNs) are large, multidomain actin cross-linking proteins with diverse functions. Besides regulating the actin cytoskeleton, they serve as important links between the extracellular matrix and the cytoskeleton by binding cell surface receptors, functioning as scaffolds for signaling proteins, and binding several other cytoskeletal proteins that regulate cell adhesion dynamics. Structurally, FLNs are formed of an amino terminal actin-binding domain followed by 24 immunoglobulin-like domains (IgFLNs). Recent studies have demonstrated that myosin-mediated contractile forces can reveal hidden protein binding sites in the domain pairs IgFLNa18-19 and 20-21, enabling FLNs to transduce mechanical signals in cells. The atomic structures of these mechanosensor domain pairs in the resting state are known, as well as the structures of individual IgFLN21 with ligand peptides. However, little experimental data is available on how interacting protein binding deforms the domain pair structures. Here, using small-angle x-ray scattering-based modelling, x-ray crystallography, and NMR, we show that the adaptor protein migfilin-derived peptide-bound structure of IgFLNa20-21 is flexible and adopts distinctive conformations depending on the presence or absence of the interacting peptide. The conformational changes reported here may be common for all peptides and may play a role in the mechanosensor function of the site.

Novel hydrazine molecules as tools to understand the flexibility of vascular adhesion protein-1 ligand-binding site: toward more selective inhibitors.

Vascular adhesion protein-1 (VAP-1) belongs to a family of amine oxidases. It plays a role in leukocyte trafficking and in amine compound metabolism. VAP-1 is linked to various diseases, such as Alzheimer's disease, psoriasis, depression, diabetes, and obesity. Accordingly, selective inhibitors of VAP-1 could potentially be used to treat those diseases. In this study, eight novel VAP-1 hydrazine derivatives were synthesized and their VAP-1 and monoamine oxidase (MAO) inhibition ability was determined in vitro. MD simulations of VAP-1 with these new molecules reveal that the VAP-1 ligand-binding pocket is flexible and capable of fitting substantially larger ligands than was previously believed. The increase in the size of the VAP-1 ligands, together with the methylation of the secondary nitrogen atom of the hydrazine moiety, improves the VAP-1 selectivity over MAO.

Synthesis, in vitro activity, and three-dimensional quantitative structure-activity relationship of novel hydrazine inhibitors of human vascular adhesion protein-1.

Vascular adhesion protein-1 (VAP-1) belongs to the semicarbazide-sensitive amine oxidases (SSAOs) that convert amines into aldehydes. SSAOs are distinct from the mammalian monoamine oxidases (MAOs), but their substrate specificities are partly overlapping. VAP-1 has been proposed as a target for anti-inflammatory drug therapy because of its role in leukocyte adhesion to endothelium. Here, we describe the synthesis and in vitro activities of novel series of VAP-1 selective inhibitors. In addition, the molecular dynamics simulations performed for VAP-1 reveal that the movements of Met211, Ser496, and especially Leu469 can enlarge the ligand-binding pocket, allowing larger ligands than those seen in the crystal structures to bind. Combining the data from molecular dynamics simulations, docking, and in vitro measurements, the three-dimensional quantitative structure-activity relationship (3D QSAR) models for VAP-1 (q(2)(LOO): 0.636; r(2): 0.828) and MAOs (q(2)(LOO): 0.749, r(2): 0.840) were built and employed in the development of selective VAP-1 inhibitors.

The regulation mechanism for the auto-inhibition of binding of human filamin A to integrin.

The ability of adhesion receptors to transmit biochemical signals and mechanical force across cell membranes depends on interactions with the actin cytoskeleton. Human filamins are large actin cross-linking proteins that connect integrins to the cytoskeleton. Filamin binding to the cytoplasmic tail of beta integrins has been shown to prevent integrin activation in cells, which is important for controlling cell adhesion and migration. The molecular-level mechanism for filamin binding to integrin has been unclear, however, as it was recently demonstrated that filamin undergoes intramolecular auto-inhibition of integrin binding. In this study, using steered molecular dynamics simulations, we found that mechanical force applied to filamin can expose cryptic integrin binding sites. The forces required for this are considerably lower than those for filamin immunoglobulin domain unfolding. The mechanical-force-induced unfolding of filamin and exposure of integrin binding sites occur through stable intermediates where integrin binding is possible. Accordingly, our results support filamin's role as a mechanotransducer, since force-induced conformational changes allow binding of integrin and other transmembrane and intracellular proteins. This observed force-induced conformational change can also be one of possible mechanisms involved in the regulation of integrin activation.

Lennard-Jones Parameters for B3LYP/CHARMM27 QM/MM Modeling of Nucleic Acid Bases.

Combined quantum mechanics/molecular mechanics (QM/MM) methods allow computations on chemical events in large molecular systems. Here, we have tested the suitability of the standard CHARMM27 forcefield Lennard-Jones van der Waals (vdW) parameters for the treatment of nucleic acid bases in QM/MM calculations at the B3LYP/6-311+G(d,p)-CHARMM27 level. Alternative parameters were also tested by comparing the QM/MM hydrogen bond lengths and interaction energies with full QM [B3LYP/6-311+G(d,p)] results. The optimization of vdW parameters for nucleic acid bases is challenging because of the likelihood of multiple hydrogen bonds between the nucleic acid base and a water molecule. Two sets of optimized atomic vdW parameters for polar hydrogen, carbonyl carbon, and aromatic nitrogen atoms for nucleic acid bases are reported: base-dependent and base-independent. The results indicate that, for QM/MM investigations of nucleic acids, the standard forcefield vdW parameters may not be appropriate for atoms treated by QM. QM/MM interaction energies calculated with standard CHARMM27 parameters are found to be too large, by around 3 kcal/mol. This is because of overestimation of electrostatic interactions. Interaction energies closer to the full QM results are found using the optimized vdW parameters developed here. The optimized vdW parameters [developed by reference to B3LYP/6-311+G(d,p) results] were also tested at the B3LYP/6-31G(d) QM/MM level and were found to be transferable to the lower level. The optimized parameters also model the interaction energies of charged nucleic acid bases and deprotonation energies reasonably well.

Subtype selectivity and flexibility of ionotropic glutamate receptors upon antagonist ligand binding.

The binding modes of a set of known ionotropic glutamate receptor antagonist-ligands have been studied using homology modeling, molecular docking, molecular dynamics (MD) simulations and ab initio quantum mechanical calculations. The core structure of the studied ligands is the decahydroisoquinoline ring, which has a carboxylic acid group at position three and different negatively-charged substituents (R) at position six. The binding affinities of these molecules have been reported earlier. From the current study, the carboxylate group of the decahydroisoquinoline ring hydrogen bonds with Arg485, the amino group with Pro478 and Thr480, and the negatively charged substituent R interacts with the positively charged N-terminus of helix-F. The subtype selectivity of these ligands seems to be strongly dependent on the amino acid at position 650 (GluR2: leucine, GluR5: valine), which affects the conformation of the ligand and ligand-receptor interactions, but depends considerably on the size of the R-group of the ligand. In addition, the MD simulations also revealed that the relative positions of the S1 and S2 domains can alter significantly showing different "closure" and "rotational movements" depending on the antagonist-ligand that is bound. Accordingly, molecular docking of antagonist ligands into static crystal structures cannot sufficiently explain ligand binding and subtype selectivity.

Comment on "Molecular dynamics DFT:B3LYP study of guanosinetriphosphate conversion into guanosinemonophosphate upon Mg2+ chelation of alpha and beta phosphate oxygens of the triphosphate tail" by Alexander A. Tulub, Phys. Chem. Chem. Phys., 2006, 8, 2187.

A recent paper in this journal uses molecular dynamics methods to study hydrolysis of guanosine triphosphate (GTP). The author reports that cleavage of the molecule occurs in less than 5 ps, and leads to a number of fragments including a free oxygen atom and a reduced magnesium ion. This conclusion is not in agreement with the known biochemistry and chemical reactivity of GTP or with previous computational studies of its hydrolysis reaction.

Determination of binding site residues responsible for the subunit selectivity of novel marine-derived compounds on kainate receptors.

Dysiherbaine (DH) and related molecules are high-affinity, subunit-selective kainate receptor (KAR) ligands originally isolated from a marine sponge. To elucidate why DH, an agonist, and MSVIII-19, a competitive antagonist, bind selectively to glutamate receptor (GluR) 5 but not to the KA2 KAR subunit, we used molecular dynamics simulations to generate binding models that were tested experimentally in radioligand binding and electrophysiological assays. Three candidate sites, Val685, Leu735, and Ser741 in GluR5, corresponding to Ile669, Phe719, and Met725 in KA2, were predicted to underlie the distinct binding profiles of the marine toxins. Single or multiple reciprocal mutations introduced into the receptor subunits produced a variety of effects on binding affinity. Most notably, mutation of Met725 to serine in KA2 increased the affinity of DH by 350-fold; in contrast, mutation of one or more of the residues in GluR5 did not markedly alter DH binding. MSVIII-19 affinity for the KA2 subunit was significantly increased in multiple site mutants, and reciprocal mutations in the GluR5 subunit produced substantial (700-fold) reductions in MSVIII-19 affinity. Physiological characterization of the double- and triple-mutant subunits demonstrated altered functional behavior consistent with the changes in binding affinity. The results provide experimental support for the importance of these three ligand binding domain (LBD) residues and suggest steric hindrance in the KA2 subunit LBD is largely responsible for the very low affinity for the two compounds. In this study, we identified the molecular basis for subunit selectivity of these marine-derived molecules on KARs, which could facilitate the rational design of selective ligands with distinct pharmacological profiles.