Conformer-dependent proton-transfer reactions of ubiquitin ions

The conformations of ubiquitin ions before and after being exposed to proton transfer reagents have been studied by using ion mobility/mass spectrometry techniques. Ions were produced by electrospray ionization and exposed to acetone, acetophenone, n-butylamine, and 7-methyl-1,5,7-triazabicyclo[4.4.0]dec-5-ene. Under the conditions employed, the +4 to +13 charge states were formed and a variety of conformations, which we have characterized as compact, partially folded, and elongated, have been observed. The low charge state ions have cross sections that are similar to those calculated for the crystal conformation. High charge states favor unfolded conformations. The ion mobility distributions recorded after ions have been exposed to each base show that the lowest charge state that is formed during proton-transfer reactions favors a compact conformation. More open conformations are observed for the higher charge states that remain after reaction. The results show that for a given charge state, the apparent gas-phase acidities of the different conformations are ordered as compact < partially folded < elongated.     

Evidence for unfolding and refolding of gas-phase cytochrome <Emphasis Type="Italic">c</Emphasis> ions in a Paul trap

The folding pathways of gas-phase cytochrome c ions produced by electrospray ionization have been studied by an ion trapping/ion mobility technique that allows conformations to be examined over extended timescales (10 ms to 10 s). The results show that the +9 charge state emerges from solution as a compact structure and then rapidly unfolds into several substantially more open structures, a transition that requires 30-60 ms; over substantially longer timescales (250 ms to 10 s) elongated states appear to refold into an array of folded structures. The new folded states are less compact than those that are apparent during the initial unfolding. Apparently, unfolding to highly open conformations is a key step that must occur before +9 ions can sample more compact states that are stable at longer times.     

The substituted alkyne 3-heptyne is eclipsed

Although butane exists in staggered anti and gauche conformations, when the ethyl groups are separated by a C[triple bond]C triple bond (3-hexyne), the stable conformation changes to eclipsed, (1)C(2v). Using rotational microwave spectroscopy, we have studied another example, 3-heptyne, the C[triple bond]C elongated analogue of pentane. The most stable conformer of pentane has anti-anti (AA) conformations about the central C-C bonds (C(2v)) and the next most stable has a gauche dihedral angle (GA, C1). This microwave study determines that the extended analogue of the AA form is not staggered about the C[triple bond]C axis but eclipsed (Cs). Also, the elongated analogue of the GA conformer is also not staggered but nearly eclipsed. The conformations of low-polarity substituted acetylenes is determined by dispersion attractions between the end groups. A microwave study of the AA and GA conformers of pentane is also reported.     

Helical molecular programming: folding of oligopyridine-dicarboxamides into molecular single helices

Molecular strands composed of alternating 2,6-diaminopyridine and 2,6-pyridinedicarbonyl units have been designed to self-organize into single stranded helical structures upon forming intramolecular hydrogen bonds. Pentameric strands 11, 12, and 14, heptameric strands 1 and 20, and undecameric strand 15 have been synthesized using stepwise convergent strategies. Single helical conformations have been characterized in the solid state by single crystal X-ray diffraction analysis for four of these compounds. Helices from pentameric strands 12 and 14 extend over one turn, and helices from heptameric 20 and undecameric 15 species extend to one and a half and two and a half turns, respectively. Intramolecular hydrogen bonds are responsible for the strong bending of the strands. 1H NMR shifts both in polar and nonpolar organic solvents indicate intramolecular overlap between the peripheral aromatic groups. Thus, helical conformations also predominate in solution. Molecular stochastic dynamic simulations of strand folding starting from a high energy extended linear conformer show a rapid (600 ps at 300 K) conversion into a stable helical conformation.     

Entropic stabilization of isolated beta-sheets

Temperature-dependent electric deflection measurements have been performed for a series of unsolvated alanine-based peptides (Ac-WA(n)-NH(2), where Ac = acetyl, W = tryptophan, A = alanine, and n = 3, 5, 10, 13, and 15). The measurements are interpreted using Monte Carlo simulations performed with a parallel tempering algorithm. Despite alanine's high helix propensity in solution, the results suggest that unsolvated Ac-WA(n)-NH(2) peptides with n > 10 adopt beta-sheet conformations at room temperature. Previous studies have shown that protonated alanine-based peptides adopt helical or globular conformations in the gas phase, depending on the location of the charge. Thus, the charge more than anything else controls the structure.     

A distance geometry heuristic for expanding the range of geometries sampled during conformational search

A recent study of crystal structures of protein-ligand complexes has shown that bioactive conformations tend to be more extended than random ones (Diller and Merz, J. Comput. Aid. Mol. Des. 2002, 16, 105-112). Existing conformational sampling techniques produce molecular conformations with a distribution of geometric sizes that may not cover that of the bioactive conformations. Here, we describe a simple heuristic for biasing the conformational search toward more extended or compact conformations, while maintaining excellent sampling. The method uses a boosting strategy to generate a series of conformations, each of which is at least as extended (or compact) as the previous one. We demonstrate that this method significantly expands the range of geometric sizes generated during the search and thus increases the efficiency of sampling bioactive conformations.     

Duplex formation and the onset of helicity in poly d(CG)n oligonucleotides in a solvent-free environment

The gas-phase conformations of a series of cytosine/guanine DNA duplexes were examined by ion mobility and molecular dynamics methods. Deprotonated duplex ions were formed by electrospray ionization, and their collision cross sections measured in helium were compared to calculated cross sections of theoretical models generated by molecular dynamics. The 4-mer (dCGCG) and 6-mer (dCGCGCG) duplexes were found to have globular conformations. Globular and helical structures were observed for the 8-mer (dCGCGCGCG) duplex, with the globular form being the more favored conformer. For the 10-mer (dCGCGCGCGCG), 14-mer (dCGCGCGCGCGCGCG), and 18-mer (dCGCGCGCGCGCGCGCGCG) duplexes, only helical structures were observed in the ion mobility measurements. Theory predicts that the helical structures are less stable than the globular forms in the gas phase and should collapse into the globular form given enough time. However, molecular dynamics simulations at 300 K indicate the helical structures are stable in aqueous solution and will retain their conformations for a limited time in the gas phase. The presence of helical structures in the ion mobility experiments indicates that the duplexes retain "solution structures" in the gas phase on the millisecond time scale.     

ortho-Phenylene oligomers with terminal push-pull substitution

ortho-Phenylenes are an emerging class of helical oligomers and polymers. We have synthesized a series of push-pull-substituted o-phenylene oligomers (dimethylamino/nitro) up to the octamer. Conformational analysis of the hexamer using a combination of low-temperature NMR spectroscopy and ab initio predictions of (1)H NMR chemical shifts indicates that, like other o-phenylenes, they exist as compact helices in solution. However, the substituents are found to have a significant effect on their conformational behavior: the nitro-functionalized terminus is 3-fold more likely to twist out of the helix. Protonation of the dimethylamino group favors the helical conformer. UV/vis spectroscopy indicates that the direct charge-transfer interaction between the push-pull substituents attenuates quickly compared to other conjugated systems, with no significant charge-transfer band for oligomers longer than the trimer. On protonation of the dimethylamino group, significant bathochromic shifts with increasing oligomer length are observed: the effective conjugation length is 9 repeat units, more than twice that of the parent oligomer. This behavior may be rationalized through examination of the frontier molecular orbitals of these compounds, which exhibit greater delocalization after protonation, as shown by DFT calculations.     

Solution Dependence of the Collisional Activation of Ubiquitin [M + 7H]7+ Ions

The solution dependence of gas-phase unfolding for ubiquitin [M + 7H]⁷⁺ ions has been studied by ion mobility spectrometry-mass spectrometry (IMS-MS). Different acidic water:methanol solutions are used to favor the native (N), more helical (A), or unfolded (U) solution states of ubiquitin. Unfolding of gas-phase ubiquitin ions is achieved by collisional heating and newly formed structures are examined by IMS. With an activation voltage of 100 V, a selected distribution of compact structures unfolds, forming three resolvable elongated states (E1-E3). The relative populations of these elongated structures depend strongly on the solution composition. Activation of compact ions from aqueous solutions known to favor N-state ubiquitin produces mostly the E1 type elongated state, whereas activation of compact ions from methanol containing solutions that populate A-state ubiquitin favors the E3 elongated state. Presumably, this difference arises because of differences in precursor ion structures emerging from solution. Thus, it appears that information about solution populations can be retained after ionization, selection, and activation to produce the elongated states. These data as well as others are discussed. Figure

?     

Systematic evaluation of the conformational properties of aliphatic omega-methoxy methyl esters

[structure: see text] A systematic conformational study of omega-methoxy methyl esters, CH(3)O-(CH2)n-COO-CH3 with n = 3 and 4, has been performed using quantum mechanical calculations at the MP2 level. Calculations have been carried out in both gas phase and chloroform solution, a polarizable continuum solvation model being used to represent the latter. Results have been compared with those recently obtained for the analogues omega-hydroxy acids, HO-(CH2)n-COOH with n = 3 and 4. The compounds with n = 3 clearly favor coiled conformations, the population expected for extended and semiextended conformations being very low. However, for compounds with n = 4 the minimum energy extended and semiextended structures become considerably more stable. The overall results indicate that the conformational preferences of the central aliphatic segment of omega-methoxy methyl esters and omega-hydroxy acids are not influenced by the formation of intramolecular O-H...O hydrogen bonds.     

Syntheses and self-assembly behaviors of the azobenzenyl modified beta-cyclodextrins isomers

A couple of analogues of azobenzenyl beta-cyclodextrins 1 and 2 with self-locked and self-unlocked conformations have been synthesized via the Huisgen cycloaddition from the same reactants, but in distinct reaction conditions (i.e., the hydrothermal synthesis and the "click" reaction, respectively), their conformations were sufficiently proved by X-ray crystal structural analysis, molecular modeling study, and 2D NMR spectroscopy, and their self-assembly behaviors in aqueous solution were also investigated by NMR spectroscopy. Interestingly, the self-locked conformer 1, which could be regarded as a new type of [1]rotaxane, self-assembled to a novel bimolecular capsule, where its azobenzene substituent was included in both its own cavity and the counterpart's cavity, in aqueous solution and in the solid state. In contrast, by adjusting the conformation of 1 to a self-unlocked one, the resulting conformer 2 was found to self-assemble to a linear supramolecule. These studies have shown the stronger impact of the reaction condition changing in cyclodextrin's modification products and will provide a new access to control the structure of supramolecular assemblies by tuning the conformation of building blocks.     

Importance of chirality and reduced flexibility of protein side chains: a study with square and tetrahedral lattice models

Side chains of amino acid residues are the determining factor that distinguishes proteins from other unstable chain polymers. In simple models they are often represented implicitly (e.g., by spin states) or simplified as one atom. Here we study side chain effects using two-dimensional square lattice and three-dimensional tetrahedral lattice models, with explicitly constructed side chains formed by two atoms of different chirality and flexibility. We distinguish effects due to chirality and effects due to side chain flexibilities, since residues in proteins are L residues, and their side chains adopt different rotameric states. For short chains, we enumerate exhaustively all possible conformations. For long chains, we sample effectively rare events such as compact conformations and obtain complete pictures of ensemble properties of conformations of these models at all compactness region. This is made possible by using sequential Monte Carlo techniques based on chain growth method. Our results show that both chirality and reduced side chain flexibility lower the folding entropy significantly for globally compact conformations, suggesting that they are important properties of residues to ensure fast folding and stable native structure. This corresponds well with our finding that natural amino acid residues have reduced effective flexibility, as evidenced by statistical analysis of rotamer libraries and side chain rotatable bonds. We further develop a method calculating the exact side chain entropy for a given backbone structure. We show that simple rotamer counting underestimates side chain entropy significantly for both extended and near maximally compact conformations. We find that side chain entropy does not always correlate well with main chain packing. With explicit side chains, extended backbones do not have the largest side chain entropy. Among compact backbones with maximum side chain entropy, helical structures emerge as the dominating configurations. Our results suggest that side chain entropy may be an important factor contributing to the formation of alpha helices for compact conformations.     

Conformational sampling of bioactive molecules: a comparative study

The necessity to generate conformations that sample the entire conformational space accessible to a given molecule is ubiquitous in the field of computer-aided drug design. Protein-ligand docking, 3D database searching, and 3D QSAR are three commonly used techniques that depend critically upon the quality and diversity of the generated conformers. Although there are a wide range of conformational search algorithms available, the extent to which they sample conformational space is often unclear. To address this question, we conducted a robust comparison of the search algorithms implemented in several widely used molecular modeling packages, including Catalyst, Macromodel, Omega, MOE, and Rubicon as well as our own method, stochastic proximity embedding (SPE). We found that SPE used in conjunction with conformational boosting, a heuristic for biasing conformational search toward more extended or compact geometries, along with Catalyst, are significantly more effective in sampling the full range of conformational space compared to the other methods, which show distinct preferences for either more extended or more compact geometries.     

Solid-state molecular folding and supramolecular structures of triptycene-derived secondary dicarboxamides

The synthesis and X-ray crystal structures of triptycene-derived secondary dicarboxamides 1 and 4-7 and reference compounds 2, 3, and 8 are reported. For comparison, molecular conformations of 1-8 in the gas phase and those of 1 and 3-6 in CD2Cl2 investigated by AM1 modeling and 1H NMR spectroscopy, respectively, are also included. The solid-state conformations of 1 and 5-8 are folded and compact, resulting from the cooperative effects of intramolecular amide-amide hydrogen bonding and edge-to-face arene-arene interactions between the triptycene and the N-acetylsulfanilyl groups. The sulfonyl ester groups are also essential in the folding of 1 and 5-8 and function as structural turn units. In contrast, the conformations of 2-4 are unfolded due to the lack of one of these three essentials. The extended triptycene ring systems in 6 and 7 provide an arene-arene contact mode that is different from that for 1 and 5. While AM1 calculations suggest that the two possible arene-arene contact modes in 6 and 7 have similar conformational energies, the one observed in the solid state is also favored in solutions. To achieve a more regular shape for compact crystal packing, the bulky triptycene groups tend to pack in pairs. As a result, the intermolecular amide-amide hydrogen bonding is perturbed and modified with the participation of either the sulfonyl groups or the methanol solvent molecules, leading to various hydrogen-bonding motifs for these triptycene diamides.     

Temperature-dependent H/D exchange of compact and elongated cytochrome c ions in the gas phase

Isotopic exchange reactions of compact and elongated conformations of gaseous cytochrome c ions (+5 and +9 states) with D2O have been measured as a function of temperature (from 300 to approximately 440 K) using ion mobility techniques. Rate constants for those sites that exchange at high temperatures (>400 K) are about an order of magnitude smaller than rate constants for sites that exchange at 300 K. Although the exchange rates decrease, the maximum exchange levels for rapidly exchanging sites increase with temperature. At 300 K, exchange levels of 53 +/- 3 and 63 +/- 3 are measured for the compact and elongated states, respectively. From 300 to 335 K, the exchange levels increase slightly to approximately 60 to 70 hydrogens. Above 335 K, the levels increase to a value of approximately 200 for the +5 state and approximately 190 for the +9 state, near the maximum possible levels, 200 and 204 for these respective charge states. Molecular dynamics simulations have been carried out on structures having calculated cross sections that are near the experimental values in order to explore the exchange process. Overall, it appears that charge site and exchange site proximities are important factors in the exchange profiles for the elongated +9 ion and the compact +5 ion.     

Transfer of structural elements from compact to extended states in unsolvated ubiquitin

Multidimensional ion mobility spectrometry techniques (IMS-IMS and IMS-IMS-IMS) combined with mass spectrometry are used to study structural transitions of ubiquitin ions in the gas phase. It is possible to select and activate narrow distributions of compact and partially folded conformation types and examine new distributions of structures that are formed. Different compact conformations unfold, producing a range of new partially folded states and three resolvable peaks associated with elongated conformers. Under gentle activation conditions, the final populations of the three elongated forms depend on the initial structures of the selected ions. This requires that some memory of the compact state (most likely secondary structure) is preserved along the unfolding pathway. Activation of selected, partially folded intermediates (formed from specific compact states) leads to elongated state populations that are consistent with the initial selected compact form-evidence that intermediates not only retain elements of initial structure but also are capable of transmitting structure to final states.     

Flaws in foldamers: conformational uniformity and signal decay in achiral helical peptide oligomers

Although foldamers, by definition, are extended molecular structures with a well-defined conformation, minor conformers must be populated at least to some extent in solution. We present a quantitative analysis of these minor conformers for a series of helical oligomers built from achiral but helicogenic α-amino acids. By measuring the chain length dependence or chain position dependence of NMR or CD quantities that measure screw-sense preference in a helical oligomer, we quantify values for the decay constant of a conformational signal as it passes through the molecular structure. This conformational signal is a perturbation of the racemic mixture of M and P helices that such oligomers typically adopt by the inclusion of an N or C terminal chiral inducer. We show that decay constants may be very low (<1% signal loss per residue) in non-polar solvents, and we evaluate the increase in decay constant that results in polar solvents, at higher temperatures, and with more conformationally flexible residues such as Gly. Decay constants are independent of whether the signal originates from the N or the C terminus. By interpreting the decay constant in terms of the probability with which conformations containing a screw-sense reversal are populated, we quantify the populations of these alternative minor conformers within the overall ensemble of secondary structures adopted by the foldamer. We deduce helical persistence lengths for Aib polymers that allow us to show that in a non-polar solvent a peptide helix, even in the absence of chiral residues, may continue with the same screw sense for approximately 200 residues.     

Role of entropy in increased rates of intramolecular reactions

The rates of intramolecular condensation of a series of monoesters of dicarboxylic acids have been shown to be highly dependent on the nature of the intervening groups. To understand the origin of this effect, we estimated DeltaS(NAC,S), the entropy difference between the ensemble of accessible ground state conformers and a single ground state conformer having transition-state-like geometry. DeltaS(NAC,S) differs from the activation entropy for the reaction by DeltaS(TS,NAC), the difference in vibrational entropy between the selected ground state conformer and the transition state. The estimated values of DeltaS(NAC,S) correlate well (R(2) = 0.96 and 0.73 using dielectric constant values of 80 and 1, respectively) with experimentally determined reaction rate constants. Normal-mode analysis performed on minimized ground state conformations of each molecule suggests that the change in vibrational entropy makes only a small contribution to the total activation entropy. These results indicate that the conformational entropy difference between the transition and the ground states contributes significantly to the free energy of activation.     

The energy landscape of unsolvated peptides: the role of context in the stability of alanine/glycine helices

Ion mobility measurements have been used to examine the conformations present for unsolvated Ac-(AG)(7)A+H(+) and (AG)(7)A+H(+) peptides (Ac = acetyl, A = alanine, and G = glycine) over a broad temperature range (100-410 K). The results are compared to those recently reported for Ac-A(4)G(7)A(4)+H(+) and A(4)G(7)A(4)+H(+), which have the same compositions but different sequences. Ac-(AG)(7)A+H(+) shows less conformational diversity than Ac-A(4)G(7)A(4)+H(+); it is much less helical than Ac-A(4)G(7)A(4)+H(+) at the upper end of the temperature range studied, and at low temperatures, one of the two Ac-A(4)G(7)A(4)+H(+) features assigned to helical conformations is missing for Ac-(AG)(7)A+H(+). Molecular dynamics simulations suggest that the different conformational preferences are not due to differences in the stabilities of the helical states, but differences in the nonhelical states: it appears that Ac-(AG)(7)A+H(+) is more flexible and able to adopt lower energy globular conformations (compact random looking three-dimensional structures) than Ac-A(4)G(7)A(4)+H(+). The helix to globule transition that occurs for Ac-(AG)(7)A+H(+) at around 250-350 K is not a direct (two-state) process, but a creeping transition that takes place through at least one and probably several intermediates.     

Conformational behavior of peptides containing residues of 3-azetidinesulfonic (3AzeS) and 4-piperidinemethanesulfonic (4PiMS) acids

The conformational behavior of four model peptides containing residues of 3-azetidinesulfonic (3AzeS) and 4-piperidinemethane sulfonic (4PiMS) acids was studied in both a crystalline state and in solution using X-ray, NMR, and IR experiments. It was found that in the crystalline state, both of the models di- and tripeptides studied adopted extended conformations and demonstrated considerable conformational flexibility. In solution, it is likely that a flexible ensemble of conformations is adopted, including extended structures and more compact ones without persistent hydrogen bonds. One of the most interesting features of the peptides was the axial chirality observed due to the slow rotation around the amide bond formed by the endocyclic nitrogen atoms of the non-chiral 3AzeS and 4PiMS residues. It was shown for one of the derivatives that the configuration of the chiral axis had an impact on the conformation of the neighboring amino acid residue.