Dissection of the conformational cycle of the multidrug/lipidA ABC exporter MsbA

Recent crystal structures of the multidrug ATP‐binding cassette (ABC) exporters Sav1866 from Staphylococcus aureus, MsbA from Escherichia coli, Vibrio cholera, and Salmonella typhimurium, and mouse ABCB1a suggest a common alternating access mechanism for export. However, the molecular framework underlying this mechanism is critically dependent on assumed conformational relationships between nonidentical crystal structures and therefore requires biochemical verification. The structures of homodimeric MsbA reveal a pair of glutamate residues (E208 and E208′) in the intracellular domains of its two half‐transporters, close to the nucleotide‐binding domains (NBDs), which are in close proximity of each other in the outward‐facing state but not in the inward‐facing state. Using intermolecular cysteine crosslinking between E208C and E208C′ in E. coli MsbA, we demonstrate that the NBDs dissociate in nucleotide‐free conditions and come close on ATP binding and ADP·vanadate trapping. Interestingly, ADP alone separates the half‐transporters like a nucleotide‐free state, presumably for the following catalytic cycle. Our data fill persistent gaps in current studies on the conformational dynamics of a variety of ABC exporters. Based on a single biochemical method, the findings describe a conformational cycle for a single ABC exporter at major checkpoints of the ATPase reaction under experimental conditions, where the exporter is transport active. Proteins 2010. © 2010 Wiley‐Liss, Inc.     

Gain of local structure in an amphipathic peptide does not require a specific tertiary framework

In this work, we studied how an amphipathic peptide of the surface of the globular protein thioredoxin, TRX94‐108, acquires a native‐like structure when it becomes involved in an apolar interaction network. We designed peptide variants where the tendency to form α‐helical conformation is modulated by replacing each of the leucine amino acid residues by an alanine. The induction of structure caused by sodium dodecyl sulfate (SDS) binding was studied by capillary zone electrophoresis, circular dichroism, DOSY‐NMR, and molecular dynamics simulations (MDS). In addition, we analyzed the strength of the interaction between a C18 RP‐HPLC matrix and the peptides. The results presented here reveal that (a) critical elements in the sequence of the wild‐type peptide stabilize a SDS/peptide supramolecular cluster; (b) the hydrophobic nature of the interaction between SDS molecules and the peptide constrains the ensemble of conformations; (c) nonspecific apolar surfaces are sufficient to stabilize peptide secondary structure. Remarkably, MDS shed light on a contact network formed by a limited number of SDS molecules that serves as a structural scaffold preserving the helical conformation of this module. This mechanism might prevail when a peptide with low helical propensity is involved in structure consolidation. We suggest that folding of peptides sharing this feature does not require a preformed tightly‐packed protein core. Thus, the formation of specific tertiary interactions would be the consequence of peptide folding and not its cause. In this scenario, folding might be thought of as a process that includes unspecific rounds of structure stabilization guiding the protein to the native state. Proteins 2010. © 2010 Wiley‐Liss, Inc.     

The effect of calcium on the conformation of cobalamin transporter BtuB

BtuB is a β‐barrel membrane protein that facilitates transport of cobalamin (vitamin B12) from the extracellular medium across the outer membrane of Escherichia coli. It is thought that binding of B12 to BtuB alters the conformation of its periplasm‐exposed N‐terminal residues (the TonB box), which enables subsequent binding of a TonB protein and leads to eventual uptake of B12 into the cytoplasm. Structural studies determined the location of the B12 binding site at the top of the BtuB's β‐barrel, surrounded by extracellular loops. However, the structure of the loops was found to depend on the method used to obtain the protein crystals, which—among other factors—differed in calcium concentration. Experimentally, calcium concentration was found to modulate the binding of the B12 substrate to BtuB. In this study, we investigate the effect of calcium ions on the conformation of the extracellular loops of BtuB and their possible role in B12 binding. Using all‐atom molecular dynamics, we simulate conformational fluctuations of several X‐ray structures of BtuB in the presence and absence of calcium ions. These simulations demonstrate that calcium ions can stabilize the conformation of loops 3–4, 5–6, and 15–16, and thereby prevent occlusion of the binding site. Furthermore, binding of calcium ions to extracellular loops of BtuB was found to enhance correlated motions in the BtuB structure, which is expected to promote signal transduction. Finally, we characterize conformation dynamics of the TonB box in different X‐ray structures and find an interesting correlation between the stability of the TonB box structure and calcium binding. Proteins 2010. © 2009 Wiley‐Liss, Inc.     

Population dynamics of the West Indian sweetpotato weevil Euscepes postfasciatus (Fairmaire): a simulation analysis

The West Indian sweetpotato weevil Euscepes postfasciatus (Fairmaire) is a major pest of the sweet potato Ipomoea batatas (L.) Lam. and this weevil is a target of an eradication program using the Sterile Insect Technique in Okinawa Prefecture, Japan. Understanding the population ecology is essential in the planning of an eradication program; hence, a host‐plant infestation survey and light trap survey have been conducted to monitor the population dynamics of the weevil on Kume Island (Okinawa Prefecture), which is the target area of the trial weevil eradication project. Seasonal tendencies of weevil density were found in these field surveys, but the tendency found in the host‐plant infestation survey was not seen every year, and the effectiveness of the light trap is somewhat suspect. To confirm the reliability of the tendency observed in these field surveys, the present study attempted to explain the tendency by a seasonal temperature change using a temperature‐based model of weevil population dynamics. The seasonal changes of weevil density differed according to host plants and host‐plant fields. The seasonal changes of weevil density inside the host plant Ipomoea indica and outside the host plants in I. indica fields were consistent with those predicted by the model. However, those inside the host plant Ipomoea pes‐caprae in the host‐plant infestation survey were contrary to the predicted ones, and those observed outside host plants in I. pes‐caprae fields by the light trap survey were not in good agreement with the predicted ones. It was concluded that the seasonal change of the weevil density observed in I. indica and I. indica fields can be explained by a seasonal temperature change, but factors other than seasonal temperature change are needed to explain those in I. pes‐caprae and I. pes‐caprae fields.     

Investigating the disorder–order transition of calmodulin binding domain upon binding calmodulin using molecular dynamics simulation

We have studied the conformational transition of the calmodulin binding domains (CBD) in several calmodulin‐binding kinases, in which CBD changes from the disordered state to the ordered state when binding with calmodulin (CaM). Targeted molecular dynamics simulation was used to investigate the binding process of CaM and CBD of CaM‐dependent kinase I (CaMKI–CBD). The results show that CaMKI–CBD began to form an α‐helix and the interaction free energy between CaM and CaMKI–CBD increased once CaM fully encompassed CaMKI–CBD. Two series of CaM/CBD complex systems, including the complexes of CaM with the initially disordered and the final ordered CBD, were constructed to study the interaction using molecular dynamics simulations. Our analyses suggest that the VDW interaction plays a dominant role in CaM/CBD binding and is a key factor in the disorder–order transition of CBD. Additionally, the entropy effect is not in favor of the formation of the CaM/CBD complex, which is consistent with the experimental evidence. Based on the results, it appears that the CBD conformational change from a non‐compact extended structure to compact α‐helix is critical in gaining a favorable VDW interaction and interaction free energy. Copyright © 2009 John Wiley & Sons, Ltd.     

Contribution of the trifluoroacetyl group in the thermodynamics of antigen–antibody binding

We analyzed the binding of the 7C8 antibody to the chloramphenicol phosphonate antigens—one containing a trifluoroacetyl group (CP‐F) and the other containing an acetyl group (CP‐H)—by using isothermal titration calorimetry (ITC). The thermodynamic difference due to the substitution of F by H was evaluated using free energy calculations based on molecular dynamics (MD) simulations. We have previously shown that another antibody, namely, 6D9, binds more weakly to CP‐H than to CP‐F, mainly due to the different hydration free energies of the dissociated state and not due to the unfavorable hydrophobic interactions with the antibody in the bound state. Unlike in the binding of the trifluoroacetyl group with 6D9, in its binding with 7C8, it is exposed to the solvent, as seen in the crystal structure of the complex of 7C8 with CP‐F. The thermodynamic analysis performed in this study showed that the binding affinity of 7C8 for CP‐H is similar to that for CP‐F, but this binding to CP‐H is accompanied with less favorable enthalpy and more favorable entropy changes. The free energy calculations indicated that, upon the substitution of F by H, enthalpy and entropy changes in the associated and dissociated states were decreased, but the magnitude of enthalpy and entropy changes in the dissociated state was larger than that in the associated state. The differences in binding free energy, enthalpy, and entropy changes determined by the free energy calculations for the substitution of F by H are in good agreement with the experimental results. Copyright © 2009 John Wiley & Sons, Ltd.     

Protein folding pathways and state transitions described by classical equations of motion of an elastic network model

Protein topology defined by the matrix of residue contacts has proved to be a fruitful basis for the study of protein dynamics. The widely implemented coarse-grained elastic network model of backbone fluctuations has been used to describe crystallographic temperature factors, allosteric couplings, and some aspects of the folding pathway. In the present study, we develop a model of protein dynamics based on the classical equations of motion of a damped network model (DNM) that describes the folding path from a completely unfolded state to the native conformation through a single-well potential derived purely from the native conformation. The kinetic energy gained through the collapse of the protein chain is dissipated through a friction term in the equations of motion that models the water bath. This approach is completely general and sufficiently fast that it can be applied to large proteins. Folding pathways for various proteins of different classes are described and shown to correlate with experimental observations and molecular dynamics and Monte Carlo simulations. Allosteric transitions between alternative protein structures are also modeled within the DNM through an asymmetric double-well potential.     

Solvent-induced lid opening in lipases: A molecular dynamics study

In most lipases, a mobile lid covers the substrate binding site. In this closed structure, the lipase is assumed to be inactive. Upon activation of the lipase by contact with a hydrophobic solvent or at a hydrophobic interface, the lid opens. In its open structure, the substrate binding site is accessible and the lipase is active. The molecular mechanism of this interfacial activation was studied for three lipases (from Candida rugosa, Rhizomucor miehei, and Thermomyces lanuginosa) by multiple molecular dynamics simulations for 25 ns without applying restraints or external forces. As initial structures of the simulations, the closed and open structures of the lipases were used. Both the closed and the open structure were simulated in water and in an organic solvent, toluene. In simulations of the closed lipases in water, no conformational transition was observed. However, in three independent simulations of the closed lipases in toluene the lid gradually opened. Thus, pathways of the conformational transitions were investigated and possible kinetic bottlenecks were suggested. The open structures in toluene were stable, but in water the lid of all three lipases moved towards the closed structure and partially unfolded. Thus, in all three lipases opening and closing was driven by the solvent and independent of a bound substrate molecule.     

Osmolytes modulate conformational exchange in solvent‐exposed regions of membrane proteins

Site‐directed spin labeling (SDSL) was used to investigate local structure and conformational exchange in two bacterial outer‐membrane TonB‐dependent transporters, BtuB and FecA. Protecting osmolytes, such as polyethylene glycols (PEGs) are known to modulate a substrate‐dependent conformational equilibrium in the energy coupling motif (Ton box) of BtuB. Here, we demonstrate that a segment that is N‐terminal to the Ton box in BtuB, is in conformational exchange between ordered and disordered states with or without substrate. Protecting osmolytes shift this equilibrium to favor the more ordered, folded state. However, a segment of BtuB that is C‐terminal to the Ton box that is not solvent exposed is insensitive to PEGs. Protecting osmolytes also modulate a conformational equilibrium in the Ton box of FecA, with larger molecular weight PEGs producing the largest shifts in the conformational free energy. These data indicate that solvent‐exposed regions of these transporters undergo conformational exchange and that regions of these transporters that are involved in protein–protein interactions sample multiple conformational substates. The sensitivity to solute provides an explanation for differences seen between two high‐resolution structures of BtuB, which each likely represent one conformation from a subset of states that are normally sampled by the protein. This work also illustrates how SDSL and osmolytes may be used to characterize and quantitate conformational equilibria in membrane proteins.