A statistical mechanical model for hydrogen exchange in globular proteins.

We develop a statistical mechanical theory for the mechanism of hydrogen exchange in globular proteins. Using the HP lattice model, we explore how the solvent accessibilities of chain monomers vary as proteins fluctuate from their stable native conformations. The model explains why hydrogen exchange appears to involve two mechanisms under different conditions of protein stability; (1) a “global unfolding” mechanism by which all protons exchange at a similar rate, approaching that of the denatured protein, and (2) a “stable-state” mechanism by which protons exchange at rates that can differ by many orders of magnitude. There has been some controversy about the stable-state mechanism: does exchange take place inside the protein by solvent penetration, or outside the protein by the local unfolding of a subregion? The present model indicates that the stable-state mechanism of exchange occurs through an ensemble of conformations, some of which may bear very little resemblance to the native structure. Although most fluctuations are small-amplitude motions involving solvent penetration or local unfolding, other fluctuations (the conformational distant relatives) can involve much larger transient excursions to completely different chain folds.     

Comparison of theoretical and experimental particle diffusion data within human airway casts

The dose delivered to airway cells is a critical factor whether one is addressing the therapeutic (i.e., positive) effects of inhaled pharmacologic agents or the toxic (i.e., negative) effects of pollutants. In this study, theoretical models describing particle deposition have been compared with experimental data from the literature. In the simulations, airways can be either roughor smooth-walled to be consistent with human lungs which can be either lined by cartilaginous rings (i.e., rough) or muscle (i.e., smooth). Particle motion for rough-walled airways within generations I=1–6 is calculated using the formula proposed by Martonen et al. (1). For smooth-walled airways within generations I=7–10, particle motion is calculated using the formula proposed by Martonen et al. (2). Theoretical predictions of particle deposition efficiencies are not only in agreement with the overall best fit empirical correlation presented by Cohen and Asgharian (3) over a wide range of dimensionless diffusion parameters, but also match individual experimental measurements (only available in I=1–6) with regard to effects of the parameters of particle size, flow rate, and airway dimensions. The mean difference in the ratio of experimental-to-theoretical particle diffusion values is 0.9 for a flow rate of 18 L/min and 1.1 for a flow rate of 34 L/min (i.e., the difference is only about 10%) within the upper airways of the casts (airway generations I=1–6), the mean difference for the whole casts was much greater. This may be attributed to the uncertainty of flow conditions in the peripheral airways as a result of the trimmmed nature of the casts. Overall, the findings suggest that the model can be a valuable component of aerosol therapy and risk assessment protocols, especially to address effects of enhanced deposition of pharmacologic drugs and radionuclides at sites within the human tracheobronchial tree.     

Dynamics of hydrogen-deuterium exchange in Chlamydomonas centrin

Chlamydomonas reinhardtii centrin is a 169-amino acid residue calcium binding protein belonging to the EF-hand protein superfamily. Centrin is associated with the microtubule organizing center (MTOC) in all eukaryotes, and in Chlamydomonas, centrin is a component of the flagellar basal body apparatus. Recombinant full-length centrin, calmodulin, and terminal domain fragments [Ccen-N (residues 1-94) and Ccen-C (residues 99-169)] were used to examine hydrogen-deuterium (H --> D) exchange dynamics using combined attenuated total reflectance (ATR) Fourier transform-infrared (FT-IR) spectroscopy, curve fit, and two-dimensional correlation analysis. Analysis of the Ccen-N and Ccen-C fragments allowed separation of domain specific solvent exchange events and together with analysis of the full-length proteins provides novel insight into domain accessibility to the aqueous environment and the internal dynamics of the protein.     

A theoretical approach to the binding of amphipathic molecules to globular proteins.

1. A theory based on a multiple equilibrium model with stoicheiometric binding constants has been formulated. It is applicable to the interaction of amphipathic molecules with charged macromolecules. 2. The theory has been applied to the binding of sodium n-dodecyl sulphate to ribonuclease A, beta-lactoglobulin and bovine serum albumin. 3. Over the ranges of surfactant concentration where binding is non-co-operative and co-operative, the experimental data can be satisfactorily fitted with energy and co-operatively parameters which are of comparable magnitude for the three globular proteins. 4. The results imply that the energy of interaction of sodium n-dodecyl sulphate with the proteins is equivalent to the formation of approximately four CH2 hydrophobic bonds.     

Topology of globular proteins

This paper inquires whether it is reasonable to expect the native structure of proteins to be “knotted”. To this end, some topological properties of polypeptides containing disulfide bridges are discussed using notions from mathematical knot theory and graph theory. The probability of occurrence of knots in random cyclic polymers is calculated as a function of chain length by elementary Monte Carlo methods. The implications of this for protein renaturation and for determining the tertiary structure of proteins are discussed.     

Thermodynamics of globular proteins

The analysis of temperature-induced unfolding of proteins in aqueous solutions was performed. Based on the data of thermodynamic parameters of protein unfolding and using the method of semi-empirical calculations of hydration parameters at reference temperature 298 K, we obtained numerical values of enthalpy, free energy, and entropy which characterize the unfolding of proteins in the ‘gas phase’. It was shown that specific values of the energy of weak intramolecular bonds (?H^int), conformational free energy (?G^conf) and entropy (?S^conf) are the same for proteins with molecular weight 7–25 kDa. Using the energy value (?H^int) and the proposed approach for estimation of the conformational entropy of native protein (S^NC), numerical values of the absolute free energy (G^NC) were obtained.     

Potentiometric titration studies on globular proteins

A power series method was applied to solve the Poisson-Boltzmann equation for the spherical polyelectrolyte model and numerical calculation with an electronic computer was performed to obtain surface electric potential on rigid globular proteins. Deviation from the ideal linear relationship in Linderstrom-Lang's plot was found to become noticeable as the surface charge density and the radius of protein increases and ionic strength decreases. The calculated surface potential was compared with potentiometric titration data of several proteins whose radii have been analyzed. Assuming the radius of the counterions to be equal to about 1.0 Å, the data for phenolic groups in ribonuclease and for carboxyl groups in conalbumin were interpreted. Reversible intramolecular transformation was found for α-lactalbumin by comparing the present results with the potentiometric titration data for carboxyl groups. The molecular size of each protein was discussed.     

Effects of net charge on the hydrogen-deuterium exchange parameters of lysozyme.

Possible effects of changes in net charge on protein hydrogen exchange rates were investigated by desalting hen egg-white lysozyme, which allowed its net charge to increase with decreasing pH in the acid region. Chloride ion-binding ratios, expressed as ratios of free to total Cl^?, were measured with a chloride-specific electrode at pH 5 on a 2.4% solution of a five-time-desalted product. This ratio was used to show a 97% reduction of the 11% Cl^? present in a commercial lysozyme preparation upon three passes of the enzyme through a column of ion-retardation resin. Net charges on the purified product were assigned from a combination of electrophoretic mobility and proton titration data gathered under minimal ionic strength conditions. The net charge on the desalted product increased by 1.64 units between pH 5.0 and 3.0. Hydrogendeuterium exchange studies on the purified lysozyme in D₂O were obtained using the near-infrared region of a Cary 14R spectrophotometer. The rate-pD profile for k₂, the rate constant for the intermediate class of exchanging hydrogens, showed a decrease in the apparent pD of minimum exchange rate of 0.3 units, when compared to that obtained earlier in 0.2 m added NaCl. However, the rate of exchange at pD minimum and the number of hydrogens in the class remained largely unaffected. A similar shift was observed for the rate-pD profile of the class 1 hydrogens. Thus, the effect of an increase in net positive charge is to shift the rate-pD profile to a lower pD. Moreover, the effect extended to the interior peptide hydrogens of this globular protein. Consequently, the exchange rates of all the observable hydrogens are altered by the net charge changes, and the effect appeared uniform. The shift can be accounted for quantitatively by applying electrostatic interaction terms to the acid and base catalytic constants characterizing the exchange process. The calculated electrostatic interaction factors in minimal salt and 0.2 m added NaCl were found to be 29 and 18% lower, respectively, than those obtained theoretically. Therefore, under conditions where changes in net charge may occur for a globular protein, the effect on hydrogen exchange rates can be estimated fairly well theoretically, especially at moderate ionic strengths.     

On the mechanism of hydrogen-deuterium exchange in bacteriorhodopsin.

Continuous-flow resonance Raman experiments carried out in bacteriorhodopsin show that the exchange of a deuteron on the Schiff base with a proton takes place in times shorter than 3 ms. Exchange mechanisms based on a base-catalyzed deprotonation followed by reprotonation of the Schiff base are excluded. A mechanism is suggested in which a water molecule interacts directly with the Schiff base deuteron in a concerted exchange mechanism. It appears that in the dark, the binding site is more accessible to neutral water molecules than to charged protons.     

Theory for the folding and stability of globular proteins

Using lattice statistical mechanics, we develop theory to account for the folding of a heteropolymer molecule such as a protein to the globular and soluble state. Folding is assumed to be driven by the association of solvophobic monomers to avoid solvent and opposed by the chain configurational entropy. Theory predicts a phase transition as a function of temperature or solvent character. Molecules that are too short or too long or that have too few solvophobic residues are predicted not to fold. Globular molecules should have a largely solvophobic core, but there is an entropic tendency for some residues to be "out of place", particularly in small molecules. For long chains, molecules comprised of globular domains are predicted to be thermodynamically more stable than spherical molecules. The number of accessible conformations in the globular state is calculated to be an exceedingly small fraction of the number available to the random coil. Previous estimates of this number, which have motivated kinetic theories of folding, err by many tens of orders of magnitude.     

Conformational stability of globular proteins

The conformational stability of ribonuclease T1 has been measured as a function of the variables of most interest to biochemists: temperature, pH, salt concentration, disulfide-bond content and amino acid sequence. The results provide insight into the forces that stabilize globular proteins.     

Nanofibers made of globular proteins

Strong nanofibers composed entirely of a model globular protein, namely, bovine serum albumin (BSA), were produced by electrospinning directly from a BSA solution without the use of chemical cross-linkers. Control of the spinnability and the mechanical properties of the produced nanofibers was achieved by manipulating the protein conformation, protein aggregation, and intra/intermolecular disulfide bonds exchange. In this manner, a low-viscosity globular protein solution could be modified into a polymer-like spinnable solution and easily spun into fibers whose mechanical properties were as good as those of natural fibers made of fibrous protein. We demonstrate here that newly formed disulfide bonds (intra/intermolecular) have a dominant role in both the formation of the nanofibers and in providing them with superior mechanical properties. Our approach to engineer proteins into biocompatible fibrous structures may be used in a wide range of biomedical applications such as suturing, wound dressing, and wound closure.