Controlling protein molecular dynamics:how to accelerate folding while preserving the native state

The dynamics of peptides and proteins generated by classical molecular dynamics (MD) is described by using a Markov model. The model is built by clustering the trajectory into conformational states and estimating transition probabilities between the states. Assuming that it is possible to influence the dynamics of the system by varying simulation parameters, we show how to use the Markov model to determine the parameter values that preserve the folded state of the protein and at the same time, reduce the folding time in the simulation. We investigate this by applying the method to two systems. The first system is an imaginary peptide described by given transition probabilities with a total folding time of 1 micros. We find that only small changes in the transition probabilities are needed to accelerate (or decelerate) the folding. This implies that folding times for slowly folding peptides and proteins calculated using MD cannot be meaningfully compared to experimental results. The second system is a four residue peptide valine-proline-alanine-leucine in water. We control the dynamics of the transitions by varying the temperature and the atom masses. The simulation results show that it is possible to find the combinations of parameter values that accelerate the dynamics and at the same time preserve the native state of the peptide. A method for accelerating larger systems without performing simulations for the whole folding process is outlined.

Why are MD simulated protein folding times wrong?

The question of significant deviations of protein folding times simulated using molecular dynamics from experimental values is investigated. It is shown that in the framework of Markov State Model (MSM) describing the conformational dynamics of peptides and proteins, the folding time is very sensitive to the simulation model parameters, such as forcefield and temperature. Using two peptides as examples, we show that the deviations in the folding times can reach an order of magnitude for modest variations of the molecular model. We, therefore, conclude that the folding rate values obtained in molecular dynamics simulations have to be treated with care.

Dynamic complexity of chaotic transitions in high-dimensional classical dynamics:Leu-Enkephalin folding

Leu-Enkephalin in explicit water is simulated using classical molecular dynamics. A ß-turn transition is investigated by calculating the topological complexity (in the "computational mechanics" framework [J. P. Crutchfield and K. Young, Phys. Rev. Lett., 63, 105 (1989)]) of the dynamics of both the peptide and the neighbouring water molecules. The complexity of the atomic trajectories of the (relatively short) simulations used in this study reflect the degree of phase space mixing in the system. It is demonstrated that the dynamic complexity of the hydrogen atoms of the peptide and almost all of the hydrogens of the neighbouring waters exhibit a minimum precisely at the moment of the ß-turn transition. This indicates the appearance of simplified periodic patterns in the atomic motion, which could correspond to high-dimensional tori in the phase space. It is hypothesized that this behaviour is the manifestation of the effect described in the approach to molecular transitions by Komatsuzaki and Berry [T. Komatsuzaki and R.S. Berry, Adv. Chem. Phys., 123, 79 (2002)], where a "quasi-regular" dynamics at the transition is suggested. Therefore, for the first time, the less chaotic character of the folding transition in a realistic molecular system is demonstrated. © Springer-Verlag Berlin Heidelberg 2006.

Using directional curvatures to visualize folding patterns of the GTM projection manifolds

In data visualization, characterizing local geometric properties of non-linear projection manifolds provides the user with valuable additional information that can influence further steps in the data analysis. We take advantage of the smooth character of GTM projection manifold and analytically calculate its local directional curvatures. Curvature plots are useful for detecting regions where geometry is distorted, for changing the amount of regularization in non-linear projection manifolds, and for choosing regions of interest when constructing detailed lower-level visualization plots.

Functional split and crosslinking of the membrane domain of the β subunit of proton-translocating transhydrogenase from Escherichia coli

Proton pumping nicotinamide nucleotide transhydrogenase from Escherichia coli contains an α subunit with the NAD(H)-binding domain I and a β subunit with the NADP(H)-binding domain III. The membrane domain (domain II) harbors the proton channel and is made up of the hydrophobic parts of the α and β subunits. The interface in domain II between the α and the β subunits has previously been investigated by cross-linking loops connecting the four transmembrane helices in the α subunit and loops connecting the nine transmembrane helices in the β subunit. However, to investigate the organization of the nine transmembrane helices in the β subunit, a split was introduced by creating a stop codon in the loop connecting transmembrane helices 9 and 10 by a single mutagenesis step, utilizing an existing downstream start codon. The resulting enzyme was composed of the wild-type α subunit and the two new peptides β1 and β2. As compared to other split membrane proteins, the new transhydrogenase was remarkably active and catalyzed activities for the reduction of 3-acetylpyridine-NAD + by NADPH, the cyclic reduction of 3-acetylpyridine-NAD + by NADH (mediated by bound NADP(H)), and proton pumping, amounting to about 50-107% of the corresponding wild-type activities. These high activities suggest that the α subunit was normally folded, followed by a concerted folding of β1 + β2. Cross-linking of a βS105C-βS237C double cysteine mutant in the functional split cysteine-free background, followed by SDS-PAGE analysis, showed that helices 9, 13, and 14 were in close proximity. This is the first time that cross-linking between helices in the same β subunit has been demonstrated.

Mutations of beta-arrestin 2 that limit self-association also interfere with interactions with the beta2-adrenoceptor and the ERK1/2 MAPKs:implications for beta2-adrenoceptor signalling via the ERK1/2 MAPKs

FRET (fluorescence resonance energy transfer) and co-immunoprecipitation studies confirmed the capacity of beta-arrestin 2 to self-associate. Amino acids potentially involved in direct protein-protein interaction were identified via combinations of spot-immobilized peptide arrays and mapping of surface exposure. Among potential key amino acids, Lys(285), Arg(286) and Lys(295) are part of a continuous surface epitope located in the polar core between the N- and C-terminal domains. Introduction of K285A/R286A mutations into beta-arrestin 2-eCFP (where eCFP is enhanced cyan fluorescent protein) and beta-arrestin 2-eYFP (where eYFP is enhanced yellow fluorescent protein) constructs substantially reduced FRET, whereas introduction of a K295A mutation had a more limited effect. Neither of these mutants was able to promote beta2-adrenoceptor-mediated phosphorylation of the ERK1/2 (extracellular-signal-regulated kinase 1/2) MAPKs (mitogen-activated protein kinases). Both beta-arrestin 2 mutants displayed limited capacity to co-immunoprecipitate ERK1/2 and further spot-immobilized peptide arrays indicated each of Lys(285), Arg(286) and particularly Lys(295) to be important for this interaction. Direct interactions between beta-arrestin 2 and the beta2-adrenoceptor were also compromised by both K285A/R286A and K295A mutations of beta-arrestin 2. These were not non-specific effects linked to improper folding of beta-arrestin 2 as limited proteolysis was unable to distinguish the K285A/R286A or K295A mutants from wild-type beta-arrestin 2, and the interaction of beta-arrestin 2 with JNK3 (c-Jun N-terminal kinase 3) was unaffected by the K285A/R286A or L295A mutations. These results suggest that amino acids important for self-association of beta-arrestin 2 also play an important role in the interaction with both the beta2-adrenoceptor and the ERK1/2 MAPKs. Regulation of beta-arrestin 2 self-association may therefore control beta-arrestin 2-mediated beta2-adrenoceptor-ERK1/2 MAPK signalling.

Undue influence and the family home

This book provides a comprehensive analysis of the doctrine of undue influence in the context of the family home and fully incorporates the recent House of Lords ruling in Royal Bank of Scotland v Etridge (No 2) (2001). It is aimed predominantly at the legal practitioner, but will also act as a useful source of reference for academics and students of contract, land law and equity and trusts. Emphasis is placed on claims brought by spouses (usually the wife) seeking to set aside a charge over the matrimonial home made in favour of a lending institution. The role of lenders in this context is also examined in depth, as is the part played by the solicitor acting on behalf of the parties. Apart from providing an exposition of the doctrine and its key elements, the book also gives a broader outlook by examining the Commonwealth experience (notably in Australia, Canada and New Zealand) and suggesting an underlying concept of unconscionability as governing undue influence claims. There is also a separate chapter on remedies, as well as an appendix containing a number of draft pleadings for use by the legal practitioner. In the foreword, the Honourable Mr Justice Neuberger writes: 'Pawlowski and Brown are to be congratulated for having produced a book ...as comprehensive and user-friendly as this volume. Not only have they considered the effect of the authorities in a clear and logical way, but they have also highlighted problems which have yet to be resolved and questions which have yet to be answered ...one of the hallmarks of a good legal book.'

Mapping the yeast host cell response to recombinant membrane protein production:relieving the biological bottlenecks

Understanding the structures and functions of membrane proteins is an active area of research within bioscience. Membrane proteins are key players in essential cellular processes such as the uptake of nutrients, the export of waste products, and the way in which cells communicate with their environment. It is therefore not surprising that membrane proteins are targeted by over half of all prescription drugs. Since most membrane proteins are not abundant in their native membranes, it is necessary to produce them in recombinant host cells to enable further structural and functional studies. Unfortunately, achieving the required yields of functional recombinant membrane proteins is still a bottleneck in contemporary bioscience. This has highlighted the need for defined and rational optimization strategies based upon experimental observation rather than relying on trial and error. We have published a transcriptome and subsequent genetic analysis that has identified genes implicated in high-yielding yeast cells. These results have highlighted a role for alterations to a cell's protein synthetic capacity in the production of high yields of recombinant membrane protein: paradoxically, reduced protein synthesis favors higher yields. These results highlight a potential bottleneck at the protein folding or translocation stage of protein production.

Multiple roles of charged amino acids in cytoplasmic loop 7 for expression and function of the multidrug and organic anion transporter MRP1 (ABCC1)

Multidrug resistance protein MRP1 mediates the ATP-dependent efflux of many chemotherapeutic agents and organic anions. MRP1 has two nucleotide binding sites (NBSs) and three membrane spanning domains (MSDs) containing 17 transmembrane helices linked by extracellular and cytoplasmic loops (CL). Homology models suggest that CL7 (amino acids 1141-1195) is in a position where it could participate in signaling between the MSDs and NBSs during the transport process. We have individually replaced eight charged residues in CL7 with Ala, and in some cases, an amino acid with the same charge, and then investigated the effects on MRP1 expression, transport activity, and nucleotide and substrate interactions. A triple mutant in which Glu(1169), Glu(1170), and Glu(1172) were all replaced with Ala was also examined. The properties of R1173A and E1184A were comparable with those of wild-type MRP1, whereas the remaining mutants were either poorly expressed (R1166A, D1183A) or exhibited reduced transport of one or more organic anions (E1144A, D1179A, K1181A, (1169)AAQA). Same charge mutant D1183E was also not expressed, whereas expression and activity of R1166K were similar to wild-type MRP1. The moderate substrate-selective changes in transport activity displayed by mutants E1144A, D1179A, K1181A, and (1169)AAQA were accompanied by changes in orthovanadate-induced trapping of [alpha-(32)P]azidoADP by NBS2 indicating changes in ATP hydrolysis or release of ADP. In the case of E1144A, estradiol glucuronide no longer inhibited trapping of azidoADP. Together, our results demonstrate the extreme sensitivity of CL7 to mutation, consistent with its critical and complex dual role in both the proper folding and transport activity of MRP1.

Self-contained encrypted image folding

The recently introduced approach for Encrypted Image Folding is generalized to make it self-contained. The goal is achieved by enlarging the folded image so as to embed all the necessary information for the image recovery. The need for extra size is somewhat compensated by considering a transformation with higher folding capacity. Numerical examples show that the size of the resulting cipher image may be significantly smaller than the plain text one. The implementation of the approach is further extended to deal also with color images.

Sparsity and "something else": an approach to encrypted image folding

A property of sparse representations in relation to their capacity for information storage is discussed. It is shown that this feature can be used for an application that we term Encrypted Image Folding. The proposed procedure is realizable through any suitable transformation. In particular, in this paper we illustrate the approach by recourse to the Discrete Cosine Transform and a combination of redundant Cosine and Dirac dictionaries. The main advantage of the proposed technique is that both storage and encryption can be achieved simultaneously using simple processing steps.

Mapping interaction sites within the N-terminus of the calcitonin gene-related peptide receptor; the role of residues 23-60 of the calcitonin receptor-like receptor

The calcitonin receptor-like receptor (CLR) acts as a receptor for the calcitonin gene-related peptide (CGRP) but in order to recognize CGRP, it must form a complex with an accessory protein, receptor activity modifying protein 1 (RAMP1). Identifying the protein/protein and protein/ligand interfaces in this unusual complex would aid drug design. The role of the extreme N-terminus of CLR (Glu23-Ala60) was examined by an alanine scan and the results were interpreted with the help of a molecular model. The potency of CGRP at stimulating cAMP production was reduced at Leu41Ala, Gln45Ala, Cys48Ala and Tyr49Ala; furthermore, CGRP-induced receptor internalization at all of these receptors was also impaired. Ile32Ala, Gly35Ala and Thr37Ala all increased CGRP potency. CGRP specific binding was abolished at Leu41Ala, Ala44Leu, Cys48Ala and Tyr49Ala. There was significant impairment of cell surface expression of Gln45Ala, Cys48Ala and Tyr49Ala. Cys48 takes part in a highly conserved disulfide bond and is probably needed for correct folding of CLR. The model suggests that Gln45 and Tyr49 mediate their effects by interacting with RAMP1 whereas Leu41 and Ala44 are likely to be involved in binding CGRP. Ile32, Gly35 and Thr37 form a separate cluster of residues which modulate CGRP binding. The results from this study may be applicable to other family B GPCRs which can associate with RAMPs.

Hijacked then lost in translation:the plight of the recombinant host cell in membrane protein structural biology projects

Membrane protein structural biology is critically dependent upon the supply of high-quality protein. Over the last few years, the value of crystallising biochemically characterised, recombinant targets that incorporate stabilising mutations has been established. Nonetheless, obtaining sufficient yields of many recombinant membrane proteins is still a major challenge. Solutions are now emerging based on an improved understanding of recombinant host cells; as a 'cell factory' each cell is tasked with managing limited resources to simultaneously balance its own growth demands with those imposed by an expression plasmid. This review examines emerging insights into the role of translation and protein folding in defining high-yielding recombinant membrane protein production in a range of host cells.

The impact of ageing and gender on visual mental imagery processes:a study of performance on tasks from the Complete Visual Mental Imagery Battery (CVMIB)

In this study we aim to evaluate the impact of ageing and gender on different visual mental imagery processes. Two hundred and fifty-one participants (130 women and 121 men; age range = 18–77 years) were given an extensive neuropsychological battery including tasks probing the generation, maintenance, inspection, and transformation of visual mental images (Complete Visual Mental Imagery Battery, CVMIB). Our results show that all mental imagery processes with the exception of the maintenance are affected by ageing, suggesting that other deficits, such as working memory deficits, could account for this effect. However, the analysis of the transformation process, investigated in terms of mental rotation and mental folding skills, shows a steeper decline in mental rotation, suggesting that age could affect rigid transformations of objects and spare non-rigid transformations. Our study also adds to previous ones in showing gender differences favoring men across the lifespan in the transformation process, and, interestingly, it shows a steeper decline in men than in women in inspecting mental images, which could partially account for the mixed results about the effect of ageing on this specific process. We also discuss the possibility to introduce the CVMIB in clinical assessment in the context of theoretical models of mental imagery.