The rigid connecting loop stabilizes hairpin folding of the two helices of the ATP synthase subunit c

Mitochondrial ATP synthase c-subunit leak channel triggers cell death upon loss of its F1 subcomplex

Mitochondrial ATP synthase is vital not only for cellular energy production but also for energy dissipation and cell death. ATP synthase c-ring was suggested to house the leak channel of mitochondrial permeability transition (mPT), which activates during excitotoxic ischemic insult. In this present study, we purified human c-ring from both eukaryotic and prokaryotic hosts to biophysically characterize its channel activity. We show that purified c-ring forms a large multi-conductance, voltage-gated ion channel that is inhibited by the addition of ATP synthase F₁ subcomplex. In contrast, dissociation of F₁ from F_O occurs during excitotoxic neuronal death suggesting that the F₁ constitutes the gate of the channel. mPT is known to dissipate the osmotic gradient across the inner membrane during cell death. We show that ATP synthase c-subunit knock down (KD) prevents the osmotic change in response to high calcium and eliminates large conductance, Ca²⁺ and CsA sensitive channel activity of mPT. These findings elucidate the gating mechanism of the ATP synthase c-subunit leak channel (ACLC) and suggest how ACLC opening is regulated by cell stress in a CypD-dependent manner.

Sequence and conformational preferences at termini of α-helices in membrane proteins: role of the helix environment

α-Helices are amongst the most common secondary structural elements seen in membrane proteins and are packed in the form of helix bundles. These α-helices encounter varying external environments (hydrophobic, hydrophilic) that may influence the sequence preferences at their N and C-termini. The role of the external environment in stabilization of the helix termini in membrane proteins is still unknown. Here we analyze α-helices in a high-resolution dataset of integral α-helical membrane proteins and establish that their sequence and conformational preferences differ from those in globular proteins. We specifically examine these preferences at the N and C-termini in helices initiating/terminating inside the membrane core as well as in linkers connecting these transmembrane helices. We find that the sequence preferences and structural motifs at capping (Ncap and Ccap) and near-helical (N' and C') positions are influenced by a combination of features including the membrane environment and the innate helix initiation and termination property of residues forming structural motifs. We also find that a large number of helix termini which do not form any particular capping motif are stabilized by formation of hydrogen bonds and hydrophobic interactions contributed from the neighboring helices in the membrane protein. We further validate the sequence preferences obtained from our analysis with data from an ultradeep sequencing study that identifies evolutionarily conserved amino acids in the rat neurotensin receptor. The results from our analysis provide insights for the secondary structure prediction, modeling and design of membrane proteins.

Residues in the polar loop of subunit c in Escherichia coli ATP synthase function in gating proton transport to the cytoplasm

Rotary catalysis in F1F0 ATP synthase is powered by proton translocation through the membrane-embedded F0 sector. Proton binding and release occur in the middle of the membrane at Asp-61 on the second transmembrane helix (TMH) of subunit c, which folds in a hairpin-like structure with two TMHs. Previously, the aqueous accessibility of Cys substitutions in the transmembrane regions of subunit c was probed by testing the inhibitory effects of Ag(+) or Cd(2+) on function, which revealed extensive aqueous access in the region around Asp-61 and on the half of TMH2 extending to the cytoplasm. In the current study, we surveyed the Ag(+) and Cd(2+) sensitivity of Cys substitutions in the loop of the helical hairpin and used a variety of assays to categorize the mechanisms by which Ag(+) or Cd(2+) chelation with the Cys thiolates caused inhibition. We identified two distinct metal-sensitive regions in the cytoplasmic loop where function was inhibited by different mechanisms. Metal binding to Cys substitutions in the N-terminal half of the loop resulted in an uncoupling of F1 from F0 with release of F1 from the membrane. In contrast, substitutions in the C-terminal half of the loop retained membrane-bound F1 after metal treatment. In several of these cases, inhibition was shown to be due to blockage of passive H(+) translocation through F0 as assayed with F0 reconstituted into liposomes. The results suggest that the C-terminal domain of the cytoplasmic loop may function in gating H(+) translocation to the cytoplasm.

Conformational sampling in template-free protein loop structure modeling: an overview

Accurately modeling protein loops is an important step to predict three-dimensional structures as well as to understand functions of many proteins. Because of their high flexibility, modeling the three-dimensional structures of loops is difficult and is usually treated as a "mini protein folding problem" under geometric constraints. In the past decade, there has been remarkable progress in template-free loop structure modeling due to advances of computational methods as well as stably increasing number of known structures available in PDB. This mini review provides an overview on the recent computational approaches for loop structure modeling. In particular, we focus on the approaches of sampling loop conformation space, which is a critical step to obtain high resolution models in template-free methods. We review the potential energy functions for loop modeling, loop buildup mechanisms to satisfy geometric constraints, and loop conformation sampling algorithms. The recent loop modeling results are also summarized.

Production and initial structural characterization of the TM4TM5 helix-loop-helix domain of the translocator protein

Mainly present in the mitochondria, the translocator protein, TSPO, previously known as the peripheral benzodiazepine receptor, is a small essential membrane protein, involved in the translocation of cholesterol across mitochondrial membranes, a rate determining step in steroids biosynthesis. We previously reported the structure of five fragments encompassing the five putative transmembrane helices and showed that each of these fragments constitutes an autonomous folding unit. To further characterize the structural determinants responsible for helix-helix association of this membrane protein, we now investigate the folding of double transmembrane domains in various detergent micelles. Herein, we present the successful biosynthesis of a double transmembrane domain encompassing the last two C-terminal helices (TM4TM5). For optimal production of this domain in Escherichia coli, the evaluation of various peptide constructs, including TM4TM5 fused to different purification tags or to solubilizing proteins, was necessary. The protocol of production of TM4TM5 with more than 95% purity is reported. This domain was further characterized using circular dichroism and solution state NMR. Far-UV circular dichroism studies indicate that the secondary structure of TM4TM5 is highly helical when solubilized in various detergent micelles including n-dodecyl-β-d-maltoside, n-octyl-β-d-glucoside, n-dodecylphosphocholine, 1,2-dihexanoyl-sn-glycero-3-phosphocholine (DHPC), and 1-palmitoyl-2-hydroxy-sn-glycero-3-phospho-(1'-rac-glycerol). In addition, the solubilization conditions of the domain were optimized for NMR experiments, and preliminary analysis indicates that TM4TM5 adopts a stable tertiary fold within the TM4TM5-DHPC complex.

Sampling multiple scoring functions can improve protein loop structure prediction accuracy

Accurately predicting loop structures is important for understanding functions of many proteins. In order to obtain loop models with high accuracy, efficiently sampling the loop conformation space to discover reasonable structures is a critical step. In loop conformation sampling, coarse-grain energy (scoring) functions coupling with reduced protein representations are often used to reduce the number of degrees of freedom as well as sampling computational time. However, due to implicitly considering many factors by reduced representations, the coarse-grain scoring functions may have potential insensitivity and inaccuracy, which can mislead the sampling process and consequently ignore important loop conformations. In this paper, we present a new computational sampling approach to obtain reasonable loop backbone models, so-called the Pareto optimal sampling (POS) method. The rationale of the POS method is to sample the function space of multiple, carefully selected scoring functions to discover an ensemble of diversified structures yielding Pareto optimality to all sampled conformations. The POS method can efficiently tolerate insensitivity and inaccuracy in individual scoring functions and thereby lead to significant accuracy improvement in loop structure prediction. We apply the POS method to a set of 4-12-residue loop targets using a function space composed of backbone-only Rosetta and distance-scale finite ideal-gas reference (DFIRE) and a triplet backbone dihedral potential developed in our lab. Our computational results show that in 501 out of 502 targets, the model sets generated by POS contain structure models are within subangstrom resolution. Moreover, the top-ranked models have a root mean square deviation (rmsd) less than 1 A in 96.8, 84.1, and 72.2% of the short (4-6 residues), medium (7-9 residues), and long (10-12 residues) targets, respectively, when the all-atom models are generated by local optimization from the backbone models and are ranked by our recently developed Pareto optimal consensus (POC) method. Similar sampling effectiveness can also be found in a set of 13-residue loop targets.

Improving predicted protein loop structure ranking using a Pareto-optimality consensus method

Background

Accurate protein loop structure models are important to understand functions of many proteins. Identifying the native or near-native models by distinguishing them from the misfolded ones is a critical step in protein loop structure prediction.

Results

We have developed a Pareto Optimal Consensus (POC) method, which is a consensus model ranking approach to integrate multiple knowledge- or physics-based scoring functions. The procedure of identifying the models of best quality in a model set includes: 1) identifying the models at the Pareto optimal front with respect to a set of scoring functions, and 2) ranking them based on the fuzzy dominance relationship to the rest of the models. We apply the POC method to a large number of decoy sets for loops of 4- to 12-residue in length using a functional space composed of several carefully-selected scoring functions: Rosetta, DOPE, DDFIRE, OPLS-AA, and a triplet backbone dihedral potential developed in our lab. Our computational results show that the sets of Pareto-optimal decoys, which are typically composed of ~20% or less of the overall decoys in a set, have a good coverage of the best or near-best decoys in more than 99% of the loop targets. Compared to the individual scoring function yielding best selection accuracy in the decoy sets, the POC method yields 23%, 37%, and 64% less false positives in distinguishing the native conformation, indentifying a near-native model (RMSD < 0.5A from the native) as top-ranked, and selecting at least one near-native model in the top-5-ranked models, respectively. Similar effectiveness of the POC method is also found in the decoy sets from membrane protein loops. Furthermore, the POC method outperforms the other popularly-used consensus strategies in model ranking, such as rank-by-number, rank-by-rank, rank-by-vote, and regression-based methods.

Conclusions

By integrating multiple knowledge- and physics-based scoring functions based on Pareto optimality and fuzzy dominance, the POC method is effective in distinguishing the best loop models from the other ones within a loop model set.

F1F0-ATP synthases of alkaliphilic bacteria: lessons from their adaptations

This review focuses on the ATP synthases of alkaliphilic bacteria and, in particular, those that successfully overcome the bioenergetic challenges of achieving robust H+-coupled ATP synthesis at external pH values>10. At such pH values the protonmotive force, which is posited to provide the energetic driving force for ATP synthesis, is too low to account for the ATP synthesis observed. The protonmotive force is lowered at a very high pH by the need to maintain a cytoplasmic pH well below the pH outside, which results in an energetically adverse pH gradient. Several anticipated solutions to this bioenergetic conundrum have been ruled out. Although the transmembrane sodium motive force is high under alkaline conditions, respiratory alkaliphilic bacteria do not use Na+- instead of H+-coupled ATP synthases. Nor do they offset the adverse pH gradient with a compensatory increase in the transmembrane electrical potential component of the protonmotive force. Moreover, studies of ATP synthase rotors indicate that alkaliphiles cannot fully resolve the energetic problem by using an ATP synthase with a large number of c-subunits in the synthase rotor ring. Increased attention now focuses on delocalized gradients near the membrane surface and H+ transfers to ATP synthases via membrane-associated microcircuits between the H+ pumping complexes and synthases. Microcircuits likely depend upon proximity of pumps and synthases, specific membrane properties and specific adaptations of the participating enzyme complexes. ATP synthesis in alkaliphiles depends upon alkaliphile-specific adaptations of the ATP synthase and there is also evidence for alkaliphile-specific adaptations of respiratory chain components.

Probing the rotor subunit interface of the ATP synthase from Ilyobacter tartaricus

The interaction between the c(11)ring and the gammaepsilon complex, forming the rotor of the Ilyobacter tartaricus ATP synthase, was probed by surface plasmon resonance spectroscopy and in vitro reconstitution analysis. The results provide, for the first time, a direct and quantitative assessment of the stability of the rotor. The data indicated very tight binding between the c(11)ring and the gammaepsilon complex, with an apparent K(d) value of approximately 7.4nm. The rotor assembly was primarily dependent on the interaction of the cring with the gammasubunit, and binding of the cring to the free epsilon subunit was not observed. Mutagenesis of selected conserved amino acid residues of all three rotor components (cR45, cQ46, gammaE204, gammaF203 and epsilonH38) severely affected rotor assembly. The interaction kinetics between the gammaepsilon complex and c(11)ring mutants suggested that the assembly of the c(11)gammaepsiloncomplex was governed by interactions of low and high affinity. Low-affinity binding was observed between the polar loops of the cring subunits and the bottom part of the gamma subunit. High-affinity interactions, involving the two residues gammaE204 and epsilonH38, stabilized the holo-c(11)gammaepsilon complex. NMR experiments indicated the acquisition of conformational order in otherwise flexible C- and N-terminal regions of the gamma subunit on rotor assembly. The results of this study suggest that docking of the central stalk of the F(1)complex to the rotor ring of F(o) to form tight, but reversible, contacts provides an explanation for the relative ease of dissociation and reconstitution of F(1)F(o)complexes.