The role of extra-membranous inter-helical loops in helix-helix interactions

Membrane Interactions of Cy3 and Cy5 Fluorophores and Their Effects on Membrane-Protein Dynamics

Organic fluorophores, such as Cy3 and Cy5, have been widely used as chemical labels to probe the structure and dynamics of membrane proteins. Although a number of previous studies have reported on the possibility of some of the water-soluble fluorophores to interact with lipid bilayers, detailed fluorophore-lipid interactions and, more importantly, the potential effect of such interactions on the natural dynamics of the labeled membrane proteins have not been well studied. We have performed a large set of all-atom molecular dynamics simulations employing the highly mobile membrane mimetic model to describe spontaneous partitioning of the fluorophores into lipid bilayers with different lipid compositions. Spontaneous membrane partitioning of Cy3 and Cy5 fluorophores captured in these simulations proceeds in two steps. Electrostatic interaction between the fluorophores and the lipid headgroups facilitates the initial, fast membrane association of the fluorophores, followed by slow insertion of hydrophobic moieties into the lipid bilayer core. After the conversion of the resulting membrane-bound systems to full-membrane representations, biased-exchange umbrella sampling simulations are performed for free energy calculations, revealing a higher energy barrier for partitioning into negatively charged (phosphatidylserine or phosphatidylcholine) membranes than purely zwitterionic (phosphatidylcholine or phosphatidylethanolamine) ones. Furthermore, the potential effect of fluorophore-lipid interactions on membrane proteins has been examined by covalently linking Cy5 to single- and multipass transmembrane helical proteins. Equilibrium simulations show strong position-dependent effects of Cy5-tagging on the structure and natural dynamics of membrane proteins. Interactions between the tagged protein and Cy5 were also observed. Our results suggest that fluorophore-lipid interactions can affect the structure and dynamics of membrane proteins to various extents, especially in systems with higher structural flexibility.

The soluble loop BC region guides, but not dictates, the assembly of the transmembrane cytochrome b6

Studying folding and assembly of naturally occurring α-helical transmembrane proteins can inspire the design of membrane proteins with defined functions. Thus far, most studies have focused on the role of membrane-integrated protein regions. However, to fully understand folding pathways and stabilization of α–helical membrane proteins, it is vital to also include the role of soluble loops. We have analyzed the impact of interhelical loops on folding, assembly and stability of the heme-containing four-helix bundle transmembrane protein cytochrome b₆ that is involved in charge transfer across biomembranes. Cytochrome b₆ consists of two transmembrane helical hairpins that sandwich two heme molecules. Our analyses strongly suggest that the loop connecting the helical hairpins is not crucial for positioning the two protein “halves” for proper folding and assembly of the holo-protein. Furthermore, proteolytic removal of any of the remaining two loops, which connect the two transmembrane helices of a hairpin structure, appears to also not crucially effect folding and assembly. Overall, the transmembrane four-helix bundle appears to be mainly stabilized via interhelical interactions in the transmembrane regions, while the soluble loop regions guide assembly and stabilize the holo-protein. The results of this study might steer future strategies aiming at designing heme-binding four-helix bundle structures, involved in transmembrane charge transfer reactions.

Molecular evolution of a widely-adopted taxonomic marker (COI) across the animal tree of life

DNA barcodes are widely used for identification and discovery of species. While such use draws on information at the DNA level, the current amassment of ca. 4.7 million COI barcodes also offers a unique resource for exploring functional constraints on DNA evolution. Here, we explore amino acid variation in a crosscut of the entire animal kingdom. Patterns of DNA variation were linked to functional constraints at the level of the amino acid sequence in functionally important parts of the enzyme. Six amino acid sites show variation with possible effects on enzyme function. Overall, patterns of amino acid variation suggest convergent or parallel evolution at the protein level connected to the transition into a parasitic life style. Denser sampling of two diverse insect taxa revealed that the beetles (Coleoptera) show more amino acid variation than the butterflies and moths (Lepidoptera), indicating fundamental difference in patterns of molecular evolution in COI. Several amino acid sites were found to be under notably strong purifying selection in Lepidoptera as compared to Coleoptera. Overall, these findings demonstrate the utility of the global DNA barcode library to extend far beyond identification and taxonomy, and will hopefully be followed by a multitude of work.

Distinct second extracellular loop structures of the brain cannabinoid CB(1) receptor: implication in ligand binding and receptor function

The G-protein-coupled receptor (GPCR) second extracellular loop (E2) is known to play an important role in receptor structure and function. The brain cannabinoid (CB(1)) receptor is unique in that it lacks the interloop E2 disulfide linkage to the transmembrane (TM) helical bundle, a characteristic of many GPCRs. Recent mutation studies of the CB(1) receptor, however, suggest the presence of an alternative intraloop disulfide bond between two E2 Cys residues. Considering the oxidation state of these Cys residues, we determine the molecular structures of the 17-residue E2 in the dithiol form (E2(dithiol)) and in the disulfide form (E2(disulfide)) of the CB(1) receptor in a fully hydrated 1-palmitoyl-2-oleoyl-sn-glycero-3-phosphocholine bilayer, using a combination of simulated annealing and molecular dynamics simulation approaches. We characterize the CB(1) receptor models with these two E2 forms, CB(1)(E2(dithiol)) and CB(1)(E2(disulfide)), by analyzing interaction energy, contact number, core crevice, and cross correlation. The results show that the distinct E2 structures interact differently with the TM helical bundle and uniquely modify the TM helical topology, suggesting that E2 of the CB(1) receptor plays a critical role in stabilizing receptor structure, regulating ligand binding, and ultimately modulating receptor activation. Further studies on the role of E2 of the CB(1) receptor are warranted, particularly comparisons of the ligand-bound form with the present ligand-free form.

Conformational and dynamics simulation study of antimicrobial peptide hedistin-heterogeneity of its helix-turn-helix motif

Hedistin is an antimicrobial peptide isolated from the coelomocytes of Nereis diversicolor, possessing activity against a large spectrum of bacteria including the methicillin resistant Staphylococcus aureus and Vibrio alginolyticus. The three-dimensional structure of hedistin in both aqueous solution and deuterated dodecylphosphocholine (DPC) micelles was examined using circular dichroism (CD) and nuclear magnetic resonance (NMR) techniques. And, the early events of the antibacterial process of hedistin were simulated using palmitoyl-oleoyl-phophatidylcholine (POPC) lipid bilayers and molecular dynamics (MD) simulation methods. Hedistin lacks secondary structure in aqueous solution, however, in DPC micelles, it features with a heterogeneous helix-turn-helix moiety and exhibits obvious amphipathic nature. The turn region (residues Val9-Thr12) in the moiety is a four-residue hinge, lying in between the first N-terminal alpha-helix (residues Leu5-Lys8) and the second alpha-helix (residues Val13-Ala17) regions and causing an approximately 120 degrees angle between the axes of the two helices. The segmental and nonlinear nature of hedistin structure is referred to as the heterogeneity of its helix-turn-helix motif which was found to be corresponding to a kind of discrete dynamics behavior, herein coined as its dynamical heterogeneity, at the early stage (0-50 ns) of the MD simulations. That is, the first helix segment, prior to (at 310 K) or following (at 363 K) the second helix, binds to the lipid head-group region and subsequently permeates into the hydrophobic lipid tail region, and the hinge is the last portion entering the lipid environment. This result implies that hedistin may adopt a "carpet" model action when disrupting bacterial membrane.

Transmembrane helical domain of the cannabinoid CB1 receptor

Brain cannabinoid (CB(1)) receptors are G-protein coupled receptors and belong to the rhodopsin-like subfamily. A homology model of the inactive state of the CB(1) receptor was constructed using the x-ray structure of beta(2)-adrenergic receptor (beta(2)AR) as the template. We used 105 ns duration molecular-dynamics simulations of the CB(1) receptor embedded in a 1-palmitoyl-2-oleoyl-sn-glycero-3-phosphocholine (POPC) bilayer to gain some insight into the structure and function of the CB(1) receptor. As judged from the root mean-square deviations combined with the detailed structural analyses, the helical bundle of the CB(1) receptor appears to be fully converged in 50 ns of the simulation. The results reveal that the helical bundle structure of the CB(1) receptor maintains a topology quite similar to the x-ray structures of G-protein coupled receptors overall. It is also revealed that the CB(1) receptor is stabilized by the formation of extensive, water-mediated H-bond networks, aromatic stacking interactions, and receptor-lipid interactions within the helical core region. It is likely that these interactions, which are often specific to functional motifs, including the S(N)LAxAD, D(E)RY, CWxP, and NPxxY motifs, are the molecular constraints imposed on the inactive state of the CB(1) receptor. It appears that disruption of these specific interactions is necessary to release the molecular constraints to achieve a conformational change of the receptor suitable for G-protein activation.

Stable interactions between the transmembrane domains of the adenosine A2A receptor

G-protein-coupled receptors (GPCRs) must properly insert and fold in the membrane to adopt a stable native structure and become biologically active. The interactions between transmembrane (TM) helices are believed to play a major role in these processes. Previous studies in our group showed that specific interactions between TM helices occur, leading to an increase in helical content, especially in weakly helical TM domains, suggesting that helix-helix interactions in addition to helix-lipid interactions facilitate helix formation. They also demonstrated that TM peptides interact in a similar fashion in micelles and lipid vesicles, as they exhibit relatively similar thermal stability and alpha-helicity inserted in SDS micelles to that observed in liposomes. In this study, we perform an analysis of pairwise interactions between peptides corresponding to the seven TM domains of the human A(2A) receptor (A(2A)R). We used a combination of Förster resonance energy transfer (FRET) measurement and circular dichroism (CD) spectroscopy to detect and analyze these interactions in detergent micelles. We found that strong and specific interactions occur in only seven of the 28 possible peptide pairs. Furthermore, not all interactions, identified by FRET, lead to a change in helicity. Our results identify stabilizing contacts that are likely related to the stability of the receptor and that are consistent with what is known about the three-dimensional structure and stability of rhodopsin and the beta(2) adrenergic receptor.

Binding mode prediction of conformationally restricted anandamide analogs within the CB1 receptor

BACKGROUND

CB1 cannabinoid receptors are G-protein coupled receptors for endocannabinoids including anandamide and 2-arachidonoylglycerol. Because these arachidonic acid metabolites possess a 20-carbon polyene chain as the alkyl terminal moiety, they are highly flexible with the potential to adopt multiple biologically relevant conformations, particularly those in a bent form. To better understand the molecular interactions associated with binding and steric trigger mechanisms of receptor activation, a series of conformationally-restricted anandamide analogs having a wide range of affinity and efficacy were evaluated.

RESULTS

A CB1 receptor model was constructed to include the extracellular loops, particularly extracellular loop 2 which possesses an internal disulfide linkage. Using both Glide (Schrödinger) and Affinity (Accelrys) docking programs, binding conformations of six anandamide analogs were identified that conform to rules applicable to the potent, efficacious and stereoselective non-classical cannabinoid CP55244. Calculated binding energies of the optimum structures from both procedures correlated well with the reported binding affinity values. The most potent and efficacious of the ligands adopted conformations characterized by interactions with both the helix-3 lysine and hydrophobic residues that interact with CP55244. The other five compounds formed fewer or less energetically favorable interactions with these critical residues. The flexibility of the tested anandamide analogs, measured by torsion angles around the benzene as well as the stretch between side chain moieties, could contribute to the differences in ability to interact with the CB1 receptor.

CONCLUSION

Analyses of multiple poses of conformationally-restricted anandamide analogs permitted identification of favored amino acid interactions within the CB1 receptor binding pocket. A ligand possessing both high affinity and cannabinoid agonist efficacy was able to interact with both polar and hydrophobic interaction sites utilized by the potent and efficacious non-classical cannabinoid CP55940. In contrast, other analogs characterized by reduced affinity or efficacy exhibited less favorable interactions with those key residues.

Molecular dynamics simulations suggest a mechanism for translocation of the HIV-1 TAT peptide across lipid membranes

The recombinant HIV-1 Tat protein contains a small region corresponding to residues (47)YGRKKRRQRR(57)R, which is capable of translocating cargoes of different molecular sizes, such as proteins, DNA, RNA, or drugs, across the cell membrane in an apparently energy-independent manner. The pathway that these peptides follow for entry into the cell has been the subject of strong controversy for the last decade. This peptide is highly basic and hydrophilic. Therefore, a central question that any candidate mechanism has to answer is how this highly hydrophilic peptide is able to cross the hydrophobic barrier imposed by the cell membrane. We propose a mechanism for the spontaneous translocation of the Tat peptides across a lipid membrane. This mechanism involves strong interactions between the Tat peptides and the phosphate groups on both sides of the lipid bilayer, the insertion of charged side chains that nucleate the formation of a transient pore, followed by the translocation of the Tat peptides by diffusing on the pore surface. This mechanism explains how key ingredients, such as the cooperativity among the peptides, the large positive charge, and specifically the arginine amino acids, contribute to the uptake. The proposed mechanism also illustrates the importance of membrane fluctuations. Indeed, mechanisms that involve large fluctuations of the membrane structure, such as transient pores and the insertion of charged amino acid side chains, may be common and perhaps central to the functions of many membrane protein functions.

The influence of protonation states on the dynamics of the NhaA antiporter from Escherichia coli

The crystal structure of NhaA Na(+)/H(+) antiporter of Escherichia coli has provided a basis to explore the mechanism of Na(+) and H(+) exchange and its regulation by pH. However, the dynamics and nature of the pH-induced changes in the proteins remained unknown. Using molecular mechanics methods, we studied the dynamic behavior of the hydrogen-bonded network in NhaA on shifting the pH from 4 to 8. The helical regions preserved the general architecture of NhaA throughout the pH change. In contrast, large conformational drifts occurred at pH 8 in the loop regions, and an increased flexibility of helix IVp was observed on the pH shift. A remarkable pH-induced conformational reorganization was found: at acidic pH helix X is slightly curved, whereas at alkaline pH, it is kinked around residue Lys(300). The barrier that exists between the cytoplasmic and periplasmic funnels at low pH is removed, and the two funnels are bridged by hydrogen bonds between water molecules and residues located in the TMSs IV/XI assembly and helix X at alkaline pH. In the variant Gly(338)Ser that lost pH control, a hydrogen-bonded chain between Ser(338) and Lys(300) was found to block the pH-induced conformational reorganization of helix X.