The structures of COPI-coated vesicles reveal alternate coatomer conformations and interactions

Transport between compartments of eukaryotic cells is mediated by coated vesicles. The archetypal protein coats COPI, COPII, and clathrin are conserved from yeast to human. Structural studies of COPII and clathrin coats assembled in vitro without membranes suggest that coat components assemble regular cages with the same set of interactions between components. Detailed three-dimensional structures of coated membrane vesicles have not been obtained. Here, we solved the structures of individual COPI-coated membrane vesicles by cryoelectron tomography and subtomogram averaging of in vitro reconstituted budding reactions. The coat protein complex, coatomer, was observed to adopt alternative conformations to change the number of other coatomers with which it interacts and to form vesicles with variable sizes and shapes. This represents a fundamentally different basis for vesicle coat assembly.

Vibrational spectroscopic studies of wool

Fourier transfonn (FT) Raman, Raman microspectroscopy and Fourier transform infrared (FTIR) spectroscopy have been used for the structural analysis and characterisation of untreated and chemically treated wool fibres. For FT -Raman spectroscopy novel methods of sample presentation have been developed and optimised for the analysis of wool. No significant fluorescence was observed and the spectra could be obtained routinely. The stability of wool keratin to the laser source was investigated and the visual and spectroscopic signs of sample damage were established. Wool keratin was found to be extremely robust with no signs of sample degradation observed for laser powers of up to 600 m W and for exposure times of up to seven and half hours. Due to improvements in band resolution and signal-to-noise ratio, several previously unobserved spectral features have become apparent. The assignment of the Raman active vibrational modes of wool have been reviewed and updated to include these features. The infrared spectroscopic techniques of attenuated total reflectance (ATR) and photoacoustic (P A) have been used to examine shrinkproofed and mothproofed wool samples. Shrinkproofing is an oxidative chemical treatment used to selectively modifY the surface of a wool fibre. Mothproofing is a chemical treatment applied to wool for the prevention of insect attack. The ability of PAS and A TR to vary the penetration depth by varying certain instrumental parameters was used to obtain spectra of the near surface regions of these chemically treated samples. These spectra were compared with those taken with a greater penetration depth, which therefore represent more of the bulk wool sample. The PA and ATR spectra demonstrated that oxidation was restricted to the near-surface layer of wool. Extensive curve fitting of ATR spectra of untreated wool indicated that cuticle was composed of a mixed protein conformation, but was predominately that of an a.-helix. The cortex was proposed to be a mixture of both a.helical and ~-pleated sheet protein conformations. These findings were supported by PAS depth profiling results. Raman microspectroscopy was used in an extensive investigation of the molecular structure of the wool fibre. This included determining the orientation of certain functional groups within the wool fibre and the symmetry of particular vibrations. The orientation ofbonds within the wool fibre was investigated by orientating the wool fibre axis parallel and then perpendicular to the plane of polarisation of the electric vector of the incident radiation. It was experimentally determined that the majority of C=O and N-H bonds of the peptide bond of wool lie parallel to the fibre axis. Additionally, a number of the important vibrations associated with the a-helix were also found to lie parallel to the fibre axis. Further investigation into the molecular structure of wool involved determining what effect stretching the wool fibre had on bond orientation. Raman spectra of stretched and unstretched wool fibres indicated that extension altered the orientation ofthe aromatic rings, the CH2 and CH3 groups of the amino acids. Curve fitting results revealed that extension resulted in significant destruction of the a-helix structure a substantial increase in the P-pleated sheet structure. Finally, depolarisation ratios were calculated for Raman spectra. The vibrations associated with the aromatic rings of amino acids had very low ratios which indicated that the vibrations were highly symmetrical.

Rational design of serine protease inhibitors

Proteases regulate a spectrum of diverse physiological processes, and dysregulation of proteolytic activity drives a plethora of pathological conditions. Understanding protease function is essential to appreciating many aspects of normal physiology and progression of disease. Consequently, development of potent and specific inhibitors of proteolytic enzymes is vital to provide tools for the dissection of protease function in biological systems and for the treatment of diseases linked to aberrant proteolytic activity. The studies in this thesis describe the rational design of potent inhibitors of three proteases that are implicated in disease development. Additionally, key features of the interaction of proteases and their cognate inhibitors or substrates are analysed and a series of rational inhibitor design principles are expounded and tested. Rational design of protease inhibitors relies on a comprehensive understanding of protease structure and biochemistry. Analysis of known protease cleavage sites in proteins and peptides is a commonly used source of such information. However, model peptide substrate and protein sequences have widely differing levels of backbone constraint and hence can adopt highly divergent structures when binding to a protease’s active site. This may result in identical sequences in peptides and proteins having different conformations and diverse spatial distribution of amino acid functionalities. Regardless of this, protein and peptide cleavage sites are often regarded as being equivalent. One of the key findings in the following studies is a definitive demonstration of the lack of equivalence between these two classes of substrate and invalidation of the common practice of using the sequences of model peptide substrates to predict cleavage of proteins in vivo. Another important feature for protease substrate recognition is subsite cooperativity. This type of cooperativity is commonly referred to as protease or substrate binding subsite cooperativity and is distinct from allosteric cooperativity, where binding of a molecule distant from the protease active site affects the binding affinity of a substrate. Subsite cooperativity may be intramolecular where neighbouring residues in substrates are interacting, affecting the scissile bond’s susceptibility to protease cleavage. Subsite cooperativity can also be intermolecular where a particular residue’s contribution to binding affinity changes depending on the identity of neighbouring amino acids. Although numerous studies have identified subsite cooperativity effects, these findings are frequently ignored in investigations probing subsite selectivity by screening against diverse combinatorial libraries of peptides (positional scanning synthetic combinatorial library; PS-SCL). This strategy for determining cleavage specificity relies on the averaged rates of hydrolysis for an uncharacterised ensemble of peptide sequences, as opposed to the defined rate of hydrolysis of a known specific substrate. Further, since PS-SCL screens probe the preference of the various protease subsites independently, this method is inherently unable to detect subsite cooperativity. However, mean hydrolysis rates from PS-SCL screens are often interpreted as being comparable to those produced by single peptide cleavages. Before this study no large systematic evaluation had been made to determine the level of correlation between protease selectivity as predicted by screening against a library of combinatorial peptides and cleavage of individual peptides. This subject is specifically explored in the studies described here. In order to establish whether PS-SCL screens could accurately determine the substrate preferences of proteases, a systematic comparison of data from PS-SCLs with libraries containing individually synthesised peptides (sparse matrix library; SML) was carried out. These SML libraries were designed to include all possible sequence combinations of the residues that were suggested to be preferred by a protease using the PS-SCL method. SML screening against the three serine proteases kallikrein 4 (KLK4), kallikrein 14 (KLK14) and plasmin revealed highly preferred peptide substrates that could not have been deduced by PS-SCL screening alone. Comparing protease subsite preference profiles from screens of the two types of peptide libraries showed that the most preferred substrates were not detected by PS SCL screening as a consequence of intermolecular cooperativity being negated by the very nature of PS SCL screening. Sequences that are highly favoured as result of intermolecular cooperativity achieve optimal protease subsite occupancy, and thereby interact with very specific determinants of the protease. Identifying these substrate sequences is important since they may be used to produce potent and selective inhibitors of protolytic enzymes. This study found that highly favoured substrate sequences that relied on intermolecular cooperativity allowed for the production of potent inhibitors of KLK4, KLK14 and plasmin. Peptide aldehydes based on preferred plasmin sequences produced high affinity transition state analogue inhibitors for this protease. The most potent of these maintained specificity over plasma kallikrein (known to have a very similar substrate preference to plasmin). Furthermore, the efficiency of this inhibitor in blocking fibrinolysis in vitro was comparable to aprotinin, which previously saw clinical use to reduce perioperative bleeding. One substrate sequence particularly favoured by KLK4 was substituted into the 14 amino acid, circular sunflower trypsin inhibitor (SFTI). This resulted in a highly potent and selective inhibitor (SFTI-FCQR) which attenuated protease activated receptor signalling by KLK4 in vitro. Moreover, SFTI-FCQR and paclitaxel synergistically reduced growth of ovarian cancer cells in vitro, making this inhibitor a lead compound for further therapeutic development. Similar incorporation of a preferred KLK14 amino acid sequence into the SFTI scaffold produced a potent inhibitor for this protease. However, the conformationally constrained SFTI backbone enforced a different intramolecular cooperativity, which masked a KLK14 specific determinant. As a consequence, the level of selectivity achievable was lower than that found for the KLK4 inhibitor. Standard mechanism inhibitors such as SFTI rely on a stable acyl-enzyme intermediate for high affinity binding. This is achieved by a conformationally constrained canonical binding loop that allows for reformation of the scissile peptide bond after cleavage. Amino acid substitutions within the inhibitor to target a particular protease may compromise structural determinants that support the rigidity of the binding loop and thereby prevent the engineered inhibitor reaching its full potential. An in silico analysis was carried out to examine the potential for further improvements to the potency and selectivity of the SFTI-based KLK4 and KLK14 inhibitors. Molecular dynamics simulations suggested that the substitutions within SFTI required to target KLK4 and KLK14 had compromised the intramolecular hydrogen bond network of the inhibitor and caused a concomitant loss of binding loop stability. Furthermore in silico amino acid substitution revealed a consistent correlation between a higher frequency of formation and the number of internal hydrogen bonds of SFTI-variants and lower inhibition constants. These predictions allowed for the production of second generation inhibitors with enhanced binding affinity toward both targets and highlight the importance of considering intramolecular cooperativity effects when engineering proteins or circular peptides to target proteases. The findings from this study show that although PS-SCLs are a useful tool for high throughput screening of approximate protease preference, later refinement by SML screening is needed to reveal optimal subsite occupancy due to cooperativity in substrate recognition. This investigation has also demonstrated the importance of maintaining structural determinants of backbone constraint and conformation when engineering standard mechanism inhibitors for new targets. Combined these results show that backbone conformation and amino acid cooperativity have more prominent roles than previously appreciated in determining substrate/inhibitor specificity and binding affinity. The three key inhibitors designed during this investigation are now being developed as lead compounds for cancer chemotherapy, control of fibrinolysis and cosmeceutical applications. These compounds form the basis of a portfolio of intellectual property which will be further developed in the coming years.

Homo- and heteronuclear meso,meso-(E)-Ethene-1,2-diyl-Linked Diporphyrins : preparation, X-ray crystal structure, electronic absorption and emission spectra and density functional theory calculations

Homo-and heteronuclear meso,meso-(E)-ethene-1,2-diyl-linked diporphyrins have been prepared by the Suzuki coupling of porphyrinylboronates and iodovinylporphyrins. Combinations comprising 5,10,15-triphenylporphyrin (TriPP) on both ends of the ethene-1,2-diyl bridge M 210 (M 2=H 2/Ni, Ni 2, Ni/Zn, H 4, H 2Zn, Zn 2) and 5,15-bis(3,5-di-tert-butylphenyl)porphyrinato-nickel(II) on one end and H 2, Ni, and ZnTriPP on the other (M 211), enable the first studies of this class of compounds possessing intrinsic polarity. The compounds were characterized by electronic absorption and steady state emission spectra, 1H NMR spectra, and for the Ni 2 bis(TriPP) complex Ni 210, single crystal X-ray structure determination. The crystal structure shows ruffled distortions of the porphyrin rings, typical of Ni II porphyrins, and the (E)-C 2H 2 bridge makes a dihedral angle of 50° with the mean planes of the macrocycles. The result is a stepped parallel arrangement of the porphyrin rings. The dihedral angles in the solid state reflect the interplay of steric and electronic effects of the bridge on interporphyrin communication. The emission spectra in particular, suggest energy transfer across the bridge is fast in conformations in which the bridge is nearly coplanar with the rings. Comparisons of the fluorescence behaviour of H 410 and H 2Ni10 show strong quenching of the free base fluorescence when the complex is excited at the lower energy component of the Soret band, a feature associated in the literature with more planar conformations. TDDFT calculations on the gas-phase optimized geometry of Ni 210 reproduce the features of the experimental electronic absorption spectrum within 0.1 eV. © 2012 WILEY-VCH Verlag GmbH & Co. KGaA, Weinheim.

Bis(benzylaminium) 4,5-dichlorobenzene-1,2-dicarboxylate monohydrate

In the structure of the title salt 2C7H10N+.C8H2Cl2O4(2-) .H2O the two benzylaminium anions have different conformations, one being essentially planar the other having the side-chain rotated out of the benzene plane (minimum ring to side-chain C-C-C-N torsion angles = -3.6(6) and 50.1(5)\%, respectively). In the 4,5-dichlorophthalate dianion the carboxyl groups make ihedral angles of 23.0(2) and 76.5(2)\% with the benzene ring. Aminium N-H...O and water O-H...O hydrogen-bonding associations with carboxyl O-atom acceptors give a two-dimensional duplex sheet structure which extends along the (011) plane. Weak pi-pi interactions are also present between the benzene ring and one of the cation rings [minimum ring centroid separation = 3.749(3)Ang.

Structure of the immature retroviral capsid at 8 Å resolution by cryo-electron microscopy

The assembly of retroviruses such as HIV-1 is driven by oligomerization of their major structural protein, Gag. Gag is a multidomain polyprotein including three conserved folded domains: MA (matrix), CA (capsid) and NC (nucleocapsid)(1). Assembly of an infectious virion proceeds in two stages(2). In the first stage, Gag oligomerization into a hexameric protein lattice leads to the formation of an incomplete, roughly spherical protein shell that buds through the plasma membrane of the infected cell to release an enveloped immature virus particle. In the second stage, cleavage of Gag by the viral protease leads to rearrangement of the particle interior, converting the non-infectious immature virus particle into a mature infectious virion. The immature Gag shell acts as the pivotal intermediate in assembly and is a potential target for anti-retroviral drugs both in inhibiting virus assembly and in disrupting virus maturation(3). However, detailed structural information on the immature Gag shell has not previously been available. For this reason it is unclear what protein conformations and interfaces mediate the interactions between domains and therefore the assembly of retrovirus particles, and what structural transitions are associated with retrovirus maturation. Here we solve the structure of the immature retroviral Gag shell from Mason-Pfizer monkey virus by combining cryo-electron microscopy and tomography. The 8-angstrom resolution structure permits the derivation of a pseudo-atomic model of CA in the immature retrovirus, which defines the protein interfaces mediating retrovirus assembly. We show that transition of an immature retrovirus into its mature infectious form involves marked rotations and translations of CA domains, that the roles of the amino-terminal and carboxy-terminal domains of CA in assembling the immature and mature hexameric lattices are exchanged, and that the CA interactions that stabilize the immature and mature viruses are almost completely distinct.

Electrochemical detection of dopamine and cytochrome c at a nanostructured gold electrode

A nanostructured gold surface consisting of closely packed outwardly growing spikes is investigated for the electrochemical detection of dopamine and cytochrome c. A significant electrocatalytic effect for the electrooxidation of both dopamine and ascorbic acid at the nanostructured electrode was found due to the presence of surface active sites which allowed the detection of dopamine in the presence of excess ascorbic acid to be achieved by differential pulse voltammetry. By simple modification with a layer of Nafion, the enhanced electrocatalytic properties of the nanostructured surface was maintained while increasing the selectivity of dopamine detection in the presence of interfering species such as excess ascorbic and uric acids. Also, upon modification of the nanostructured surface with a monolayer of cysteine, the electrochemical response of immobilised cytochrome c in two distinct conformations was observed. This opens up the possibility of using such a nanostructured surface for the characterisation of other biomolecules and in bio-electroanalytical applications.

Differential Regulation of Clusterin Isoforms by the Androgen Receptor

Clusterin (CLU) was initially reported as an androgen-repressed gene which is now shown to be an androgen-regulated ATP-independent cytoprotective molecular chaperone. CLU binds to a wide variety of client proteins to potently inhibit stress-induced protein aggregation and chaperone or stabilise conformations of proteins at times of cell stress. CLU is an enigmatic protein, being ascribed both pro- and anti-apoptotic roles. Recent evidence has shown that both secreted (sCLU) and nuclear (nCLU) isoforms can be produced, and that protein function is dependent on the sub-cellular localisation. We and others have shown that sCLU is cytoprotective, while nCLU is pro-apoptotic. It now seems likely that the apparently dichotomous functions of CLU result from the expression of different but related CLU isoforms and splice variants, and that cell survival depends in part on the relative expression of pro- versus anti-apoptotic CLU proteins. In cancer cells, increased sCLU expression is associated with increased resistance to apoptotic triggers and treatment resistance. CLU is a stress-induced protein upregulated after apoptotic triggers like androgen ablation and chemotherapy. Treatment strategies targeting stress-associated increases in sCLU expression enhance treatment-induced apoptosis and delay the emergence of androgen independence. Differential regulation of CLU isoforms and splice variants by androgens may be a pathway whereby cancer cells develop treatment resistance and evade apoptosis.

Sulfur-containing amide-based [2]rotaxanes and molecular shuttles

We report the synthesis, structure and properties of [2]rotaxanes prepared by the assembly of benzylic amide macrocycles around a series of amide and sulfide-/sulfoxide-/sulfone-containing threads. The efficacy of rotaxane formation is related to the hydrogen bond accepting properties of the various sulfur-containing functional groups in the thread, with the highest yields (up to 63% with a rigid vinyl spacer in the template site) obtained for sulfoxide rotaxanes. X-Ray crystallography of a sulfoxide rotaxane, 5, shows that the macrocycle adopts a highly symmetrical chair-like conformation in the solid state, with short hydrogen bonds between the macrocycle isophthalamide NH-protons and the amide carbonyl and sulfoxide S-O of the thread. In contrast, in the X-ray crystal structures of the analogous sulfide (4) and sulfone (6) rotaxanes the macrocycle adopts boat-like conformations with long intercomponent NH…O=SO and NH…S hydrogen bonds (in addition to several intercomponent amide-amide hydrogen bonds). Taking advantage of the different hydrogen bonding modes of the sulfur-based functional groups, a switchable molecular shuttle was prepared in which the oxidation level of sulfur determines the position of the macrocycle on the thread.

Crystal structures and hydrogen bonding in the anhydrous tryptaminium salts of the isomeric (2,4-dichlorophenoxy)acetic and (3,5-dichlorophenoxy)acetic acids

The anhydrous salts of 1H-indole-3-ethanamine (tryptamine) with isomeric (2,4-dichlorophenoxy)acetic acid (2,4-D) and (3,5-dichlorophenoxy)acetic (3,5-D), C10H13N2+ (C8H5Cl2O3)-, [(I) and (II), respectively] have been determined and their one-dimensional hydrogen-bonded polymeric structures are described. In the crystal of (I),the aminium H-atoms are involved in three separate inter-species N-H...O hydrogen-bonding interactions, two with carboxyl O-atom acceptors and the third in an asymmetric three-centre bidentate carboxyl O,O' chelate [graph set R2/1(4)]. The indole H-atom forms an N-H...O~carboxyl~ hydrogen bond, extending the chain structure along the b axial direction. In (II), two of the three aminium H-atoms are also involved in N-H...O(carboxyl) hydrogen bonds similar to (I) but with the third, a three-centre asymmetric interaction with carboxyl and phenoxy O-atoms is found [graph set R2/1(5)]. The chain polymeric extension is also along b. There are no pi--pi ring interactions in either of the structures. The aminium side chain conformations differ significantly between the two structures, reflecting the conformational ambivalence of the tryptaminium cation, as found also in the benzoate salts.

Experimental and theoretical characterization of 5,10- diminoporphodimethenes: Dearomatized porphyrinoids from palladium-catalyzed hydrazinations of 5,10-diarylporphyrins

A new dearomatized porphyrinoid, 5,10-diiminoporphodimethene (5,10-DIPD), has been prepared by palladium-catalyzed hydrazination of 5,10-dibromo-15,20-bis(3,5-di-tert-butylphenyl)porphyrin and its nickel(II) complex, by using ethyl and 4-methoxybenzyl carbazates. The oxidative dearomatization of the porphyrin ring occurs in high yield. Further oxidation with 2,3-dichloro-5,6-dicyanobenzoquinone forms the corresponding 5,10-bis(azocarboxylates), thereby restoring the porphyrin aromaticity. The UV/visible spectra of the NiII DIPDs exhibit remarkable redshifts of the lowest-energy bands to 780 nm, and differential pulse voltammetry reveals a contracted electrochemical HOMO–LUMO gap of 1.44 V. Density functional theory (DFT) was used to calculate the optimized geometries and frontier molecular orbitals of model 5,10-DIPD Ni7c and 5,10-bis(azocarboxylate) Ni8c. The conformations of the carbamate groups and the configurations of the CNZ unit were considered in conjunction with the NOESY spectra, to generate the global minimum geometry and two other structures with slightly higher energies. In the absence of solution data regarding conformations, ten possible local minimum conformations were considered for Ni8c. Partition of the porphyrin macrocycle into tri- and monopyrrole fragments in Ni7c and the inclusion of terminal conjugating functional groups generate unique frontier molecular orbital distributions and a HOMO–LUMO transition with a strong element of charge transfer from the monopyrrole ring. Time-dependent DFT calculations were performed for the three lowest-energy structures of Ni7c and Ni8c, and weighting according to their energies allowed the prediction of the electronic spectra. The calculations reproduce the lower-energy regions of the spectra and the overall forms of the spectra with high accuracy, but agreement is not as good in the Soret region below 450 nm.

Crystal structures of the endoplasmic reticulum aminopeptidase-1 (ERAP1) reveal the molecular basis for N-terminal peptide trimming

Endoplasmatic reticulum aminopeptidase 1 (ERAP1) is a multifunctional enzyme involved in trimming of peptides to an optimal length for presentation by major histocompatibility complex (MHC) class I molecules. Polymorphisms in ERAP1 have been associated with chronic inflammatory diseases, including ankylosing spondylitis (AS) and psoriasis, and subsequent in vitro enzyme studies suggest distinct catalytic properties of ERAP1 variants. To understand structure-activity relationships of this enzyme we determined crystal structures in open and closed states of human ERAP1, which provide the first snapshots along a catalytic path. ERAP1 is a zinc-metallopeptidase with typical H-E-X-X-H-(X)18-E zinc binding and G-A-M-E-N motifs characteristic for members of the gluzincin protease family. The structures reveal extensive domain movements, including an active site closure as well as three different open conformations, thus providing insights into the catalytic cycle. A K 528R mutant strongly associated with AS in GWAS studies shows significantly altered peptide processing characteristics, which are possibly related to impaired interdomain interactions.

Computational analyses of the catalytic and heparin-binding sites and their interactions with glycosaminoglycans in glycoside hydrolase family 79 endo-d-glucuronidase (heparanase)

Mammalian heparanase is an endo-β-glucuronidase associated with cell invasion in cancer metastasis, angiogenesis and inflammation. Heparanase cleaves heparan sulfate proteoglycans in the extracellular matrix and basement membrane, releasing heparin/heparan sulfate oligosaccharides of appreciable size. This in turn causes the release of growth factors, which accelerate tumor growth and metastasis. Heparanase has two glycosaminoglycan-binding domains; however, no three-dimensional structure information is available for human heparanase that can provide insights into how the two domains interact to degrade heparin fragments. We have constructed a new homology model of heparanase that takes into account the most recent structural and bioinformatics data available. Heparin analogs and glycosaminoglycan mimetics were computationally docked into the active site with energetically stable ring conformations and their interaction energies were compared. The resulting docked structures were used to propose a model for substrates and conformer selectivity based on the dimensions of the active site. The docking of substrates and inhibitors indicates the existence of a large binding site extending at least two saccharide units beyond the cleavage site (toward the nonreducing end) and at least three saccharides toward the reducing end (toward heparin-binding site 2). The docking of substrates suggests that heparanase recognizes the N-sulfated and O-sulfated glucosamines at subsite +1 and glucuronic acid at the cleavage site, whereas in the absence of 6-O-sulfation in glucosamine, glucuronic acid is docked at subsite +2. These findings will help us to focus on the rational design of heparanase-inhibiting molecules for anticancer drug development by targeting the two heparin/heparan sulfate recognition domains.

Can current force fields reproduce ring puckering in 2-O-sulfo-α-l-iduronic acid? A molecular dynamics simulation study

The monosaccharide 2-O-sulfo-α-l-iduronic acid (IdoA2S) is one of the major components of glycosaminoglycans. The ability of molecular mechanics force fields to reproduce ring-puckering conformational equilibrium is important for the successful prediction of the free energies of interaction of these carbohydrates with proteins. Here we report unconstrained molecular dynamics simulations of IdoA2S monosaccharide that were carried out to investigate the ability of commonly used force fields to reproduce its ring conformational flexibility in aqueous solution. In particular, the distribution of ring conformer populations of IdoA2S was determined. The GROMOS96 force field with the SPC/E water potential can predict successfully the dominant skew-boat to chair conformational transition of the IdoA2S monosaccharide in aqueous solution. On the other hand, the GLYCAM06 force field with the TIP3P water potential sampled transitional conformations between the boat and chair forms. Simulations using the GROMOS96 force field showed no pseudorotational equilibrium fluctuations and hence no inter-conversion between the boat and twist boat ring conformers. Calculations of theoretical proton NMR coupling constants showed that the GROMOS96 force field can predict the skew-boat to chair conformational ratio in good agreement with the experiment, whereas GLYCAM06 shows worse agreement. The omega rotamer distribution about the C5–C6 bond was predicted by both force fields to have torsions around 10°, 190°, and 360°.

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.