Translational dynamics of antifreeze glycoprotein in supercooled water

Structure and dynamics of biomolecules in supercooled water assume a particular and distinct importance in the case of antifreeze glycoproteins (AFGPs), which function at sub-zero temperatures. To investigate whether any large-scale structural digressions in the supercooled state are correlated to the function of AFGPs, self-diffusion behavior of the AFGP8, the smallest AFGP is monitored as a function of temperature from 243 to 303 K using nuclear magnetic resonance (NMR) spectroscopy. The experimental results are compared with the hydrodynamic calculations using the viscosity of water at the same temperature range. In order to evaluate results on AFGP8, the smallest AFGP, constituting approximately two-thirds of the total AFGP fraction in fish blood serum, similar experimental and computational calculations were also performed on a set of globular proteins. These results show that even though the general trend of translational dynamics of AFGP is similar to that of the other globular proteins, AFGP8 appears to be more hydrated (approximately 30% increase in the bead radius) than the others over the temperature range studied. These results also suggest that local conformational changes such as segmental librations or hydrogen bond dynamics that are closer to the protein surface are more likely the determining dynamic factors for the function of AFGPs rather than any large-scale structural rearrangements.     

Blood glycoprotein from antarctic fish possible conformational origin of antifreeze activity

High resolution ¹H NMR and circular dichroism (CD) measurements have been performed on aqueous solutions of antarctic fish antifreeze glycoprotein. The carbohydrate contribution to the observed CD spectrum has been estimated from closely analogous model compounds. The residual peptide contribution cannot be interpreted of the known spectral behaviour of α-helix, β-sheet and random coil. Instead it resembles the CD spectrum of β-structure in position, magnitude and spectral form, but is of opposite sign, indicating a specific but unusual peptide conformation, which we suggest may be stabilized by non-bonded interactions between the peptide backbone and the carbohydrate sidechains. Previous evidence which supports this interpretation is reviewed. NMR and CD measurements between ?2 and +30°C are consistent with conformational stability throughout the biologically relevant temperature range. The mechanism of the antifreeze activity is discussed in terms of the spatial and orientational correlations of sugar hydroxy groups and water in the liquid and solid states. The implication of an ordered peptide structure is explained by the comparison of the antifreeze glycoprotein with synthetic water-soluble polymers which also exhibit limited antifreeze properties.     

Blood glycoprotein from antarctic fish. Possible conformational origin of antifreeze activity

Two purine nucleoside phosphorylases (purine-nucleoside:orthophosphate ribosyltransferase, EC 2.4.2.1) were purified from vegetative Bacillus subtilis cells. One enzyme, inosine-guanosine phosphorylase, showed great similarity to the homologous enzyme of Bacillus cereus. It appeared to be a tetramer of molecular weight 95 000. The other enzyme, adenosine phosphorylase, was specific for adenosine and deoxyadenosine. The molecular weight of the native enzyme was 153 000 +/- 10% and the molecular weight of the subunits was 25 500 +/- 5%. This indicates a hexameric structure. The adenosine phosphorylase was inactivated by 10(-3) M p-chloromercuribenzoate and protected against this inactivation by phosphate, adenosine and ribose 1-phosphate.     

Synthesis and conformational studies of proline-containing tripeptides

Using tripeptides of the type Boc-Pro-X-Gly-OEt and Boc-X-Pro-Gly-OEt where X = Val, Leu, Ile, Nle we have studied the influence of the X residue on the stability of folded conformations, most probably the β-turn, in these peptides. In addition, the substitution Gly→β-Ala was also investigated. Our c.d. and i.r. studies show significant changes in β-turn stability depending on the nature and the position of X and on specific solute-solvent interactions.     

Allostery and conformational free energy changes in human tryptophanyl‐tRNA synthetase from essential dynamics and structure networks

The interdependence of the concept of allostery and enzymatic catalysis, and they being guided by conformational mobility is gaining increased prominence. However, to gain a molecular level understanding of allostery and hence of enzymatic catalysis, it is of utter importance that the networks of amino acids participating in allostery be deciphered. Our lab has been exploring the methods of network analysis combined with molecular dynamics simulations to understand allostery at molecular level. Earlier we had outlined methods to obtain communication paths and then to map the rigid/flexible regions of proteins through network parameters like the shortest correlated paths, cliques, and communities. In this article, we advance the methodology to estimate the conformational populations in terms of cliques/communities formed by interactions including the side‐chains and then to compute the ligand‐induced population shift. Finally, we obtain the free‐energy landscape of the protein in equilibrium, characterizing the free‐energy minima accessed by the protein complexes. We have chosen human tryptophanyl‐tRNA synthetase (hTrpRS), a protein responsible for charging tryptophan to its cognate tRNA during protein biosynthesis for this investigation. This is a multidomain protein exhibiting excellent allosteric communication. Our approach has provided valuable structural as well as functional insights into the protein. The methodology adopted here is highly generalized to illuminate the linkage between protein structure networks and conformational mobility involved in the allosteric mechanism in any protein with known structure. Proteins 2010. © 2009 Wiley‐Liss, Inc.     

Aggregation of antifreeze glycoprotein fraction 8 and its effect on antifreeze activity

Antifreeze glycoproteins (AFGPs) have many potential applications ranging from the cryopreservation and hypothermic storage of tissues and organs to the preservation of various frozen food products. Since supplying native AFGP for these applications is a labor-intensive and costly process, the rational design and synthesis of functional AFGP analogues is a very attractive alternative. While structure-function studies have implicated specific structural motifs as essential for antifreeze activity in AFGP, the relationship between solution conformation and antifreeze activity is poorly understood. Toward this end, we have analyzed AFGP8 in aqueous solutions using dynamic light scattering (DLS) and circular dichroism (CD). Our results indicate that AFGP8 forms discrete aggregates in solution. These aggregates are predominantly composed of dimers that form at solution concentrations greater than 20 mM. CD spectroscopy indicates that the preferred solution conformation of AFGP8 is consistent with that of random coil. However, significant beta-sheet and alpha-helix character is observed in more concentrated solutions, indicating that these glycopeptides are highly flexible in solution. Aggregation appears to have a minimal effect on the overall solution conformation. Thermal hysteresis (TH) activity of the aggregated solutions is much higher than that of less concentrated solutions that do not form aggregates. While cooperative functioning between lower and higher molecular weight AFGPs has been reported, this is the first instance where cooperative functioning in lower molecular weight AFGPs has been observed.     

Prion glycoprotein: structure, dynamics, and roles for the sugars

The prion protein contains two N-linked glycosylation sites and a glycosylphosphatidylinositol (GPI) anchor. The large size of the N-linked sugars, together with their dynamic properties, enables them to shield two orthogonal faces of the protein almost completely. Thus, the sugars can protect large regions of the protein surface from proteases and from nonspecific protein-protein interactions. Immunoprecipitation of prion protein with calnexin suggests that in the ER the oligosaccharides may provide a route for protein folding via the calnexin pathway. Major questions relate to the relevance of the glycoform distribution (as defined by glycan site occupancy) to strain type and disease transmission. Glycan analysis has shown that prion protein contains at least 52 different sugars, that these consist of a subset of brain sugars, and that there is site specific glycan processing. PrP(Sc) from the brains of Syrian hamsters contains the same set of glycans as PrP(C), but a higher proportion of tri- and tetra-antennary sugars. This may be attributed to a decrease in the activity of GnTIII. The GPI anchor, which is modified with sialic acid, may allow the prion protein to be mobile in the lipid bilayer. Potentially, this provides a possible means for translocating the prions from one cell to another.     

Structure and thermodynamics of RNA-protein binding: using molecular dynamics and free energy analyses to calculate the free energies of binding and conformational change

An adaptive binding mechanism, requiring large conformational rearrangements, occurs commonly with many RNA-protein associations. To explore this process of reorganization, we have investigated the conformational change upon spliceosomal U1A-RNA binding with molecular dynamics (MD) simulations and free energy analyses. We computed the energetic cost of conformational change in U1A-hairpin and U1A-internal loop binding using a hybrid of molecular mechanics and continuum solvent methods. Encouragingly, in all four free energy comparisons (two slightly different proteins, two different RNAs), the free macromolecule was more stable than the bound form by the physically reasonable value of approximately 10 kcal/mol. We calculated the absolute binding free energies for both complexes to be in the same range as that found experimentally.     

Crystal structure,conformational analysis,and molecular dynamics of tetra-O-methyl-(+) -catechin

The structure of tetra-O-methyl- (+) -catechin has been determined in the crystalline state. Two independent molecules, denoted structure A and structure B, exist in the unit cell. Crystals are triclinic, space group P1, a = 4.8125(2) Å, b = 12.9148(8) Å, c = 13.8862(11) Å, α = 86.962(6) °, β = 89.120(5)°, γ = 88.044(5)°, Z = 2, D_c = 1.336 g cm^?3, R = 0.033 for 6830 observations. The heterocyclic rings of the crystal structures are compared to previous results for 8-bromotetra-O-methyl-(+)-catechin, penta-O-acetyl-(+)-catechin, and (?) -epicatechin. One of the two molecules has a heterocyclic ring conformation similar to that observed previously for (?)-epicatechin, and the other has a heterocyclic ring conformation similar to one predicted earlier in a theoretical analysis of dimers of (+)-catechin and (?) -epicatechin. Both structure A and structure B in the crystal have heterocyclic ring conformations that place the dimethoxyphenyl substituent at C(2) in the equatorial position. However, this heterocyclic ring conformation does not explain the proton nmr coupling constant measured in solution. Molecular dynamics simulations show an equatorial ? axial interconversion of the heterocyclic ring, which can explain the nmr results. © 1993 John Wiley & Sons, Inc.     

Minimum energy conformations of proline-containing helices.

Proline occurs frequently in transmembrane alpha-helices of transport and receptor proteins even though statistical surveys demonstrate the overwhelming preference of this residue for a non-alpha-helical, hydrophilic environment. As a result, membrane-buried proline has been proposed to be functionally important, with function arising from structural discontinuity or destabilization of the helix. Destabilization may occur by Pro-mediated conformational transitions between discrete states, and may be manifested in membrane protein systems through reversible processes such as channel opening and closing or signal transduction. In this study, computer modeling of a model transmembrane alpha-helix, (Ala)8-Leu-Pro-Phe-(Ala)8, in a medium of low polarity (dielectric = 2), is used to examine the occurrence and energetic accessibility of Pro-mediated conformational interconversions. Leu psi and chi 1, Pro psi, and Phe phi and chi 1 torsion angles were assigned random values so that a data base of 200 conformations for each of the cis and trans states was generated. The conformations were minimized and low-energy structures organized into families. This analysis demonstrated that the most populated lowest energy family is the Trans-I conformation, corresponding to proline in a kinked alpha-helix. Two additional trans structures, Trans-II and Trans-III, as well as a cis conformation, Cis-I, are also energetically competitive. Interconversions between the trans states could thus be mediated by changes at a single torsion angle, accompanied by minor local hydrogen-bonding rearrangements. This work substantiates that membrane-buried proline can provide the basis for conformational transitions between discrete alpha-helix-based structures in a nonpolar environment.     

Murine interleukin-3: structure, dynamics, and conformational heterogeneity in solution

Interleukin-3 (IL-3), a cytokine produced primarily by activated T-cells during immune responses, is a crucial regulator of allergic inflammation. The three-dimensional structure of murine IL-3 (mIL-3) has remained elusive owing to its poor solubility and strong tendency toward aggregation under solution conditions typically used for structural studies. Here we describe the solution properties and structure of mIL-3 determined by NMR using an engineered construct of mIL-3 (mIL-3(33-156)). mIL-3 adopts a four-helical bundle fold, typical of proteins belonging to the short-chain cytokine family, and features a core of highly conserved hydrophobic residues. While significant line broadening and peak disappearance were observed in NMR spectra at higher temperatures, there was no evidence for temperature-dependent changes of the oligomeric state of mIL-3(33-156). Further analysis of the temperature dependence of amide (1)H chemical shifts and backbone (15)N relaxation parameters, including (15)N relaxation dispersion, revealed the presence of significant conformational exchange and local conformational heterogeneity. Residues recently shown by mutagenesis to play key roles in β(IL-3) receptor recognition and activation, which are located within the α(A) and α(C) helices and aligned on one face of the mIL-3(33-156) structure, are relatively rigid. In contrast, pronounced conformational heterogeneity was observed for a cluster of residues located on the opposite side of mIL-3, which corresponds spatially to sites in the related cytokines human IL-3, IL-5, and GM-CSF that are known to mediate interactions with their respective α-receptor subunits. Such conformational heterogeneity may facilitate the interaction of mIL-3 with each of two naturally occurring mIL-3Rα isoforms, leading to structurally distinct high-affinity complexes.     

A serendipitous discovery of antifreeze protein-specific activity in C-linked antifreeze glycoprotein analogs

Structurally diverse carbon-linked (C-linked) analogs of antifreeze glycoprotein (AFGP) have been prepared via linear or convergent solid phase synthesis. These analogs range in molecular weight from approx 1.5–4.1 KDa and do not possess the β-d-galactose-1,3-α-d-N-acetylgalactosamine carbohydrate moiety or the l-threonine-l-alanine-l-alanine polypeptide backbone native to the AFGP wild-type. Despite these dramatic structural modifications, the 2.7-KDa and 4.1-KDa analogs possess antifreeze protein-specific activity as determined by recrystallization-inhibition (RI) and thermal hysteresis (TH) assays. These analogs are weaker than the wild-type in their activity, but nanoliter osmometry indicates that these compounds are binding to ice and affecting a localized freezing point depression. This is the first example of a C-linked AFGP analog that possesses TH and RI activity and suggests that the rational design and synthesis of chemically and biologically stable AFGP analogs is a feasible and worthwhile endeavor. Given the low degree of TH activity, these compounds may prove useful for the protection of cells during freezing and thawing cycles.     

A kinetic description of antifreeze glycoprotein activity

The antifreeze glycoproteins (AFGP) of polar fish have the ability to depress the freezing temperature of water approximately 500 times the amount expected based on the number of AFGP molecules in solution; yet AFGP solutions have a purely colligative melting point depression. The difference of solution melting and freezing temperatures is the antifreeze activity of AFGP. One characteristic of AFGP activity that requires further examination is the effect of concentration on antifreeze activity, especially whether the activity saturates at high concentrations or the measured activity increases ad infinitum. This study first surveys the activity of the various antifreeze components from both Pagothenia borchgrevinki and the Arg-containing antifreeze glycoprotein from Eleginus gracilis (EgAF). It was found that all AFGP components examined have a plateau in activity at high concentration, but the actual value of the plateau activity differs between the different length AFGP components and between AFGP and EgAF. While the low molecular weight components of both AFGP and EgAF lose activity at deep supercooling, at high concentration activity is restored. The activity data is then shown to fit a reversible kinetic model of AFGP activity, and the coefficients obtained are used to compare the activity differences between AFGP components and between AFGP and EgAF. The model is also shown to describe the activity of the antifreeze protein of the fish Pseudopleuronectes americanus and the thermal hysteresis protein of the insect, Tenebrio molitor.     

The free energy landscapes governing conformational changes in a glutamate receptor ligand-binding domain

Ionotropic glutamate receptors are ligand-gated transmembrane ion channels activated by the binding of glutamate. The free energy landscapes governing the opening/closing of the GluR2 S1S2 ligand-binding domain in the apo, DNQX-, and glutamate-bound forms are computed by using all-atom molecular dynamics simulations with explicit solvent, in conjunction with an umbrella sampling strategy. The apo S1S2 easily accesses low-energy conformations that are more open than observed in X-ray crystal structures. A free energy of 9-12 kcal/mol becomes available upon glutamate binding for driving conformational changes in S1S2 associated with receptor activation. Small-angle X-ray scattering profiles calculated from computed ensemble averages agree better with experimental results than profiles calculated from static X-ray crystal structures. Water molecules in the cleft may contribute to stabilizing the apo S1S2 in open conformations. Free energy landscapes were also computed for the glutamate-bound T686A and T686S S1S2 mutants, and the results elaborate on findings from experimental functional studies.     

Molecular dynamics and free energy analysis of neuraminidase-ligand interactions

We report molecular dynamics calculations of neuraminidase in complex with an inhibitor, 4-amino-2-deoxy-2,3-didehydro-N-acetylneuraminic acid (N-DANA), with subsequent free energy analysis of binding by using a combined molecular mechanics/continuum solvent model approach. A dynamical model of the complex containing an ionized Glu119 amino acid residue is found to be consistent with experimental data. Computational analysis indicates a major van der Waals component to the inhibitor-neuraminidase binding free energy. Based on the N-DANA/neuraminidase molecular dynamics trajectory, a perturbation methodology was used to predict the binding affinity of related neuraminidase inhibitors by using a force field/Poisson-Boltzmann potential. This approach, incorporating conformational search/local minimization schemes with distance-dependent dielectric or generalized Born solvent models, correctly identifies the most potent neuraminidase inhibitor. Mutation of the key ligand four-substituent to a hydrogen atom indicates no favorable binding free energy contribution of a hydroxyl group; conversely, cationic substituents form favorable electrostatic interactions with neuraminidase. Prospects for further development of the method as an analysis and rational design tool are discussed.     

Conformation of the antifreeze glycoprotein of polar fish

High-field proton and 13C NMR spectroscopy has been used to test and refine the recent proposal, based on vacuum uv circular dichroism results, of a threefold left-handed helical conformation for antifreeze glycoprotein (AFGP). Partial assignment of the protons of the glycotripeptide repeating unit has been made by comparison with spectra of model compounds, by selective decoupling, and by measurements of nuclear Overhauser effect (nOe). At 40 degrees C, AFGP fraction 8 (Mr 2600) shows 2-Hz linewidths which broaden at lower temperature. Neither 1H nor 13C chemical shifts depend strongly on temperature, suggesting no abrupt conformational transition. The nOe between alanine alpha and beta protons vary with temperature and with field strength, from small positive enhancements at 50 degrees C and 80 MHz to large negative effects at 3 degrees C and 300 MHz, indicating a substantial change of rotational correlation time with temperature. The higher-molecular-weight fraction 1-4 shows negative nOe at all temperatures. The CD spectra of fraction 1-4 show bands characteristic of the polyproline II structure at both 3 and 50 degrees C, while those bands in fraction 8 are weaker at 50 than 3 degrees C. The 1H nOe, the 13C T1, and CD data are interpreted as indicating that AFGP fraction 8 is an extended "rod-like" conformation at low temperature which becomes a flexible coil at high temperature, while fraction 1-4 is a flexible rod with sufficient segmental mobility to eliminate any long-range order.     

Spruce budworm antifreeze protein: changes in structure and dynamics at low temperature

Antifreeze proteins (AFPs) prevent the growth of ice, and are used by some organisms that live in sub-zero environments for protection against freezing. All AFPs are thought to function by an adsorption inhibition process. In order to elucidate the ice-binding mechanism, the structures of several AFPs have been determined, and have been shown to consist of different folds. Recently, the first structures of the highly active insect AFPs have been characterized. These proteins have a beta-helix structure, which adds yet another fold to the AFP family. The 90-residue spruce budworm (Choristoneura fumiferana) AFP consists of a beta-helix with 15 residues per coil. The structure contains two ranks of aligned threonine residues (known as the TXT motif), which were shown by mutagenesis experiments to be located in the ice-binding face. In our previous NMR study of this AFP at 30 degrees C, we found that the TXT face was not optimally defined because of the broadening of NMR resonances potentially due to weak oligomerization. We present here a structure of spruce budworm AFP determined at 5 degrees C, where this broadening is reduced. In addition, the 1H-15N NMR dynamics of the protein were examined at 30 degrees C and 5 degrees C. The results show that the spruce budworm AFP is more structured at 5 degrees C, and support the general observation that AFPs become more rigid as the temperature is lowered.     

Membrane glycoprotein and surface free energy changes in hypoxic fibroblast cells

Hypoxia affects the biochemistry of mammalian cells and thus alters their sensitivity to subsequent chemo- and radiotherapy. When V79 Chinese hamster lung fibroblasts were grown under conditions of extreme hypoxia (less than 10 ppm O2) there was a significant shift in the membrane glycoprotein composition. Scanning electron microscopy revealed altered cell surface morphology including loss of pseudopodial projections. Experiments to determine changes in interfacial free energy of these cells using equilibrium two phase systems of poly(ethylene glycol) (PEG) and dextran were carried out. Test fluid droplets of the denser dextran-rich phase were formed on layers of cells in the PEG-rich phase as the bathing medium, and the contact angles the droplets made with the cell layers were measured from photomicrographs. The contact angles on cells in the plateau phase increased significantly with time of exposure to hypoxia, from 25 degrees (zero time) to 35 degrees (6 h) to 60 degrees (9 h). Contact angles on cells in the exponential phase increased from 80 degrees (zero time) to 150 degrees after 20 h of hypoxia. It appears that the altered contact angles reflect changes in cell surface hydrophobicity that may, in part, reflect alterations in the membrane glycoprotein composition.     

Characterization of proline-containing right-handed alpha-helix by molecular dynamics studies

Proline residues play a special role in shaping the secondary and tertiary structures of proteins. Many of these aspects have been studied in great detail. Current interest lies in elucidating the structure of right-handed alpha-helical fragments which contain proline in the middle of the helix. Such structures play an important role in membrane proteins and in the tight packing of globular proteins. Analysis of several crystal structures and energy minimization using flexible geometry have elucidated the nature of the bend produced by proline in the right-handed alpha-helical structure. Molecular dynamics (MD) simulation studies are ideally suited to characterize rigidity or flexibility in different parts of the molecule and can also give an idea of various conformations of the molecule which can exist at a given temperature. Hence, MD studies on Ace-(Ala)6-Pro-(Ala)3-NHMe have been carried out for 100 ps after equilibration and the resulting trajectories have been analyzed. Information regarding the average values, r.m.s. fluctuations of internal parameters and the time spent in different conformations are discussed. Energy minimization has been carried out on selected MD simulated points in order to analyze the characteristics of different conformations.