Decreased expression of surfactant Protein-C and CD74 in alveolar epithelial cells during influenza virus A(H1N1)pdm09 and H3N2 infection.

The primary replication site of Influenza A virus (IAV) is type II alveolar epithelial cells (AECII), which are central to normal lung function and present important immune functions. Surfactant components are synthesized primarily by AECII, which play a crucial role in host defense against infection. The aim of this study was to analyze if the impact of influenza infection is differential between A(H1N1)pdm09 and A/Victoria/3/75 (H3N2) on costimulatory molecules and ProSP-C expression in AECII from BALB/c mice infected and A549 cell line infected with both strains. Pandemic A(H1N1)pdm09 and A/Victoria/3/75 (H3N2) were used to infect BALB/c mice and the A549 cell line. We evaluated the surface expression of co-stimulatory molecules (CD45/CD31/CD74/ProSP-C) in AECII and A549 cell lines. Our results showed a significant decrease in ProSP-C+ CD31- CD45- and CD74+ CD31- CD45- expression in AECII and A549 cell line with the virus strain A(H1N1)pdm09 versus A/Victoria/3/75 (H3N2) and controls (non-infection conditions). Our findings indicate that changes in the expression of ProSP-C in AECII and A549 cell lines in infection conditions could result in dysfunction leading to decreased lung compliance, increased work of breathing and increased susceptibility to injury.

Hydrophobic hydration is an important source of elasticity in elastin-based biopolymers.

Molecular dynamics simulations with explicit waters have been employed to investigate the dominant source of elastin's elasticity. An elastin-like peptide, (VPGVG)(18), was pulled and released in molecular dynamics simulations, at 10 and 42 degrees C, lasting several nanoseconds, which is consistent with the experimentally determined dielectric and NMR relaxation time scales. At elastin's physiological temperature and degree of extension, the simulations indicate that the orientational entropy of waters hydrating hydrophobic groups decreases during pulling of the molecule, but it increases upon release. In contrast, the main-chain fluctuations and other measures of mobility suggest that elastin's backbone is more dynamic in the extended than released state. These results and the agreement between the simulations with various experimental observations suggest that hydrophobic hydration is an important source of the entropy-based elasticity of elastin. Moreover, elastin tends to reorder itself to form a hydrophobic globule when it was held in its extended state, indicating that the hydrophobic effect also contributes in the holding process. On the whole, our simulations support the hydrophobic mechanism of elasticity and provide a framework for description of the molecular basis of this phenomenon.

Mapping the early steps in the pH-induced conformational conversion of the prion protein.

Under certain conditions, the prion protein (PrP) undergoes a conformational change from the normal cellular isoform, PrP(C), to PrP(Sc), an infectious isoform capable of causing neurodegenerative diseases in many mammals. Conversion can be triggered by low pH, and in vivo this appears to take place in an endocytic pathway and/or caveolae-like domains. It has thus far been impossible to characterize the conformational change at high resolution by experimental methods. Therefore, to investigate the effect of acidic pH on PrP conformation, we have performed 10-ns molecular dynamics simulations of PrP(C) in water at neutral and low pH. The core of the protein is well maintained at neutral pH. At low pH, however, the protein is more dynamic, and the sheet-like structure increases both by lengthening of the native beta-sheet and by addition of a portion of the N terminus to widen the sheet by another two strands. The side chain of Met-129, a polymorphic codon in humans associated with variant Creutzfeldt-Jakob disease, pulls the N terminus into the sheet. Neutralization of Asp-178 at low pH removes interactions that inhibit conversion, which is consistent with the Asp-178-Asn mutation causing human prion diseases.

The molecular basis for the inverse temperature transition of elastin.

Elastin undergoes an "inverse temperature transition" such that it becomes more ordered as the temperature increases. To investigate the molecular basis for this behavior, molecular dynamics simulations were conducted above and below the transition temperature. Simulations of a 90-residue elastin peptide, (VPGVG)(18), with explicit water molecules were performed at seven different temperatures between 7 and 42 degrees C, for a total of 80 ns. Beginning from an idealized beta-spiral structure, hydrophobic collapse was observed over a narrow temperature range in the simulations. Moreover, simulations above and below elastin's transition temperature indicate that elastin has more turns and distorted beta-structure at higher temperatures. Water was critical to the inverse temperature transition and elastin-associated water molecules can be divided into three categories: those closely associated with beta II turns; those that form hydrogen bonds with the main-chain groups; and those hydrating the hydrophobic side-chains. Water-swollen, monomeric elastin above the transition temperature is best described as a compact amorphous structure with distorted beta-strands, fluctuating turns, buried hydrophobic residues, and main-chain polar atoms that participate in hydrogen bonds with water. Below the transition temperature, elastin is expanded with approximately 40 % local beta-spiral structure. Overall the simulations are in agreement with experiment and therefore appear to provide an atomic-level description of the conformational properties of elastin monomers and the basis for their elastomeric properties.

Characterization of the unfolding pathway of the cell-cycle protein p13suc1 by molecular dynamics simulations: implications for domain swapping.

BACKGROUND: The p13suc1 gene product is a member of the cks (cyclin-dependent protein kinase subunit) protein family and has been implicated in regulation of the cell cycle. Various crystal structures of suc1 are available, including a globular, monomeric form and a beta-strand exchanged dimer. It has been suggested that conversions between these forms, and perhaps others, may be important in the regulation of the cell cycle. RESULTS: We have undertaken molecular dynamics simulations of protein unfolding to investigate the conformational properties of suc1. Unfolding transition states were identified for each of four simulations. These states contain some native secondary structure, primarily helix alpha1 and the core of the beta sheet. The hydrophobic core is loosely packed. Further unfolding leads to an intermediate state that is slightly more expanded than the transition state, but with considerably fewer nonlocal, tertiary packing contacts and less secondary structure. The helices are fluctuating but partially formed in the denatured state and beta2 and beta4 remain associated. CONCLUSIONS: It appears that suc1 folds by a nucleation-condensation mechanism, similar to that observed for two-state folding proteins. However, suc1 forms an intermediate during unfolding and contains considerable residual structure in the denatured state. The stability of the beta2-beta4 residual structure is surprising, because beta4 is the strand involved in domain swapping. This stability suggests that the domain-swapping event, if physiologically relevant, may require the assistance of additional factors in vivo or occur early in the folding process.

Staphylococcal protein A: unfolding pathways, unfolded states, and differences between the B and E domains.

We present multiple native and denaturation simulations of the B and E domains of the three-helix bundle protein A, totaling 60 ns. The C-terminal helix (H3) consistently denatures later than either of the other two helices and contains residual helical structure in the denatured state. These results are consistent with experiments suggesting that the isolated H3 fragment is more stable than H1 and H2 and that H3 forms early in folding. Interestingly, the denatured state of the B domain is much more compact than that of the E domain. This sequence-dependent effect on the dimensions of the denatured state and the lack of correlation with structure suggest that the radius of gyration can be a misleading reaction coordinate for unfolding/folding. Various unfolding and refolding events are observed in the denaturation simulations. In some cases, the transitions are facilitated through interactions with other portions of the protein-contact-assisted helix formation. In the native simulations, the E domain is very stable: after 6 ns, the C(alpha) root-mean-square deviation from the starting structure is less than 1.4 A. In contrast, the native state of the B domain deviates more and its inter-helical angles fluctuate. In apparent contrast, we note that the B domain is thermodynamically more stable than the E domain. The simulations suggest that the increased stability of the B domain may be due to heightened mobility, and therefore entropy, in the native state and decreased mobility/entropy in the more compact denatured state.

Molecular dynamics simulations of hydrophobic collapse of ubiquitin.

Nine nonnative conformations of ubiquitin, generated during two different thermal denaturation trajectories, were simulated under nearly native conditions (62 degrees C). The simulations included all protein and solvent atoms explicitly, and simulation times ranged from 1-2.4 ns. The starting structures had alpha-carbon root-mean-square deviations (RMSDs) from the crystal structure of 4-12 A and radii of gyration as high as 1.3 times that of the native state. In all but one case, the protein collapsed when the temperature was lowered and sampled conformations as compact as those reached in a control simulation beginning from the crystal structure. In contrast, the protein did not collapse when simulated in a 60% methanol:water mixture. The behavior of the protein depended on the starting structure: during simulation of the most native-like starting structures (<5 A RMSD to the crystal structure) the RMSD decreased, the number of native hydrogen bonds increased, and the secondary and tertiary structure increased. Intermediate starting structures (5-10 A RMSD) collapsed to the radius of gyration of the control simulation, hydrophobic residues were preferentially buried, and the protein acquired some native contacts. However, the protein did not refold. The least native starting structures (10-12 A RMSD) did not collapse as completely as the more native-like structures; instead, they experienced large fluctuations in radius of gyration and went through cycles of expansion and collapse, with improved burial of hydrophobic residues in successive collapsed states.

Theoretical studies of sequence effects on the conformational properties of a fragment of the prion protein: implications for scrapie formation.

BACKGROUND: Prion diseases are neurodegenerative disorders that appear to be due to a conformational change, involving the conversion of alpha-helices in the normal, cellular isoform of the prion protein (PrPC) to beta-structure in the infectious scrapie form (PrPSc). One form of Gerstmann-Sträussler-Scheinker syndrome (GSS), an inherited prion disease, is caused by mutation of Ala117 of PrPC to Val. We therefore set out to evaluate the effects of this mutation on the stability of the PrPC form. RESULTS: We have performed molecular dynamics simulations of a portion of the PrPC sequence (residues 109-122, termed H1) that is proposed to figure prominently in the conversion of PrPC to PrPSc. In particular, beginning with H1 in the alpha-helical state, the conformational consequences of sequence changes at position 117 were investigated for six hydrophobic mutations. Of these, only the Val mutation was helix-destabilizing. Portions of this mutant peptide adopted and retained an extended conformation during a 2 ns simulation of the peptide in water. CONCLUSIONS: The conformational transitions and structures observed in the simulation of the mutant peptide with Val at position 117 provide insight into the possible early steps in the conversion of PrPC to PrPSc.

Molecular dynamics simulations of protein unfolding and limited refolding: characterization of partially unfolded states of ubiquitin in 60% methanol and in water.

Extensive experimental data are available on the native, partially and fully unfolded states of ubiquitin. Two and three-dimensional NMR experiments of a partially unfolded form of the protein in 60% methanol indicate that approximately one-half of the molecule contains disrupted but native-like structure while the other half is unstructured and/or contains non-native structure. In contrast, the interpretation of hydrogen-exchange data have led to the conclusion that this state is native-like. Thus, there are discrepancies between the experimental studies, or interpretations based on the data. We compare the results of molecular dynamics simulations, under varying conditions, with the experimental results. The simulations extend past 0.5 ns and include explicit solvent molecules: either pure water or 60% methanol. To begin with, ubiquitin was thermally denatured in water (at 498 K). Two particular structures, or "aliquots", during the unfolding process were selected for further study (60 and 198 ps). These structures were then simulated separately in water and 60% methanol at a lower and experimentally meaningful temperature (335 K). The conformations generated from the structure extracted later in the simulation contained significant amounts of non-native structure in the presence of methanol while satisfying both the NMR and hydrogen exchange data. In fact, clearly non-native regions of the structure yielded the desired protection from hydrogen exchange. In contrast, an earlier, more native-like, intermediate did not do as well at predicting the hydrogen-exchange behavior and was inconsistent with the NMR data. These data suggest that the results and interpretations using the different experimental techniques can be reconciled by a single state. This finding also brings into question the practice of interpreting protection to hydrogen exchange in terms of native secondary and tertiary structure, especially when one has weak patterns and low protection factors. When the partially unfolded states were placed in pure water, the protein collapsed and began to refold. Therefore, the desired solvent-dependent properties were observed: the partially unfolded conformations with increased exposure of hydrophobic residues remained expanded in methanol but collapsed in water as the non-polar groups minimized their exposure to solvent.

The three states of globular proteins: acid denaturation.

We describe statistical mechanical theory that aims to predict protein stabilities as a function of temperature, pH, and salt concentration, from the physical properties of the constituent amino acids: (1) the number of nonpolar groups, (2) the chain length, (3) the temperature-dependent free energy of transfer, (4) the pKa's (including those in the native state) and their temperature dependencies. We calculate here the phase diagrams for apomyoglobin and hypothetical variant proteins. The theory captures essential features of protein stability including myoglobin's Tm vs pH as measured by P. L. Privalov [(1979) Advances in Protein Chemistry, Vol. 33, pp. 167-241] and its ionic strength vs pH phase diagram as measured by Y. Goto and A. L. Fink [(1990) Journal of Molecular Biology, Vol. 214, pp. 803-805]. The main predictions here are the following: (1) There are three stable states, corresponding to native (N), compact denatured (C), and highly unfolded (U), with transitions between them. (2) In agreement with experiments, the compact denatured state is predicted to have enthalpy closer to U than N because even though there is considerable hydrophobic "clustering" in C, this nevertheless represents a major loss of hydrophobic contacts relative to configurations (N) that have a hydrophobic "core." (3) C becomes more prominent in the phase diagram with increasing nonpolar content or decreasing chain length, perhaps thus accounting for (a) why lysozyme and alpha-lactalbumin differ in their denatured states, and (b) why shortened Staph nuclease molecules are compact. (4) Of major importance for protein calorimetry is Privalov's observation that the enthalpy of folding, delta H (T, pH) is independent of pH. The theory accounts for this through the prediction that the main electrostatic contribution to stability is not enthalpic; the main contribution is the entropy, mainly due to the different distributions of protons and small ions in the native and denatured states.

Solvent denaturation and stabilization of globular proteins.

Statistical thermodynamic theory has recently been developed to account for the stabilities of globular proteins. Here we extend that work to predict the effects of solvents on protein stability. Folding is assumed to be driven by solvophobic interactions and opposed by conformational entropy. The solvent dependence of the solvophobic interactions is taken from transfer experiments of Nozaki and Tanford on amino acids into aqueous solutions of urea or guanidine hydrochloride (GuHCl). On the basis of the assumption of two pathways involving collapse and formation of a core, the theory predicts that increasing denaturant should lead to a two-state denaturation transition (i.e., there is a stable state along each path separated by a free energy barrier). The denaturation midpoint is predicted to occur at higher concentrations of urea than of GuHCl. At neutral pH, the radius of the solvent-denatured state should be much smaller than for a random-flight chain and increase with either denaturant concentration or number of polar residues in the chain. A question of interest is whether free energies of folding should depend linearly on denaturant, as is often assumed. The free energy is predicted to be linear for urea but to have some small curvature for GuHCl. Predicted slopes and exposed areas of the unfolded states are found to be in generally good agreement with experiments. We also discuss stabilizing solvents and compare thermal with solvent denaturation.

Protein stability: electrostatics and compact denatured states.

Globular proteins can be denatured by changing pH and ionic strength. Much recent evidence has led to the surprising conclusion that there are two acid-denatured states: one highly unfolded and the other more compact, sometimes called the "molten globule." Here we describe a molecular theory for electrostatic stability of globular proteins based on the properties of the constituent amino acids: oil/water partition coefficients, pK values of the titratable groups, and their temperature dependences. Predicted denaturation temperatures vs. pH are in good agreement with experiments of other workers on myoglobin. The theory also predicts two populations of denatured species, one open and the other more compact, with densities in the range found experimentally for molten globular states. In addition, it predicts a phase diagram (stability vs. pH, ionic strength) in good agreement with experiments of Goto and Fink [Goto, Y. & Fink, A. L. (1989) Biochemistry 28, 945-952; and Goto, Y. & Fink, A. L. (1990) J. Mol. Biol. 214, 803-805]. The well-known salt destabilization of myoglobin has been generally considered evidence for ion pairing, but the present theory, based on smeared charge repulsion, explains the salt destabilization at low pH without ion pairing. In addition, for myoglobin the theory predicts salt stabilization at high pH, as observed for beta-lactamase by Goto and Fink.

Thermal stabilities of globular proteins.

Statistical thermodynamic theory has recently been developed to account for the stabilities of globular proteins. Here we extend that work to predict the dependence on temperature. Folding is assumed to be driven by solvophobic interactions and opposed by the conformational entropy. The temperature dependence of the solvophobic interaction is taken from the transfer experiments on amino acids by Tanford and Nozaki and on model solutes by Gill and Wadsö. One long-standing puzzle has been why proteins denature upon heating, since the solvophobic force to fold strengthens with increasing temperature. This is resolved by the theory, which predicts two first-order phase transitions. "Cold denaturation" is driven principally by the weakening of the solvophobic interaction, but normal denaturation is driven principally by the gain of conformational entropy of the chain. Predictions of the thermodynamic state functions are in reasonable agreement with the calorimetric experiments of Privalov and Khechinashvili. Comparison of the theory with experiments suggests that there may be an additional enthalpic driving force toward folding which is not due to the solvophobic interactions.