Ceramide-1-phosphate transfer protein promotes sphingolipid reorientation needed for binding during membrane interaction

Lipid transfer proteins acquire and release their lipid cargoes by interacting transiently with source and destination biomembranes. In the GlycoLipid Transfer Protein (GLTP) superfamily, the two-layer all-α-helical GLTP-fold defines proteins that specifically target sphingolipids (SLs) containing either sugar or phosphate headgroups via their conserved but evolutionarily-modified SL recognitions centers. Despite comprehensive structural insights provided by X-ray crystallography, the conformational dynamics associated with membrane interaction and SL uptake/release by GLTP superfamily members have remained unknown. Herein, we report insights gained from molecular dynamics (MD) simulations into the conformational dynamics that enable ceramide-1-phosphate transfer proteins (CPTPs) to acquire and deliver ceramide-1-phosphate (C1P) during interaction with 1-palmitoyl-2-oleoyl phosphatidylcholine bilayers. The focus on CPTP reflects this protein’s involvement in regulating pro-inflammatory eicosanoid production and autophagy-dependent inflammasome assembly that drives interleukin (IL-1β and IL-18) production and release by surveillance cells. We found that membrane penetration by CPTP involved α-6 helix and the α-2 helix N-terminal region, was confined to one bilayer leaflet, and was relatively shallow. Large-scale dynamic conformational changes were minimal for CPTP during membrane interaction or C1P uptake except for the α-3/α-4 helices connecting loop, which is located near the membrane interface and interacts with certain phosphoinositide headgroups. Apart from functioning as a shallow membrane-docking element, α-6 helix was found to adeptly reorient membrane lipids to help guide C1P hydrocarbon chain insertion into the interior hydrophobic pocket of the SL binding site.These findings support a proposed ‘hydrocarbon chain-first’ mechanism for C1P uptake, in contrast to the ‘lipid polar headgroup-first’ uptake used by most lipid-transfer proteins.     

Combining protein sequence,structure, and dynamics: A novel approach for functional evolution analysis of PAS domain superfamily

PAS domains are widespread in archaea, bacteria, and eukaryota, and play important roles in various functions. In this study, we aim to explore functional evolutionary relationship among proteins in the PAS domain superfamily in view of the sequence‐structure‐dynamics‐function relationship. We collected protein sequences and crystal structure data from RCSB Protein Data Bank of the PAS domain superfamily belonging to three biological functions (nucleotide binding, photoreceptor activity, and transferase activity). Protein sequences were aligned and then used to select sequence‐conserved residues and build phylogenetic tree. Three‐dimensional structure alignment was also applied to obtain structure‐conserved residues. The protein dynamics were analyzed using elastic network model (ENM) and validated by molecular dynamics (MD) simulation. The result showed that the proteins with same function could be grouped by sequence similarity, and proteins in different functional groups displayed statistically significant difference in their vibrational patterns. Interestingly, in all three functional groups, conserved amino acid residues identified by sequence and structure conservation analysis generally have a lower fluctuation than other residues. In addition, the fluctuation of conserved residues in each biological function group was strongly correlated with the corresponding biological function. This research suggested a direct connection in which the protein sequences were related to various functions through structural dynamics. This is a new attempt to delineate functional evolution of proteins using the integrated information of sequence, structure, and dynamics.     

Computational studies on mutant protein stability: The correlation between surface thermal expansion and protein stability

Thermal stability of mutant proteins has been investigated using temperature dependent molecular dynamics (MD) simulations in vacuo. The numerical modeling was aimed at mimicking protein expansion upon heating. After the conditions for an expanding protein accessible surface area were established for T4 lysozyme and barnase wild-type proteins, MD simulations were carried out under the same conditions using the crystal structures of several mutant proteins. The computed thermal expansion of the accessible surface area of mutant proteins was found to be strongly correlated with their experimentally measured stabilities. A similar, albeit weaker, correlation was observed for model mutant proteins. This opens the possibility of obtaining stability information directly from protein structure.     

Theoretical Model of the Three-dimensional Structure of a Disease Resistance Gene Homolog Encoding Resistance Protein in Vigna mungo

Abstract

Plant disease resistance (R) genes, the key players of innate immunity system in plants encode ‘R’ proteins. ‘R’ protein recognizes product of avirulance gene from the pathogen and activate downstream signaling responses leading to disease resistance. No three dimensional (3D) structural information of any ‘R’ proteins is available as yet. We have reported a ‘R’ gene homolog, the 'VMYR1′, encoding ‘R’ protein in Vigna mungo. Here, we describe the homology modeling of the 'VMYR1′ protein. The model was created by using the 3D structure of an ATP-binding cassette transporter protein from Vibrio cholerae as a template. The strategy for homology modeling was based on the high structural conservation in the superfamily of P-loop containing nucleoside triphosphate hydrolase in which target and template proteins belong. This is the first report of theoretical model structure of any ‘R’ proteins.     

The binding mechanism of a novel nicotinamide isostere inhibiting with TNKSs: a molecular dynamic simulation and binding free energy calculation

Tankyrases (TNKSs), a member of human poly (ADP-ribose) polymerase (PARP) protein superfamily, plays a key role in regulation of cell proliferation. Among the representative proteins of the PARPs family, it is found that the inhibitors have high selectivity for Tankyrase1 (TNKS1). The specific binding modes are investigated between the TNKS1 protein and nicotinamide isostere (ISX) which functions as an inhibitor of TNKS1. The stabilities of ISX-TNKS1 and AVA939-TNKS1 complexes are estimated by molecular dynamics (MD) simulations and free energy calculations; a good agreement with experimental results is reached. On the basis of the calculated results of MD simulations, we found that the inhibitors influence the conformational flexibility of TNKS1 and the XAV939 binding drive the peptide Ile1228-Gly1229-Gly1230 to form a helical structure while the ISX binding drive the peptide to form a turn structure. Moreover, the formed important hydrogen bonds of Tyr1203 residue with XVA939 and WAT1551 with ISX enhance stabilities of the complexes, and the electrostatic interactions in XAV939-TNKS1 and van der Waals interactions in ISX-TNKS1 system are main driving forces for affinity. According to the results of the decomposition of binding free energy, it is obvious that the residues Try1224 and Lys1220 make the most favorable contributions to the binding in, respectively, ISX and XAV939 complexes. Taken together, the obtained results are useful for studying the binding mechanisms of TNKSs and inhibitors and for designing potent inhibitors.     

Polymorphism of the epidermal growth factor receptor extracellular ligand binding domain: the dimer interface depends on domain stabilization

Epidermal growth factor receptors (EGFRs) and their cytoplasmic tyrosine kinases play important roles in cell proliferation and signaling. The EGFR extracellular domain (sEGFR) forms a dimer upon the binding of ligands, such as epidermal growth factor (EGF) and transforming growth factor α (TGFα). In this study, multiple molecular dynamics (MD) simulations of the 2:2 EGF·sEGFR3-512 dimer and the 2:2 TGFα·sEGFR3-512 dimer were performed in solvent and crystal environments. The simulations of systems comprising up to half a million atoms reveal part of the structural dynamics of which sEGFR dimers are capable. The solvent simulations consistently exhibited a prominent conformational relaxation from the initial crystal structures on the nanosecond time scale, leading to symmetry breaking and more extensive contacts between the two sEGFR monomers. In the crystal control simulation, this symmetry breaking and compaction was largely suppressed by crystal packing contacts. The simulations also provided evidence that the disordered domain IV of sEGFR may act as a stabilizing spacer in the dimer. Thus, the simulations suggest that the sEGFR dimer can take diverse configurations in solvent environments. These biologically relevant conformations of the EGFR signal transduction network can be controlled by contacts among the structural domains of sEGFR and its ligands.     

Carbon monoxide binding to the heme group at the dimeric interface modulates structure and copper accessibility in the Cu,Zn superoxide dismutase from Haemophilus ducreyi: in silico and in vitro evidences

X-ray absorption near-edge structure (XANES) spectroscopy and molecular dynamics (MD) simulations have been jointly applied to the study of the Cu,Zn superoxide dismutase from Haemophilus ducreyi (HdSOD) in interaction with the carbon monoxide molecule. The configurational flexibility of the Fe(II)-heme group, intercalated between the two subunits, has been sampled by MD simulations and included in the XANES data analysis without optimization in the structural parameter space. Our results provide an interpretation of the observed discrepancy in the Fe-heme distances as detected by extended X-ray absorption fine structure (EXAFS) spectroscopy and the classical XANES analysis, in which the structural parameters are optimized in a unique structure. Moreover, binding of the CO molecule to the heme induces a long range effect on the Cu,Zn active site, as evidenced by both MD simulations and in vitro experiments. MD simulation of the CO bound system, in fact, highlighted a structural rearrangement of the protein-protein hydrogen bond network in the region of the Cu,Zn active site, correlated with an increase in water accessibility at short distance from the copper atom. In line, in vitro experiments evidenced an increase of copper accessibility to a chelating agent when the CO molecule binds to the heme group, as compared to a heme deprived HdSOD. Altogether, our results support the hypothesis that the HdSOD is a heme-sensor protein, in which binding to small gaseous molecules modulates the enzyme superoxide activity as an adaptive response to the bacterial environment.     

Coupled intra- and interdomain dynamics support domain cross-talk in Pin1

The functional mechanisms of multidomain proteins often exploit interdomain interactions, or “cross-talk.” An example is human Pin1, an essential mitotic regulator consisting of a Trp–Trp (WW) domain flexibly tethered to a peptidyl-prolyl isomerase (PPIase) domain, resulting in interdomain interactions important for Pin1 function. Substrate binding to the WW domain alters its transient contacts with the PPIase domain via means that are only partially understood. Accordingly, we have investigated Pin1 interdomain interactions using NMR paramagnetic relaxation enhancement (PRE) and molecular dynamics (MD) simulations. The PREs show that apo-Pin1 samples interdomain contacts beyond the range suggested by previous structural studies. They further show that substrate binding to the WW domain simultaneously alters interdomain separation and the internal conformation of the WW domain. A 4.5-μs all-atom MD simulation of apo-Pin1 suggests that the fluctuations of interdomain distances are correlated with fluctuations of WW domain interresidue contacts involved in substrate binding. Thus, the interdomain/WW domain conformations sampled by apo-Pin1 may already include a range of conformations appropriate for binding Pin1''s numerous substrates. The proposed coupling between intra-/interdomain conformational fluctuations is a consequence of the dynamic modular architecture of Pin1. Such modular architecture is common among cell-cycle proteins; thus, the WW–PPIase domain cross-talk mechanisms of Pin1 may be relevant for their mechanisms as well.     

A comparative analysis of the equilibrium dynamics of a designed protein inferred from NMR,X‐ray,and computations

A detailed analysis of high‐resolution structural data and computationally predicted dynamics was carried out for a designed sugar‐binding protein. The mean‐square deviations in the positions of residues derived from nuclear magnetic resonance (NMR) models and those inferred from X‐ray crystallographic B‐factors for two different crystal forms were compared with the predictions based on the Gaussian Network Model (GNM) and the results from molecular dynamics (MD) simulations. GNM systematically yielded a higher correlation than MD, with experimental data, suggesting that the lack of atomistic details in the coarse‐grained GNM is more than compensated for by the mathematically exact evaluation of fluctuations using the native contacts topology. Evidence is provided that particular loop motions are curtailed by intermolecular contacts in the crystal environment causing a discrepancy between theory and experiments. Interestingly, the information conveyed by X‐ray crystallography becomes more consistent with NMR models and computational predictions when ensembles of X‐ray models are considered. Less precise (broadly distributed) ensembles indeed appear to describe the accessible conformational space under native state conditions better than B‐factors. Our results highlight the importance of using multiple conformations obtained by alternative experimental methods, and analyzing results from both coarse‐grained models and atomic simulations, for accurate assessment of motions accessible to proteins under native state conditions. Proteins 2009. © 2009 Wiley‐Liss, Inc.     

Structural insight into GRIP1-PDZ6 in Alzheimer’s disease: study from protein expression data to molecular dynamics simulations

Protein–protein interaction domain, PDZ, plays a critical role in efficient synaptic transmission in brain. Dysfunction of synaptic transmission is thought to be the underlying basis of many neuropsychiatric and neurodegenerative disorders including Alzheimer’s disease (AD). In this study, Glutamate Receptor Interacting Protein1 (GRIP1) was identified as one of the most important differentially expressed, topologically significant proteins in the protein–protein interaction network. To date, very few studies have analyzed the detailed structural basis of PDZ-mediated protein interaction of GRIP1. In order to gain better understanding of structural and dynamic basis of these interactions, we employed molecular dynamics (MD) simulations of GRIP1-PDZ6 dimer bound with Liprin-alpha and GRIP1-PDZ6 dimer alone each with 100 ns simulations. The analyses of MD simulations of Liprin-alpha bound GRIP1-PDZ6 dimer show considerable conformational differences than that of peptide-free dimer in terms of SASA, hydrogen bonding patterns, and along principal component 1 (PC1). Our study also furnishes insight into the structural attunement of the PDZ6 domains of Liprin-alpha bound GRIP1 that is attributed by significant shift of the Liprin-alpha recognition helix in the simulated peptide-bound dimer compared to the crystal structure and simulated peptide-free dimer. It is evident that PDZ6 domains of peptide-bound dimer show differential movements along PC1 than that of peptide-free dimers. Thus, Liprin-alpha also serves an important role in conferring conformational changes along the dimeric interface of the peptide-bound dimer. Results reported here provide information that may lead to novel therapeutic approaches in AD.     

Assessing and refining molecular dynamics simulations of proteins with nuclear magnetic resonance data

The sophistication of the force fields, algorithms and hardware used for molecular dynamics (MD) simulations of proteins is continuously increasing. No matter how advanced the methodology, however, it is essential to evaluate the appropriateness of the structures sampled in a simulation by comparison with quantitative experimental data. Solution nuclear magnetic resonance (NMR) data are particularly useful for checking the quality of protein simulations, as they provide both structural and dynamic information on a variety of temporal and spatial scales. Here, various features and implications of using NMR data to validate and bias MD simulations are outlined, including an overview of the different types of NMR data that report directly on structural properties and of relevant simulation techniques. The focus throughout is on how to properly account for conformational averaging, particularly within the context of the assumptions inherent in the relationships that link NMR data to structural properties.     

Functional conformational changes of endo-1,4-xylanase II from Trichoderma reesei: A molecular dynamics study

Recent crystallographic studies have revealed a range of structural changes in the three-dimensional structure of endo-1,4-xylanase (XYNII) from Trichoderma reesei. The observed conformational changes can be described as snapshots of an open-close movement of the active site of XYNII. These structures were further analyzed in this study. In addition, a total of four 1 ns molecular dynamics (MD) simulations were performed representing different states of the enzyme. A comparison of the global and local changes found in the X-ray structures and the MD runs suggested that the simulations reproduced a similar kind of active site opening and closing as predicted by the crystal structures. The open-close movement was characterized by the use of distance difference matrixes and the Hingefind program (Wriggers and Schulten, Proteins 29:1–14, 1997) to be a ‘hinge-bending’ motion involving two large rigidly-moving regions and an extended hinge. This conformational feature is probably inherent to this molecular architecture and probably plays a role in the function of XYNII. Proteins 31:434–444, 1998. © 1998 Wiley-Liss, Inc.     

The conformational feasibility for the formation of reaching dimer in ASV and HIV integrase: a molecular dynamics study

Retroviral integrases are reported to form alternate dimer assemblies like the core–core dimer and reaching dimer. The core–core dimer is stabilized predominantly by an extensive interface between two catalytic core domains. The reaching dimer is stabilized by N-terminal domains that reach to form intermolecular interfaces with the other subunit’s core and C-terminal domains (CTD), as well as CTD–CTD interactions. In this study, molecular dynamics (MD), Brownian dynamics (BD) simulations, and free energy analyses, were performed to elucidate determinants for the stability of the reaching dimer forms of full-length Avian Sarcoma Virus (ASV) and Human Immunodeficiency Virus (HIV) IN, and to examine the role of the C-tails (the last ~16–18 residues at the C-termini) in their structural dynamics. The dynamics of an HIV reaching dimer derived from small angle X-ray scattering and protein crosslinking data, was compared with the dynamics of a core–core dimer model derived from combining the crystal structures of two-domain fragments. The results showed that the core domains in the ASV reaching dimer express free dynamics, whereas those in the HIV reaching dimer are highly stable. BD simulations suggest a higher rate of association for the HIV core–core dimer than the reaching dimer. The predicted stability of these dimers was therefore ranked in the following order: ASV reaching dimer < HIV reaching dimer < composite core–core dimer. Analyses of MD trajectories have suggested residues that are critical for intermolecular contacts in each reaching dimer. Tests of these predictions and insights gained from these analyses could reveal a potential pathway for the association and dissociation of full-length IN multimers.     

Molecular Dynamics Study on Protein and it's Water Structure at High Pressure

Abstract

We have performed NPT molecular dynamics simulations (Langevin Piston Method) on two types of solvated proteins-‘denaturation-unfavorable’ protein (insulin) and ‘denaturation-favorable protein’ (ribonuclease A) at high pressure (from 1 bar up to 20 kbar). The method is based on the extended system formalism introduced by Andersen, where the deterministic equations of motion for the piston degree of freedom are replaced by Langevin equation. We report the structural changes of proteins (ribonuclease A and insulin) and water molecules through radius of gyration, solvent accessible surface area, hydrogen bond pattern, and the topology of water clusters connected by the hydrogen bonded circular network. The solvent accessibility of ribonuclease A is mainly decreased by hydrophilic residues rather than hydrophobic residues under high pressure. From the results of hydrogen bond analysis, we have found that α-helix is more stable than β-sheet under high pressure. In addition, from the analysis of the water cluster, we have observed that for ribonuclease A, 5-membered ring structure is more favorable than 6-membered ring at higher pressure. However, for insulin, the ratio of 5 to 6-ring is constant over the pressure ranges for which we have performed MD simulation. This indicates that the water structure around insulin does not change under high pressure.     

细胞膜蛋白与细胞骨架蛋白相互作用研究进展

细胞膜蛋白与胞浆骨架蛋白的相互作用对于维持细胞正常形态 ,细胞粘附与信号传导有重要作用。含有 4 .1 JEF结构域的蛋白 4 .1超家族与含有PDZ结构域的MAGUK蛋白家族能结合多种膜蛋白胞内区与胞浆蛋白 ,在膜蛋白与胞浆蛋白之间建立联系 ,对于细胞、细胞 -细胞间连接的正常结构与功能的维持有着重要作用。     

Evolutionarily conserved regions and hydrophobic contacts at the superfamily level: The case of the fold-type I, pyridoxal-5'-phosphate-dependent enzymes

The wealth of biological information provided by structural and genomic projects opens new prospects of understanding life and evolution at the molecular level. In this work, it is shown how computational approaches can be exploited to pinpoint protein structural features that remain invariant upon long evolutionary periods in the fold-type I, PLP-dependent enzymes. A nonredundant set of 23 superposed crystallographic structures belonging to this superfamily was built. Members of this family typically display high-structural conservation despite low-sequence identity. For each structure, a multiple-sequence alignment of orthologous sequences was obtained, and the 23 alignments were merged using the structural information to obtain a comprehensive multiple alignment of 921 sequences of fold-type I enzymes. The structurally conserved regions (SCRs), the evolutionarily conserved residues, and the conserved hydrophobic contacts (CHCs) were extracted from this data set, using both sequence and structural information. The results of this study identified a structural pattern of hydrophobic contacts shared by all of the superfamily members of fold-type I enzymes and involved in native interactions. This profile highlights the presence of a nucleus for this fold, in which residues participating in the most conserved native interactions exhibit preferential evolutionary conservation, that correlates significantly (r = 0.70) with the extent of mean hydrophobic contact value of their apolar fraction.     

Subtle functional collective motions in pancreatic-like ribonucleases: from ribonuclease A to angiogenin

The analysis of the dynamic behavior of enzymes is fundamental to structural biology. A direct relationship between protein flexibility and biological function has been shown for bovine pancreatic ribonuclease (RNase A) (Rasmussen et al., Nature 1992;357:423-424). More recently, crystallographic studies have shown that functional motions in RNase A involve the enzyme beta-sheet regions that move concertedly on substrate binding and release (Vitagliano et al., Proteins 2002;46:97-104). These motions have been shown to correspond to intrinsic dynamic properties of the native enzyme by molecular dynamics (MD) simulations. To unveil the occurrence of these collective motions in other members of pancreatic-like superfamily, we carried out MD simulations on human angiogenin (Ang). Essential dynamics (ED) analyses performed on the trajectories reveal that Ang exhibits collective motions similar to RNase A, despite the limited sequence identity (33%) of the two proteins. Furthermore, we show that these collective motions are also present in ensembles of experimentally determined structures of both Ang and RNase A. Finally, these subtle concerted beta-sheet motions were also observed for other two members of the pancreatic-like superfamily by comparing the ligand-bound and ligand-free structures of these enzymes. Taken together, these findings suggest that pancreatic-like ribonucleases share an evolutionary conserved dynamic behavior consisting of subtle beta-sheet motions, which are essential for substrate binding and release.