Partial molar volumes and adiabatic compressibilities of unfolded protein states

Volumetric Properties of Four-Stranded DNA Structures

Four-stranded non-canonical DNA structures including G-quadruplexes and i-motifs have been found in the genome and are thought to be involved in regulation of biological function. These structures have been implicated in telomere biology, genomic instability, and regulation of transcription and translation events. To gain an understanding of the molecular determinants underlying the biological role of four-stranded DNA structures, their biophysical properties have been extensively studied. The limited libraries on volume, expansibility, and compressibility accumulated to date have begun to provide insights into the molecular origins of helix-to-coil and helix-to-helix conformational transitions involving four-stranded DNA structures. In this article, we review the recent progress in volumetric investigations of G-quadruplexes and i-motifs, emphasizing how such data can be used to characterize intra-and intermolecular interactions, including solvation. We describe how volumetric data can be interpreted at the molecular level to yield a better understanding of the role that solute-solvent interactions play in modulating the stability and recognition events of nucleic acids. Taken together, volumetric studies facilitate unveiling the molecular determinants of biological events involving biopolymers, including G-quadruplexes and i-motifs, by providing one more piece to the thermodynamic puzzle describing the energetics of cellular processes in vitro and, by extension, in vivo.

Pressure Unfolding of Proteins: New Insights into the Role of Bound Water

High pressures can be detrimental for protein stability, resulting in unfolding and loss of function. This phenomenon occurs because the unfolding transition is accompanied by a decrease in volume, which is typically attributed to the elimination of cavities that are present within the native state as a result of packing defects. We present a novel computational approach that enables the study of pressure unfolding in atomistically detailed protein models in implicit solvent. We include the effect of pressure using a transfer free energy term that allows us to decouple the effect of protein residues and bound water molecules on the volume change upon unfolding. We discuss molecular dynamics simulations results using this protocol for two model proteins, Trp-cage and staphylococcal nuclease (SNase). We find that the volume reduction of bound water is the key energetic term that drives protein denaturation under the effect of pressure, for both Trp-cage and SNase. However, we note differences in unfolding mechanisms between the smaller Trp-cage and the larger SNase protein. Indeed, the unfolding of SNase, but not Trp-cage, is seen to be further accompanied by a reduction in the volume of internal cavities. Our results indicate that, for small peptides, like Trp-cage, pressure denaturation is driven by the increase in solvent accessibility upon unfolding, and the subsequent increase in the number of bound water molecules. For larger proteins, like SNase, the cavities within the native fold act as weak spots, determining the overall resistance to pressure denaturation. Our simulations display a striking agreement with the pressure-unfolding profile experimentally obtained for SNase and represent a promising approach for a computationally efficient and accurate exploration of pressure-induced denaturation of proteins.

Toward Biotherapeutics Formulation Composition Engineering using Site-Identification by Ligand Competitive Saturation (SILCS)

Formulation of protein-based therapeutics employ advanced formulation and analytical technologies for screening various parameters such as buffer, pH, and excipients. At a molecular level, physico-chemical properties of a protein formulation depend on self-interaction between protein molecules, protein-solvent and protein-excipient interactions. This work describes a novel in silico approach, SILCS-Biologics, for structure-based modeling of protein formulations. SILCS Biologics is based on the Site-Identification by Ligand Competitive Saturation (SILCS) technology and enables modeling of interactions among different components of a formulation at an atomistic level while accounting for protein flexibility. It predicts potential hotspot regions on the protein surface for protein-protein and protein-excipient interactions. Here we apply SILCS-Biologics on a Fab domain of a monoclonal antibody (mAbN) to model Fab-Fab interactions and interactions with three amino acid excipients, namely, arginine HCl, proline and lysine HCl. Experiments on 100 mg/ml formulations of mAbN showed that arginine increased, lysine reduced, and proline did not impact viscosity. We use SILCS-Biologics modeling to explore a structure-based hypothesis for the viscosity modulating effect of these excipients. Current efforts are aimed at further validation of this novel computational framework and expanding the scope to model full mAb and other protein therapeutics.

Monitoring of lysozyme thermal denaturation by volumetric measurements and nanoDSF technique in the presence of N-butylurea

The results of thermal studies of denaturation of hen egg white lysozyme (HEWL) in water and an aqueous solution of N-butylurea (BU) are presented. High-precision densimetric measurements were used to characterize and analyze the changes of the specific volume, v, during temperature elevation. The temperature of the midpoint of protein denaturation was also determined by nanoDSF technique (differential scanning fluorimetry). The densities of lysozyme solutions were measured at temperatures ranging from 298.15 to 353.15 K with an interval of 5 K at atmospheric pressure (0.1 MPa). The concentration of the protein covered the range from 2 to 20 mg per 1 ml of the solution. The optimal range of the concentration for the densimetric measurements was roughly estimated. In the transition region, the structural changes of the protein are accompanied by the biggest increase of ν values with temperature. Our measurements show that this effect can be monitored from volumetric data without precise determination of protein concentration. The results prove that a two-state model of denaturation could be used for data interpretation. Contrary to common misconception, the volumetric measurements suggest that the denatured protein does not necessarily need to be in a fully extended state. In this way, the 'protein volume paradox' could be explained. The surface area of the protein remains unchanged and thus the increase of the specific volume of the protein is relatively small. Additionally, the self-stabilizing effect of the protein in BU solution was reported. For the HEWL in pure water, this phenomenon was not observed.

On empirical decomposition of volumetric data

Volumetric characterization of proteins and their recognition events has been instrumental in providing information on the role of intra- and intermolecular interactions, including hydration, in stabilizing biomolecules. The credibility of molecular models and interpretation schemes used to rationalize experimental data are essential for the validity of microscopic insights derived from volumetric results. Current empirical schemes used to interpret volumetric data suffer from a lack of theoretical and computational substantiation. In this contribution, we take advantage age of recent MD simulations of proteins in solution coupled with Voronoi-Delaunay tessellation of simulated structures that have provided an exceptional level of structural detail on the nature of protein-water interfaces. We use these structural insights to re-evaluate empirical frameworks used for interpretation of volumetric data. An important issue in this respect is the actual dividing surface between water and protein atoms that is used in volumetric studies when the solute and solvent are treated as hard spheres enclosed within their respective van der Waals surfaces. In one development, using Voronoi tessellation of MD simulated protein-water systems the dividing surface has been defined as the points equidistant from the water and protein atoms. The interstitial void volume between the solute and the dividing surface corresponds to thermal volume envisaged by Scaled Particle Theory. In this communication, we explicitly account for the contributions of thermal volume to the partial molar volume, compressibility, and expansibility of proteins and re-examine and redefine the intrinsic and hydration volumetric contributions. We discuss the implications of our results for protein transitions and association events.

Molecular Determinants of Temperature Dependence of Protein Volume Change upon Unfolding

Pressure is a well-known environmental stressor that can either stabilize or destabilize proteins. The volumetric change upon protein unfolding determines the effect of pressure on protein stability, where negative volume changes destabilized proteins at high pressures. High temperature often accompanies high pressure, for example, in the ocean depths near hydrothermal vents or near faults, so it is important to study the effect of temperature on the volumetric change upon unfolding. We previously detailed the magnitude and sign of the molecular determinants of volumetric change, allowing us to quantitatively predict the volumetric change upon protein unfolding. Here, we present a comprehensive analysis of the temperature dependence of the volumetric components of proteins, showing that hydration volume is the primary component that defines expansivities of the native and unfolded states and void volume only contributes slightly to the folded state expansivity.

Molecular determinant of the effects of hydrostatic pressure on protein folding stability

Hydrostatic pressure is an important environmental variable that plays an essential role in biological adaptation for many extremophilic organisms (for example, piezophiles). Increase in hydrostatic pressure, much like increase in temperature, perturbs the thermodynamic equilibrium between native and unfolded states of proteins. Experimentally, it has been observed that increase in hydrostatic pressure can both increase and decrease protein stability. These observations suggest that volume changes upon protein unfolding can be both positive and negative. The molecular details of this difference in sign of volume changes have been puzzling the field for the past 50 years. Here we present a comprehensive thermodynamic model that provides in-depth analysis of the contribution of various molecular determinants to the volume changes upon protein unfolding. Comparison with experimental data shows that the model allows quantitative predictions of volume changes upon protein unfolding, thus paving the way to proteome-wide computational comparison of proteins from different extremophilic organisms.

Proteins can be both stabilized and destabilized by pressure. Here the authors analyse the factors contributing to both negative and positive protein volume change upon denaturation, and shed light on the molecular determinants allowing proteins to be stable at high pressures.

Volume and Compressibility of Proteins

The partial specific (or molar) volume, expansibility, and compressibility of a protein are fundamental thermodynamic quantities for characterizing its structure in solution. We review the definitions, measurements, and implications of these volumetric quantities in relation to protein structural biology. The partial specific volumes under constant molality (isomolal) and chemical potential (isopotential) conditions of the cosolvent (multicomponent systems) are explained in terms of preferential solvent interactions relevant to the solubility and stability of proteins. The partial expansibility is briefly discussed in terms of the effects of temperature on protein-solvent interactions (hydration) and internal packing defects (cavities). We discuss the compressibility-structure-function relationships of proteins based on analyses of the correlations between the partial adiabatic compressibilities and the structures or functions of various globular proteins (including mutants), focusing on the roles of the internal cavities in structural fluctuations. The volume and compressibility changes associated with various conformational transitions are also discussed in terms of the changes in hydration and cavities in order to elucidate the nonnative structures and the transition mechanisms, especially those associated with pressure denaturation.

Driving Forces in Pressure-Induced Protein Transitions

The molecular mechanisms underlying pressure-induced protein denaturation can be analyzed based on the pressure-dependent differences in the apparent volume occupied by amino acids inside the protein and when exposed to water in an unfolded conformation. This chapter presents a volumetric analysis of the peptide group and the 20 naturally occurring amino acid side chains in the interior of the native state, the micelle-like interior of the pressure-induced denatured state, and in the unfolded conformation modeled by low-molecular analogs of proteins. The transfer of a peptide group from the protein interior to water becomes increasingly favorable as pressure increases. This observation classifies solvation of peptide groups as a major driving force in pressure-induced protein denaturation. Polar side chains do not appear to exhibit significant pressure-dependent changes in their preference for the protein interior or solvent. The transfer of nonpolar side chains from the protein interior to water becomes more unfavorable as pressure increases. An inference can be drawn that a sizeable population of nonpolar side chains remains buried inside a solvent-inaccessible core of the pressure-induced denatured state. At elevated pressures this core, owing to the absence of structural constraints, may become packed almost as tightly as the interior of the native state. The presence and partial disappearance of large intraglobular voids is another driving force facilitating pressure-induced protein denaturation. Volumetric data presented here have implications for the kinetics of protein folding and shed light on the nature of the folding transition state ensembles.

Profiling of substrate specificities of 3C-like proteases from group 1, 2a, 2b, and 3 coronaviruses

Background

Coronaviruses (CoVs) can be classified into alphacoronavirus (group 1), betacoronavirus (group 2), and gammacoronavirus (group 3) based on diversity of the protein sequences. Their 3C-like protease (3CL^pro), which catalyzes the proteolytic processing of the polyproteins for viral replication, is a potential target for anti-coronaviral infection.

Methodology/Principal Findings

Here, we profiled the substrate specificities of 3CL^pro from human CoV NL63 (group 1), human CoV OC43 (group 2a), severe acute respiratory syndrome coronavirus (SARS-CoV) (group 2b) and infectious bronchitis virus (IBV) (group 3), by measuring their activity against a substrate library of 19×8 of variants with single substitutions at P5 to P3' positions. The results were correlated with structural properties like side chain volume, hydrophobicity, and secondary structure propensities of substituting residues. All 3CL^pro prefer Gln at P1 position, Leu at P2 position, basic residues at P3 position, small hydrophobic residues at P4 position, and small residues at P1' and P2' positions. Despite 3CL^pro from different groups of CoVs share many similarities in substrate specificities, differences in substrate specificities were observed at P4 positions, with IBV 3CL^pro prefers P4-Pro and SARS-CoV 3CL^pro prefers P4-Val. By combining the most favorable residues at P3 to P5 positions, we identified super-active substrate sequences ‘VARLQ↓SGF’ that can be cleaved efficiently by all 3CL^pro with relative activity of 1.7 to 3.2, and ‘VPRLQ↓SGF’ that can be cleaved specifically by IBV 3CL^pro with relative activity of 4.3.

Conclusions/Significance

The comprehensive substrate specificities of 3CL^pro from each of the group 1, 2a, 2b, and 3 CoVs have been profiled in this study, which may provide insights into a rational design of broad-spectrum peptidomimetic inhibitors targeting the proteases.

Profiling of substrate specificity of SARS-CoV 3CL

BACKGROUND

The 3C-like protease (3CL(pro)) of severe acute respiratory syndrome-coronavirus is required for autoprocessing of the polyprotein, and is a potential target for treating coronaviral infection.

METHODOLOGY/PRINCIPAL FINDINGS

To obtain a thorough understanding of substrate specificity of the protease, a substrate library of 198 variants was created by performing saturation mutagenesis on the autocleavage sequence at P5 to P3' positions. The substrate sequences were inserted between cyan and yellow fluorescent proteins so that the cleavage rates were monitored by in vitro fluorescence resonance energy transfer. The relative cleavage rate for different substrate sequences was correlated with various structural properties. P5 and P3 positions prefer residues with high β-sheet propensity; P4 prefers small hydrophobic residues; P2 prefers hydrophobic residues without β-branch. Gln is the best residue at P1 position, but observable cleavage can be detected with His and Met substitutions. P1' position prefers small residues, while P2' and P3' positions have no strong preference on residue substitutions. Noteworthy, solvent exposed sites such as P5, P3 and P3' positions favour positively charged residues over negatively charged one, suggesting that electrostatic interactions may play a role in catalysis. A super-active substrate, which combined the preferred residues at P5 to P1 positions, was found to have 2.8 fold higher activity than the wild-type sequence.

CONCLUSIONS/SIGNIFICANCE

Our results demonstrated a strong structure-activity relationship between the 3CL(pro) and its substrate. The substrate specificity profiled in this study may provide insights into a rational design of peptidomimetic inhibitors.

Molecular determinants of the pKa values of Asp and Glu residues in staphylococcal nuclease

Prior computational studies of the acid-unfolding behavior of staphylococcal nuclease (SNase) suggest that the pK(a) values of its carboxylic groups are difficult to reproduce with electrostatics calculations with continuum methods. To examine the molecular determinants of the pK(a) values of carboxylic groups in SNase, the pK(a) values of all 20 Asp and Glu residues were measured with multidimensional and multinuclear NMR spectroscopy in an acid insensitive variant of SNase. The crystal structure of the protein was obtained to describe the microenvironments of the carboxylic groups. Fourteen Asp and Glu residues titrate with relatively normal pK(a) values that are depressed by less than 1.1 units relative to the normal pK(a) of Asp and Glu in water. Only six residues have pK(a) values shifted by more than 1.5 units. Asp-21 has an unusually high pK(a) of 6.5, which is probably the result of interactions with other carboxylic groups at the active site. The most perturbed pK(a) values appear to be governed by hydrogen bonding and not by Coulomb interactions. The pK(a) values calculated with standard continuum electrostatics methods applied to static structures are more depressed than the measured values because Coulomb effects are exaggerated in the calculations. The problems persist even when the protein is treated with the dielectric constant of water. This can be interpreted to imply that structural relaxation is an important determinant of the pK(a) values; however, no major pH-sensitive conformational reorganization of the backbone was detected using NMR spectroscopy.

Use of pressure perturbation calorimetry to characterize the volumetric properties of proteins

Pressure perturbation calorimetry (PPC) is a new technique that makes possible to study the volumetric changes that occur upon protein unfolding. Here, we summarize the thermodynamic foundation of the method and introduce a two-state model for the analysis of the unfolding data monitored by PPC. Several examples of data analysis illustrating potential pitfalls and solutions are discussed.