Changes in dynamics upon oligomerization regulate substrate binding and allostery in amino acid kinase family members

Oligomerization is a functional requirement for many proteins. The interfacial interactions and the overall packing geometry of the individual monomers are viewed as important determinants of the thermodynamic stability and allosteric regulation of oligomers. The present study focuses on the role of the interfacial interactions and overall contact topology in the dynamic features acquired in the oligomeric state. To this aim, the collective dynamics of enzymes belonging to the amino acid kinase family both in dimeric and hexameric forms are examined by means of an elastic network model, and the softest collective motions (i.e., lowest frequency or global modes of motions) favored by the overall architecture are analyzed. Notably, the lowest-frequency modes accessible to the individual subunits in the absence of multimerization are conserved to a large extent in the oligomer, suggesting that the oligomer takes advantage of the intrinsic dynamics of the individual monomers. At the same time, oligomerization stiffens the interfacial regions of the monomers and confers new cooperative modes that exploit the rigid-body translational and rotational degrees of freedom of the intact monomers. The present study sheds light on the mechanism of cooperative inhibition of hexameric N-acetyl-L-glutamate kinase by arginine and on the allosteric regulation of UMP kinases. It also highlights the significance of the particular quaternary design in selectively determining the oligomer dynamics congruent with required ligand-binding and allosteric activities.     

Intermonomer electron transfer in the bc1 complex dimer is controlled by the energized state and by impaired electron transfer between low and high potential hemes

The cytochrome bc(1) complex (commonly called Complex III) is the central enzyme of respiratory and photosynthetic electron transfer chains. X-ray structures have revealed the bc(1) complex to be a dimer, and show that the distance between low potential (b(L)) and high potential (b(H)) hemes, is similar to the distance between low potential hemes in different monomers. This suggests that electron transfer between monomers should occur at the level of the b(L) hemes. Here, we show that although the rate constant for b(L)-->b(L) electron transfer is substantial, it is slow compared to the forward rate from b(L) to b(H), and the intermonomer transfer only occurs after equilibration within the first monomer. The effective rate of intermonomer transfer is about 2-orders of magnitude slower than the direct intermonomer electron transfer.     

Association/Dissociation of the Nucleotide-binding Domains of the ATP-binding Cassette Protein MsbA Measured during Continuous Hydrolysis

In ATP-binding cassette proteins, the two nucleotide-binding domains (NBDs) work as dimers to bind and hydrolyze ATP, but the molecular mechanism of nucleotide hydrolysis is controversial. It is still unresolved whether hydrolysis leads to dissociation of the ATP-induced dimers or partial opening of the dimers such that the NBDs remain in contact during the hydrolysis cycle. We studied the bacterial lipid flippase MsbA by luminescence resonance energy transfer (LRET). The LRET signal between optical probes reacted with single-cysteine mutants was employed to follow NBD association/dissociation in real time. The intermonomer distances calculated from LRET data indicate that the NBDs separate completely following ATP hydrolysis, even in the presence of mm MgATP, and that the dissociation occurs following each hydrolysis cycle. The results support association/dissociation, as opposed to constant contact models, for the mode of operation of ATP-binding cassette proteins.     

Dissociation of ATP-binding cassette nucleotide-binding domain dimers into monomers during the hydrolysis cycle

ATP-binding cassette (ABC) proteins have two nucleotide-binding domains (NBDs) that work as dimers to bind and hydrolyze ATP, but the molecular mechanism of nucleotide hydrolysis is controversial. In particular, it is still unresolved whether hydrolysis leads to dissociation of the ATP-induced dimers or opening of the dimers, with the NBDs remaining in contact during the hydrolysis cycle. We studied a prototypical ABC NBD, the Methanococcus jannaschii MJ0796, using spectroscopic techniques. We show that fluorescence from a tryptophan positioned at the dimer interface and luminescence resonance energy transfer between probes reacted with single-cysteine mutants can be used to follow NBD association/dissociation in real time. The intermonomer distances calculated from luminescence resonance energy transfer data indicate that the NBDs separate completely following ATP hydrolysis, instead of opening. The results support ABC protein NBD association/dissociation, as opposed to constant-contact models.     

Intermonomer electron transfer between the low-potential b hemes of cytochrome bc₁

Cytochrome (cyt) bc(1) is a structural dimer with its monomers consisting of the Fe-S protein, cyt b, and cyt c(1) subunits. Its three-dimensional architecture depicts it as a symmetrical homodimer, but the mobility of the head domain of the Fe-S protein indicates that the functional enzyme exists in asymmetrical heterodimeric conformations. Here, we report a new genetic system for studying intra- and intermonomer interactions within the cyt bc(1) using the facultative phototrophic bacterium Rhodobacter capsulatus. The system involves two different sets of independently expressed cyt bc(1) structural genes carried by two plasmids that are coharbored by a cell without its endogenous enzyme. Our results indicate that coexpressed cyt bc(1) subunits were matured, assorted, and assembled in vivo into homo- and heterodimeric enzymes that can bear different mutations in each monomer. Using the system, the occurrence of intermonomer electron transfer between the low-potential b hemes of cyt bc(1) was probed by choosing mutations that perturb electron transfer at the hydroquinone oxidation (Q(o)) and quinone reduction (Q(i)) sites of the enzyme. The data demonstrate that active heterodimeric variants, formed of monomers carrying mutations that abolish only one of the two (Q(o) or Q(i)) active sites of each monomer, are produced, and they support photosynthetic growth of R. capsulatus. Detailed analyses of the physicochemical properties of membranes of these mutants, as well as purified homo- and heterodimeric cyt bc(1) preparations, demonstrated that efficient and productive electron transfer occurs between the low-potential b(L) hemes of the monomers in a heterodimeric enzyme. Overall findings are discussed with respect to intra- and intermonomer interactions that take place during the catalytic turnover of cyt bc(1).     

Recruiting monomer for dimer formation: resolving the antagonistic mechanisms of novel immune check point inhibitors against Programmed Death Ligand-1 in cancer immunotherapy

The design of small molecule antagonists against Programmed Death Ligand-1 (PD-L1) has been the recent highlight of the immune checkpoint blockade therapy. This interventive approach has been potentiated by the development of BMS compounds; BMS-1001 and BMS-1166, which exert their therapeutic activities by inducing dimerisation of PD-L1; a molecular mechanism that has remained unclear. For the first time, we resolve the dynamical events that underlie the antagonistic mechanisms of BMS-1001 and BMS-1166 when bound to PD-L1 using an all-atom molecular dynamics (MD) simulations approach and free binding energy Molecular Mechanics/Poisson-Boltzmann Surface Area (MM/PBSA) calculations. Time-scale dynamical findings revealed that upon binding a PD-L1 monomer, the BMS-compounds gradually facilitated the ‘inbound’ motion of another PD-L1 monomer in the same conformational phase space up till dimer formation. Moreover, the non-liganded PD-L1 monomer exhibited the highest structural flexibility and atomistic motions relative to the BMS-liganded monomer as revealed by post-MD trajectory analyses using root mean square deviation (RMSD) and root mean square fluctuations (RMSF) parameters. Trajectory investigations into ligand motions also revealed that the BMS compounds exhibited mechanistic transitions from the monomeric binding site (monomer A) where they were initially bound, to the second monomeric site (monomer B) where they were strongly bound, followed by eventual high-affinity interactions at the tunnel-like binding cleft formed upon the dimerisation of both PD-L1 monomers. These findings present a model that describes the mechanism by which the BMS compounds induce PD-L1 dimerisation and could further enhance the design of highly selective and novel monomeric recruiters of PD-L1 in cancer immunotherapy.     

Steric interactions of valines 1, 5, and 7 in [valine 5, D-alanine 8] gramicidin A channels.

When the central valine residues 6, 7, and 8 of gramicidin A (gA) are shifted by one position, the resulting [Val(5), D-Ala(8)]gA forms right-handed channels with a single-channel conductance and average duration somewhat less than gA channels. The reduction in channel duration has been attributed to steric conflict between the side chains of Val(1) and Val(5) in opposing monomers (Koeppe, R. E. II, D. V. Greathouse, A. Jude, G. Saberwal, L. L. Providence, and O. S. Andersen. 1994. J. Biol. Chem. 269:12567-12576). To investigate the orientations and motions of valines in [Val(5), D-Ala(8)]gA, we have incorporated (2)H labels at Val 1, 5, or 7 and recorded (2)H-NMR spectra of oriented and nonoriented samples in hydrated dimyristoylphosphatidylcholine. Spectra of nonoriented samples at 4 degrees C reveal powder patterns that indicate rapid side chain "hopping" for Val(5), and an intermediate rate of hopping for Val(1) and Val(7) that is somewhat slower than in gA. Oriented samples of deuterated Val(1) and Val(7) show large changes in the methyl and C(beta)-(2)H quadrupolar splittings (Deltanu(q)) when Ala(5) in native gA is changed to Val(5). Three or more peaks for the Val(1) methyls with Deltanu(q) values that vary with the echo delay, together with an intermediate spectrum for nonoriented samples at 4 degrees C, suggest unusual side chain dynamics for Val(1) in [Val(5), D-Ala(8)]gA. These results are consistent with a steric conflict that has been introduced between the two opposing monomers. In contrast, the acylation of gA has little influence on the side chain dynamics of Val(1), regardless of the identity of residue 5.     

Common mechanism of pore opening shared by five different potassium channels

A fundamental question associated with the function of ion channels is the conformational changes that allow for reversibly opening/occluding the pore through which the cations permeate. The recently elucidated crystal structures of potassium channels reveal similar structural motifs at their pore-forming regions, suggesting that they share a common gating mechanism. The validity of this hypothesis is explored by analyzing the collective dynamics of five known K(+) channel structures. Normal-mode analysis using the Gaussian network model strikingly reveals that all five structures display the same intrinsic motions at their pore-forming region despite the differences in their sequences, structures, and activation mechanisms. Superposition of the most cooperative mode profiles shows that the identified common mechanism is a global corkscrew-like counterrotation of the extracellular and cytoplasmic (CP) regions, leading to the opening of the CP end of the pore. A second cooperative mode shared by all five K(+) channels is the extension of the extracellular and/or CP ends via alternating anticorrelated fluctuations of pairs of diagonally opposite monomers. Residues acting as hinges/anchors in both modes are highly conserved across the members of the family of K(+) channel proteins, consistent with their presently disclosed critical mechanical role in pore gating.     

Proton transfer in gramicidin channels is modulated by the thickness of monoglyceride bilayers

The thickness of monoglyceride planar bilayers has significant effects on the transfer of protons in both native gramicidin A (gA) and in covalently linked SS- and RR-dioxolane-linked gA proteins. Planar bilayers with various thicknesses were formed from an appropriate combination of monoglyceride with various fatty acid lengths and solvent. Bilayer thicknesses ranged from 25 A (monoolein in squalene) to 54 A (monoeicosenoin in decane). Single-channel conductances to protons (g(H)) were measured in the concentration range of 10-5000 mM HCl. In native gA as well as in RR channels, the shape of the log(g(H))-log([H(+)]) relationships was nonlinear and remained basically unaltered in monoglyceride bilayers with various thicknesses. For both native gA and RR channels, g(H) values were systematically and significantly larger in thin than in thick bilayers. By contrast, the shape of the log(g(H))-log([H(+)]) relationships in the SS channel was linear (with a slope considerably smaller than 1) in thick (>37 A) bilayers. However, in thin (<37 A) bilayers these plots became nonlinear and g(H) values approached those obtained in native gA channels. The linearization of the log-log plots in the SS channel in thick bilayers is a consequence of a dramatic increase (instead of a decrease as in native gA and RR channels) of g(H) in these bilayers in [H(+)] <1 M. The gating characteristics of the various gA channels as a function of bilayer thickness followed the same pattern as described previously. It was noticed, however, that in the thickest monoglyceride bilayer used in this study, both the SS- and RR-dioxolane-linked channels opened in a mode of bursting activity instead of remaining in the open state as in thin bilayers. It is proposed that the thickness of monoglyceride bilayers modulates proton transfer in native gA channels by a combination of factors including the access resistances of channels to H(+), and fluctuations in both the structure of the lipid bilayer and in the distance between gA monomers. The differential effects of relatively thick monoglyceride bilayers on proton transfer in both dioxolane-linked gA channels must relate to distinct interactions between the bilayers and the SS and RR dioxolanes.     

On the origin of closing flickers in gramicidin channels: a new hypothesis

The submillisecond closing events (flickers) and the single channel conductances to protons (g(H)) were studied in native gramicidin A (gA) and in the SS and RR diastereoisomers of dioxolane-linked gA channels in planar bilayers. Bilayers were formed from glycerylmonooleate (GMO) in various solvents. In GMO/decane (thick) bilayers, the largest flicker frequency occurred in the SS channel (39 s(-1)), followed by the RR (4 s(-1)) and native gA channels (3 s(-1)). These frequencies were attenuated in GMO/squalene (thin) bilayers by 100-, 30-, and 70-fold in the SS, RR, and native gA channels, respectively. In thin bilayers, the average burst duration of native gA channels was 30-fold longer than in thick bilayers. The RR dioxolane-linked gA dimer "inactivated" in GMO/decane but not in squalene-containing bilayers. The mean closed time of flickers (approximately 0.12 ms) was essentially the same in various gA channels. In thin bilayers, g(H) values were larger by approximately 10% (SS), 30% (RR), and 20% (native gA) in relation to thick bilayers. It is concluded that flickers are not related to pre-dissociation or dissociation states of gA monomers, and do not seem to be caused by intrinsic conformational changes of channel proteins. It is proposed that flickers are caused by undulations of the bilayer that obliterate the openings of gA channels. Differences between flicker frequencies in various gA channels are likely to result from variations in channel geometries at the bilayer/channel interface. The smaller g(H) in thick bilayers suggests that the deformation of these bilayers around the gA channel creates a diffusional pathway next to the mouths of the channel that is longer and more restrictive than in thin GMO bilayers. A possible molecular interpretation for these effects is attempted.     

Curcumin is a modulator of bilayer material properties

Curcumin (1,7-bis(4-hydroxy-3-methoxyphenyl)-1,6-heptadiene-3,5-dione) is the major bioactive compound in turmeric (Curcuma longa) with antioxidant, antiinflammatory, anticarcinogenic, and antimutagenic effects. At low muM concentrations, curcumin modulates many structurally and functionally unrelated proteins, including membrane proteins. Because the cell membranes' lipid bilayer serves as a gate-keeper and regulator of many cell functions, we explored whether curcumin modifies general bilayer properties using channels formed by gramicidin A (gA). gA channels form when two monomers from opposing monolayers associate to form a conducting dimer with a hydrophobic length that is less than the bilayer hydrophobic thickness; gA channel formation thus causes a local bilayer thinning. The energetic cost of this bilayer deformation alters the gA monomer <--> dimer equilibrium, which makes the channels' appearance rate and lifetime sensitive to changes in bilayer material properties, and the gA channels become probes for changes in bilayer properties. Curcumin decreases bilayer stiffness, increasing both gA channel lifetimes and appearance rates, meaning that the energetic cost of the gA-induced bilayer deformation is reduced. These results show that curcumin may exert some of its effects on a diverse range of membrane proteins through a bilayer-mediated mechanism.     

The membrane interface dictates different anchor roles for "inner pair" and "outer pair" tryptophan indole rings in gramicidin A channels

We investigated the effects of substituting two of the four tryptophans (the "inner pair" Trp(9) and Trp(11) or the "outer pair" Trp(13) and Trp(15)) in gramicidin A (gA) channels. The conformational preferences of the doubly substituted gA analogues were assessed using circular dichroism spectroscopy and size-exclusion chromatography, which show that the inner tryptophans 9 and 11 are critical for the gA's conformational preference in lipid bilayer membranes. [Phe(13,15)]gA largely retains the single-stranded helical channel structure, whereas [Phe(9,11)]gA exists primarily as double-stranded conformers. Within this context, the (2)H NMR spectra from labeled tryptophans were used to examine the changes in average indole ring orientations, induced by the Phe substitutions and by the shift in conformational preference. Using a method for deuterium labeling of already synthesized gAs, we introduced deuterium selectively onto positions C2 and C5 of the remaining tryptophan indole rings in the substituted gA analogues for solid-state (2)H NMR spectroscopy. The (least possible) changes in orientation and overall motion of each indole ring were estimated from the experimental spectra. Regardless of the mixture of backbone folds, the indole ring orientations observed in the analogues are similar to those found previously for gA channels. Both Phe-substituted analogues form single-stranded channels, as judged from the formation of heterodimeric channels with the native gA. [Phe(13,15)]gA channels have Na(+) currents that are ~50% and lifetimes that are ~80% of those of native gA channels. The double-stranded conformer(s) of [Phe(9,11)]gA do not form detectable channels. The minor single-stranded population of [Phe(9,11)]gA forms channels with Na(+) currents that are ~25% and single-channel lifetimes that are ~300% of those of native gA channels. Our results suggest that Trp(9) and Trp(11), when "reaching" for the interface, tend to drive both monomer folding (to "open" a channel) and dimer dissociation (to "close" a channel). Furthermore, the dipoles of Trp(9) and Trp(11) are relatively more important for the single-channel conductance than are the dipoles of Trp(13) and Trp(15).     

Coarse-grained modeling of the actin filament derived from atomistic-scale simulations

A coarse-grained (CG) procedure that incorporates the information obtained from all-atom molecular dynamics (MD) simulations is presented and applied to actin filaments (F-actin). This procedure matches the averaged values and fluctuations of the effective internal coordinates that are used to define a CG model to the values extracted from atomistic MD simulations. The fluctuations of effective internal coordinates in a CG model are computed via normal-mode analysis (NMA), and the computed fluctuations are matched with the atomistic MD results in a self-consistent manner. Each actin monomer (G-actin) is coarse-grained into four sites, and each site corresponds to one of the subdomains of G-actin. The potential energy of a CG G-actin contains three bonds, two angles, and one dihedral angle; effective harmonic bonds are used to describe the intermonomer interactions in a CG F-actin. The persistence length of a CG F-actin was found to be sensitive to the cut-off distance of assigning intermonomer bonds. Effective harmonic bonds for a monomer with its third nearest neighboring monomers are found to be necessary to reproduce the values of persistence length obtained from all-atom MD simulations. Compared to the elastic network model, incorporating the information of internal coordinate fluctuations enhances the accuracy and robustness for a CG model to describe the shapes of low-frequency vibrational modes. Combining the fluctuation-matching CG procedure and NMA, the achievable time- and length scales of modeling actin filaments can be greatly enhanced. In particular, a method is described to compute the force-extension curve using the CG model developed in this work and NMA. It was found that F-actin is easily buckled under compressive deformation, and a writhing mode is developed as a result. In addition to the bending and twisting modes, this novel writhing mode of F-actin could also play important roles in the interactions of F-actin with actin-binding proteins and in the force-generation process via polymerization.     

Structural Basis of pH Dependence of Neoculin,a Sweet Taste-Modifying Protein

Among proteins utilized as sweeteners, neoculin and miraculin are taste-modifying proteins that exhibit pH-dependent sweetness. Several experiments on neoculin have shown that His11 of neoculin is responsible for pH dependence. We investigated the molecular mechanism of the pH dependence of neoculin by molecular dynamics (MD) calculations. The MD calculations for the dimeric structures of neoculin and His11 mutants showed no significant structural changes for each monomer at neutral and acidic pH levels. The dimeric structure of neoculin dissociated to form isolated monomers under acidic conditions but was maintained at neutral pH. The dimeric structure of the His11Ala mutant, which is sweet at both neutral and acidic pH, showed dissociation at both pH 3 and 7. The His11 residue is located at the interface of the dimer in close proximity to the Asp91 residue of the other monomer. The MD calculations for His11Phe and His11Tyr mutants demonstrated the stability of the dimeric structures at neutral pH and the dissociation of the dimers to isolated monomers. The dissociation of the dimer caused a flexible backbone at the surface that was different from the dimeric interface at the point where the other monomer interacts to form an oligomeric structure. Further MD calculations on the tetrameric structure of neoculin suggested that the flexible backbone contributed to further dissociation of other monomers under acidic conditions. These results suggest that His11 plays a role in the formation of oligomeric structures at pH 7 and that the isolated monomer of neoculin at acidic pH is responsible for sweetness.     

Functionally relevant conformational dynamics of water-soluble proteins

A study of the functional-relevant dynamics of three typical water-soluble proteins, including Calmodulin, Src-tyrosine kinase, and a repressor of the trip operon, has been reported. The application of state-of-art methods of structural bioinformatics allowed us to identify the dynamics seen in the X-ray structures of the investigated proteins associated with their specific biological functions. In addition, technique of normal mode analysis reveals the most probable directions of the functionally relevant motions for all proteins. Importantly, the overall type of the motions observed in the lowest-frequency modes was very similar to the motions seen from the analysis of the X-ray data of the examined macromolecules. Thus, it was shown that the large-scale, as well as local, conformational motions of the proteins might already be predetermined at the level of their tertiary structures. In particular, the determining factor might be the specific fold of the α-helices. Thus, the functionally relevant in vivo dynamics of the investigated proteins might be evolutionarily formed by natural selection at the level of the spatial topology.     

Inherent Dynamics of the Acid-Sensing Ion Channel 1 Correlates with the Gating Mechanism

The acid-sensing ion channel 1 (ASIC1) is a key receptor for extracellular protons. Although numerous structural and functional studies have been performed on this channel, the structural dynamics underlying the gating mechanism remains unknown. We used normal mode analysis, mutagenesis, and electrophysiological methods to explore the relationship between the inherent dynamics of ASIC1 and its gating mechanism. Here we show that a series of collective motions among the domains and subdomains of ASIC1 correlate with its acid-sensing function. The normal mode analysis result reveals that the intrinsic rotation of the extracellular domain and the collective motions between the thumb and finger induced by proton binding drive the receptor to experience a deformation from the extracellular domain to the transmembrane domain, triggering the channel pore to undergo “twist-to-open” motions. The movements in the transmembrane domain indicate that the likely position of the channel gate is around Leu440. These motion modes are compatible with a wide body of our complementary mutations and electrophysiological data. This study provides the dynamic fundamentals of ASIC1 gating.     

Peptide-induced membrane elastic deformations decelerate gramicidin dimer-monomer equilibration

Gramicidin A (gA) is a hydrophobic pentadecapeptide readily incorporating into a planar bilayer lipid membrane (BLM), thereby inducing a large macroscopic current across the BLM. This current results from ion-channel formation due to head-to-head transbilayer dimerization of gA monomers with rapidly established monomer-dimer equilibrium. Any disturbance of the equilibrium, e.g., by sensitized photoinactivation of a portion of gA monomers, causes relaxation toward a new equilibrium state. According to previous studies, the characteristic relaxation time of the gA-mediated electric current decreases as the current increases upon elevating the gA concentration in the membrane. Here, we report data on the current relaxation kinetics for gA analogs with N-terminal valine replaced by glycine or tyrosine. Surprisingly, the relaxation time increased rather than decreased upon elevation of the total membrane conductance induced by these gA analogs, thus contradicting the classical kinetic scheme. We developed a general theoretical model that accounts for lateral interaction of monomers and dimers mediated by membrane elastic deformations. The modified kinetic scheme of the gramicidin dimerization predicts the reverse dependence of the relaxation time on membrane conductance for gA analogs, with a decreased dimerization constant that is in a good agreement with our experimental data. The equilibration process may be also modulated by incorporation of other peptides (“impurities”) into the lipid membrane.     

Fluorescent intercalator displacement replacement (FIDR) assay: determination of relative thermodynamic and kinetic parameters in triplex formation--a case study using triplex-forming LNAs

Triplex forming oligonucleotides (TFOs) are the most commonly used approach for site-specific targeting of double stranded DNA (dsDNA). Important parameters describing triplex formation include equilibrium binding constants (K(eq)) and association/dissociation rate constants (k(on) and k(off)). The 'fluorescent intercalator displacement replacement' (FIDR) assay is introduced herein as an operationally simple approach toward determination of these parameters for triplexes involving TC-motif TFOs. Briefly described, relative rate constants are determined from fluorescence intensity changes upon: (i) TFO-mediated displacement of pre-intercalated and fluorescent ethidium from dsDNA targets (triplex association) and (ii) Watson-Crick complement-mediated displacement of the TFO and replacement with ethidium (triplex dissociation). The assay is used to characterize triplexes between purine-rich dsDNA targets and TC-motif TFOs modified with six different locked nucleic acid (LNA) monomers, i.e. conventional and C5-alkynyl-functionalized LNA and α-L-LNA pyrimidine monomers. All of the studied monomers increase triplex stability by decreasing the triplex dissociation rate. LNA-modified TFOs form more stable triplexes than α-L-LNA-modified counterparts owing to slower triplex dissociation. Triplexes modified with C5-(3-aminopropyn-1-yl)-LNA-U monomer Z are particularly stable. The study demonstrates that three affinity-enhancing features can be combined into one high-affinity TFO monomer: conformational restriction of the sugar ring, expansion of the pyrimidine π-stacking surface and introduction of an exocyclic amine.     

Formation of non-beta 6.3-helical gramicidin channels between sequence-substituted gramicidin analogues.

Using the linear gramicidins as an example, we have previously shown how the statistical properties of heterodimeric (hybrid) channels (formed between the parent [Val1]gramicidin A (gA) and a sequence-altered analogue) can be used to assess whether the analogue forms channels that are structurally equivalent to the parent channels (Durkin, J. T., R. E. Koeppe II, and O. S. Andersen. 1990. J. Mol. Biol. 211:221-234). Generally, the gramicidins are tolerant of amino acid sequence alterations. We report here an exception. The optically reversed analogue, gramicidin M- (gM-) (Heitz, F., G. Spach, and Y. Trudelle. 1982. Biophys. J. 40:87-89), forms channels that are the mirror-image of [Val1]gA channels; gM- should thus form no hybrid channels with analogues having the same helix sense as [Val1]gA. Surprisingly, however, gM- forms hybrid channels with the shortened analogues des-Val1-[Ala2]gA and des-Val1-gC, but these channels differ fundamentally from the parent channels: (a) the appearance rate of these heterodimers is only approximately 1/10 of that predicted from the random assortment of monomers into conducting dimers, indicating the existence of an energy barrier to their formation (e.g., monomer refolding into a new channel-forming conformation); and (b), once formed, the hybrid channels are stabilized approximately 1,000-fold relative to the parent channels. The increased stability suggests a structure that is joined by many hydrogen bonds, such as one of the double-stranded helical dimers shown to be adopted by gramicidins in organic solvents (Veatch, W. R., E. T. Fossel, and E. R. Blout. 1974. Biochemistry. 13:5249-5256).