Analysis of Molecular Dynamics Simulations of Protein Folding

Microsecond long Molecular Dynamics (MD) trajectories of biomolecular processes are now possible due to advances in computer technology. Soon, trajectories long enough to probe dynamics over many milliseconds will become available. Since these timescales match the physiological timescales over which many small proteins fold, all atom MD simulations of protein folding are now becoming popular. To distill features of such large folding trajectories, we must develop methods that can both compress trajectory data to enable visualization, and that can yield themselves to further analysis, such as the finding of collective coordinates and reduction of the dynamics. Conventionally, clustering has been the most popular MD trajectory analysis technique, followed by principal component analysis (PCA). Simple clustering used in MD trajectory analysis suffers from various serious drawbacks, namely, (i) it is not data driven, (ii) it is unstable to noise and change in cutoff parameters, and (iii) since it does not take into account interrelationships amongst data points, the separation of data into clusters can often be artificial. Usually, partitions generated by clustering techniques are validated visually, but such validation is not possible for MD trajectories of protein folding, as the underlying structural transitions are not well understood. Rigorous cluster validation techniques may be adapted, but it is more crucial to reduce the dimensions in which MD trajectories reside, while still preserving their salient features. PCA has often been used for dimension reduction and while it is computationally inexpensive, being a linear method, it does not achieve good data compression. In this thesis, I propose a different method, a nonmetric multidimensional scaling (nMDS) technique, which achieves superior data compression by virtue of being nonlinear, and also provides a clear insight into the structural processes underlying MD trajectories. I illustrate the capabilities of nMDS by analyzing three complete villin headpiece folding and six norleucine mutant (NLE) folding trajectories simulated by Freddolino and Schulten [1]. Using these trajectories, I make comparisons between nMDS, PCA and clustering to demonstrate the superiority of nMDS. The three villin headpiece trajectories showed great structural heterogeneity. Apart from a few trivial features like early formation of secondary structure, no commonalities between trajectories were found. There were no units of residues or atoms found moving in concert across the trajectories. A flipping transition, corresponding to the flipping of helix 1 relative to the plane formed by helices 2 and 3 was observed towards the end of the folding process in all trajectories, when nearly all native contacts had been formed. However, the transition occurred through a different series of steps in all trajectories, indicating that it may not be a common transition in villin folding. The trajectories showed competition between local structure formation/hydrophobic collapse and global structure formation in all trajectories. Our analysis on the NLE trajectories confirms the notion that a tight hydrophobic core inhibits correct 3-D rearrangement. Only one of the six NLE trajectories folded, and it showed no flipping transition. All the other trajectories get trapped in hydrophobically collapsed states. The NLE residues were found to be buried deeply into the core, compared to the corresponding lysines in the villin headpiece, thereby making the core tighter and harder to undo for 3-D rearrangement. Our results suggest that the NLE may not be a fast folder as experiments suggest. The tightness of the hydrophobic core may be a very important factor in the folding of larger proteins. It is likely that chaperones like GroEL act to undo the tight hydrophobic core of proteins, after most secondary structure elements have been formed, so that global rearrangement is easier. I conclude by presenting facts about chaperone-protein complexes and propose further directions for the study of protein folding.

Protein encapsulation and release from PEO-b-polyphosphoester templated calcium carbonate particles

International audience

Biochemical studies of postsynaptic density signaling proteins with a focus on synaptic GTPase activating protein and PDZ domains

Memory storage in the brain involves adjustment of the strength of existing synapses and formation of new neural networks. A key process underlying memory formation is synaptic plasticity, the ability of excitatory synapses to strengthen or weaken their connections in response to patterns of activity between their connected neurons. Synaptic plasticity is governed by the precise pattern of Ca²⁺ influx through postsynaptic N-methyl-D-aspartate-type glutamate receptors (NMDARs), which can lead to the activation of the small GTPases Ras and Rap. Differential activation of Ras and Rap acts to modulate synaptic strength by promoting the insertion or removal of 2-amino-3-(3-hydroxy-5-methyl-isoxazol-4-yl)propanoic acid receptors (AMPARs) from the synapse. Synaptic GTPase activating protein (synGAP) regulates AMPAR levels by catalyzing the inactivation of GTP-bound (active) Ras or Rap. synGAP is positioned in close proximity to the cytoplasmic tail regions of the NMDAR through its association with the PDZ domains of PSD-95. SynGAP’s activity is regulated by the prominent postsynaptic protein kinase, Ca²⁺/calmodulin-dependent protein kinase II (CaMKII) and cyclin-dependent kinase 5 (CDK5), a known binding partner of CaMKII. Modulation of synGAP’s activity by phosphorylation may alter the ratio of active Ras to Rap in spines, thus pushing the spine towards the insertion or removal of AMPARs, subsequently strengthening or weakening the synapse. To date, all biochemical studies of the regulation of synGAP activity by protein kinases have utilized impure preparations of membrane bound synGAP. Here we have clarified the effects of phosphorylation of synGAP on its Ras and Rap GAP activities by preparing and utilizing purified, soluble recombinant synGAP, Ras, Rap, CaMKII, CDK5, PLK2, and CaM. Using mass spectrometry, we have confirmed the presence of previously identified CaMKII and CDK5 sites in synGAP, and have identified novel sites of phosphorylation by CaMKII, CDK5, and PLK2. We have shown that the net effect of phosphorylation of synGAP by CaMKII, CDK5, and PLK2 is an increase in its GAP activity toward HRas and Rap1. In contrast, there is no effect on its GAP activity toward Rap2. Additionally, by assaying the GAP activity of phosphomimetic synGAP mutants, we have been able to hypothesize the effects of CDK5 phosphorylation at specific sites in synGAP. In the course of this work, we also found, unexpectedly, that synGAP is itself a Ca²⁺/CaM binding protein. While Ca²⁺/CaM binding does not directly affect synGAP activity, it causes a conformational change in synGAP that increases the rate of its phosphorylation and exposes additional phosphorylation sites that are inaccessible in the absence of Ca²⁺/CaM.

The postsynaptic density (PSD) is an electron-dense region in excitatory postsynaptic neurons that contains a high concentration of glutamate receptors, cytoskeletal proteins, and associated signaling enzymes. Within the PSD, three major classes of scaffolding molecules function to organize signaling enzymes and glutamate receptors. PDZ domains present in the Shank and PSD-95 scaffolds families serve to physically link AMPARs and NMDARs to signaling molecules in the PSD. Because of the specificity and high affinity of PDZ domains for their ligands, I reasoned that these interacting pairs could provide the core components of an affinity chromatography system, including affinity resins, affinity tags, and elution agents. I show that affinity columns containing the PDZ domains of PSD-95 can be used to purify active PDZ domain-binding proteins to very high purity in a single step. Five heterologously expressed neuronal proteins containing endogenous PDZ domain ligands (NMDAR GluN2B subunit Tail, synGAP, neuronal nitric oxide synthase PDZ domain, cysteine rich interactor of PDZ three and cypin) were purified using PDZ domain resin, with synthetic peptides having the sequences of cognate PDZ domain ligands used as elution agents. I also show that conjugation of PDZ domain-related affinity tags to Proteins Of Interest (POIs) that do not contain endogenous PDZ domains or ligands does not alter protein activity and enables purification of the POIs on PDZ domain-related affinity resins.

Taurine Supplementation Increases K(atp) Channel Protein Content, Improving Ca2+ Handling And Insulin Secretion In Islets From Malnourished Mice Fed On A High-fat Diet.

Pancreatic β-cells are highly sensitive to suboptimal or excess nutrients, as occurs in protein-malnutrition and obesity. Taurine (Tau) improves insulin secretion in response to nutrients and depolarizing agents. Here, we assessed the expression and function of Cav and KATP channels in islets from malnourished mice fed on a high-fat diet (HFD) and supplemented with Tau. Weaned mice received a normal (C) or a low-protein diet (R) for 6 weeks. Half of each group were fed a HFD for 8 weeks without (CH, RH) or with 5% Tau since weaning (CHT, RHT). Isolated islets from R mice showed lower insulin release with glucose and depolarizing stimuli. In CH islets, insulin secretion was increased and this was associated with enhanced KATP inhibition and Cav activity. RH islets secreted less insulin at high K(+) concentration and showed enhanced KATP activity. Tau supplementation normalized K(+)-induced secretion and enhanced glucose-induced Ca(2+) influx in RHT islets. R islets presented lower Ca(2+) influx in response to tolbutamide, and higher protein content and activity of the Kir6.2 subunit of the KATP. Tau increased the protein content of the α1.2 subunit of the Cav channels and the SNARE proteins SNAP-25 and Synt-1 in CHT islets, whereas in RHT, Kir6.2 and Synt-1 proteins were increased. In conclusion, impaired islet function in R islets is related to higher content and activity of the KATP channels. Tau treatment enhanced RHT islet secretory capacity by improving the protein expression and inhibition of the KATP channels and enhancing Synt-1 islet content.

A Putative Otu Domain-containing Protein 1 Deubiquitinating Enzyme Is Differentially Expressed In Thyroid Cancer And Identifies Less-aggressive Tumours.

This study aimed to identify novel biomarkers for thyroid carcinoma diagnosis and prognosis. We have constructed a human single-chain variable fragment (scFv) antibody library that was selected against tumour thyroid cells using the BRASIL method (biopanning and rapid analysis of selective interactive ligands) and phage display technology. One highly reactive clone, scFv-C1, with specific binding to papillary thyroid tumour proteins was confirmed by ELISA, which was further tested against a tissue microarray that comprised of 229 thyroid tissues, including: 110 carcinomas (38 papillary thyroid carcinomas (PTCs), 42 follicular carcinomas, 30 follicular variants of PTC), 18 normal thyroid tissues, 49 nodular goitres (NG) and 52 follicular adenomas. The scFv-C1 was able to distinguish carcinomas from benign lesions (P=0.0001) and reacted preferentially against T1 and T2 tumour stages (P=0.0108). We have further identified an OTU domain-containing protein 1, DUBA-7 deubiquitinating enzyme as the scFv-binding antigen using two-dimensional polyacrylamide gel electrophoresis and mass spectrometry. The strategy of screening and identifying a cell-surface-binding antibody against thyroid tissues was highly effective and resulted in a useful biomarker that recognises malignancy among thyroid nodules and may help identify lower-risk cases that can benefit from less-aggressive management.

Core-valence Correlation Effects On Ir Calculations: The Bf3 And Bcl3 Cases.

The first theoretical results of core-valence correlation effects are presented for the infrared wavenumbers and intensities of the BF3 and BCl3 molecules, using (double- and triple-zeta) Dunning core-valence basis sets at the CCSD(T) level. The results are compared with those calculated in the frozen core approximation with standard Dunning basis sets at the same correlation level and with the experimental values. The general conclusion is that the effect of core-valence correlation is, for infrared wavenumbers and intensities, smaller than the effect of adding augmented diffuse functions to the basis set, e.g., cc-pVTZ to aug-cc-pVTZ. Moreover, the trends observed in the data are mainly related to the augmented functions rather than the core-valence functions added to the basis set. The results obtained here confirm previous studies pointing out the large descrepancy between the theoretical and experimental intensities of the stretching mode for BCl3.

Iis--integrated Interactome System: A Web-based Platform For The Annotation, Analysis And Visualization Of Protein-metabolite-gene-drug Interactions By Integrating A Variety Of Data Sources And Tools.

High-throughput screening of physical, genetic and chemical-genetic interactions brings important perspectives in the Systems Biology field, as the analysis of these interactions provides new insights into protein/gene function, cellular metabolic variations and the validation of therapeutic targets and drug design. However, such analysis depends on a pipeline connecting different tools that can automatically integrate data from diverse sources and result in a more comprehensive dataset that can be properly interpreted. We describe here the Integrated Interactome System (IIS), an integrative platform with a web-based interface for the annotation, analysis and visualization of the interaction profiles of proteins/genes, metabolites and drugs of interest. IIS works in four connected modules: (i) Submission module, which receives raw data derived from Sanger sequencing (e.g. two-hybrid system); (ii) Search module, which enables the user to search for the processed reads to be assembled into contigs/singlets, or for lists of proteins/genes, metabolites and drugs of interest, and add them to the project; (iii) Annotation module, which assigns annotations from several databases for the contigs/singlets or lists of proteins/genes, generating tables with automatic annotation that can be manually curated; and (iv) Interactome module, which maps the contigs/singlets or the uploaded lists to entries in our integrated database, building networks that gather novel identified interactions, protein and metabolite expression/concentration levels, subcellular localization and computed topological metrics, GO biological processes and KEGG pathways enrichment. This module generates a XGMML file that can be imported into Cytoscape or be visualized directly on the web. We have developed IIS by the integration of diverse databases following the need of appropriate tools for a systematic analysis of physical, genetic and chemical-genetic interactions. IIS was validated with yeast two-hybrid, proteomics and metabolomics datasets, but it is also extendable to other datasets. IIS is freely available online at: http://www.lge.ibi.unicamp.br/lnbio/IIS/.

The Effect Of Celastrol, A Triterpene With Antitumorigenic Activity, On Conformational And Functional Aspects Of The Human 90kda Heat Shock Protein Hsp90α, A Chaperone Implicated In The Stabilization Of The Tumor Phenotype.

Hsp90 is a molecular chaperone essential for cell viability in eukaryotes that is associated with the maturation of proteins involved in important cell functions and implicated in the stabilization of the tumor phenotype of various cancers, making this chaperone a notably interesting therapeutic target. Celastrol is a plant-derived pentacyclic triterpenoid compound with potent antioxidant, anti-inflammatory and anticancer activities; however, celastrol's action mode is still elusive. In this work, we investigated the effect of celastrol on the conformational and functional aspects of Hsp90α. Interestingly, celastrol appeared to target Hsp90α directly as the compound induced the oligomerization of the chaperone via the C-terminal domain as demonstrated by experiments using a deletion mutant. The nature of the oligomers was investigated by biophysical tools demonstrating that a two-fold excess of celastrol induced the formation of a decameric Hsp90α bound throughout the C-terminal domain. When bound, celastrol destabilized the C-terminal domain. Surprisingly, standard chaperone functional investigations demonstrated that neither the in vitro chaperone activity of protecting against aggregation nor the ability to bind a TPR co-chaperone, which binds to the C-terminus of Hsp90α, were affected by celastrol. Celastrol interferes with specific biological functions of Hsp90α. Our results suggest a model in which celastrol binds directly to the C-terminal domain of Hsp90α causing oligomerization. However, the ability to protect against protein aggregation (supported by our results) and to bind to TPR co-chaperones are not affected by celastrol. Therefore celastrol may act primarily by inducing specific oligomerization that affects some, but not all, of the functions of Hsp90α. To the best of our knowledge, this study is the first work to use multiple probes to investigate the effect that celastrol has on the stability and oligomerization of Hsp90α and on the binding of this chaperone to Tom70. This work provides a novel mechanism by which celastrol binds Hsp90α.

Treadmill Training Increases Sirt-1 And Pgc-1 α Protein Levels And Ampk Phosphorylation In Quadriceps Of Middle-aged Rats In An Intensity-dependent Manner.

The present study investigated the effects of running at 0.8 or 1.2 km/h on inflammatory proteins (i.e., protein levels of TNF- α , IL-1 β , and NF- κ B) and metabolic proteins (i.e., protein levels of SIRT-1 and PGC-1 α , and AMPK phosphorylation) in quadriceps of rats. Male Wistar rats at 3 (young) and 18 months (middle-aged rats) of age were divided into nonexercised (NE) and exercised at 0.8 or 1.2 km/h. The rats were trained on treadmill, 50 min per day, 5 days per week, during 8 weeks. Forty-eight hours after the last training session, muscles were removed, homogenized, and analyzed using biochemical and western blot techniques. Our results showed that: (a) running at 0.8 km/h decreased the inflammatory proteins and increased the metabolic proteins compared with NE rats; (b) these responses were lower for the inflammatory proteins and higher for the metabolic proteins in young rats compared with middle-aged rats; (c) running at 1.2 km/h decreased the inflammatory proteins and increased the metabolic proteins compared with 0.8 km/h; (d) these responses were similar between young and middle-aged rats when trained at 1.2 km. In summary, the age-related increases in inflammatory proteins, and the age-related declines in metabolic proteins can be reversed and largely improved by treadmill training.

Reduced Insulin Clearance And Lower Insulin-degrading Enzyme Expression In The Liver Might Contribute To The Thrifty Phenotype Of Protein-restricted Mice.

Nutrient restriction during the early stages of life usually leads to alterations in glucose homeostasis, mainly insulin secretion and sensitivity, increasing the risk of metabolic disorders in adulthood. Despite growing evidence regarding the importance of insulin clearance during glucose homeostasis in health and disease, no information exists about this process in malnourished animals. Thus, in the present study, we aimed to determine the effect of a nutrient-restricted diet on insulin clearance using a model in which 30-d-old C57BL/6 mice were exposed to a protein-restricted diet for 14 weeks. After this period, we evaluated many metabolic variables and extracted pancreatic islet, liver, gastrocnemius muscle (GCK) and white adipose tissue samples from the control (normal-protein diet) and restricted (low-protein diet, LP) mice. Insulin concentrations were determined using RIA and protein expression and phosphorylation by Western blot analysis. The LP mice exhibited lower body weight, glycaemia, and insulinaemia, increased glucose tolerance and altered insulin dynamics after the glucose challenge. The improved glucose tolerance could partially be explained by an increase in insulin sensitivity through the phosphorylation of the insulin receptor/protein kinase B and AMP-activated protein kinase/acetyl-CoA carboxylase in the liver, whereas the changes in insulin dynamics could be attributed to reduced insulin secretion coupled with reduced insulin clearance and lower insulin-degrading enzyme (IDE) expression in the liver and GCK. In summary, protein-restricted mice not only produce and secrete less insulin, but also remove and degrade less insulin. This phenomenon has the double benefit of sparing insulin while prolonging and potentiating its effects, probably due to the lower expression of IDE in the liver, possibly with long-term consequences.

Islet Neogenesis-associated Protein (ingap): The Role Of Its Endogenous Production As A Positive Modulator Of Insulin Secretion.

Islet neogenesis-associated protein (INGAP) is a peptide found in pancreatic exocrine-, duct- and islet- non-β-cells from normal hamsters. Its increase induced by either its exogenous administration or by the overexpression of its gene enhances β-cell secretory function and increases β-cell mass by a combination of stimulation of cell replication and islet neogenesis and reduction of β-cell apoptosis. We studied the potential modulatory role of endogenous INGAP in insulin secretion using two different experimental approaches. Hamster islets transfected with INGAP-small interfering RNA (INGAP-siRNA) were used to study glucose-stimulated insulin secretion (GSIS). In parallel, freshly isolated islets were incubated with high glucose and the same concentration of either a specific anti-INGAP rabbit serum or normal rabbit serum. INGAP-siRNA transfected islets reduced their INGAP mRNA and protein content by 35.1% and 47.2%, respectively whereas GSIS decreased by 25.8%. GSIS by transfected islets attained levels comparable to those recorded in control islets when INGAP pentadecapeptide (INGAP-PP) was added to the culture medium. INGAP antibody in the medium decreased significantly GSIS in a dose-dependent manner. These results indicate that endogenous INGAP plays a physiological positive modulatory role in insulin secretion, supporting its possible use in the treatment of prediabetes and Type 2 diabetes.

Taurine Supplementation Reduces Blood Pressure And Prevents Endothelial Dysfunction And Oxidative Stress In Post-weaning Protein-restricted Rats.

Taurine is a sulfur-containing amino acid that exerts protective effects on vascular function and structure in several models of cardiovascular diseases through its antioxidant and anti-inflammatory properties. Early protein malnutrition reprograms the cardiovascular system and is linked to hypertension in adulthood. This study assessed the effects of taurine supplementation in vascular alterations induced by protein restriction in post-weaning rats. Weaned male Wistar rats were fed normal- (12%, NP) or low-protein (6%, LP) diets for 90 days. Half of the NP and LP rats concomitantly received 2.5% taurine supplementation in the drinking water (NPT and LPT, respectively). LP rats showed elevated systolic, diastolic and mean arterial blood pressure versus NP rats; taurine supplementation partially prevented this increase. There was a reduced relaxation response to acetylcholine in isolated thoracic aortic rings from the LP group that was reversed by superoxide dismutase (SOD) or apocynin incubation. Protein expression of p47phox NADPH oxidase subunit was enhanced, whereas extracellular (EC)-SOD and endothelial nitric oxide synthase phosphorylation at Ser 1177 (p-eNOS) were reduced in aortas from LP rats. Furthermore, ROS production was enhanced while acetylcholine-induced NO release was reduced in aortas from the LP group. Taurine supplementation improved the relaxation response to acetylcholine and eNOS-derived NO production, increased EC-SOD and p-eNOS protein expression, as well as reduced ROS generation and p47phox expression in the aortas from LPT rats. LP rats showed an increased aortic wall/lumen ratio and taurine prevented this remodeling through a reduction in wall media thickness. Our data indicate a protective role of taurine supplementation on the high blood pressure, endothelial dysfunction and vascular remodeling induced by post-weaning protein restriction. The beneficial vascular effect of taurine was associated with restoration of vascular redox homeostasis and improvement of NO bioavailability.

Histologic Correlation Of Expression Of Ki-67 In Squamous Cell Carcinoma Of The Glottis According To The Degree Of Cell Differentiation.

Squamous cell carcinoma is the most common neoplasm of the larynx and glottis, and its prognosis depends on the size of the lesion, level of local invasion, cervical lymphatic spread, and presence of distant metastases. Ki-67 (MKI67) is a protein present in the core, whose function is related to cell proliferation. To evaluate the expression of marker Ki-67 in squamous cell carcinoma of the larynx and glottis and its correlation to pathological findings. Experimental study with immunohistochemistry analysis of Ki-67, calculating the percentage of the cell proliferation index in glottic squamous cell carcinomas. Sixteen cases were analyzed, with six well-differentiated and 10 poorly/moderately differentiated tumors. There was a correlation between cell proliferation index and degree of cell differentiation, with higher proliferation in poorly/moderately differentiated tumors. The cell proliferation index, as measured by Ki-67, may be useful in the characterization of histological degree in glottic squamous cell tumors.