Search NASASearch

SEARCH · Search NASA

Results for “Protein Structure, Secondary”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 19 records

Computer simulation of ion channel gating: the M(2) channel of influenza A virus in a lipid bilayer

The transmembrane fragment of the influenza virus M(2) protein forms a homotetrameric channel that transports protons. In this paper, we use molecular dynamics simulations to help elucidate the mechanism of channel gating by four histidines that occlude the channel lumen in the closed state. We test two competing hypotheses. In the "shuttle" mechanism, the delta nitrogen atom on the extracellular side of one histidine is protonated by the incoming proton, and, subsequently, the proton on the epsilon nitrogen atom is released on the opposite side. In the "water-wire" mechanism, the gate opens because of electrostatic repulsion between four simultaneously biprotonated histidines. This allows for proton transport along the water wire that penetrates the gate. For each system, composed of the channel embedded in a hydrated phospholipid bilayer, a 1.3-ns trajectory was obtained. It is found that the states involved in the shuttle mechanism, which contain either single-protonated histidines or a mixture of single-protonated histidines plus one biprotonated residue, are stable during the simulations. Furthermore, the orientations and dynamics of water molecules near the gate are conducive to proton transfer. In contrast, the fully biprotonated state is not stable. Additional simulations show that if only two histidines are biprotonated, the channel deforms but the gate remains closed. These results support the shuttle mechanism but not the gate-opening mechanism of proton gating in M(2).

NASA Center ARC

Early events in the folding of an amphipathic peptide: A multinanosecond molecular dynamics study

Folding of the capped LQQLLQQLLQL peptide is investigated at the water-hexane interface by molecular dynamics simulations for 161.5 ns. Initially placed in the aqueous phase as a beta-strand, the peptide rapidly adsorbs to the interface, where it adopts an amphipathic conformation. The marginal presence of nonamphipathic structures throughout the complete trajectory indicates that the corresponding conformations are strongly disfavored at the interface. It is further suggestive that folding in an interfacial environment proceeds through a pathway of successive amphipathic intermediates. The energetic and entropic penalties involved in the conformational changes along this pathway markedly increase the folding time scales of LQQLLQQLLQL, explaining why the alpha-helix, the hypothesized lowest free energy structure for a sequence with a hydrophobic periodicity of 3.6, has not been reached yet. The formation of a type I beta-turn at the end of the simulation confirms the importance of such motifs as initiation sites allowing the peptide to coalesce towards a secondary structure. Proteins 1999;36:383-399. Copyright 1999 Wiley-Liss, Inc.

NASA Center ARC

Helicity, membrane incorporation, orientation and thermal stability of the large conductance mechanosensitive ion channel from E. coli

In this report, we present structural studies on the large conductance mechanosensitive ion channel (MscL) from E. coli in detergent micelles and lipid vesicles. Both transmission Fourier transform infrared spectroscopy and circular dichroism (CD) spectra indicate that the protein is highly helical in detergents as well as liposomes. The secondary structure of the proteins was shown to be highly resistant towards denaturation (25-95 degrees C) based on an ellipticity thermal profile. Amide H+/D+ exchange was shown to be extensive (ca. 66%), implying that two thirds of the protein are water accessible. MscL, reconstituted in oriented lipid bilayers, was shown to possess a net bilayer orientation using dichroic ratios measured by attenuated total-reflection Fourier transform infrared spectroscopy. Here, we present and discuss this initial set of structural data on this new family of ion-channel proteins.

NASA Discipline Cell Biology

Multistep modeling of protein structure: application to bungarotoxin

Modelling of bungarotoxin in atomic details is presented in this article. The model-building procedure utilizes the low-resolution crystal coordinates of the c-alpha atoms of bungarotoxin, sequence homology within the neurotoxin family, as well as high-resolution x-ray diffraction data of cobratoxin and erabutoxin. Our model-building procedure involves: (a) principles of comparative modelling, (b) embedding procedures of distance geometry, and (c) use of molecular mechanics for optimizing packing. The model is not only consistent with the c-alpha coordinates of crystal structure, but also agrees with solution conformational features of the triple-stranded beta sheet as observed by NOE measurements.

NASA Discipline Exobiology

Proteolytic dissection of Zab, the Z-DNA-binding domain of human ADAR1

Zalpha is a peptide motif that binds to Z-DNA with high affinity. This motif binds to alternating dC-dG sequences stabilized in the Z-conformation by means of bromination or supercoiling, but not to B-DNA. Zalpha is part of the N-terminal region of double-stranded RNA adenosine deaminase (ADAR1), a candidate enzyme for nuclear pre-mRNA editing in mammals. Zalpha is conserved in ADAR1 from many species; in each case, there is a second similar motif, Zbeta, separated from Zalpha by a more divergent linker. To investigate the structure-function relationship of Zalpha, its domain structure was studied by limited proteolysis. Proteolytic profiles indicated that Zalpha is part of a domain, Zab, of 229 amino acids (residues 133-361 in human ADAR1). This domain contains both Zalpha and Zbeta as well as a tandem repeat of a 49-amino acid linker module. Prolonged proteolysis revealed a minimal core domain of 77 amino acids (positions 133-209), containing only Zalpha, which is sufficient to bind left-handed Z-DNA; however, the substrate binding is strikingly different from that of Zab. The second motif, Zbeta, retains its structural integrity only in the context of Zab and does not bind Z-DNA as a separate entity. These results suggest that Zalpha and Zbeta act as a single bipartite domain. In the presence of substrate DNA, Zab becomes more resistant to proteases, suggesting that it adopts a more rigid structure when bound to its substrate, possibly with conformational changes in parts of the protein.

Non-NASA Center

Calcium-stimulated autophosphorylation site of plant chimeric calcium/calmodulin-dependent protein kinase

The existence of two molecular switches regulating plant chimeric Ca(2+)/calmodulin-dependent protein kinase (CCaMK), namely the C-terminal visinin-like domain acting as Ca(2+)-sensitive molecular switch and calmodulin binding domain acting as Ca(2+)-stimulated autophosphorylation-sensitive molecular switch, has been described (Sathyanarayanan, P. V., Cremo, C. R., and Poovaiah, B. W. (2000) J. Biol. Chem. 275, 30417-30422). Here we report the identification of Ca(2+)-stimulated autophosphorylation site of CCaMK by matrix-assisted laser desorption ionization time of flight-mass spectrometry. Thr(267) was confirmed as the Ca(2+)-stimulated autophosphorylation site by post-source decay experiments and by site-directed mutagenesis. The purified T267A mutant form of CCaMK did not show Ca(2+)-stimulated autophosphorylation, autophosphorylation-dependent variable calmodulin affinity, or Ca(2+)/calmodulin stimulation of kinase activity. Sequence comparison of CCaMK from monocotyledonous plant (lily) and dicotyledonous plant (tobacco) suggests that the autophosphorylation site is conserved. This is the first identification of a phosphorylation site specifically responding to activation by second messenger system (Ca(2+) messenger system) in plants. Homology modeling of the kinase and calmodulin binding domain of CCaMK with the crystal structure of calcium/calmodulin-dependent protein kinase 1 suggests that the Ca(2+)-stimulated autophosphorylation site is located on the surface of the kinase and far from the catalytic site. Analysis of Ca(2+)-stimulated autophosphorylation with increasing concentration of CCaMK indicates the possibility that the Ca(2+)-stimulated phosphorylation occurs by an intermolecular mechanism.

Non-NASA Center

Crystal structures of pinoresinol-lariciresinol and phenylcoumaran benzylic ether reductases and their relationship to isoflavone reductases

Despite the importance of plant lignans and isoflavonoids in human health protection (e.g. for both treatment and prevention of onset of various cancers) as well as in plant biology (e.g. in defense functions and in heartwood development), systematic studies on the enzymes involved in their biosynthesis have only recently begun. In this investigation, three NADPH-dependent aromatic alcohol reductases were comprehensively studied, namely pinoresinol-lariciresinol reductase (PLR), phenylcoumaran benzylic ether reductase (PCBER), and isoflavone reductase (IFR), which are involved in central steps to the various important bioactive lignans and isoflavonoids. Of particular interest was in determining how differing regio- and enantiospecificities are achieved with the different enzymes, despite each apparently going through similar enone intermediates. Initially, the three-dimensional x-ray crystal structures of both PLR_Tp1 and PCBER_Pt1 were solved and refined to 2.5 and 2.2 A resolutions, respectively. Not only do they share high gene sequence similarity, but their structures are similar, having a continuous alpha/beta NADPH-binding domain and a smaller substrate-binding domain. IFR (whose crystal structure is not yet obtained) was also compared (modeled) with PLR and PCBER and was deduced to have the same overall basic structure. The basis for the distinct enantio-specific and regio-specific reactions of PCBER, PLR, and IFR, as well as the reaction mechanism and participating residues involved (as identified by site-directed mutagenesis), are discussed.

NASA Discipline Plant Biology

IR spectroscopic characteristics of cell cycle and cell death probed by synchrotron radiation based Fourier transform IR spectromicroscopy

Synchrotron radiation based Fourier transform IR (SR-FTIR) spectromicroscopy allows the study of individual living cells with a high signal to noise ratio. Here we report the use of the SR-FTIR technique to investigate changes in IR spectral features from individual human lung fibroblast (IMR-90) cells in vitro at different points in their cell cycle. Clear changes are observed in the spectral regions corresponding to proteins, DNA, and RNA as a cell changes from the G(1)-phase to the S-phase and finally into mitosis. These spectral changes include markers for the changing secondary structure of proteins in the cell, as well as variations in DNA/RNA content and packing as the cell cycle progresses. We also observe spectral features that indicate that occasional cells are undergoing various steps in the process of cell death. The dying or dead cell has a shift in the protein amide I and II bands corresponding to changing protein morphologies, and a significant increase in the intensity of an ester carbonyl C===O peak at 1743 cm(-1) is observed. Copyright John Wiley & Sons, Inc. Biopolymers (Biospectroscopy) 57: 329-335, 2000.

Non-NASA Center

Secondary structure in solution of two anti-HIV-1 hammerhead ribozymes as investigated by two-dimensional 1H 500 MHz NMR spectroscopy in water

Two hammerhead chimeric RNA/DNA ribozymes (HRz) were synthesized in pure form. Both were 30 nucleotides long, and the sequences were such that they could be targeted to cleave the HIV-1 gag RNA. Named HRz-W and HRz-M, the former had its invariable core region conserved, the latter had a uridine in the invariable region replaced by a guanine. Their secodary structures were determined by 2D NOESY 1H 500 MHz NMR spectroscopy in 90% water and 10% D2(0), following the imino protons. The data show that both HRz-M and HRz-W form identical secondary structures with stem regions consisting of continuous stacks of AT and GT pairs. An energy minimized computer model of this stem region is provided. The results suggest that the loss of catalytic activity that is known to result when an invariant core residue is replaced is not related to the secondary structure of the ribozymes in the absence of substrate.

NASA Discipline Exobiology

Structural models of the MscL gating mechanism

Three-dimensional structural models of the mechanosensitive channel of large conductance, MscL, from the bacteria Mycobacterium tuberculosis and Escherichia coli were developed for closed, intermediate, and open conformations. The modeling began with the crystal structure of M. tuberculosis MscL, a homopentamer with two transmembrane alpha-helices, M1 and M2, per subunit. The first 12 N-terminal residues, not resolved in the crystal structure, were modeled as an amphipathic alpha-helix, called S1. A bundle of five parallel S1 helices are postulated to form a cytoplasmic gate. As membrane tension induces expansion, the tilts of M1 and M2 are postulated to increase as they move away from the axis of the pore. Substantial expansion is postulated to occur before the increased stress in the S1 to M1 linkers pulls the S1 bundle apart. During the opening transition, the S1 helices and C-terminus amphipathic alpha-helices, S3, are postulated to dock parallel to the membrane surface on the perimeter of the complex. The proposed gating mechanism reveals critical spatial relationships between the expandable transmembrane barrel formed by M1 and M2, the gate formed by S1 helices, and "strings" that link S1s to M1s. These models are consistent with numerous experimental results and modeling criteria.

NASA Discipline Cell Biology

Constructing high complexity synthetic libraries of long ORFs using in vitro selection

We present a method that can significantly increase the complexity of protein libraries used for in vitro or in vivo protein selection experiments. Protein libraries are often encoded by chemically synthesized DNA, in which part of the open reading frame is randomized. There are, however, major obstacles associated with the chemical synthesis of long open reading frames, especially those containing random segments. Insertions and deletions that occur during chemical synthesis cause frameshifts, and stop codons in the random region will cause premature termination. These problems can together greatly reduce the number of full-length synthetic genes in the library. We describe a strategy in which smaller segments of the synthetic open reading frame are selected in vitro using mRNA display for the absence of frameshifts and stop codons. These smaller segments are then ligated together to form combinatorial libraries of long uninterrupted open reading frames. This process can increase the number of full-length open reading frames in libraries by up to two orders of magnitude, resulting in protein libraries with complexities of greater than 10(13). We have used this methodology to generate three types of displayed protein library: a completely random sequence library, a library of concatemerized oligopeptide cassettes with a propensity for forming amphipathic alpha-helical or beta-strand structures, and a library based on one of the most common enzymatic scaffolds, the alpha/beta (TIM) barrel. Copyright 2000 Academic Press.

Non-NASA Center

Classification and evolution of EF-hand proteins

Forty-five distinct subfamilies of EF-hand proteins have been identified. They contain from two to eight EF-hands that are recognizable by amino acid sequence as being statistically similar to other EF-hand domains. All proteins within one subfamily are congruent to one another, i.e. the dendrogram computed from one of the EF-hand domains is similar, within statistical error, to the dendrogram computed from another(s) domain. Thirteen subfamilies--including Calmodulin, Troponin C, Essential light chain, Regulatory light chain--referred to collectively as CTER, are congruent with one another. They appear to have evolved from a single ur-domain by two cycles of gene duplication and fusion. The subfamilies of CTER subsequently evolved by gene duplications and speciations. The remaining 32 subfamilies do not show such general patterns of congruence; however, some--such as S100, intestinal calcium binding protein (calbindin 9 kd), and trichohylin--do not form congruent clusters of subfamilies. Nearly all of the domains 1, 3, 5, and 7 are most similar to other ODD domains. Correspondingly the EVEN numbered domains of all 45 subfamilies most closely resemble EVEN domains of other subfamilies. Many sequence and chemical characteristics do not show systemic trends by subfamily or species of host organisms; such homoplasy is widespread. Eighteen of the subfamilies are heterochimeric; in addition to multiple EF-hands they contain domains of other evolutionary origins.

Non-NASA Center

Single residue substitutions that change the gating properties of a mechanosensitive channel in Escherichia coli

MscL is a channel that opens a large pore in the Escherichia coli cytoplasmic membrane in response to mechanical stress. Previously, we highly enriched the MscL protein by using patch clamp as a functional assay and cloned the corresponding gene. The predicted protein contains a largely hydrophobic core spanning two-thirds of the molecule and a more hydrophilic carboxyl terminal tail. Because MscL had no homology to characterized proteins, it was impossible to predict functional regions of the protein by simple inspection. Here, by mutagenesis, we have searched for functionally important regions of this molecule. We show that a short deletion from the amino terminus (3 amino acids), and a larger deletion of 27 amino acids from the carboxyl terminus of this protein, had little if any effect in channel properties. We have thus narrowed the search of the core mechanosensitive mechanism to 106 residues of this 136-amino acid protein. In contrast, single residue substitutions of a lysine in the putative first transmembrane domain or a glutamine in the periplasmic loop caused pronounced shifts in the mechano-sensitivity curves and/or large changes in the kinetics of channel gating, suggesting that the conformational structure in these regions is critical for normal mechanosensitive channel gating.

Non-NASA Center

Positive selection moments identify potential functional residues in human olfactory receptors

Correlated mutation analysis and molecular models of olfactory receptors have provided evidence that residues in the transmembrane domains form a binding pocket for odor ligands. As an independent test of these results, we have calculated positive selection moments for the alpha-helical sixth transmembrane domain (TM6) of human olfactory receptors. The moments can be used to identify residues that have been preferentially affected by positive selection and are thus likely to interact with odor ligands. The results suggest that residue 622, which is commonly a serine or threonine, could form critical H-bonds. In some receptors a dual-serine subsite, formed by residues 622 and 625, could bind hydroxyl determinants on odor ligands. The potential importance of these residues is further supported by site-directed mutagenesis in the beta-adrenergic receptor. The findings should be of practical value for future physiological studies, binding assays, and site-directed mutagenesis.

Non-NASA Center

Approximate scaling properties of RNA free energy landscapes

RNA free energy landscapes are analysed by means of "time-series" that are obtained from random walks restricted to excursion sets. The power spectra, the scaling of the jump size distribution, and the scaling of the curve length measured with different yard stick lengths are used to describe the structure of these "time series". Although they are stationary by construction, we find that their local behavior is consistent with both AR(1) and self-affine processes. Random walks confined to excursion sets (i.e., with the restriction that the fitness value exceeds a certain threshold at each step) exhibit essentially the same statistics as free random walks. We find that an AR(1) time series is in general approximately self-affine on timescales up to approximately the correlation length. We present an empirical relation between the correlation parameter rho of the AR(1) model and the exponents characterizing self-affinity.

Signal Processing, Computer-Assisted

Finding the global minimum: a fuzzy end elimination implementation

The 'fuzzy end elimination theorem' (FEE) is a mathematically proven theorem that identifies rotameric states in proteins which are incompatible with the global minimum energy conformation. While implementing the FEE we noticed two different aspects that directly affected the final results at convergence. First, the identification of a single dead-ending rotameric state can trigger a 'domino effect' that initiates the identification of additional rotameric states which become dead-ending. A recursive check for dead-ending rotameric states is therefore necessary every time a dead-ending rotameric state is identified. It is shown that, if the recursive check is omitted, it is possible to miss the identification of some dead-ending rotameric states causing a premature termination of the elimination process. Second, we examined the effects of removing dead-ending rotameric states from further considerations at different moments of time. Two different methods of rotameric state removal were examined for an order dependence. In one case, each rotamer found to be incompatible with the global minimum energy conformation was removed immediately following its identification. In the other, dead-ending rotamers were marked for deletion but retained during the search, so that they influenced the evaluation of other rotameric states. When the search was completed, all marked rotamers were removed simultaneously. In addition, to expand further the usefulness of the FEE, a novel method is presented that allows for further reduction in the remaining set of conformations at the FEE convergence. In this method, called a tree-based search, each dead-ending pair of rotamers which does not lead to the direct removal of either rotameric state is used to reduce significantly the number of remaining conformations. In the future this method can also be expanded to triplet and quadruplet sets of rotameric states. We tested our implementation of the FEE by exhaustively searching ten protein segments and found that the FEE identified the global minimum every time. For each segment, the global minimum was exhaustively searched in two different environments: (i) the segments were extracted from the protein and exhaustively searched in the absence of the surrounding residues; (ii) the segments were exhaustively searched in the presence of the remaining residues fixed at crystal structure conformations. We also evaluated the performance of the method for accurately predicting side chain conformations. We examined the influence of factors such as type and accuracy of backbone template used, and the restrictions imposed by the choice of potential function, parameterization and rotamer database. Conclusions are drawn on these results and future prospects are given.

NASA Program Exobiology

Alpha-helical hydrophobic polypeptides form proton-selective channels in lipid bilayers

Proton translocation is important in membrane-mediated processes such as ATP-dependent proton pumps, ATP synthesis, bacteriorhodopsin, and cytochrome oxidase function. The fundamental mechanism, however, is poorly understood. To test the theoretical possibility that bundles of hydrophobic alpha-helices could provide a low energy pathway for ion translocation through the lipid bilayer, polyamino acids were incorporated into extruded liposomes and planar lipid membranes, and proton translocation was measured. Liposomes with incorporated long-chain poly-L-alanine or poly-L-leucine were found to have proton permeability coefficients 5 to 7 times greater than control liposomes, whereas short-chain polyamino acids had relatively little effect. Potassium permeability was not increased markedly by any of the polyamino acids tested. Analytical thin layer chromatography measurements of lipid content and a fluorescamine assay for amino acids showed that there were approximately 135 polyleucine or 65 polyalanine molecules associated with each liposome. Fourier transform infrared spectroscopy indicated that a major fraction of the long-chain hydrophobic peptides existed in an alpha-helical conformation. Single-channel recording in both 0.1 N HCl and 0.1 M KCl was also used to determine whether proton-conducting channels formed in planar lipid membranes (phosphatidylcholine/phosphatidylethanolamine, 1:1). Poly-L-leucine and poly-L-alanine in HCl caused a 10- to 30-fold increase in frequency of conductive events compared to that seen in KCl or by the other polyamino acids in either solution. This finding correlates well with the liposome observations in which these two polyamino acids caused the largest increase in membrane proton permeability but had little effect on potassium permeability. Poly-L-leucine was considerably more conductive than poly-L-alanine due primarily to larger event amplitudes and, to a lesser extent, a higher event frequency. Poly-L-leucine caused two populations of conductive events, one in the 0.1-0.5 pA range, and one in the 1.0-5.0 pA range, whereas nearly all events caused by poly-L-alanine were in the 0.1-0.5 pA range at an applied voltage of +60 mV. The channel-like activity appeared to switch between conductive and nonconductive states, with most open-times in the range of 50-200 ms. We conclude that hydrophobic polyamino acids produce proton-conducting defects in lipid bilayers that may be used to model functional proton channels in biological membranes.

NASA Discipline Exobiology