Dynamics of ATP-binding cassette contribute to allosteric control, nucleotide binding and energy transduction in ABC transporters

Combining theoretical and experimental data to decipher CFTR 3D structures and functions

Cryo-electron microscopy (cryo-EM) has recently provided invaluable experimental data about the full-length cystic fibrosis transmembrane conductance regulator (CFTR) 3D structure. However, this experimental information deals with inactive states of the channel, either in an apo, quiescent conformation, in which nucleotide-binding domains (NBDs) are widely separated or in an ATP-bound, yet closed conformation. Here, we show that 3D structure models of the open and closed forms of the channel, now further supported by metadynamics simulations and by comparison with the cryo-EM data, could be used to gain some insights into critical features of the conformational transition toward active CFTR forms. These critical elements lie within membrane-spanning domains but also within NBD1 and the N-terminal extension, in which conformational plasticity is predicted to occur to help the interaction with filamin, one of the CFTR cellular partners.

Solution NMR studies of periplasmic binding proteins and their interaction partners

Periplasmic binding proteins (PBPs) are a crucial part of ATP-binding cassette import systems in Gram-negative bacteria. Central to their function is the ability to undergo a large-scale conformational rearrangement from open-unliganded to closed-liganded, which signals the presence of substrate and starts its translocation. Over the years, PBPs have been extensively studied not only owing to their essential role in nutrient uptake but also because they serve as excellent models for both practical applications (e.g., biosensor technology) and basic research (e.g., allosteric mechanisms). Although much of our knowledge at atomic level has been inferred from the detailed, static pictures afforded by crystallographic studies, nuclear magnetic resonance (NMR) has been able to fill certain gaps in such body of work, particularly with regard to dynamic processes. Here, we review NMR studies on PBPs, and their unique insights on conformation, dynamics, energetics, substrate binding, and interactions with related transport proteins. Based on the analysis of recent paramagnetic NMR results, as well as crystallographic and functional observations, we propose a mechanism that could explain the ability of certain PBPs to achieve a closed conformation in absence of ligand while others seem to remain open until ligand-mediated closure.

Robust patterning of gene expression based on internal coordinate system of cells

Cell-to-cell communication in multicellular organisms is established through the transmission of various kinds of chemical substances such as proteins. It is well known that gene expression triggered by a chemical substance in individuals has stable spatial patterns despite the individual differences in concentration patterns of the chemical substance. This fact reveals an important property of multicellular organisms called "robustness", which allows the organisms to generate their forms while maintaining proportion. Robustness has been conventionally accounted for by the stability of solutions of dynamical equations that represent a specific interaction network of chemical substances. However, any biological system is composed of autonomous elements. In general, an autonomous element does not merely accept information on the chemical substance from the environment; instead, it accepts the information based on its own criteria for reaction. Therefore, this phenomenon needs to be considered from the viewpoint of cells. Such a viewpoint is expected to allow the consideration of the autonomy of cells in multicellular organisms. This study aims to explain theoretically the robust patterning of gene expression from the viewpoint of cells. For this purpose, we introduced a new operator for transforming a state variable of a chemical substance from an external coordinate system to an internal coordinate system of each cell, which describes the observation of the chemical substance by cells. We then applied this operator to the simplest reaction-diffusion model of the chemical substance to investigate observation effects by cells. Our mathematical analysis of this extended model indicates that the robust patterning of gene expression against individual differences in concentration pattern of the chemical substance can be explained from the viewpoint of cells if there is a regulation field that compensates for the difference between cells seen in the observation results. This result provides a new insight into the investigation of the mechanism of robust patterning in biological systems composed of individual elements.

Successful development and use of a thermodynamic stability screen for optimizing the yield of nucleotide binding domains

ATP sensitive potassium (KATP) channels consist of four copies of a pore-forming inward rectifying potassium channel (Kir6.1 or Kir6.2) and four copies of a sulfonylurea receptor (SUR1, SUR2A, or SUR2B). SUR proteins are members of the ATP-binding cassette superfamily of proteins. Binding of ATP to the Kir6.x subunit mediates channel inhibition, whereas MgATP binding and hydrolysis at the SUR NBDs results in channel opening. Mutations in SUR1 and SUR2A NBDs cause diseases of insulin secretion and cardiac disorders, respectively, underlying the importance of studying the NBDs. Although purification of SUR2A NBD1 in a soluble form is possible, the lack of long-term sample stability of the protein in a concentrated form has precluded detailed studies of the protein aimed at gaining a molecular-level understanding of how SUR mutations cause disease. Here we use a convenient and cost-effective thermodynamic screening method to probe stabilizing conditions for SUR2A NBD1. Results from the screen are used to alter the purification protocol to allow for significantly increased yields of the purified protein. In addition, the screen provides strategies for long-term storage of NBD1 and generating NBD1 samples at high concentrations suitable for NMR studies. NMR spectra of NBD1 with MgAMP-PNP are of higher quality compared to using MgATP, indicating that MgAMP-PNP be used as the ligand in future NMR studies. The screen presented here can be expanded to using different additives and can be employed to enhance purification yields, sample life times, and storage of other low stability nucleotide binding domains, such as GTPases.

Chemical exchange in biomacromolecules: past, present, and future

The perspective reviews quantitative investigations of chemical exchange phenomena in proteins and other biological macromolecules using NMR spectroscopy, particularly relaxation dispersion methods. The emphasis is on techniques and applications that quantify the populations, interconversion kinetics, and structural features of sparsely populated conformational states in equilibrium with a highly populated ground state. Applications to folding, molecular recognition, catalysis, and allostery by proteins and nucleic acids are highlighted.

ATP binding to two sites is necessary for dimerization of nucleotide-binding domains of ABC proteins

ATP binding cassette (ABC) transporters have a functional unit formed by two transmembrane domains and two nucleotide binding domains (NBDs). ATP-bound NBDs dimerize in a head-to-tail arrangement, with two nucleotides sandwiched at the dimer interface. Both NBDs contribute residues to each of the two nucleotide-binding sites (NBSs) in the dimer. In previous studies, we showed that the prototypical NBD MJ0796 from Methanocaldococcus jannaschii forms ATP-bound dimers that dissociate completely following hydrolysis of one of the two bound ATP molecules. Since hydrolysis of ATP at one NBS is sufficient to drive dimer dissociation, it is unclear why all ABC proteins contain two NBSs. Here, we used luminescence resonance energy transfer (LRET) to study ATP-induced formation of NBD homodimers containing two NBSs competent for ATP binding, and NBD heterodimers with one active NBS and one binding-defective NBS. The results showed that binding of two ATP molecules is necessary for NBD dimerization. We conclude that ATP hydrolysis at one nucleotide-binding site drives NBD dissociation, but two binding sites are required to form the ATP-sandwich NBD dimer necessary for hydrolysis.

Backbone resonance assignments for a homolog of the essential ribosome biogenesis factor Fap7 from P. horikoshii in its nucleotide-free and -bound forms

The protein "factor activating Pos9 (Skn7)", Fap7, is an essential protein in yeast and plays an important role in the biogenesis of the small ribosomal subunit. In eukaryotes, the final processing step of the small ribosomal subunit RNA is the endonucleolytic cleavage of 20S pre-rRNA at cleavage site D yielding mature 18S rRNA. Depletion of Fap7 in yeast leads to a dramatic accumulation of 20S pre-rRNA and a concomitant decrease in 18S rRNA in the cytoplasm. In addition, these cells contain higher levels of 60S, but decreased numbers of 40S ribosomal subunits. Fap7 contains a P-loop like motif placing it in a class with NTPases and kinases and a role for it as an adenylate kinase has been suggested. Up to now both the structure of Fap7 and its detailed function during ribosome biogenesis remain elusive. Here, we present the backbone NMR assignments of a Fap7 homolog from the thermophilic archaeon Pyrococcus horikoshii in its nucleotide free form and bound to the adenylate kinase inhibitor AP5A.

Allosteric coupling between the intracellular coupling helix 4 and regulatory sites of the first nucleotide-binding domain of CFTR

Cystic fibrosis is caused by mutations in CFTR (cystic fibrosis transmembrane conductance regulator), leading to folding and processing defects and to chloride channel gating misfunction. CFTR is regulated by ATP binding to its cytoplasmic nucleotide-binding domains, NBD1 and NBD2, and by phosphorylation of the NBD1 regulatory insert (RI) and the regulatory extension (RE)/R region. These regulatory effects are transmitted to the rest of the channel via NBD interactions with intracellular domain coupling helices (CL), particularly CL4. Using a sensitive method for detecting inter-residue correlations between chemical shift changes in NMR spectra, an allosteric network was revealed within NBD1, with a construct lacking RI. The CL4-binding site couples to the RI-deletion site and the C-terminal residues of NBD1 that precede the R region in full-length CFTR. Titration of CL4 peptide into NBD1 perturbs the conformational ensemble in these sites with similar titration patterns observed in F508del, the major CF-causing mutant, and in suppressor mutants F494N, V510D and Q637R NBD1, as well as in a CL4-NBD1 fusion construct. Reciprocally, the C-terminal mutation, Q637R, perturbs dynamics in these three sites. This allosteric network suggests a mechanism synthesizing diverse regulatory NBD1 interactions and provides biophysical evidence for the allosteric coupling required for CFTR function.

Backbone NMR resonance assignments of the nucleotide binding domain of the ABC multidrug transporter LmrA from Lactococcus lactis in its ADP-bound state

LmrA from Lactococcus lactis is a multidrug transporter and a member of the ATP binding cassette (ABC) transporter family. ABC transporters consist of a transmembrane domain (TMD) and a nucleotide binding domain (NBD). The NBD contains the highly conserved signature motifs of this transporter superfamily. In the case of LmrA, the TMD and the NBD are expressed as a single polypeptide. LmrA catalyzes the extrusion of hydrophobic compounds including antibiotics from the cell membrane at the expense of ATP hydrolysis. ATP binds to the NBD, where binding and hydrolysis induce conformational changes that lead to the extrusion of the substrate via the TMD. Here, we report the (1)H, (13)C and (15)N backbone chemical shift assignments of the isolated 263 amino acid containing NBD of LmrA in its ADP bound state.

Kinetics of the association/dissociation cycle of an ATP-binding cassette nucleotide-binding domain

Most ATP binding cassette (ABC) proteins are pumps that transport substrates across biological membranes using the energy of ATP hydrolysis. Functional ABC proteins have two nucleotide-binding domains (NBDs) that bind and hydrolyze ATP, but the molecular mechanism of nucleotide hydrolysis is unresolved. This is due in part to the limited kinetic information on NBD association and dissociation. Here, we show dimerization of a catalytically active NBD and follow in real time the association and dissociation of NBDs from the changes in fluorescence emission of a tryptophan strategically located at the center of the dimer interface. Spectroscopic and structural studies demonstrated that the tryptophan can be used as dimerization probe, and we showed that under hydrolysis conditions (millimolar MgATP), not only the dimer dissociation rate increases, but also the dimerization rate. Neither dimer formation or dissociation are clearly favored, and the end result is a dynamic equilibrium where the concentrations of monomer and dimer are very similar. We proposed that based on their variable rates of hydrolysis, the rate-limiting step of the hydrolysis cycle may differ among full-length ABC proteins.

Conformational changes induced by ATP-hydrolysis in an ABC transporter: a molecular dynamics study of the Sav1866 exporter

ATP-Binding Cassette (ABC) transporters are ubiquitous membrane proteins that use energy from ATP binding or/and hydrolysis to actively transport allocrites across membranes. In this study, we identify ATP-hydrolysis induced conformational changes in a complete ABC exporter (Sav1866) from Staphylococcus aureaus, using molecular dynamics (MD) simulations. By performing MD simulations on the ATP and ADP+IP bound states, we identify the conformational consequences of hydrolysis, showing that the major rearrangements are not restricted to the NBDs, but extend to the transmembrane domains (TMDs) external regions. For the first time, to our knowledge, we see, within the context of a complete transporter, NBD dimer opening in the ADP+IP state in contrast with all ATP-bound states. This opening results from the dissociation of the ABC signature motif from the nucleotide. In addition, in both states, we observe the opening of a gate entrance in the intracellular loop region leading to the exposure of the TMDs internal cavity to the cytoplasm. To see if this opening was large enough to allow allocrite transport, the adiabatic energy profile for doxorubicin passage was determined. For both states, this profile, although an approximation, is overall downhill from the cytoplasmatic to the extracellular side, and the local energy barriers along the TMDs are relatively small, evidencing the exporter nature of Sav1866. The major difference between states is an energy barrier located in the cytoplasmic gate region, which becomes reduced upon hydrolysis, suggesting that allocrite passage is facilitated, and evidencing a possible molecular mechanism for the active transport in these proteins.

NMR spectroscopy to study the dynamics and interactions of CFTR

The cystic fibrosis transmembrane conductance regulator (CFTR) is a multi-domain membrane chloride channel whose activity is regulated by ATP at two nucleotide-binding domains (NBD1 and NBD2) and by phosphorylation of the regulatory (R) region. The NBDs and the R region have functionally relevant motions that are critical for channel gating. Nuclear magnetic resonance (NMR) spectroscopy is a highly useful technique for obtaining information on the structure and interactions of CFTR and is extremely powerful for probing dynamics. NMR approaches for studying CFTR are reviewed, using our previous NBD1 and the R region results to provide examples. These NMR data are yielding insights into the dynamic properties and interactions that facilitate normal CFTR regulation as well as pathological effects of mutations, including the most common disease mutant, deletion of F508 in NBD1.

Arginine kinase: joint crystallographic and NMR RDC analyses link substrate-associated motions to intrinsic flexibility

The phosphagen kinase family, including creatine and arginine kinases (AKs), catalyzes the reversible transfer of a "high-energy" phosphate between ATP and a phosphoguanidino substrate. They have become a model for the study of both substrate-induced conformational change and intrinsic protein dynamics. Prior crystallographic studies indicated large substrate-induced domain rotations, but differences among a recent set of AK structures were interpreted as a plastic deformation. Here, the structure of Limulus substrate-free AK is refined against high-resolution crystallographic data and compared quantitatively with NMR chemical shifts and residual dipolar couplings (RDCs). This demonstrates the feasibility of this type of RDC analysis of proteins that are large by NMR standards (42 kDa) and illuminates the solution structure, free from crystal-packing constraints. Detailed comparison of the 1.7 Å resolution substrate-free crystal structure against the 1.7 Å transition-state analog complex shows large substrate-induced domain motions that can be broken down into movements of smaller quasi-rigid bodies. The solution-state structure of substrate-free AK is most consistent with an equilibrium of substrate-free and substrate-bound structures, with the substrate-free form dominating, but with varying displacements of the quasi-rigid groups. Rigid-group rotations evident from the crystal structures are about axes previously associated with intrinsic millisecond dynamics using NMR relaxation dispersion. Thus, "substrate-induced" motions are along modes that are intrinsically flexible in the substrate-free enzyme and likely involve some degree of conformational selection.

The Calcium-Dependent Interaction of S100B with Its Protein Targets

S100B is a calcium signaling protein that is a member of the S100 protein family. An important feature of S100B and most other S100 proteins (S100s) is that they often bind Ca²⁺ ions relatively weakly in the absence of a protein target; upon binding their target proteins, Ca²⁺-binding then increases by as much as from 200- to 400-fold. This manuscript reviews the structural basis and physiological significance of increased Ca²⁺-binding affinity in the presence of protein targets. New information regarding redundancy among family members and the structural domains that mediate the interaction of S100B, and other S100s, with their targets is also presented. It is the diversity among individual S100s, the protein targets that they interact with, and the Ca²⁺ dependency of these protein-protein interactions that allow S100s to transduce changes in [Ca²⁺]_{intracellular} levels into spatially and temporally unique biological responses.

Insights into the molecular mechanism of an ABC transporter: conformational changes in the NBD dimer of MJ0796

Despite the rapid advances in the study of ABC transporters, many fundamental questions linked to ATP binding/hydrolysis and its relation to the transport cycle remain unanswered. In particular, it is still neither clear nor consensual how the ATP energy is used by the nucleotide binding domains (NBDs) to produce mechanical work and drive the substrate translocation. The major conformational changes in the NBDs following ATP hydrolysis during the transport cycle and the role played by the conserved family motifs in harnessing the energy associated with nucleotide hydrolysis are yet unknown. Additionally, the way energy is transmitted from the catalytic to the membrane domains, in order to drive substrate translocation, is also a fundamental question that remains unanswered. Due to the high structure similarities of the NBD architecture throughout the whole ABC family, it is likely that the mechanism of ATP binding, hydrolysis, and communication with the transmembrane domains is similar in all family members, independently of the nature of the transported substrate. In this work, we focused our attention on the consequences of ATP hydrolysis in the NBDs, especially on the structural changes that occur during this process. For that, we use molecular dynamics simulation techniques taking as a starting point the X-ray structure of the MJ0796 dimer from Methanococcus jannaschii. Several potential intermediate states of the ATP hydrolytic cycle are investigated, each consisting of different combinations of nucleotide-bound forms. The results obtained allowed us to identify the conformational rearrangements induced by hydrolysis on the catalytic subunits, as well as the residues involved in this reorganization. The major changes are localized at specific regions of the protein, namely, involving segments 11-19 and 93-124. Additionally, our results together with the knowledge of complete ABC transporter X-ray structures suggest a possible NBD:TMD signal transmission interface.

Structure and dynamics of a processive Brownian motor: the translating ribosome

There is mounting evidence indicating that protein synthesis is driven and regulated by mechanisms that direct stochastic, large-scale conformational fluctuations of the translational apparatus. This mechanistic paradigm implies that a free-energy landscape governs the conformational states that are accessible to and sampled by the translating ribosome. This scenario presents interdependent opportunities and challenges for structural and dynamic studies of protein synthesis. Indeed, the synergism between cryogenic electron microscopic and X-ray crystallographic structural studies, on the one hand, and single-molecule fluorescence resonance energy transfer (smFRET) dynamic studies, on the other, is emerging as a powerful means for investigating the complex free-energy landscape of the translating ribosome and uncovering the mechanisms that direct the stochastic conformational fluctuations of the translational machinery. In this review, we highlight the principal insights obtained from cryogenic electron microscopic, X-ray crystallographic, and smFRET studies of the elongation stage of protein synthesis and outline the emerging themes, questions, and challenges that lie ahead in mechanistic studies of translation.

Intrinsic domain and loop dynamics commensurate with catalytic turnover in an induced-fit enzyme

Arginine kinase catalyzes reversible phosphoryl transfer between ATP and arginine, buffering cellular ATP concentrations. Structures of substrate-free and -bound enzyme have highlighted a range of conformational changes thought to occur during the catalytic cycle. Here, NMR is used to characterize the intrinsic backbone dynamics over multiple timescales. Relaxation dispersion indicates rigid-body motion of the N-terminal domain and flexible dynamics in the I182-G209 loop, both at millisecond rates commensurate with k(cat), implying that either might be rate limiting upon catalysis. Lipari-Szabo analysis indicates backbone flexibility on the nanosecond timescale in the V308-V322 loop, while the rest of the enzyme is more rigid in this timescale. Thus, intrinsic dynamics are most prominent in regions that have been independently implicated in conformational changes. Substrate-free enzyme may sample an ensemble of different conformations, of which a subset is selected upon substrate binding, with critical active site residues appropriately configured for binding and catalysis.

Millisecond timescale fluctuations in dihydrofolate reductase are exquisitely sensitive to the bound ligands

Enzyme catalysis can be described as progress over a multi-dimensional energy landscape where ensembles of interconverting conformational substates channel the enzyme through its catalytic cycle. We applied NMR relaxation dispersion to investigate the role of bound ligands in modulating the dynamics and energy landscape of Escherichia coli dihydrofolate reductase to obtain insights into the mechanism by which the enzyme efficiently samples functional conformations as it traverses its reaction pathway. Although the structural differences between the occluded substrate binary complexes and product ternary complexes are very small, there are substantial differences in protein dynamics. Backbone fluctuations on the micros-ms timescale in the cofactor binding cleft are similar for the substrate and product binary complexes, but fluctuations on this timescale in the active site loops are observed only for complexes with substrate or substrate analog and are not observed for the binary product complex. The dynamics in the substrate and product binary complexes are governed by quite different kinetic and thermodynamic parameters. Analogous dynamic differences in the E:THF:NADPH and E:THF:NADP(+) product ternary complexes are difficult to rationalize from ground-state structures. For both of these complexes, the nicotinamide ring resides outside the active site pocket in the ground state. However, they differ in the structure, energetics, and dynamics of accessible higher energy substates where the nicotinamide ring transiently occupies the active site. Overall, our results suggest that dynamics in dihydrofolate reductase are exquisitely "tuned" for every intermediate in the catalytic cycle; structural fluctuations efficiently channel the enzyme through functionally relevant conformational space.

Structure and dynamics of NBD1 from CFTR characterized using crystallography and hydrogen/deuterium exchange mass spectrometry

The DeltaF508 mutation in nucleotide-binding domain 1 (NBD1) of the cystic fibrosis transmembrane conductance regulator (CFTR) is the predominant cause of cystic fibrosis. Previous biophysical studies on human F508 and DeltaF508 domains showed only local structural changes restricted to residues 509-511 and only minor differences in folding rate and stability. These results were remarkable because DeltaF508 was widely assumed to perturb domain folding based on the fact that it prevents trafficking of CFTR out of the endoplasmic reticulum. However, the previously reported crystal structures did not come from matched F508 and DeltaF508 constructs, and the DeltaF508 structure contained additional mutations that were required to obtain sufficient protein solubility. In this article, we present additional biophysical studies of NBD1 designed to address these ambiguities. Mass spectral measurements of backbone amide (1)H/(2)H exchange rates in matched F508 and DeltaF508 constructs reveal that DeltaF508 increases backbone dynamics at residues 509-511 and the adjacent protein segments but not elsewhere in NBD1. These measurements also confirm a high level of flexibility in the protein segments exhibiting variable conformations in the crystal structures. We additionally present crystal structures of a broader set of human NBD1 constructs, including one harboring the native F508 residue and others harboring the DeltaF508 mutation in the presence of fewer and different solubilizing mutations. The only consistent conformational difference is observed at residues 509-511. The side chain of residue V510 in this loop is mostly buried in all non-DeltaF508 structures but completely solvent exposed in all DeltaF508 structures. These results reinforce the importance of the perturbation DeltaF508 causes in the surface topography of NBD1 in a region likely to mediate contact with the transmembrane domains of CFTR. However, they also suggest that increased exposure of the 509-511 loop and increased dynamics in its vicinity could promote aggregation in vitro and aberrant intermolecular interactions that impede trafficking in vivo.

NMR evidence for differential phosphorylation-dependent interactions in WT and DeltaF508 CFTR

The most common cystic fibrosis (CF)-causing mutation in the cystic fibrosis transmembrane conductance regulator (CFTR) is deletion of Phe508 (DeltaF508) in the first of two nucleotide-binding domains (NBDs). Nucleotide binding and hydrolysis at the NBDs and phosphorylation of the regulatory (R) region are required for gating of CFTR chloride channel activity. We report NMR studies of wild-type and DeltaF508 murine CFTR NBD1 with the C-terminal regulatory extension (RE), which contains residues of the R region. Interactions of the wild-type NBD1 core with the phosphoregulatory regions, the regulatory insertion (RI) and RE, are disrupted upon phosphorylation, exposing a potential binding site for the first coupling helix of the N-terminal intracellular domain (ICD). Phosphorylation of DeltaF508 NBD1 does not as effectively disrupt interactions with the phosphoregulatory regions, which, along with other structural differences, leads to decreased binding of the first coupling helix. These results provide a structural basis by which phosphorylation of CFTR may affect the channel gating of full-length CFTR and expand our understanding of the molecular basis of the DeltaF508 defect.

NMR and EPR studies of membrane transporters

In order to fulfill their function, membrane transport proteins have to cycle through a number of conformational and/or energetic states. Thus, understanding the role of conformational dynamics seems to be the key for elucidation of the functional mechanism of these proteins. However, membrane proteins in general are often difficult to express heterologously and in sufficient amounts for structural studies. It is especially challenging to trap a stable energy minimum, e.g., for crystallographic analysis. Furthermore, crystallization is often only possible by subjecting the protein to conditions that do not resemble its native environment and crystals can only be snapshots of selected conformational states. Nuclear magnetic resonance (NMR) and electron paramagnetic resonance (EPR) spectroscopy are complementary methods that offer unique possibilities for studying membrane proteins in their natural membrane environment and for investigating functional conformational changes, lipid interactions, substrate-lipid and substrate-protein interactions, oligomerization states and overall dynamics of membrane transporters. Here, we review recent progress in the field including studies from primary and secondary active transporters.

Solution structure of S100A1 bound to the CapZ peptide (TRTK12)

As is typical for S100-target protein interactions, a Ca 2+-dependent conformational change in S100A1 is required to bind to a 12-residue peptide (TRTK12) derived from the actin-capping protein CapZ. In addition, the Ca 2+-binding affinity of S100A1 is found to be tightened (greater than threefold) when TRTK12 is bound. To examine the biophysical basis for these observations, we determined the solution NMR structure of TRTK12 in a complex with Ca 2+-loaded S100A1. When bound to S100A1, TRTK12 forms an amphipathic helix (residues N6 to S12) with several favorable hydrophobic interactions observed between W7, I10, and L11 of the peptide and a well-defined hydrophobic binding pocket in S100A1 that is only present in the Ca 2+-bound state. Next, the structure of S100A1-TRTK12 was compared to that of another S100A1-target complex (i.e., S100A1-RyRP12), which illustrated how the binding pocket in Ca 2+-S100A1 can accommodate peptide targets with varying amino acid sequences. Similarities and differences were observed when the structures of S100A1-TRTK12 and S100B-TRTK12 were compared, providing insights regarding how more than one S100 protein can interact with the same peptide target. Such comparisons, including those with other S100-target and S100-drug complexes, provide the basis for designing novel small-molecule inhibitors that could be specific for blocking one or more S100-target protein interactions.

Selected-fit versus induced-fit protein binding: kinetic differences and mutational analysis

The binding of a ligand molecule to a protein is often accompanied by conformational changes of the protein. A central question is whether the ligand induces the conformational change (induced-fit), or rather selects and stabilizes a complementary conformation from a pre-existing equilibrium of ground and excited states of the protein (selected-fit). We consider here the binding kinetics in a simple four-state model of ligand-protein binding. In this model, the protein has two conformations, which can both bind the ligand. The first conformation is the ground state of the protein when the ligand is off, and the second conformation is the ground state when the ligand is bound. The induced-fit mechanism corresponds to ligand binding in the unbound ground state, and the selected-fit mechanism to ligand binding in the excited state. We find a simple, characteristic difference between the on- and off-rates in the two mechanisms if the conformational relaxation into the ground states is fast. In the case of selected-fit binding, the on-rate depends on the conformational equilibrium constant, whereas the off-rate is independent. In the case of induced-fit binding, in contrast, the off-rate depends on the conformational equilibrium, while the on-rate is independent. Whether a protein binds a ligand via selected-fit or induced-fit thus may be revealed by mutations far from the protein's binding pocket, or other "perturbations" that only affect the conformational equilibrium. In the case of selected-fit, such mutations will only change the on-rate, and in the case of induced-fit, only the off-rate.

Dynamic properties of a type II cadherin adhesive domain: implications for the mechanism of strand-swapping of classical cadherins

Cadherin-mediated cell adhesion is achieved through dimerization of cadherin N-terminal extracellular (EC1) domains presented from apposed cells. The dimer state is formed by exchange of N-terminal beta strands and insertion of conserved tryptophan indole side chains from one monomer into hydrophobic acceptor pockets of the partner molecule. The present work characterizes individual monomer and dimer states and the monomer-dimer equilibrium of the mouse Type II cadherin-8 EC1 domain using NMR spectroscopy. Limited picosecond-to-nanosecond timescale dynamics of the tryptophan indole moieties for both monomer and dimer states are consistent with well-ordered packing of the N-terminal beta strands intramolecularly and intermolecularly, respectively. However, pronounced microsecond-to-millisecond timescale dynamics of the side chains are observed for the monomer but not the dimer state, suggesting that monomers transiently sample configurations in which the indole moieties are exposed. The results suggest possible kinetic mechanisms for EC1 dimerization.

The intact CFTR protein mediates ATPase rather than adenylate kinase activity

The two NBDs (nucleotide-binding domains) of ABC (ATP-binding-cassette) proteins function in a complex to mediate ATPase activity and this activity has been linked to their regulated transport activity. A similar model has been proposed for CFTR (cystic fibrosis transmembrane conductance regulator), the chloride channel defective in cystic fibrosis, wherein ATP binding and hydrolysis regulate the channel gate. Recently, it was shown that the individual NBDs isolated from CFTR primarily mediate adenylate kinase activity, raising the possibility that this activity may also contribute to gating of the CFTR channel. However, this present study shows that whereas the isolated NBDs exhibit adenylate kinase activity, the full-length purified and reconstituted CFTR protein functions as an ATPase, arguing that the enzymatic activity of the NBDs is dependent on their molecular context and appropriate domain-domain assembly. As expected, the disease-causing mutant bearing a mutation in the ABC signature motif, CFTR-G551D, exhibited a markedly reduced ATPase activity. Furthermore, mutation of the putative catalytic base in CFTR caused a reduction in ATPase activity, with the CFTR-E1371Q mutant supporting a low level of residual activity. Neither of these mutants exhibited detectable adenylate kinase activity. Together, these findings support the concept that the molecular mechanism of action of CFTR is dependent on ATP binding and hydrolysis, and that the structure of prokaryotic ABC ATPases provide a useful template for understanding their mechanism of action.

Context shapes: Efficient complementary shape matching for protein-protein docking

We describe an efficient method for partial complementary shape matching for use in rigid protein-protein docking. The local shape features of a protein are represented using boolean data structures called Context Shapes. The relative orientations of the receptor and ligand surfaces are searched using precalculated lookup tables. Energetic quantities are derived from shape complementarity and buried surface area computations, using efficient boolean operations. Preliminary results indicate that our context shapes approach outperforms state-of-the-art geometric shape-based rigid-docking algorithms.

Nucleotide-induced structural changes in P-glycoprotein observed by electron microscopy

P-glycoprotein (Pgp) is an ATP hydrolysis driven multidrug efflux pump, which, when overexpressed in the plasma membrane of certain cancers, can lead to the failure of chemotherapy. Previously, we have presented a projection structure of nucleotide-free mouse Pgp from electron microscopic images of lipid monolayer-generated two-dimensional crystals ( Lee, J. Y., Urbatsch, I. L., Senior, A. E., and Wilkens, S. (2002) J. Biol. Chem. 277, 40125-40131 ). Here we have analyzed the structure of cysteine-free human Pgp from two-dimensional crystals that were generated with the same lipid-monolayer technique in the absence and presence of various nucleotides. The images show that human Pgp has a similar structure to the mouse protein. Furthermore, the analysis of projection structures obtained under different nucleotide conditions suggests that Pgp can exist in at least two major conformations, one of which shows a central cavity between the N- and C-terminal halves of the molecule and another in which the two halves have moved sideways, thereby closing the central cavity. Intermediate conformations were observed for some nucleotide/vanadate combinations. A low-resolution, three-dimensional model of human Pgp was calculated from tilted specimen crystallized in the presence of the non-hydrolyzable nucleotide analog, adenosine 5'-O-(thiotriphosphate). The structural analysis presented here adds to the emerging picture that multidrug ABC transporters function by switching between two major conformations in a nucleotide-dependent manner.

Nucleotide-dependent allostery within the ABC transporter ATP-binding cassette: a computational study of the MJ0796 dimer

ATP-binding cassette transporters perform energy-dependent transmembrane solute trafficking in all organisms. These proteins often mediate cellular resistance to therapeutic drugs and are involved in a range of human genetic diseases. Enzymological studies have implicated a helical subdomain within the ATP-binding cassette nucleotide-binding domain in coupling ATP hydrolysis to solute transport in the transmembrane domains. Consistent with this, structural and computational analyses have indicated that the helical subdomain undergoes nucleotide-dependent movement relative to the core of the nucleotide-binding domain fold. Here we use theoretical methods to examine the allosteric nucleotide dependence of helical subdomain transitions to further elucidate its role in interactions between the transmembrane and nucleotide-binding domains. Unrestrained 30-ns molecular dynamics simulations of the ATP-bound, ADP-bound, and apo states of the MJ0796 monomer support the idea that interaction of a conserved glutamine residue with the catalytic metal mediates the rotation of the helical subdomain in response to nucleotide binding and hydrolysis. Simulations of the nucleotide-binding domain dimer revealed that ATP hydrolysis induces a large transition of one helical subdomain, resulting in an asymmetric conformation of the dimer not observed previously. A coarse-grained elastic network analysis supports this finding, revealing the existence of corresponding dynamic modes intrinsic to the contact topology of the protein. The implications of these findings for the coupling of ATP hydrolysis to conformational changes in the transmembrane domains required for solute transport are discussed in light of recent whole transporter structures.

Protein free energy landscapes remodeled by ligand binding

Glucose/galactose binding protein (GGBP) functions in two different larger systems of proteins used by enteric bacteria for molecular recognition and signaling. Here we report on the thermodynamics of conformational equilibrium distributions of GGBP. Three fluorescence components appear at zero glucose concentration and systematically transition to three components at high glucose concentration. Fluorescence anisotropy correlations, fluorescent lifetimes, thermodynamics, computational structure minimization, and literature work were used to assign the three components as open, closed, and twisted conformations of the protein. The existence of three states at all glucose concentrations indicates that the protein continuously fluctuates about its conformational state space via thermally driven state transitions; glucose biases the populations by reorganizing the free energy profile. These results and their implications are discussed in terms of the two types of specific and nonspecific interactions GGBP has with cytoplasmic membrane proteins.

Substrate dependence of conformational changes in the RNA-binding domain of tristetraprolin assessed by fluorescence spectroscopy of tryptophan mutants

Association of tristetraprolin (TTP) with mRNAs containing selected AU-rich mRNA-destabilizing elements (AREs) initiates rapid cytoplasmic degradation of these transcripts. The RNA-binding activity of TTP is mediated by an internal tandem zinc finger domain that preferentially recognizes U-rich RNA ligands containing adjacent UUAU half-sites and is accompanied by conformational changes within the peptide. Here, we have used analogues of the TTP RNA-binding domain containing specific tryptophan substitutions to probe the Zn2+ and RNA substrate dependence of conformational events within individual zinc fingers. Fluorescence methods demonstrate that the N-terminal, but not C-terminal, zinc finger domain adopts a stably folded conformation in the presence of Zn2+. Denaturant titrations suggest that both the N- and C-terminal zinc fingers exhibit limited structural heterogeneity in the absence of RNA substrates, although this is more pronounced for the C-terminal finger. Binding to a cognate ARE substrate induced significant conformational changes within each zinc finger, which also included increased resistance to chemical denaturation. Studies with mutant ARE ligands revealed that a single UUAU half-site was sufficient to induce structural modulation of the N-terminal finger. However, RNA-dependent folding of the C-terminal zinc finger was only observed in the presence of tandem UUAU half-sites, suggesting that the conformation of this domain is linked not only to RNA substrate recognition but also to the ligand occupancy and/or conformational status of the N-terminal finger. Coupled with previous structural and thermodynamic analyses, these data provide a mechanistic framework for discrimination of RNA substrates involving ligand-dependent conformational adaptation of both zinc fingers within the TTP RNA-binding domain.

Set theory formulation of the model-free problem and the diffusion seeded model-free paradigm

Model-free analysis of NMR relaxation data, which describes the motion of individual atoms, is a problem intricately linked to the Brownian rotational diffusion of the macromolecule. The diffusion tensor parameters strongly influence the optimisation of the various model-free models and the subsequent model selection between them. Finding the optimal model of the dynamics of the system among the numerous diffusion and model-free models is hence quite complex. Using set theory, the entirety of this global problem has been encapsulated by the universal set Ll, and its resolution mathematically formulated as the universal solution Ll. Ever since the original Lipari and Szabo papers the model-free dynamics of a molecule has most often been solved by initially estimating the diffusion tensor. The model-free models which depend on the diffusion parameter values are then optimised and the best model is chosen to represent the dynamics of the residue. Finally, the global model of all diffusion and model-free parameters is optimised. These steps are repeated until convergence. For simplicity this approach to Ll will be labelled the diffusion seeded model-free paradigm. Although this technique suffers from a number of problems many have been solved. All aspects of the diffusion seeded paradigm and its consequences, together with a few alternatives to the paradigm, will be reviewed through the use of set notation.

Computer simulations of ABC transporter componentsThis paper is one of a selection of papers published in this Special Issue, entitled CSBMCB — Membrane Proteins in Health and Disease

Current computer simulation techniques provide robust tools for studying the detailed structure and functional dynamics of proteins, as well as their interaction with each other and with other biomolecules. In this minireview, we provide an illustration of recent progress and future challenges in computer modeling by discussing computational studies of ATP-binding cassette (ABC) transporters. ABC transporters have multiple components that work in a well coordinated fashion to enable active transport across membranes. The mechanism by which members of this superfamily execute transport remains largely unknown. Molecular dynamics simulations initiated from high-resolution crystal structures of several ABC transporters have proven to be useful in the investigation of the nature of conformational coupling events that may drive transport. In addition, fruitful efforts have been made to predict unknown structures of medically relevant ABC transporters, such as P-glycoprotein, using homology-based computational methods. The various techniques described here are also applicable to gaining an atomically detailed understanding of the functional mechanisms of proteins in general.

A novel ligand bound ABC transporter, LolCDE, provides insights into the molecular mechanisms underlying membrane detachment of bacterial lipoproteins

The LolCDE complex of Escherichia coli belongs to the ABC transporter superfamily and initiates the lipoprotein sorting to the outer membrane by catalysing their release from the inner membrane. LolC and/or LolE, membrane subunits, recognize lipoproteins anchored to the outer surface of the inner membrane, while LolD hydrolyses ATP on its inner surface. We report here that ligand-bound LolCDE can be purified from the inner membrane in the absence of ATP. Liganded LolCDE represents an intermediate of the release reaction and exhibits higher affinity for ATP than the unliganded form. ATP binding to LolD weakens the interaction between LolCDE and lipoproteins and causes their dissociation in a detergent solution, while lipoprotein release from membranes requires ATP hydrolysis. Liganded LolCDE thus reveals molecular events brought about through ATP binding and hydrolysis. LolCDE is the first example of an ABC transporter purified with tightly bound native substrates. A single molecule of lipoprotein is found to bind per molecule of the LolCDE complex.

A mutation in the S-switch region of the Runt domain alters the dynamics of an allosteric network responsible for CBFbeta regulation

The Runt domain is the DNA binding domain of the core binding factor (CBF) Runx subunits. The CBFs are transcription factors that play critical roles in hematopoiesis, bone, and neuron development in mammals. A common non-DNA binding CBFbeta subunit heterodimerizes with the Runt domain of the Runx proteins and allosterically regulates its affinity for DNA. Previous NMR dynamics studies suggested a model whereby CBFbeta allosterically regulates DNA binding by quenching conformational exchange in the Runt domain, particularly in the S-switch region and the betaE'-F loop. We sought to test this model, and to this end introduced all possible single amino acid substitutions into the S-switch region and the betaE'-F loop, and screened for mutations that enhanced DNA-binding. We demonstrate that one Runt domain mutant, R164N, binds both DNA and CBFbeta with higher affinity, but it is less sensitive to allosteric regulation by CBFbeta. Analysis of NMR relaxation data shows that the chemical exchange exhibited by the wild-type Runt domain is largely quenched by the R164N substitution. These data support a model in which the dynamic behavior of a network of residues connecting the CBFbeta and DNA binding sites on the Runt domain plays a critical role in the mechanism of allosteric regulation. This study provides an important functional link between dynamic behavior and protein allosteric function, consistent with results on other allosterically regulated proteins.

The dynamic energy landscape of dihydrofolate reductase catalysis

We used nuclear magnetic resonance relaxation dispersion to characterize higher energy conformational substates of Escherichia coli dihydrofolate reductase. Each intermediate in the catalytic cycle samples low-lying excited states whose conformations resemble the ground-state structures of preceding and following intermediates. Substrate and cofactor exchange occurs through these excited substates. The maximum hydride transfer and steady-state turnover rates are governed by the dynamics of transitions between ground and excited states of the intermediates. Thus, the modulation of the energy landscape by the bound ligands funnels the enzyme through its reaction cycle along a preferred kinetic path.

Protein flexibility and rigidity predicted from sequence

Structural flexibility has been associated with various biological processes such as molecular recognition and catalytic activity. In silico studies of protein flexibility have attempted to characterize and predict flexible regions based on simple principles. B-values derived from experimental data are widely used to measure residue flexibility. Here, we present the most comprehensive large-scale analysis of B-values. We used this analysis to develop a neural network-based method that predicts flexible-rigid residues from amino acid sequence. The system uses both global and local information (i.e., features from the entire protein such as secondary structure composition, protein length, and fraction of surface residues, and features from a local window of sequence-consecutive residues). The most important local feature was the evolutionary exchange profile reflecting sequence conservation in a family of related proteins. To illustrate its potential, we applied our method to 4 different case studies, each of which related our predictions to aspects of function. The first 2 were the prediction of regions that undergo conformational switches upon environmental changes (switch II region in Ras) and the prediction of surface regions, the rigidity of which is crucial for their function (tunnel in propeller folds). Both were correctly captured by our method. The third study established that residues in active sites of enzymes are predicted by our method to have unexpectedly low B-values. The final study demonstrated how well our predictions correlated with NMR order parameters to reflect motion. Our method had not been set up to address any of the tasks in those 4 case studies. Therefore, we expect that this method will assist in many attempts at inferring aspects of function.

Heteronuclear multidimensional NMR and homology modelling studies of the C-terminal nucleotide-binding domain of the human mitochondrial ABC transporter ABCB6

Human ATP-binding cassette, sub-family B, member 6 (ABCB6) is a mitochondrial ABC transporter, and presumably contributes to iron homeostasis. Aimed at understanding the structural basis for the conformational changes accompanying the substrate-transportation cycle, we have studied the C-terminal nucleotide-binding domain of ABCB6 (ABCB6-C) in both the nucleotide-free and ADP-bound states by heteronuclear multidimensional NMR and homology modelling. A non-linear sampling scheme was utilised for indirectly acquired 13C and 15N dimensions of all 3D triple-resonance NMR experiments, in order to overcome the instability and the low solubility of ABCB6-C. The backbone resonances for approximately 25% of non-proline residues, which are mostly distributed around the functionally important loops and in the Helical domain, were not observed for nucleotide-free form of ABCB6-C. From the pH, temperature and magnetic field strength dependencies of the resonance intensities, we concluded that this incompleteness in the assignments is mainly due to the exchange between multiple conformations at an intermediate rate on the NMR timescale. These localised conformational dynamics remained in ADP-bound ABCB6-C except for the loops responsible for adenine base and alpha/beta-phosphate binding. These results revealed that the localised dynamic cooperativity, which was recently proposed for a prokaryotic ABC MJ1267, also exists in a higher eukaryotic ABC, and is presumably shared by all members of the ABC family. Since the Helical domain is the putative interface to the transmembrane domain, this cooperativity may explain the coupled functions between domains in the substrate-transportation cycle.

Attempts to characterize the NBD heterodimer of MRP1: transient complex formation involves Gly771 of the ABC signature sequence but does not enhance the intrinsic ATPase activity

MRP1 (multidrug-resistance-associated protein 1; also known as ABCC1) is a member of the human ABC (ATP-binding cassette) transporter superfamily that confers cell resistance to chemotherapeutic agents. Considering the structural and functional similarities to the other ABC proteins, the interaction between its two NBDs (nucleotide-binding domains), NBD1 (N-terminal NBD) and NBD2 (C-terminal NBD), is proposed to be essential for the regulation of the ATP-binding/ATP-hydrolysis cycle of MRP1. We were interested in the ability of recombinant NBD1 and NBD2 to interact with each other and to influence ATPase activity. We purified NBD1 (Asn642-Ser871) and NBD2 (Ser1286-Val1531) as soluble monomers under native conditions. We measured extremely low intrinsic ATPase activity of NBD1 (10(-5) s(-1)) and NBD2 (6x10(-6) s(-1)) and no increase in the ATP-hydrolysis rate could be detected in an NBD1+NBD2 mixture, with concentrations up to 200 microM. Despite the fact that both monomers bind ATP, no stable NBD1.NBD2 heterodimer could be isolated by gel-filtration chromatography or native-PAGE, but we observed some significant modifications of the heteronuclear single-quantum correlation NMR spectrum of 15N-NBD1 in the presence of NBD2. This apparent NBD1.NBD2 interaction only occurred in the presence of Mg2+ and ATP. Partial sequential assignment of the NBD1 backbone resonances shows that residue Gly771 of the LSGGQ sequence is involved in NBD1.NBD2 complex formation. This is the first NMR observation of a direct interaction between the ABC signature and the opposite NBD. Our study also reveals that the NBD1.NBD2 heterodimer of MRP1 is a transient complex. This labile interaction is not sufficient to induce an ATPase co-operativity of the NBDs and suggests that other structures are required for the ATPase activation mechanism.

PROFbval: predict flexible and rigid residues in proteins

The mobility of a residue on the protein surface is closely linked to its function. The identification of extremely rigid or flexible surface residues can therefore contribute information crucial for solving the complex problem of identifying functionally important residues in proteins. Mobility is commonly measured by B-value data from high-resolution three-dimensional X-ray structures. Few methods predict B-values from sequence. Here, we present PROFbval, the first web server to predict normalized B-values from amino acid sequence. The server handles amino acid sequences (or alignments) as input and outputs normalized B-value and two-state (flexible/rigid) predictions. The server also assigns a reliability index for each prediction. For example, PROFbval correctly identifies residues in active sites on the surface of enzymes as particularly rigid.

AVAILABILITY

http://www.rostlab.org/services/profbval

CONTACT

profbval@rostlab.org

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Solution structure of the low-molecular-weight protein tyrosine phosphatase from Tritrichomonas foetus reveals a flexible phosphate binding loop

Eukaryotic low-molecular-weight protein tyrosine phosphatases (LMW PTPs) contain a conserved serine, a histidine with an elevated pKa, and an active site asparagine that together form a highly conserved hydrogen bonding network. This network stabilizes the active site phosphate binding loop for optimal substrate binding and catalysis. In the phosphatase from the bovine parasite Tritrichomonas foetus (TPTP), both the conserved serine (S37) and asparagine (N14) are present, but the conserved histidine has been replaced by a glutamine residue (Q67). Site-directed mutagenesis, kinetic, and spectroscopic experiments suggest that Q67 is located near the active site and is important for optimal catalytic activity. Kinetic experiments also suggest that S37 participates in the active site/hydrogen bonding network. Nuclear magnetic resonance spectroscopy was used to determine the three-dimensional structure of the TPTP enzyme and to further examine the roles of S37 and Q67. The backbone conformation of the TPTP phosphate binding loop is nearly superimposable with that of other tyrosine phosphatases, with N14 existing in a strained, left-handed conformation that is a hallmark of the active site hydrogen bonding network in the LMW PTPs. As expected, both S37 and Q67 are located at the active site, but in the consensus structure they are not within hydrogen bonding distance of N14. The hydrogen bond interactions that are observed in X-ray structures of LMW PTPs may in fact be transient in solution. Protein dynamics within the active site hydrogen bonding network appear to be affected by the presence of substrate or bound inhibitors such as inorganic phosphate.

Nucleotide binding to the homodimeric MJ0796 protein: a computational study of a prokaryotic ABC transporter NBD dimer

Transport by ABC proteins requires a cycle of ATP-driven conformational changes of the nucleotide binding domains (NBDs). We compare three molecular dynamics simulations of dimeric MJ0796: with ATP was present at both NBDs; with ATP at one NBD but ADP at the other; and without any bound ATP. In the simulation with ATP present at both NBDs, the dimeric protein interacts with the nucleotides in a symmetrical manner. However, if ADP is present at one binding site then both NBD-NBD and protein-ATP interactions are enhanced at the opposite site.

Role of Conformational Heterogeneity in Domain Swapping and Adapter Function of the Cks Proteins

Cks proteins are adapter molecules that coordinate the assembly of multiprotein complexes. They share the ability to domain swap by exchanging a beta-strand, beta4. Here we use NMR spectroscopy and molecular dynamics simulations to investigate the dynamic properties of human Cks1 and its response on assembly with components of the SCF(Skp2) ubiquitin ligation machinery. In the NMR experiment with the free form of Cks1, a subset of residues displayed elevated R2 values and the cross-peaks of neighboring residues were missing from the spectrum, indicating a substantial conformational exchange contribution on the microsecond to millisecond time scale. Strikingly the region of greatest conformational variability was the beta4-strand that domain swaps to form the dimer. Binding of the ligand common to all Cks proteins, Cdk2, suppressed the conformational heterogeneity. This response was specific to Cdk2 binding; in contrast, binding of Skp2, a ligand unique to human Cks1, did not alter the dynamic behavior. Short time (<5 ns) molecular dynamics simulations indicate that residues of Cks1 that form the binding site for phosphorylated ligands are considerably more flexible in the free form of Cks1 than they are in the Cdk2-Cks1 complex. A cooperative interaction between Cdk2 and Cks1 is suggested, which reduces the configurational entropy of Cks1 and therefore facilitates phosphoprotein binding. Indications of an unusual dynamic behavior of strand beta4 in the free form of Cks1 were obtained from longer time scale (50 ns) dynamics simulations. A spontaneous reversible unzipping of hydrogen bonds between beta4 and beta2 was observed, suggesting an early intermediate structure for unfolding and/or domain swapping. We propose that the dynamic properties of the beta-sheet and its modification upon ligand binding underlie the domain swapping ability and the adapter function of Cks proteins.