Simulation studies on bacteriorhodopsin bundle of transmembrane α segments

Bacteriorhodopsin (BR) is a membrane protein which pumps protons through the plasma membrane. Seven transmembrane BR helical segments are subjected to simulation studies in order to investigate the packing process of transmembrane helices. A Monte Carlo simulated annealing protocol is employed to optimize the helix bundle system. Helix packing is optimized according to a semi-empirical potential mainly composed of six components: a bilayer potential, a crossing angle potential, a helix dipole potential, a helix-helix distance potential, a helix orientation potential and a helix-helix distance restraint potential (a loop potential). Necessary parameters are derived from theoretical studies and statistical analysis of experimentally determined protein structures. The structures from the simulations are compared with the experimentally determined structures in terms of geometry. The structures generated show similar shapes to the experimentally suggested structure even without the helix-helix distance restraint potential. However, the relative locations of individual helices were reproduced only when the helix-helix distance restraint potential was used with restraint conditions. Our results suggest that transmembrane helix bundles resembling those observed experimentally may be generated by simulations using simple potentials. Received: 19 April 1999 / Revised version: 6 September 1999 / Accepted: 17 September 1999     

Modelling of transmembrane α-helix bundles

A strategy for modelling transmembrane -helix bundles has been investigated. Results concerning the rotational orientations of the helices are described and perspectives for extensions of the method are discussed.     

Characterization and sequence prediction of structural variations in α-helix

Background

Centralised resources such as GenBank and UniProt are perfect examples of the major international efforts that have been made to integrate and share biological information. However, additional data that adds value to these resources needs a simple and rapid route to public access. The Distributed Annotation System (DAS) provides an adequate environment to integrate genomic and proteomic information from multiple sources, making this information accessible to the community. DAS offers a way to distribute and access information but it does not provide domain experts with the mechanisms to participate in the curation process of the available biological entities and their annotations.

Results

We designed and developed a Collaborative Annotation System for proteins called DAS Writeback. DAS writeback is a protocol extension of DAS to provide the functionalities of adding, editing and deleting annotations. We implemented this new specification as extensions of both a DAS server and a DAS client. The architecture was designed with the involvement of the DAS community and it was improved after performing usability experiments emulating a real annotation task.

Conclusions

We demonstrate that DAS Writeback is effective, usable and will provide the appropriate environment for the creation and evolution of community protein annotation.     

Hydrophobic channels in crystals of an α-aminoisobutyric acid pentapeptide

The crystal structure of the pentapeptide p-toluene-sulfonyl-(α-aminoisobutyryl)₅-methyl ester (Tosyl-(Aib)₅-OMe) has been determined in the space group P$\text{I}$. Pentapeptide molecules are folded in the 3₁₀ helical conformation and packed together, so as to yield a hydrophobic channel with a minimim diameter of 5.2 Å. The channel contains crystallographically disordered material. This structure provides a model for channel formation by hydrophobic peptide aggregates and should prove useful in studies of alamethicin, suzukacillin and related Aib containing membrane channels. Triclinic (P$\text{I}$) crystals of the pentapeptide are obtained in the presence of LiClO₄ in aqueous methanol, whereas crystallization from methanol alone yields crystals in the space group Pbca. The conformations of the peptide in the two crystal forms are very similar and only the molecular packing is dramatically different.     

Synthesis of antiparallel 4α-helix bundle TASP by chemoselective ligation

Summary The use of chemoselective ligation methods and orthogonal protection techniques allows access to Template-Assembled Synthetic Protein (TASP) molecules exhibiting a large variety of packing topologies. This is demonstrated for the synthesis of an antiparallel 4-helical bundle TASP by condensing amphiphilic peptide blocks, containing aldehyde functions at the C- or N-terminus, to a selectively addressable topological template via oxime bond formation. The resulting antiparallel 4-helix TASP is obtained in high yield and shows a template-induced helical conformation.     

Sequence preference of α-helix N-terminal tetrapeptide

The α-helix is the most abundant secondary structure in proteins. Due to the specific i, i+4 hydrogen bond pattern, the two termini have unsatisfied hydrogen bonds, and are less constrained; in order to compensate for this, specific residues are preferred for the terminal positions. However, a naive combination of the statistically-preferred residues for each position may not result in a stable N-terminal helical sequence. In order to provide a set of preferable N-terminal peptides for α-helix design, we have studied the N-terminal tetrapeptide sequence motifs that are favorable for helix formation using statistical analysis and atomistic simulations. A set of tetrapeptide sequences including TEEE and TPEE were found to be favorable motifs. In addition to forming more hydrogen bonds in the helical conformation, the favorable motifs also tended to form more capping boxes. To empirically test our predictions, we obtained 10 peptides with different N-terminal motifs and measured their α-helical content by circular dichroism spectroscopy. The experimental results agreed qualitatively with the statistical and simulation results. Furthermore, some of the suggested preferable tetrapeptide sequences have been successfully applied in de novo protein design.     

Functional characterization of amphipathic α-helix in the osmoregulatory ABC transporter OpuA

The ATP-binding-cassette transporter OpuA from Lactococcus lactis is composed of two ATPase subunits (OpuAA) and two subunits (OpuABC) with the transmembrane domain fused to an extracellular substrate-binding protein. Of the almost 1900 homologues of OpuA known to date, a subset has an amino-terminal amphipathic helix (plus extra transmembrane segment) fused to the core of the transmembrane domain of the OpuABC subunit. FRET measurements indicate that the amphipathic α-helix is located close to the membrane surface, where its hydrophobic face interacts with the transport protein rather than the membrane lipids. Next, we determined the functional role of this accessory region by engineering the amphipathic α-helix. We analyzed the consequence of the mutations in intact cells by monitoring growth and transport of glycine betaine under normal and osmotic stress conditions. More detailed studies were performed in hybrid membrane vesicles, proteoliposomes, and bilayer nanodisks. We show that the amphipathic α-helix of OpuA is necessary for high activity of OpuA but is not critical for the biogenesis of the protein or the ionic regulation of transport.     

Simulation studies on bacteriorhodopsin α-helices

Bacteriorhodopsin (BR) is a membrane protein which pumps protons through the plasma membrane. Transmembrane BR helical segments are subjected to simulation studies in order to investigate the effect of bilayer environment in various simulation conditions. A bilayer potential is introduced to the system to mimic the lipid membrane. The structures from the simulations are compared with the experimentally determined structures in terms of geometrical properties. Electrostatic contribution to the helix packing is also investigated. The simulation results show that the packing geometry of the transmembrane helices is highly affected by the bilayer potential. The results obtained from the simulations may be used for further simulation studies and analysis in investigating transmembrane helix packing. Received: 5 July 1999 / Revised version: 5 October 1999 / Accepted: 15 October 1999     

Genomic organization of the complex α-gliadin gene loci in wheat

To better understand the molecular evolution of the large -gliadin gene family, a half-million bacterial artificial chromosome (BAC) library clones from tetraploid durum wheat, Triticum turgidum ssp. durum (2n=4x=28, genome AB), were screened for large genomic segments carrying the -gliadin genes of the Gli-2 loci on the group 6 homoeologous chromosomes. The resulting 220 positive BAC clones—each containing between one and four copies of -gliadin sequences—were fingerprinted for contig assembly to produce contiguous chromosomal regions covering the Gli-2 loci. While contigs consisting of as many as 21 BAC clones and containing up to 17 -gliadin genes were formed, many BAC clones remained as singletons. The accuracy of the order of BAC clones in the contigs was verified by Southern hybridization analysis of the BAC fingerprints using an -gliadin probe. These results indicate that -gliadin genes are not evenly dispersed in the Gli-2 locus regions. Hybridization of these BACs with probes for long terminal repeat retrotransposons was used to determine the abundance and distribution of repetitive DNA in this region. Sequencing of BAC ends indicated that 70% of the sequences were significantly similar to different classes of retrotransposons, suggesting that these elements are abundant in this region. Several mechanisms underlying the dynamic evolution of the Gli-2 loci are discussed.     

Mesoscale organization of domains in the plasma membrane – beyond the lipid raft

The plasma membrane is compartmentalized into several distinct regions or domains, which show a broad diversity in both size and lifetime. The segregation of lipids and membrane proteins is thought to be driven by the lipid composition itself, lipid–protein interactions and diffusional barriers. With regards to the lipid composition, the immiscibility of certain classes of lipids underlies the “lipid raft” concept of plasmalemmal compartmentalization. Historically, lipid rafts have been described as cholesterol and (glyco)sphingolipid-rich regions of the plasma membrane that exist as a liquid-ordered phase that are resistant to extraction with non-ionic detergents. Over the years the interest in lipid rafts grew as did the challenges with studying these nanodomains. The term lipid raft has fallen out of favor with many scientists and instead the terms “membrane raft” or “membrane nanodomain” are preferred as they connote the heterogeneity and dynamic nature of the lipid-protein assemblies. In this article, we will discuss the classical lipid raft hypothesis and its limitations. This review will also discuss alternative models of lipid-protein interactions, annular lipid shells, and larger membrane clusters. We will also discuss the mesoscale organization of plasmalemmal domains including visible structures such as clathrin-coated pits and caveolae.     

Biophysics of α-synuclein membrane interactions

Membrane proteins participate in nearly all cellular processes; however, because of experimental limitations, their characterization lags far behind that of soluble proteins. Peripheral membrane proteins are particularly challenging to study because of their inherent propensity to adopt multiple and/or transient conformations in solution and upon membrane association. In this review, we summarize useful biophysical techniques for the study of peripheral membrane proteins and their application in the characterization of the membrane interactions of the natively unfolded and Parkinson's disease (PD) related protein, α-synuclein (α-syn). We give particular focus to studies that have led to the current understanding of membrane-bound α-syn structure and the elucidation of specific membrane properties that affect α-syn-membrane binding. Finally, we discuss biophysical evidence supporting a key role for membranes and α-syn in PD pathogenesis. This article is part of a Special Issue entitled: Membrane protein structure and function.     

Orientational oscillations of the peptide groups of an α-helix

The low-energy orientational oscillations of the peptide groups of an -helix are considered and the value of the frequency is estimated to be in agreement with experiments. Approximate formulae are derived for the projection of a dipole moment on the helix axis and for the helix parameters. Within the framework of a three-chain model, the asymptotics of the soliton solution is obtained using a discrete approach.The analysis of -helix geometry exhibits two types of low-frequency oscillations of the -helix. The first one is connected with atom movements along the helix axis with the peptide groups twisting around the helix axis. Accordingly, it changes the hydrogen bond lengths between neighbouring peptide groups. In the second case, the slopes of the peptide groups to the helix axis oscillate without the helix parameters changing. Here, the energy of interactions between peptide-group dipoles is changed and, as a result, the oscillations have an optical nature. The frequency of the optical orientational oscillations is approximately 100 cm^-1.     

Simulation of the packing of idealized transmembrane α-helix bundles

The aim of this study is to investigate if the packing motifs of native transmembrane helices can be produced by simulations with simple potentials and to develop a method for the rapid generation of initial candidate models for integral membrane proteins composed of bundles of transmembrane helices. Constituent residues are mapped along the helix axis in order to maintain the amino acid sequence-dependent properties of the helix. Helix packing is optimized according to a semi-empirical potential mainly composed of four components: a bilayer potential, a crossing angle potential, a helix dipole potential and a helix-helix distance potential. A Monte Carlo simulated annealing protocol is employed to optimize the helix bundle system. Necessary parameters are derived from theoretical studies and statistical analysis of experimentally determined protein structures. Preliminary testing of the method has been conducted with idealized seven Ala₂₀ helix bundles. The structures generated show a high degree of compactness. It was observed that both bacteriorhodopsin-like and δ-endotoxin-like structures are generated in seven-helix bundle simulations, within which the composition varies dependent upon the cooling rate. The simulation method has also been employed to explore the packing of N = 4 and N = 12 transmembrane helix bundles. The results suggest that seven and 12 transmembrane helix bundles resembling those observed experimentally (e.g., bacteriorhodopsin, rhodopsin and cytochrome c oxidase subunit I) may be generated by simulations using simple potentials. Received: 16 November 1998 / Revised version: 26 March 1999 / Accepted: 8 April 1999     

Accurate prediction of the burial status of transmembrane residues of α-helix membrane protein by incorporating the structural and physicochemical features

Predicting the burial status (the residue exposure to the lipid bilayer or buried within the protein core) of transmembrane (TM) residues of α-helix membrane protein (αHMP) is of great importance for genome-wide annotation and for experimental researchers to elucidate diverse physiological processes. In this work, we developed a new computational model that can be used for predicting the burial status of TM residues of αHMP. By incorporating physicochemical scales and conservation index, an efficient prediction model using least squares support vector machine (LS-SVM) was developed. The model was developed from 43 protein chains and its prediction ability was evaluated by an independent test set of other non-redundant ten protein chains. The prediction accuracy of our method was much better than the results of the reported works. Our results demonstrate that the LS-SVM prediction model incorporating structural and physicochemical features derived from sequence information could greatly improve the prediction accuracy.     

Pressure-dependent changes in the structure of the melittin α-helix determined by NMR

A novel method is described, which uses changes in NMR chemical shifts to characterise the structural change in a protein with pressure. Melittin in methanol is a small -helical protein, and its chemical shifts change linearly and reversibly with pressure between 1 and 2000 bar. An improved relationship between structure and HN shift has been calculated, and used to drive a molecular dynamics-based calculation of the change in structure. With pressure, the helix is compressed, with the H—O distance of the NH—O=C hydrogen bonds decreased by 0.021 ± 0.039 Å, leading to an overall compression along the entire helix of about 0.4 Å, corresponding to a static compressibility of 6 ×10^–6 bar^–1. The backbone dihedral angles and are altered by no more than ± 3° for most residues with a negative correlation coefficient of –0.85 between _i and _i–1, indicating that the local conformation alters to maintain hydrogen bonds in good geometries. The method is shown to be capable of calculating structural change with high precision, and the results agree with structural changes determined using other methodologies.     

Moving Iron through ferritin protein nanocages depends on residues throughout each four α-helix bundle subunit

Eukaryotic H ferritins move iron through protein cages to form biologically required, iron mineral concentrates. The biominerals are synthesized during protein-based Fe²⁺/O₂ oxidoreduction and formation of [Fe³⁺O]_n multimers within the protein cage, en route to the cavity, at sites distributed over ∼50 Å. Recent NMR and Co²⁺-protein x-ray diffraction (XRD) studies identified the entire iron path and new metal-protein interactions: (i) lines of metal ions in 8 Fe²⁺ ion entry channels with three-way metal distribution points at channel exits and (ii) interior Fe³⁺O nucleation channels. To obtain functional information on the newly identified metal-protein interactions, we analyzed effects of amino acid substitution on formation of the earliest catalytic intermediate (diferric peroxo-A_{650 nm}) and on mineral growth (Fe³⁺O-A_{350 nm}), in A26S, V42G, D127A, E130A, and T149C. The results show that all of the residues influenced catalysis significantly (p < 0.01), with effects on four functions: (i) Fe²⁺ access/selectivity to the active sites (Glu¹³⁰), (ii) distribution of Fe²⁺ to each of the three active sites near each ion channel (Asp¹²⁷), (iii) product (diferric oxo) release into the Fe³⁺O nucleation channels (Ala²⁶), and (iv) [Fe³⁺O]_n transit through subunits (Val⁴², Thr¹⁴⁹). Synthesis of ferritin biominerals depends on residues along the entire length of H subunits from Fe²⁺ substrate entry at 3-fold cage axes at one subunit end through active sites and nucleation channels, at the other subunit end, inside the cage at 4-fold cage axes. Ferritin subunit-subunit geometry contributes to mineral order and explains the physiological impact of ferritin H and L subunits.     

Role of terminal dipole charges in aggregation of α-helix pair in the voltage gated K+ channel

The voltage sensor domain (VSD) of the potassium ion channel KvAP is comprised of four (S1–S4) α-helix proteins, which are encompassed by several charged residues. Apart from these charges, each peptide α-helix having two inherent equal and opposite terminal dipolar charges behave like a macrodipole. The activity of voltage gated ion channel is electrostatic, where all the charges (charged residues and dipolar terminal charges) interact with each other and with the transmembrane potential. There are evidences that the role of the charged residues dominate the stabilization of the conformation and the gating process of the ion channel, but the role of the terminal dipolar charges are never considered in such analysis. Here, using electrostatic theory, we have studied the role of the dipolar terminal charges in aggregation of the S3b–S4 helix pair of KvAP in the absence of any external field (V = 0). A system attains stability, when its potential energy reaches minimum values. We have shown that the presence of terminal dipole charges (1) change the total potential energy of the charges on S3b–S4, affecting the stabilization of the α-helix pair within the bilayer lipid membrane and (2) the C- and the N-termini of the α-helices favor a different dielectric medium for enhanced stability. Thus, the dipolar terminal charges play a significant role in the aggregation of the two neighboring α-helices.