Conformational changes associated with L16P and T118M mutations in the membrane-embedded PMP22 protein,consequential in CMT-1A

Peripheral myelin protein 22 (PMP22) resides in the plasma membrane and is required for myelin formation in the peripheral nervous system. Excess PMP22 mutants accumulate in the endoplasmic reticulum (ER) resulting in the inherited neuropathies of Charcot–Marie–Tooth disease. However, there was no evidence of the structure of PMP22 or how mutations affect its folding. Therefore, in this study, we combined bioinformatics and homology modeling approaches to obtain three-dimensional native and mutated PMP22 models and its anchoring to a POPC membrane, submitted to .5-μs MD simulations, to determine how the L16P and T118M mutations affect the conformational behavior of PMP22. In addition, we investigated the ability of the native and mutated species to accumulate in the ER, via interaction with RER1, by combining protein–protein docking and MD simulations, taking the conformations that were most representative of the native and mutated PMP22 systems and RER1 conformations. Principal component analysis over MD simulations revealed that L16P and T118M mutations resulted in increased structural instability compared to the native form, which is consistent with previous experimental findings of increased structural fluctuations along a loop connecting transmembrane α-helix1 and α-helix2. Docking and MD simulations coupled with the MMGBSA approach allowed the identification that the binding interface for the PMP22-RER1 complex takes place through transmembrane α-helix1 and α-helix2, with higher effective binding free energy values between the mutated PMP22 systems and RER1 than for the native PMP22, mainly through van der Waals interactions.     

Cover Image,Volume 87, Issue 2

Regulator of G protein signaling (RGS) proteins play a pivotal role in regulation of G protein-coupled receptor (GPCR) signaling and are therefore becoming an increasingly important therapeutic target. Recently discovered thiadiazolidinone (TDZD) compounds that target cysteine residues have shown different levels of specificities and potencies for the RGS4 protein, thereby suggesting intrinsic differences in dynamics of this protein upon binding of these compounds. In this work, we investigated using atomistic molecular dynamics (MD) simulations the effect of binding of several small-molecule inhibitors on perturbations and dynamical motions in RGS4. Specifically, we studied two conformational models of RGS4 in which a buried cysteine residue is solvent-exposed due to side-chain motions or due to flexibility in neighboring helices. We found that TDZD compounds with aromatic functional groups perturb the RGS4 structure more than compounds with aliphatic functional groups. Moreover, small-molecules with aromatic functional groups but lacking sulfur atoms only transiently reside within the protein and spontaneously dissociate to the solvent. We further measured inhibitory effects of TDZD compounds using a protein–protein interaction assay on a single-cysteine RGS4 protein showing trends in potencies of compounds consistent with our simulation studies. Thermodynamic analyses of RGS4 conformations in the apo-state and on binding to TDZD compounds revealed links between both conformational models of RGS4. The exposure of cysteine side-chains appears to facilitate initial binding of TDZD compounds followed by migration of the compound into a bundle of four helices, thereby causing allosteric perturbations in the RGS/Gα protein–protein interface.     

Characterization of the homodimerization interface and functional hotspots of the CXCR4 chemokine receptor

The recent crystallographic structures of the human chemokine CXC Receptor 4 (CXCR4) provide experimental evidence of a human G Protein-Coupled Receptor (GPCR) dimer in atomic detail. The CXCR4 homodimers reveal an unexpected dimerization mode involving transmembrane helices TM5 and TM6, which is examined here using all-atom molecular dynamics (MD) simulations in the physiological environment of a lipid bilayer. The bacteriophage T4 lysozyme (T4L), which was fused to the crystallized protein but absent in our simulations, is found to slightly affect the observed relative position of the protomers in the two dimers studied here, and consequently some rearrangements of the dimerization interface are proposed. In addition, the simulations provide further evidence about the role of the two stabilizing single point mutations introduced to crystallize the receptor. Finally, this work analyzes the structural and dynamic role of key residues involved both in ligand binding and in the infection process of HIV. In particular, the different side chain conformations of His113(3.39) are found to influence the dynamics of the surrounding functional hotspot region being evaluated both in the presence and in the absence of the co-crystallized ligand IT1t. The analysis reported here adds valuable knowledge for future structure-based drug design (SBDD) efforts on this pharmacological target.     

NMR resolved multiple anesthetic binding sites in the TM domains of the α4β2 nAChR

The α4β2 nicotinic acetylcholine receptor (nAChR) has significant roles in nervous system function and disease. It is also a molecular target of general anesthetics. Anesthetics inhibit the α4β2 nAChR at clinically relevant concentrations, but their binding sites in α4β2 remain unclear. The recently determined NMR structures of the α4β2 nAChR transmembrane (TM) domains provide valuable frameworks for identifying the binding sites. In this study, we performed solution NMR experiments on the α4β2 TM domains in the absence and presence of halothane and ketamine. Both anesthetics were found in an intra-subunit cavity near the extracellular end of the β2 transmembrane helices, homologous to a common anesthetic binding site observed in X-ray structures of anesthetic-bound GLIC (Nury et al., [32]). Halothane, but not ketamine, was also found in cavities adjacent to the common anesthetic site at the interface of α4 and β2. In addition, both anesthetics bound to cavities near the ion selectivity filter at the intracellular end of the TM domains. Anesthetic binding induced profound changes in protein conformational exchanges. A number of residues, close to or remote from the binding sites, showed resonance signal splitting from single to double peaks, signifying that anesthetics decreased conformation exchange rates. It was also evident that anesthetics shifted population of two conformations. Altogether, the study comprehensively resolved anesthetic binding sites in the α4β2 nAChR. Furthermore, the study provided compelling experimental evidence of anesthetic-induced changes in protein dynamics, especially near regions of the hydrophobic gate and ion selectivity filter that directly regulate channel functions.     

Solution structure of the TatB component of the twin-arginine translocation system

The twin-arginine protein transport (Tat) system translocates fully folded proteins across lipid membranes. In Escherichia coli, the Tat system comprises three essential components: TatA, TatB and TatC. The protein translocation process is proposed to initiate by signal peptide recognition and substrate binding to the TatBC complex. Upon formation of the TatBC–substrate protein complex, the TatA subunits are recruited and form the protein translocation pore. Experimental evidences suggest that TatB forms a tight complex with TatC at 1:1 molar ratio and the TatBC complex contains multiple copies of both proteins. Cross-linking experiments demonstrate that TatB functions in tetrameric units and interacts with both TatC and substrate proteins. However, structural information of the TatB protein is still lacking, and its functional mechanism remains elusive. Herein, we report the solution structure of TatB in DPC micelles determined by Nuclear Magnetic Resonance (NMR) spectroscopy. Overall, the structure shows an extended ‘L-shape’ conformation comprising four helices: a transmembrane helix (TMH) α1, an amphipathic helix (APH) α2, and two solvent exposed helices α3 and α4. The packing of TMH and APH is relatively rigid, whereas helices α3 and α4 display notably higher mobility. The observed floppiness of helices α3 and α4 allows TatB to sample a large conformational space, thus providing high structural plasticity to interact with substrate proteins of different sizes and shapes.     

Probing the structural and energetic basis of kinesin-microtubule binding using computational alanine-scanning mutagenesis

Kinesin-microtubule (MT) binding plays a critical role in facilitating and regulating the motor function of kinesins. To obtain a detailed structural and energetic picture of kinesin-MT binding, we performed large-scale computational alanine-scanning mutagenesis based on long-time molecular dynamics (MD) simulations of the kinesin-MT complex in both ADP and ATP states. First, we built three all-atom kinesin-MT models: human conventional kinesin bound to ADP and mouse KIF1A bound to ADP and ATP. Then, we performed 30 ns MD simulations followed by kinesin-MT binding free energy calculations for both the wild type and mutants obtained after substitution of each charged residue of kinesin with alanine. We found that the kinesin-MT binding free energy is dominated by van der Waals interactions and further enhanced by electrostatic interactions. The calculated mutational changes in kinesin-MT binding free energy are in excellent agreement with results of an experimental alanine-scanning study with a root-mean-square error of ~0.32 kcal/mol [Woehlke, G., et al. (1997) Cell 90, 207-216]. We identified a set of important charged residues involved in the tuning of kinesin-MT binding, which are clustered on several secondary structural elements of kinesin (including well-studied loops L7, L8, L11, and L12, helices α4, α5, and α6, and less-explored loop L2). In particular, we found several key residues that make different contributions to kinesin-MT binding in ADP and ATP states. The mutations of these residues are predicted to fine-tune the motility of kinesin by modulating the conformational transition between the ADP state and the ATP state of kinesin.     

Unveiling the Unfolding Pathway of F5F8D Disorder-Associated D81H/V100D Mutant of MCFD2 via Multiple Molecular Dynamics Simulations

Combined factor deficiency (F5F8D) is a rare autosomal recessive disorder caused by mutations in the LMAN1 or MCFD2 genes. It has been proposed that this pathogenic process occurs via a multi-step pathway involving metal loss, EF-hand-Ca21 dissociation and assembly of misfolded MCFD2-LMAN1 complex. Here, we have investigated the solution conformations of the MCFD2((D81H,V100D)) protein mutant through extensive molecular dynamics (MD) simulations. The V100D, one of the many MCFD2 mutations known to be associated to F5F8D, is difficult to be reconciled with the pathway model because it is located far from the metal sites and the MCFD2/LMAN1 interface. Consequently, an inspection of all the steps involved in D81H/V100D MCFD2 misfolding is expected to provide hints in the understanding of the molecular basis of the disease. A comparison with parallel studies carried out for the Wild-Type (WT) MCFD2 pointed out that the mutation decreases the affinity of the protein for the Ca21 ion. Multiple explicit solvents MD simulations (50_ns) performed on the two proteins revealed that in the WT protein, stable H-bond network and compact hydrophobic core region are created thus confirming a pivotal role of this region in driving the biophysical properties of the entire protein. In fact it is shown that the V100D mutation, although located far away the EF-hand domain, may induce subtle modification in the structural core of MCFD2 leading to the loosening of metal binding and to the formation of metastable intermediate states along the unfolding pathway. The native-like hydrophobic cluster formed near the V100 residue in the wild-type protein is disrupted by the negatively charged Asparagine residue. Furthermore, the presence of the D81H mutation in the EF-1 hand domain may also increase the protein unfolding rate and consequently prevent the formation of the MCFD2-LMAN1 complex. The detailed structural insights obtained from our large-scale simulations complement the clinical features and offer useful insights into the mechanism behind MCFD2 protein misfolding.     

Remodeling of HIV-1 Nef Structure by Src-Family Kinase Binding

The HIV-1 accessory protein Nef controls multiple aspects of the viral life cycle and host immune response, making it an attractive therapeutic target. Previous X-ray crystal structures of Nef in complex with key host cell binding partners have shed light on protein–protein interactions critical to Nef function. Crystal structures of Nef in complex with either the SH3 or tandem SH3–SH2 domains of Src-family kinases reveal distinct dimer conformations of Nef. However, the existence of these Nef dimer complexes in solution has not been established. Here we used hydrogen exchange mass spectrometry (HX MS) to compare the solution conformation of Nef alone and in complexes with the SH3 or the SH3–SH2 domains of the Src-family kinase Hck. HX MS revealed that interaction with the Hck SH3 or tandem SH3–SH2 domains induces protection of the Nef αB-helix from deuterium uptake, consistent with a role for αB in dimer formation. HX MS analysis of a Nef mutant (position Asp123, a site buried in the Nef:SH3 dimer but surface exposed in the Nef:SH3–SH2 complex), showed a Hck-induced conformational change in Nef relative to wild-type Nef. These results support a model in which Src-family kinase binding induces conformational changes in Nef to expose residues critical for interaction with the μ1 subunit of adaptor protein 1 and the major histocompatibility complex-1 tail, and subsequent major histocompatibility complex-1 downregulation and immune escape of HIV-infected cells required for functional interactions with downstream binding partners.     

Stabilization of the integrase‐DNA complex by Mg2+ ions and prediction of key residues for binding HIV‐1 integrase inhibitors

The HIV‐1 integrase is an attractive target for the therapeutics development against AIDS, as no host homologue of this protein has been identified. The integrase strand transfer inhibitors (INSTIs), including raltegravir, specifically target the second catalytic step of the integration process by binding to the DDE motif of the catalytic site and coordinating Mg²⁺ ions. Recent X‐ray crystallographic structures of the integrase/DNA complex from prototype foamy virus allowed to investigate the role of the different partners (integrase, DNA, Mg²⁺ ions, raltegravir) in the complex stability using molecular dynamics (MD) simulations. The presence of Mg²⁺ ions is found to be essential for the stability, whereas the simultaneous presence of raltegravir and Mg²⁺ ions has a destabilizing influence. A homology model of HIV‐1 integrase was built on the basis of the X‐ray crystallographic information, and protein marker residues for the ligand binding were detected by clustering the docking poses of known HIV‐1 integrase inhibitors on the model. Interestingly, we had already identified some of these residues to be involved in HIV‐1 resistance mutations and in the stabilization of the catalytic site during the MD simulations. Classification of protein conformations along MD simulations, as well as of ligand docking poses, was performed by using an original learning method, based on self‐organizing maps. This allows us to perform a more in‐depth investigation of the free‐energy basins populated by the complex in MD simulations on the one hand, and a straightforward classification of ligands according to their binding residues on the other hand. Proteins 2014; 82:466–478. © 2013 Wiley Periodicals, Inc.     

Extensive domain motion and electron transfer in the human electron transferring flavoprotein.medium chain Acyl-CoA dehydrogenase complex

The crystal structure of the human electron transferring flavoprotein (ETF).medium chain acyl-CoA dehydrogenase (MCAD) complex reveals a dual mode of protein-protein interaction, imparting both specificity and promiscuity in the interaction of ETF with a range of structurally distinct primary dehydrogenases. ETF partitions the functions of partner binding and electron transfer between (i) the recognition loop, which acts as a static anchor at the ETF.MCAD interface, and (ii) the highly mobile redox active FAD domain. Together, these enable the FAD domain of ETF to sample a range of conformations, some compatible with fast interprotein electron transfer. Disorders in amino acid or fatty acid catabolism can be attributed to mutations at the protein-protein interface. Crucially, complex formation triggers mobility of the FAD domain, an induced disorder that contrasts with general models of protein-protein interaction by induced fit mechanisms. The subsequent interfacial motion in the MCAD.ETF complex is the basis for the interaction of ETF with structurally diverse protein partners. Solution studies using ETF and MCAD with mutations at the protein-protein interface support this dynamic model and indicate ionic interactions between MCAD Glu(212) and ETF Arg alpha(249) are likely to transiently stabilize productive conformations of the FAD domain leading to enhanced electron transfer rates between both partners.     

Identification of interacting hot spots in the beta3 integrin stalk using comprehensive interface design

Protein-protein interfaces are usually large and complementary surfaces, but specific side chains, representing energetic "hot spots," often contribute disproportionately to binding free energy. We used a computational method, comprehensive interface design, to identify hot spots in the interface between the stalk regions of the β3 and the complementary αIIb and αv integrin subunits. Using the Rosetta alanine-scanning and design algorithms to predict destabilizing, stabilizing, and neutral mutations in the β3 region extending from residues Lys(532) through Gly(690), we predicted eight alanine mutations that would destabilize the αIIbβ3 interface as well as nine predicted to destabilize the αvβ3 interface, by at least 0.3 kcal/mol. The mutations were widely and unevenly distributed, with four between residues 552 and 563 and five between 590 and 610, but none between 565 and 589, and 611 and 655. Further, mutations destabilizing the αvβ3 and αIIbβ3 interfaces were not identical. The predictions were then tested by introducing selected mutations into the full-length integrins expressed in Chinese hamster ovary cells. Five mutations predicted to destabilize αIIb and β3 caused fibrinogen binding to αIIbβ3, whereas three of four predicted to be neutral or stabilizing did not. Conversely, a mutation predicted to destabilize αvβ3, but not αIIbβ3 (D552A), caused osteopontin binding to αvβ3, but not fibrinogen binding to αIIbβ3. These results indicate that stability of the distal stalk interface is involved in constraining integrins in stable, inactive conformations. Further, they demonstrate the ability of comprehensive interface design to identify functionally significant integrin mutations.     

Modulation of allosteric behavior through adjustment of the differential stability of the two interacting domains in E. coli cAMP receptor protein

The communication mechanism(s) responsible for the allosteric behavior of E.coli cAMP binding receptor protein, CRP, is still a subject of intense investigation. As a tool to explore the communication mechanism, the mutations at various positions in the cAMP-binding (K52N, D53H, S62F and T127L) or the DNA- binding (H159L) domain or both (K52N/H159L) were generated. The sites and specific nature of side chain substitutions were defined by earlier genetic studies, the results of which show that these mutants have a similar phenotype i.e. they are activated without exogenous cAMP. Presently, no significant changes in the structures of WT and mutant CRPs have been observed. Hence, the pressing issue is to identify a physical parameter that reflects the effects of mutations. In this study, the stability of these various CRP species in the presence of GuHCl was monitored by three spectroscopic techniques, namely, CD, tryptophan fluorescence and FT-IR which could provide data on the stability of α-helices and β-strands separately. Results of this study led to the following conclusions: 1. The α-helices can be grouped into two families with different stabilities. Mutations exert a differential effect on the stability of helices as demonstrated by a biphasic unfolding curve for the helices. 2. Regardless of the locations of mutations, the effects can be communicated to the other domain resulting in a perturbation of the stability of both domains, although the effects are more significantly expressed in the stability of the helices. 3. Although in an earlier study [Gekko, et al. Biochemistry 43 (2004) 3844] we showed that cooperativity of cAMP binding is generally correlated to the global dynamics of the protein and DNA binding affinity, in this study we found that generally there is no clear correlation between functional energetics and stability of secondary structures. Thus, results of this study imply that modulation of allostery in CRP is entropic in nature.     

Investigation of Binding Phenomenon of NSP3 and p130Cas Mutants and Their Effect on Cell Signalling

Members of the novel SH2-containing protein (NSP3) and Crk-associated substrate (p130Cas) protein families form a multi-domain signalling platforms that mediate cell signalling process. We analysed the damaging consequences of three mutations, each from NSP3 (NSP3^L469R, NSP3^L623E, NSP3^R627E) and p130Cas (p130Cas^F794R, p130Cas^L787E, p130Cas^D797R) protein with respect to their native biological partners. Mutations depicted notable loss in interaction affinity towards their corresponding biological partners. NSP3^L469R and p130Cas^D797R mutations were predicted as most prominent in docking analysis. Molecular dynamics (MD) studies were conducted to evaluate structural consequences of most prominent mutation in NSP3 and p130Cas obtained from the docking analysis. MD analysis confirmed that mutation in NSP3^L469R and p130Cas^D797R showed significant structural deviation, changes in conformations and increased flexibility, which in turn affected the binding affinity with their biological partners. Moreover, the root mean square fluctuation has indicated a rise in fluctuation of residues involved in moderate interaction acquired between the NSP3 and p130Cas. It has significantly affected the binding interaction in mutant complexes. The results obtained in this work present a detailed overview of molecular mechanisms involved in the loss of cell signalling associated with NSP3 and p130Cas protein.     

Expanding the Druggable Space of the LSD1/CoREST Epigenetic Target: New Potential Binding Regions for Drug-Like Molecules,Peptides, Protein Partners,and Chromatin

Lysine specific demethylase-1 (LSD1/KDM1A) in complex with its corepressor protein CoREST is a promising target for epigenetic drugs. No therapeutic that targets LSD1/CoREST, however, has been reported to date. Recently, extended molecular dynamics (MD) simulations indicated that LSD1/CoREST nanoscale clamp dynamics is regulated by substrate binding and highlighted key hinge points of this large-scale motion as well as the relevance of local residue dynamics. Prompted by the urgent need for new molecular probes and inhibitors to understand LSD1/CoREST interactions with small-molecules, peptides, protein partners, and chromatin, we undertake here a configurational ensemble approach to expand LSD1/CoREST druggability. The independent algorithms FTMap and SiteMap and our newly developed Druggable Site Visualizer (DSV) software tool were used to predict and inspect favorable binding sites. We find that the hinge points revealed by MD simulations at the SANT2/Tower interface, at the SWIRM/AOD interface, and at the AOD/Tower interface are new targets for the discovery of molecular probes to block association of LSD1/CoREST with chromatin or protein partners. A fourth region was also predicted from simulated configurational ensembles and was experimentally validated to have strong binding propensity. The observation that this prediction would be prevented when using only the X-ray structures available (including the X-ray structure bound to the same peptide) underscores the relevance of protein dynamics in protein interactions. A fifth region was highlighted corresponding to a small pocket on the AOD domain. This study sets the basis for future virtual screening campaigns targeting the five novel regions reported herein and for the design of LSD1/CoREST mutants to probe LSD1/CoREST binding with chromatin and various protein partners.