Molecular dynamics simulation of the phosphorylation-induced conformational changes of a tau peptide fragment

Aggregation of the microtubule associated protein tau (MAPT) within neurons of the brain is the leading cause of tauopathies such as Alzheimer's disease. MAPT is a phospho-protein that is selectively phosphorylated by a number of kinases in vivo to perform its biological function. However, it may become pathogenically hyperphosphorylated, causing aggregation into paired helical filaments and neurofibrillary tangles. The phosphorylation induced conformational change on a peptide of MAPT (htau225−250) was investigated by performing molecular dynamics simulations with different phosphorylation patterns of the peptide (pThr231 and/or pSer235) in different simulation conditions to determine the effect of ionic strength and phosphate charge. All phosphorylation patterns were found to disrupt a nascent terminal β-sheet pattern (226VAVVR230 and 244QTAPVP249), replacing it with a range of structures. The double pThr231/pSer235 phosphorylation pattern at experimental ionic strength resulted in the best agreement with NMR structural characterization, with the observation of a transient α-helix (239AKSRLQT245). PPII helical conformations were only found sporadically throughout the simulations. Proteins 2014; 82:1907–1923. © 2014 Wiley Periodicals, Inc.

A forensic investigation of single human hair fibres using FTIR-ATR spectroscopy and chemometrics

Human hair fibres are ubiquitous in nature and are found frequently at crime scenes often as a result of exchange between the perpetrator, victim and/or the surroundings according to Locard's Principle. Therefore, hair fibre evidence can provide important information for crime investigation. For human hair evidence, the current forensic methods of analysis rely on comparisons of either hair morphology by microscopic examination or nuclear and mitochondrial DNA analyses. Unfortunately in some instances the utilisation of microscopy and DNA analyses are difficult and often not feasible. This dissertation is arguably the first comprehensive investigation aimed to compare, classify and identify the single human scalp hair fibres with the aid of FTIR-ATR spectroscopy in a forensic context. Spectra were collected from the hair of 66 subjects of Asian, Caucasian and African (i.e. African-type). The fibres ranged from untreated to variously mildly and heavily cosmetically treated hairs. The collected spectra reflected the physical and chemical nature of a hair from the near-surface particularly, the cuticle layer. In total, 550 spectra were acquired and processed to construct a relatively large database. To assist with the interpretation of the complex spectra from various types of human hair, Derivative Spectroscopy and Chemometric methods such as Principal Component Analysis (PCA), Fuzzy Clustering (FC) and Multi-Criteria Decision Making (MCDM) program; Preference Ranking Organisation Method for Enrichment Evaluation (PROMETHEE) and Geometrical Analysis for Interactive Aid (GAIA); were utilised. FTIR-ATR spectroscopy had two important advantages over to previous methods: (i) sample throughput and spectral collection were significantly improved (no physical flattening or microscope manipulations), and (ii) given the recent advances in FTIR-ATR instrument portability, there is real potential to transfer this work.s findings seamlessly to on-field applications. The "raw" spectra, spectral subtractions and second derivative spectra were compared to demonstrate the subtle differences in human hair. SEM images were used as corroborative evidence to demonstrate the surface topography of hair. It indicated that the condition of the cuticle surface could be of three types: untreated, mildly treated and treated hair. Extensive studies of potential spectral band regions responsible for matching and discrimination of various types of hair samples suggested the 1690-1500 cm^-1 IR spectral region was to be preferred in comparison with the commonly used 1750-800 cm^-1. The principal reason was the presence of the highly variable spectral profiles of cystine oxidation products (1200-1000 cm^-1), which contributed significantly to spectral scatter and hence, poor hair sample matching. In the preferred 1690-1500 cm^-1 region, conformational changes in the keratin protein attributed to the α-helical to β-sheet transitions in the Amide I and Amide II vibrations and played a significant role in matching and discrimination of the spectra and hence, the hair fibre samples. For gender comparison, the Amide II band is significant for differentiation. The results illustrated that the male hair spectra exhibit a more intense β-sheet vibration in the Amide II band at approximately 1511 cm^-1 whilst the female hair spectra displayed more intense α-helical vibration at 1520-1515cm^-1. In terms of chemical composition, female hair spectra exhibit greater intensity of the amino acid tryptophan (1554 cm^-1), aspartic and glutamic acid (1577 cm^-1). It was also observed that for the separation of samples based on racial differences, untreated Caucasian hair was discriminated from Asian hair as a result of having higher levels of the amino acid cystine and cysteic acid. However, when mildly or chemically treated, Asian and Caucasian hair fibres are similar, whereas African-type hair fibres are different. In terms of the investigation's novel contribution to the field of forensic science, it has allowed for the development of a novel, multifaceted, methodical protocol where previously none had existed. The protocol is a systematic method to rapidly investigate unknown or questioned single human hair FTIR-ATR spectra from different genders and racial origin, including fibres of different cosmetic treatments. Unknown or questioned spectra are first separated on the basis of chemical treatment i.e. untreated, mildly treated or chemically treated, genders, and racial origin i.e. Asian, Caucasian and African-type. The methodology has the potential to complement the current forensic analysis methods of fibre evidence (i.e. Microscopy and DNA), providing information on the morphological, genetic and structural levels.

Effect of external high pressures on the clay mineral sodium montmorillonite intercalated with methylated octadecylammonium bromide surfactants

Raman microprobe spectra of the clay mineral Wyoming SWy-2-sodium montmorillonite intercalated with the surfactants, methyltrioctadecylammonium bromide (TOMA) dimethyldiotadecylammonium bromide (DODMA) and octadecyl-trimethylammonium bromide (ODTMA), have been measured in the CH2 stretching region at external pressures up to ~40 kbar with the aid of a diamond-anvil cell. In the case of the intercalated clays containing TOMA and DODMA, the Raman data afford evidence for gauche to trans conformational changes in the orientation of the CH2 chains in the surfactants with increasing pressure. These conformational changes are reversed completely upon the release of pressure.

The effect of various pH, ionic strength and temperature on papain hydrolysis of salivary film

Stimulated human whole saliva (WS) was used to study the dynamics of papain hydrolysis at defined pH, ionic strength and temperature with the view of reducing an acquired pellicle. A quartz crystal microbalance with dissipation (QCM-D) was used to monitor the changes in frequency due to enzyme hydrolysis of WS films and the hydrolytic parameters were calculated using an empirical model. The morphological and conformational changes of the salivary films before and after enzymatic hydrolysis were characterized by atomic force microscopy (AFM) imaging and grazing angle infrared spectroscopy (GA-FTIR) spectra, respectively. The characteristics of papain hydrolysis of WS films were pH-, ionic strength- and temperature-dependent. The WS films were partially removed by the action of enzyme, resulting thinner and smoother surfaces. The IR data suggested that hydrolysis-induced deformation did not occur onto the remnants salivary films. The processes of papain hydrolysis of WS films can be controlled by properly regulating pH, ionic strength and temperature.

Catalytic properties of 26 S and 20 S proteasomes and radiolabeling of MB1, LMP7, and C7 subunits associated with trypsin-like and chymotrypsin-like activities*

20 and 26 S proteasomes were isolated from rat liver. The procedure developed for the 26 S proteasome resulted in greatly improved yields compared with previously published methods. A comparison of the kinetic properties of 20 and 26 S proteasomes showed significant differences in the kinetic characteristics with certain substrates and differences in the effects of a protein substrate on peptidase activity. Observed differences in the kinetics of peptidylglutamyl peptide hydrolase activity suggest that the 26 S complex cannot undergo the conformational changes of 20 S proteasomes at high concentrations of the substrate benzyloxycarbonyl (Z) -Leu-Leu-Glu-b-naphthylamide. Various inhibitors that differentially affect the trypsin-like and chymotrypsin-like activities have been identified. Ala-Ala-Phe-chloromethyl (CH2Cl) inhibits chymotrypsin-like activity assayed with succinyl (Suc) -Leu-Leu-Val-Tyr-AMC, but surprisingly not hydrolysis of Ala-Ala-Phe-7-amido-4-methylcoumarin (AMC). Tyr-Gly-Arg-CH2Cl inhibits Suc-Leu-Leu-Val-Tyr-AMC hydrolysis as well as trypsinlike activity measured with t-butoxycarbonyl (Boc) -Leu-Ser-Thr-Arg-AMC, while Z-Phe-Gly-Tyr-diazomethyl (CHN2) was found to inhibit only the two chymotrypsin- like activities. Radiolabeled forms of peptidyl chloromethane and peptidyl diazomethane inhibitors, [3H]acetyl-Ala-Ala-Phe-CH2Cl, [3H]acetyland radioiodinated Tyr-Gly-Arg-CH2Cl, and Z-Phe-Gly- Tyr-(125I-CHN2), have been used to identify catalytic components associated with each of the three peptidase activities. In each case, incorporation of the label could be blocked by prior treatment of the proteasomes with known active site-directed inhibitors, calpain inhibitor 1 or 3,4-dichloroisocoumarin. Subunits of labeled proteasomes were separated either by reverse phase-HPLC and SDS-polyacrylamide gel electrophoresis or by twodimensional polyacrylamide gel electrophoresis followed by autoradiography/fluorography and immunoblotting with subunit-specific antibodies. In each case, label was found to be incorporated into subunits C7, MB1, and LMP7 but in different relative amounts depending on the inhibitor used, consistent with the observed effects on the different peptidase activities. The results strongly suggest a relationship between trypsin-like activity and chymotrypsin-like activity. They also help to relate the different subunits of the complex to the assayed multicatalytic endopeptidase activities

Catalytic components of proteasomes and the regulation of proteinase activity

The proteasome (multicatalytic proteinase complex) is a large multimeric complex which is found in the nucleus and cytoplasm of eukaryotic cells. It plays a major role in both ubiquitin-dependent and ubiquitin-independent nonlysosomal pathways of protein degradation. Proteasome subunits are encoded by members of the same gene family and can be divided into two groups based on their similarity to the c~ and /3 subunits of the simpler proteasome isolated from Thermoplasma acidophilum. Proteasomes have a cylindrical structure composed of four rings of seven subunits. The 26S form of the proteasome, which is responsible for ubiquitin-dependent proteolysis, contains additional regulatory complexes. Eukaryotic proteasomes have multiple catalytic activities which are catalysed at distinct sites. Since proteasomes are unrelated to other known proteases, there are no clues as to which are the catalytic components from sequence alignments. It has been assumed from studies with yeast mutants that /3-type subunits play a catalytic role. Using a radiolabelled peptidyl chloromethane inhibitor of rat liver proteasomes we have directly identified RC7 as a catalytic component. Interestingly, mutants in Prel, the yeast homologue of RC7, have already been reported to have defective chymotrypsin-like activity. These results taken together confirm a direct catalytic role for these/3-type subunits. Proteasome activities are sensitive to conformational changes and there are several ways in which proteasome function may be modulated in vivo. Our recent studies have shown that in animal cells at least two proteasome subunits can undergo phosphorylation, the level of which is likely to be important for determining proteasome localization, activity or ability to form larger complexes. In addition, we have isolated two isoforms of the 26S proteinase.

Modeling ion channel dynamics through reflected stochastic differential equations

Ion channels are membrane proteins that open and close at random and play a vital role in the electrical dynamics of excitable cells. The stochastic nature of the conformational changes these proteins undergo can be significant, however current stochastic modeling methodologies limit the ability to study such systems. Discrete-state Markov chain models are seen as the "gold standard," but are computationally intensive, restricting investigation of stochastic effects to the single-cell level. Continuous stochastic methods that use stochastic differential equations (SDEs) to model the system are more efficient but can lead to simulations that have no biological meaning. In this paper we show that modeling the behavior of ion channel dynamics by a reflected SDE ensures biologically realistic simulations, and we argue that this model follows from the continuous approximation of the discrete-state Markov chain model. Open channel and action potential statistics from simulations of ion channel dynamics using the reflected SDE are compared with those of a discrete-state Markov chain method. Results show that the reflected SDE simulations are in good agreement with the discrete-state approach. The reflected SDE model therefore provides a computationally efficient method to simulate ion channel dynamics while preserving the distributional properties of the discrete-state Markov chain model and also ensuring biologically realistic solutions. This framework could easily be extended to other biochemical reaction networks. © 2012 American Physical Society.

Unique residues involved in activation of the multitasking protease/chaperone HtrA from Chlamydia trachomatis

DegP, a member of the HtrA family of proteins, conducts critical bacterial protein quality control by both chaperone and proteolysis activities. The regulatory mechanisms controlling these two distinct activities, however, are unknown. DegP activation is known to involve a unique mechanism of allosteric binding, conformational changes and oligomer formation. We have uncovered a novel role for the residues at the PDZ1:protease interface in oligomer formation specifically for chaperone substrates of Chlamydia trachomatis HtrA (DegP homolog). We have demonstrated that CtHtrA proteolysis could be activated by allosteric binding and oligomer formation. The PDZ1 activator cleft was required for the activation and oligomer formation. However, unique to CtHtrA was the critical role for residues at the PDZ1:protease interface in oligomer formation when the activator was an in vitro chaperone substrate. Furthermore, a potential in vivo chaperone substrate, the major outer membrane protein (MOMP) from Chlamydia, was able to activate CtHtrA and induce oligomer formation. Therefore, we have revealed novel residues involved in the activation of CtHtrA which are likely to have important in vivo implications for outer membrane protein assembly.

Synthesis and characterisation of organo-beidellite

Today, there are growing concerns about the presence of environmental pollutants in many parts of the world. In particular, a lot of attention has been drawn to the levels of water and soil contaminants (de Paiva et al., 2008). The majority of these contaminants consist of NOCs (non-ionic organic compounds) and can enter our waterways through industrial activities, mining operations, crop and animal production, waste disposal and accidental leakage (de Paiva et al., 2008; Park et al., 2011). Therefore, there is an increased interest in the synthesis of new materials that can be used to remove potentially carcinogenic and toxic water contaminants. Smectite type organoclays are widely used in numerous applications, such as sorbent agents for environmental remediation, due to their unique properties (Jiunn-Fwu et al., 1990; Sheng et al., 1996; Zhou et al., 2007; Bektas et al., 2011; Park et al., 2011). This investigation focuses on beidellite (SBId-1), which belongs to the smectite clay family. Their properties include high cation exchange capacity (CEC), swelling properties, porous, high surface area and consequential strong adsorption/absorption capacity (Xi et al., 2007). However, swelling clays in general are not an effective sorbent agent in nature due to their hydrophilic properties. The hydrophilic properties of the clay can be changed to organophilic by intercalating a cationic surfactant. Many applications of organoclays are strongly dependent on their structural properties and hence, a better understanding of the configuration and structural change of organoclay is crucial. Organoclays were synthesised through ion exchange of 21CODTMA (MW: 392.5 g mol-1) and characterised using XRD and FTIR spectroscopy. This study investigates the structural and conformational changes of beidellite intercalated with octadecyltrimethylammonium bromide.

The M(r) 43K major capsid protein of rice ragged stunt oryzavirus is a post-translationally processed product of a M(r) 67 348 polypeptide encoded by genome segment 8

The nucleotide sequence of DNA complementary to rice ragged stunt oryzavirus (RRSV) genome segment 8 (S8) of an isolate from Thailand was determined. RRSV S8 is 1 914 bp in size and contains a single large open reading frame (ORF) spanning nucleotides 23 to 1 810 which is capable of encoding a protein of M(r) 67 348. The N-terminal amino acid sequence of a ~43K virion polypeptide matched to that inferred for an internal region of the S8 coding sequence. These data suggest that the 43K protein is encoded by S8 and is derived by a proteolytic cleavage. Predicted polypeptide sizes from this possible cleavage of S8 protein are 26K and 42K. Polyclonal antibodies raised against a maltose binding protein (MBP)-S8 fusion polypeptide (expressed in Escherichia coli) recognised four RRSV particle associated polypeptides of M(r) 67K, 46K, 43K and 26K and all except the 26K polypeptide were also highly immunoreactive to polyclonal antibodies raised against purified RRSV particles. Cleavage of the MBP-S8 fusion polypeptide with protease Factor X produced the expected 40K MBP and two polypeptides of apparent M(r) 46K and 26K. Antibodies to purified RRSV particles reacted strongly with the intact fusion protein and the 46K cleavage product but weakly to the 26K product. Furthermore, in vitro transcription and translation of the S8 coding region revealed a post-translational self cleavage of the 67K polypeptide to 46K and 26K products. These data indicate that S8 encodes a structural polypeptide, the majority of which is auto- catalytically cleaved to 26K and 46K proteins. The data also suggest that the 26K protein is the self cleaving protease and that the 46K product is further processed or undergoes stable conformational changes to a ~43K major capsid protein.

A Survey of the 2010 Quartz Crystal Microbalance Literature

In 2010 there has again been an increase in the number of papers published involving piezoelectric acoustic sensors, or quartz crystal microbalances (QCM), when compared to the last period reviewed 2006-2009. The average number of QCM publications per annum was 124 in the period 2001-2005, 223 in the period 2006-9, and 273 in 2010. There are trends towards increasing use of QCM in the study of protein adsorption to surfaces (93% increase), homeostasis (67% increase), protein-protein interactions (40% increase), and carbohydrates (43% increase). New commercial systems have been released that are driving the uptake of the technology for characterisation of binding specificities, affinities, kinetics and conformational changes associated with a molecular recognition event. This article highlights theoretical and practical aspects of the principals that underpin acoustic analysis, then reviews exemplary papers in key application areas involving small molecular weight ligands, carbohydrates, proteins, nucleic acids, viruses, bacteria, cells, and membrane interfaces.

Alteration of the C-terminal ligand specificity of the Erbin PDZ domain by allosteric mutational effects

Modulation of protein binding specificity is important for basic biology and for applied science. Here we explore how binding specificity is conveyed in PDZ (postsynaptic density protein-95/discs large/zonula occludens-1) domains, small interaction modules that recognize various proteins by binding to an extended C terminus. Our goal was to engineer variants of the Erbin PDZ domain with altered specificity for the most C-terminal position (position 0) where a Val is strongly preferred by the wild-type domain. We constructed a library of PDZ domains by randomizing residues in direct contact with position 0 and in a loop that is close to but does not contact position 0. We used phage display to select for PDZ variants that bind to 19 peptide ligands differing only at position 0. To verify that each obtained PDZ domain exhibited the correct binding specificity, we selected peptide ligands for each domain. Despite intensive efforts, we were only able to evolve Erbin PDZ domain variants with selectivity for the aliphatic C-terminal side chains Val, Ile and Leu. Interestingly, many PDZ domains with these three distinct specificities contained identical amino acids at positions that directly contact position 0 but differed in the loop that does not contact position 0. Computational modeling of the selected PDZ domains shows how slight conformational changes in the loop region propagate to the binding site and result in different binding specificities. Our results demonstrate that second-sphere residues could be crucial in determining protein binding specificity.

Molecular dynamics simulations of CXCL-8 and its interactions with a receptor peptide, heparin fragments, and sulfated linked cyclitols

CXCL-8 (Interleukin 8) is a CXC chemokine with a central role in the human immune response. We have undertaken extensive in silico analyses to elucidate the interactions of CXCL-8 with its various binding partners, which are crucial for its biological function. Sequence and structure analyses showed that residues in the thirdq β-sheet and basic residues in the heparin binding site are highly variable, while residues in the second β-sheet are highly conserved. Molecular dynamics simulations in aqueous solution of dimeric CXCL-8 have been performed with starting geometries from both X-ray and NMR structures showed shearing movements between the two antiparallel C-terminal helices. Dynamic conservation analyses of these simulations agreed with experimental data indicating that structural differences between the two structures at quaternary level arise from changes in the secondary structure of the N-terminal loop, the 310-helix, the 30s, 40s, and 50s loops and the third β-sheet, resulting in a different interhelical separation. Nevertheless, the observation of these different states indicates that CXCL-8 has the potential to undergo conformational changes, and it seems likely that this feature is relevant to the mode of binding of glycosaminoglycan (GAG) mimetics such as cyclitols. Simulations of the receptor peptide fragment−CXCL-8 complex identified several specific interactions of the receptor peptide with CXCL-8 that could be exploited in the structure-based design of competitive peptides and nonpeptidic molecules targeting CXCL-8 for combating inflammatory diseases. Simulations of the CXCL-8 dimer complexed with a 24-mer heparin fragment and of the CXCL-8−receptor peptide complex revealed that Arg60, Lys64, and Arg68 in the dimer bind to cyclitols in a horseshoe pattern, defining a region which is spatially distinct from the receptor binding site. There appears to be an optimum number of sulfates and an optimum length of alkyl spacers required for the interaction of cyclitol inhibitors with the dimeric form of CXCL-8. Calculation of the binding affinities of cyclitol inhibitors reflected satisfactorily the ranking of experimentally determined inhibitory potencies. The findings of these molecular modeling studies will help in the search for inhibitors which can modulate various CXCL-8 biological activities and serve as an excellent model system to study CXC-inhibitor interactions.

Free energy calculations of the interactions of c-Jun-based synthetic peptides with the c-Fos protein

The c-Fos–c-Jun complex forms the activator protein 1 transcription factor, a therapeutic target in the treatment of cancer. Various synthetic peptides have been designed to try to selectively disrupt the interaction between c-Fos and c-Jun at its leucine zipper domain. To evaluate the binding affinity between these synthetic peptides and c-Fos, polarizable and nonpolarizable molecular dynamics (MD) simulations were conducted, and the resulting conformations were analyzed using the molecular mechanics generalized Born surface area (MM/GBSA) method to compute free energies of binding. In contrast to empirical and semiempirical approaches, the estimation of free energies of binding using a combination of MD simulations and the MM/GBSA approach takes into account dynamical properties such as conformational changes, as well as solvation effects and hydrophobic and hydrophilic interactions. The predicted binding affinities of the series of c-Jun-based peptides targeting the c-Fos peptide show good correlation with experimental melting temperatures. This provides the basis for the rational design of peptides based on internal, van der Waals, and electrostatic interactions.

Conformational Spread in the Flagellar Motor Switch: A Model Study

The reliable response to weak biological signals requires that they be amplified with fidelity. In E. coli, the flagellar motors that control swimming can switch direction in response to very small changes in the concentration of the signaling protein CheY-P, but how this works is not well understood. A recently proposed allosteric model based on cooperative conformational spread in a ring of identical protomers seems promising as it is able to qualitatively reproduce switching, locked state behavior and Hill coefficient values measured for the rotary motor. In this paper we undertook a comprehensive simulation study to analyze the behavior of this model in detail and made predictions on three experimentally observable quantities: switch time distribution, locked state interval distribution, Hill coefficient of the switch response. We parameterized the model using experimental measurements, finding excellent agreement with published data on motor behavior. Analysis of the simulated switching dynamics revealed a mechanism for chemotactic ultrasensitivity, in which cooperativity is indispensable for realizing both coherent switching and effective amplification. These results showed how cells can combine elements of analog and digital control to produce switches that are simultaneously sensitive and reliable. © 2012 Ma et al.