Kathryn A. Scott

Hidden complexity in the mechanical properties of titin

Individual molecules of the giant protein titin span the A-bands and I-bands that make up striated muscle. The I-band region of titin is responsible for passive elasticity in such muscle, and contains tandem arrays of immunoglobulin domains. One such domain (I27) has been investigated extensively, using dynamic force spectroscopy and simulation. However, the relevance of these studies to the behaviour of the protein under physiological conditions was not established. Force studies reveal a lengthening of I27 without complete unfolding, forming a stable intermediate that has been suggested to be an important component of titin elasticity. To develop a more complete picture of the forced unfolding pathway, we use mutant titins—certain mutations allow the role of the partly unfolded intermediate to be investigated in more depth. Here we show that, under physiological forces, the partly unfolded intermediate does not contribute to mechanical strength. We also propose a unified forced unfolding model of all I27 analogues studied, and conclude that I27 can withstand higher forces in muscle than was predicted previously.

Coarse-grained MD simulations of membrane protein-bilayer self-assembly

Complete determination of a membrane protein structure requires knowledge of the protein position within the lipid bilayer. As the number of determined structures of membrane proteins increases so does the need for computational methods which predict their position in the lipid bilayer. Here we present a coarse-grained molecular dynamics approach to lipid bilayer self-assembly around membrane proteins. We demonstrate that this method can be used to predict accurately the protein position in the bilayer for membrane proteins with a range of different sizes and architectures.

Coarse-grained simulation: a high-throughput computational approach to membrane proteins

An understanding of the interactions of membrane proteins with a lipid bilayer environment is central to relating their structure to their function and stability. A high-throughput approach to prediction of membrane protein interactions with a lipid bilayer based on coarse-grained Molecular Dynamics simulations is described. This method has been used to develop a database of CG simulations (coarse-grained simulations) of membrane proteins (http://sbcb.bioch.ox.ac.uk/cgdb). Comparison of CG simulations and AT simulations (atomistic simulations) of lactose permease reveals good agreement between the two methods in terms of predicted lipid headgroup contacts. Both CG and AT simulations predict considerable local bilayer deformation by the voltage sensor domain of the potassium channel KvAP.

Dynameomics: mass annotation of protein dynamics and unfolding in water by high-throughput atomistic molecular dynamics simulations

The goal of Dynameomics is to perform atomistic molecular dynamics (MD) simulations of representative proteins from all known folds in explicit water in their native state and along their thermal unfolding pathways. Here we present 188-fold representatives and their native state simulations and analyses. These 188 targets represent 67% of all the structures in the Protein Data Bank. The behavior of several specific targets is highlighted to illustrate general properties in the full dataset and to demonstrate the role of MD in understanding protein function and stability. As an example of what can be learned from mining the Dynameomics database, we identified a protein fold with heightened localized dynamics. In one member of this fold family, the motion affects the exposure of its phosphorylation site and acts as an entropy sink to offset another portion of the protein that is relatively immobile in order to present a consistent interface for protein docking. In another member of this family, a polymorphism in the highly mobile region leads to a host of disease phenotypes. We have constructed a web site to provide access to a novel hybrid relational/multidimensional database (described in the succeeding two papers) to view and interrogate simulations of the top 30 targets: http://www.dynameomics.org. The Dynameomics database, currently the largest collection of protein simulations and protein structures in the world, should also be useful for determining the rules governing protein folding and kinetic stability, which should aid in deciphering genomic information and for protein engineering and design.

Gating at Both Ends and Breathing in the Middle: Conformational Dynamics of TolC

Drug extrusion via efflux through a tripartite complex (an inner membrane pump, an outer membrane protein, and a periplasmic protein) is a widely used mechanism in Gram-negative bacteria. The outer membrane protein (TolC in Escherichia coli; OprM in Pseudomonas aeruginosa) forms a tunnel-like pore through the periplasmic space and the outer membrane. Molecular dynamics simulations of TolC have been performed, and are compared to simulations of Y362F/R367S mutant, and to simulations of its homolog OprM. The results reveal a complex pattern of conformation dynamics in the TolC protein. Two putative gate regions, located at either end of the protein, can be distinguished. These regions are the extracellular loops and the mouth of the periplasmic domain, respectively. The periplasmic gate has been implicated in the conformational changes leading from the closed x-ray structure to a proposed open state of TolC. Between the two gates, a peristaltic motion of the periplasmic domain is observed, which may facilitate transport of the solutes from one end of the tunnel to the other. The motions observed in the atomistic simulations are also seen in coarse-grained simulations in which the protein tertiary structure is represented by an elastic network model.

Spectrin R16: Broad energy barrier or sequential transition states?

A number of models have been proposed to account for nonlinearity in the relation between observed rate constants for folding and/or unfolding and denaturant concentration. Where curvature is seen principally in the arm of a chevron plot, three explanations are proposed: a change in the ground state at increasing concentration of urea, movement of the transition state along a broad energy barrier, and a switch between two sequential transition states separated by an on‐pathway high‐energy intermediate. Here we demonstrate that the latter two models in particular can be used to describe the data for the all‐α protein spectrin R16. Further, whatever the method of analysis, the pattern of Φ‐values seen is robust; thus we would draw the same conclusions from our data set independently of the method used for analysis. While this is not a novel observation, this is the first systematic study where a comparison has been made between Φ‐values calculated using the broad and sequential transition state models.

CGDB: A database of membrane protein/lipid interactions by coarse-grained molecular dynamics simulations

Membrane protein function and stability has been shown to be dependent on the lipid environment. Recently, we developed a high-throughput computational approach for the prediction of membrane protein/lipid interactions. In the current study, we enhanced this approach with the addition of a new measure of the distortion caused by membrane proteins on a lipid bilayer. This is illustrated by considering the effect of lipid tail length and headgroup charge on the distortion caused by the integral membrane proteins MscS and FLAP, and by the voltage sensing domain from the channel KvAP. Changing the chain length of lipids alters the extent but not the pattern of distortion caused by MscS and FLAP; lipid headgroups distort in order to interact with very similar but not identical regions in these proteins for all bilayer widths investigated. Introducing anionic lipids into a DPPC bilayer containing the KvAP voltage sensor does not affect the extent of bilayer distortion.

Quinol oxidation by c-type cytochromes: structural characterization of the menaquinol binding site of NrfHA.

Membrane-bound cytochrome c quinol dehydrogenases play a crucial role in bacterial respiration by oxidizing menaquinol and transferring electrons to various periplasmic oxidoreductases. In this work, the menaquinol oxidation site of NrfH was characterized by the determination of the X-ray structure of Desulfovibrio vulgaris NrfHA nitrite reductase complex bound to 2-heptyl-4-hydroxyquinoline-N-oxide, which is shown to act as a competitive inhibitor of NrfH quinol oxidation activity. The structure, at 2.8-A resolution, reveals that the inhibitor binds close to NrfH heme 1, where it establishes polar contacts with two essential residues: Asp89, the residue occupying the heme distal ligand position, and Lys82, a strictly conserved residue. The menaquinol binding cavity is largely polar and has a wide opening to the protein surface. Coarse-grained molecular dynamics simulations suggest that the quinol binding site of NrfH and several other respiratory enzymes lie in the head group region of the membrane, which probably facilitates proton transfer to the periplasm. Although NrfH is not a multi-span membrane protein, its quinol binding site has several characteristics similar to those of quinone binding sites previously described. The data presented here provide the first characterization of the quinol binding site of the cytochrome c quinol dehydrogenase family.

Dynameomics: A Consensus View of the Protein Unfolding/Folding Transition State Ensemble across a Diverse Set of Protein Folds

The Dynameomics project aims to simulate a representative sample of all globular protein metafolds under both native and unfolding conditions. We have identified protein unfolding transition state (TS) ensembles from multiple molecular dynamics simulations of high-temperature unfolding in 183 structurally distinct proteins. These data can be used to study individual proteins and individual protein metafolds and to mine for TS structural features common across all proteins. Separating the TS structures into four different fold classes (all proteins, all-alpha, all-beta, and mixed alpha/beta and alpha +beta) resulted in no significant difference in the overall protein properties. The residues with the most contacts in the native state lost the most contacts in the TS ensemble. On average, residues beginning in an alpha-helix maintained more structure in the TS ensemble than did residues starting in beta-strands or any other conformation. The metafolds studied here represent 67% of all known protein structures, and this is, to our knowledge, the largest, most comprehensive study of the protein folding/unfolding TS ensemble to date. One might have expected broad distributions in the average global properties of the TS relative to the native state, indicating variability in the amount of structure present in the TS. Instead, the average global properties converged with low standard deviations across metafolds, suggesting that there are general rules governing the structure and properties of the TS.

Specificity of localization and phosphotransfer in the CheA proteins of Rhodobacter sphaeroides

Specificity of protein–protein interactions plays a vital role in signal transduction. The chemosensory pathway of Rhodobacter sphaeroides comprises multiple homologues of chemotaxis proteins characterized in organisms such as Escherichia coli. Three CheA homologues are essential for chemotaxis in R. sphaeroides under laboratory conditions. These CheAs are differentially localized to two chemosensory clusters, one at the cell pole and one in the cytoplasm. The polar CheA, CheA2, has the same domain structure as E. coli CheA and can phosphorylate all R. sphaeroides chemotaxis response regulators. CheA3 and CheA4 independently localize to the cytoplasmic cluster; each protein has a subset of the CheA domains, with CheA3 phosphorylating CheA4 together making a functional CheA protein. Interestingly, CheA3‐P can only phosphorylate two response regulators, CheY6 and CheB2. R. sphaeroides CheAs exhibit two interesting differences in specificity: (i) the response regulators that they phosphorylate and (ii) the chemosensory cluster to which they localize. Using a domain‐swapping approach we investigated the role of the P1 and P5 CheA domains in determining these specificities. We show that the P1 domain is sufficient to determine which response regulators will be phosphorylated in vitro while the P5 domain is sufficient to localize the CheAs to a specific chemosensory cluster.