A large compressibility change of protein induced by a single amino acid substitution.

The adiabatic compressibility (beta s) was determined, by means of the precise sound velocity and density measurements, for a series of single amino acid substituted mutant enzymes of Escherichia coli dihydrofolate reductase (DHFR) and aspartate aminotransferase (AspAT). Interestingly, the beta s values of both DHFR and AspAT were influenced markedly by the mutations at glycine-121 and valine-39, respectively, in which the magnitude of the change was proportional to the enzyme activity. This result demonstrates that the local change of the primary structure plays an important role in atomic packing and protein dynamics, which leads to the modified stability and enzymatic function. This is the first report on the compressibility of mutant proteins.     

Compressibility of protein transitions

We review the results of compressibility studies on proteins and low molecular weight compounds that model the hydration properties of these biopolymers. In particular, we present an analysis of compressibility changes accompanying conformational transitions of globular proteins. This analysis, in conjunction with experimental compressibility data on protein transitions, were used to define the changes in the hydration properties and intrinsic packing associated with native-to-molten globule, native-to-partially unfolded, and native-to-fully unfolded transitions of globular proteins. In addition, we discuss the molecular origins of predominantly positive changes in compressibility observed for pressure-induced denaturation transitions of globular proteins. Throughout this review, we emphasize the importance of compressibility data for characterizing protein transitions, while also describing how such data can be interpreted to gain insight into role that hydration and intrinsic packing play in modulating the stability of and recognition between proteins and other biologically important compounds.     

Binding of bovine pancreatic trypsin inhibitor to trypsinogen: spectroscopic and volumetric studies

We have investigated the binding of bovine pancreatic trypsin inhibitor (BPTI) to bovine trypsinogen by combining ultrasonic velocimetry, high precision densimetry, and fluorescence spectroscopy. We report the changes in volume, adiabatic compressibility, van't Hoff enthalpy, entropy, and free energy that accompany the association of the two proteins at 25 degrees C and pH 8.0. We have used the measured changes in volume and compressibility in conjunction with available structural data to characterize the binding-induced changes in the hydration properties and intrinsic packing of the two proteins. Our estimate reveals that 110 +/- 40 water molecules become released to the bulk from the hydration shells of BPTI and trypsinogen. Furthermore, we find that the intrinsic coefficient of adiabatic compressibility of the two proteins decreases by 14 +/- 2%, which is suggestive of the binding-induced rigidification of the proteins' interior. BPTI-trypsinogen association is an entropy-driven event which proceeds with an unfavorable change in enthalpy. The favorable change in entropy results from partial compensation between two predominant terms. Namely, a large favorable change in hydrational entropy slightly prevails over a close in magnitude but opposite in sign change in configurational entropy. The reduction in configurational entropy and, consequently, protein dynamics is consistent with the observed decrease in intrinsic compressibility. In general, results of this work emphasize the vital role that water plays in modulating protein recognition events.     

THz time scale structural rearrangements and binding modes in lysozyme-ligand interactions

Predicting the conformational changes in proteins that are relevant for substrate binding is an ongoing challenge in the aim of elucidating the functional states of proteins. The motions that are induced by protein-ligand interactions are governed by the protein global modes. Our measurements indicate that the detected changes in the global backbone motion of the enzyme upon binding reflect a shift from the large-scale collective dominant mode in the unbound state towards a functional twisting deformation that assists in closing the binding cleft. Correlated motion in lysozyme has been implicated in enzyme function in previous studies, but detailed characterization of the internal fluctuations that enable the protein to explore the ensemble of conformations that ultimately foster large-scale conformational change is yet unknown. For this reason, we use THz spectroscopy to investigate the picosecond time scale binding modes and collective structural rearrangements that take place in hen egg white lysozyme (HEWL) when bound by the inhibitor (NAG) ₃. These protein thermal motions correspond to fluctuations that have a role in both selecting and sampling from the available protein intrinsic conformations that communicate function. Hence, investigation of these fast, collective modes may provide knowledge about the mechanism leading to the preferred binding process in HEWL-(NAG) ₃. Specifically, in this work we find that the picosecond time scale hydrogen-bonding rearrangements taking place in the protein hydration shell with binding modify the packing density within the hydrophobic core on a local level. These localized, intramolecular contact variations within the protein core appear to facilitate the large cooperative movements within the interfacial region separating the α- and β- domain that mediate binding. The THz time-scale fluctuations identified in the protein-ligand system may also reveal a molecular mechanism for substrate recognition.     

Volumetric characterization of tri-N-acetylglucosamine binding to lysozyme

Volumetric characteristics of protein recognition events determine the direction of pressure-induced shifts in the recognition reaction, while also providing insights into the structural, dynamic, and hydration changes. We report changes in volume, ΔV, and adiabatic compressibility, ΔK(S), accompanying the binding of tri-N-acetylglucosamine [(GlcNAc)(3)] to lysozyme at 25 °C in a pH 5.5 sodium acetate buffer. We interpret our measured changes in volume and compressibility in terms of changes in hydration and dynamic properties of the protein. On the basis of our ΔV data, we find that 79 ± 44 water molecules are released to the bulk from the hydration shells of the protein and the ligand. Our ΔK(S) data suggest a 4 ± 2% decrease in the mean-square fluctuations of the intrinsic volume of the protein, <δV(M)(2)> (or a 2% decrease in δV(M)). Thus, the trisaccharide-bound state of the enzyme is less hydrated, more rigid, and less dynamic compared to the unbound state. In general, we discuss the importance of volumetric insights into the molecular origins of protein recognition events.     

Metal coordination influences substrate binding in horseradish peroxidase

To clarify the role of metal ion coordination in horseradish peroxidase C (HRPC), the effect of pressure and of an externally applied electric field on spectral holes was compared for both metal-free and Mg-mesoporphyrin-substituted horseradish peroxidase C (MP-HRP and MgMP-HRP), as affected by the binding of 2-naphthohydroxamic acid (NHA). The data are compared to earlier studies performed on the same derivatives. Results obtained for MP-HRP show the presence of a predominant MP tautomer, as well as that of another small population with different pocket field and isothermal compressibility (0.12 vs 0.24 GPa⁻¹). Binding NHA induces the formation of two new almost equal populations of MP-HRP tautomer complexes and the protein compressibility in both forms is increased to 0.50 and 0.36 GPa⁻¹. The protein structure becomes much softer than in the absence of NHA. Binding the same substrate to MgMP-HRP resulted in MgMP adopting a single conformation with no compressibility changes, while without NHA, two forms were possible. Stark effect results show charge rearrangement upon substrate binding in both cases. We propose that it is the presence of the metal that stabilizes the structure during the reorganization of the protein matrix induced by the substrate binding event. With the metal, only one conformation is adopted, without significant structural rearrangement but with charge redistribution. The dissociation constants determined for NHA binding to both derivatives and to native HRPC show that studies using mesoporphyrin and Mg-mesoporphyrin derivatives are relevant to investigating the specificity of the substrate-binding pocket in this enzyme. Received: 15 October 1999 / Revised version: 3 April 2000 / Accepted: 5 April 2000     

Softening of the packing density of horseradish peroxidase by a H-donor bound near the heme pocket

We use pressure tuning of spectral holes to estimate the compressibility of protein molecules by optical means. We found that the compressibility of mesoporphyrin-substituted horseradish peroxidase increases by a factor of three when it incorporates small aromatic H-donor molecules that bind in the vicinity of its heme pocket. Such a dramatic softening of its packing density corresponds to a jump from a compressibility range characteristic for the solid state into that characteristic for liquids.     

Nicotine and carbamylcholine binding to nicotinic acetylcholine receptors as studied in AChBP crystal structures

Nicotinic acetylcholine receptors are prototypes for the pharmaceutically important family of pentameric ligand-gated ion channels. Here we present atomic resolution structures of nicotine and carbamylcholine binding to AChBP, a water-soluble homolog of the ligand binding domain of nicotinic receptors and their family members, GABAA, GABAC, 5HT3 serotonin, and glycine receptors. Ligand binding is driven by enthalpy and is accompanied by conformational changes in the ligand binding site. Residues in the binding site contract around the ligand, with the largest movement in the C loop. As expected, the binding is characterized by substantial aromatic and hydrophobic contributions, but additionally there are close contacts between protein oxygens and positively charged groups in the ligands. The higher affinity of nicotine is due to a main chain hydrogen bond with the B loop and a closer packing of the aromatic groups. These structures will be useful tools for the development of new drugs involving nicotinic acetylcholine receptor-associated diseases.     

Effects of oxidation on changes of compressibility of bovine serum albumin

The methods of ultrasound velocity and density measurements were used to study the adiabatic compressibility of bovine serum albumin (BSA) during its oxidation by the prooxidants Cu2+ and 2,2'-azobis(2-amidinopropane) hydrochloride (AAPH). We did not find changes of compressibility of BSA in the presence of copper ions at rather high molar ratio Cu2+/BSA = 0.66 mol/mol. This can be explained by binding of the Cu2+ to the binding site of BSA and thus protecting the prooxidant action of the copper. However, AAPH-mediated oxidation of BSA resulted in an increase of its apparent specific compressibility (psik/beta0). These changes could be caused by the fragmentation of the protein.     

Molecular docking using surface complementarity

A method is described to dock a ligand into a binding site in a protein on the basis of the complementarity of the inter-molecular atomic contacts. Docking is performed by maximization of a complementarity function that is dependent on atomic contact surface area and the chemical properties of the contacting atoms. The generality and simplicity of the complementarity function ensure that a wide range of chemical structures can be handled. The ligand and the protein are treated as rigid bodies, but displacement of a small number of residues lining the ligand binding site can be taken into account. The method can assist in the design of improved ligands by indicating what changes in complementarity may occur as a result of the substitution of an atom in the ligand. The capabilities of the method are demonstrated by application to 14 protein–ligand complexes of known crystal structure. © 1996 Wiley Liss, Inc.     

Proteins will Fold Anyway!!

Abstract

In order to investigate the environmental conditions of amino acid residues in protein molecules, four kinds of packing studies (atomic, geometric, hydrophobic and hydration) were formulated and tested on two proteins; bovine pancreatic trypsin inhibitor (BPTI) and bovine pancreatic ribonuclease S (RNase S). The inter-relationship of these packings on the fluctuations of amino acid residues was analysed by comparing the packing results with the dynamical studies, such as the root-mean-square-deviation values of atomic displacements obtained from the trajectories of molecular dynamics simulation, temperature factor information from crystal structures and residue fluctuations in proteins from continuum model. These analyses yield information about the most fluctuating and most stabilizing residue sites. Comparison of the results obtained by these methods indicate a good agreement, specifying an inverse correlation between the residue packing and fluctuations. This kind of study is helpful in identifying the specific residue sites such as nucleation, receptor binding and antigenic determining sites which in a way indirectly correlates with the functional residues in protein molecules.     

Relating protein conformational changes to packing efficiency and disorder

Changes in protein conformation play key roles in facilitating various biochemical processes, ranging from signaling and phosphorylation to transport and catalysis. While various factors that drive these motions such as environmental changes and binding of small molecules are well understood, specific causative effects on the structural features of the protein due to these conformational changes have not been studied on a large scale. Here, we study protein conformational changes in relation to two key structural metrics: packing efficiency and disorder. Packing has been shown to be crucial for protein stability and function by many protein design and engineering studies. We study changes in packing efficiency during conformational changes, thus extending the analysis from a static context to a dynamic perspective and report some interesting observations. First, we study various proteins that adopt alternate conformations and find that tendencies to show motion and change in packing efficiency are correlated: residues that change their packing efficiency show larger motions. Second, our results suggest that residues that show higher changes in packing during motion are located on the changing interfaces which are formed during these conformational changes. These changing interfaces are slightly different from shear or static interfaces that have been analyzed in previous studies. Third, analysis of packing efficiency changes in the context of secondary structure shows that, as expected, residues buried in helices show the least change in packing efficiency, whereas those embedded in bends are most likely to change packing. Finally, by relating protein disorder to motions, we show that marginally disordered residues which are ordered enough to be crystallized but have sequence patterns indicative of disorder show higher dislocation and a higher change in packing than ordered ones and are located mostly on the changing interfaces. Overall, our results demonstrate that between the two conformations, the cores of the proteins remain mostly intact, whereas the interfaces display the most elasticity, both in terms of disorder and change in packing efficiency. By doing a variety of tests, we also show that our observations are robust to the solvation state of the proteins.     

Comparing Residue Clusters from Thermophilic and Mesophilic Enzymes Reveals Adaptive Mechanisms

Understanding how proteins adapt to function at high temperatures is important for deciphering the energetics that dictate protein stability and folding. While multiple principles important for thermostability have been identified, we lack a unified understanding of how internal protein structural and chemical environment determine qualitative or quantitative impact of evolutionary mutations. In this work we compare equivalent clusters of spatially neighboring residues between paired thermophilic and mesophilic homologues to evaluate adaptations under the selective pressure of high temperature. We find the residue clusters in thermophilic enzymes generally display improved atomic packing compared to mesophilic enzymes, in agreement with previous research. Unlike residue clusters from mesophilic enzymes, however, thermophilic residue clusters do not have significant cavities. In addition, anchor residues found in many clusters are highly conserved with respect to atomic packing between both thermophilic and mesophilic enzymes. Thus the improvements in atomic packing observed in thermophilic homologues are not derived from these anchor residues but from neighboring positions, which may serve to expand optimized protein core regions.     

Atomic contacts in protein structures. A detailed analysis of atomic radii, packing, and overlaps

A rigorous quantitative assessment of atomic contacts and packing in native protein structures is presented. The analysis is based on optimized atomic radii derived from a set of high-resolution protein structures and reveals that the distribution of atomic contacts and overlaps is a structural constraint in proteins, irrespective of structural or functional classification and size. Furthermore, a newly developed method for calculating packing properties is introduced and applied to sets of protein structures at different levels of resolution. The results show that limited resolution yields decreasing packing quality, which underscores the relevance of packing considerations for structure prediction, design, dynamics, and docking.     

Structural flexibility in proteins: impact of the crystal environment

MOTIVATION: In the study of the structural flexibility of proteins, crystallographic Debye-Waller factors are the most important experimental information used in the calibration and validation of computational models, such as the very successful elastic network models (ENMs). However, these models are applied to single protein molecules, whereas the experiments are performed on crystals. Moreover, the energy scale in standard ENMs is undefined and must be obtained by fitting to the same data that the ENM is trying to predict, reducing the predictive power of the model. RESULTS: We develop an elastic network model for the whole protein crystal in order to study the influence of crystal packing and lattice vibrations on the thermal fluctuations of the atom positions. We use experimental values for the compressibility of the crystal to establish the energy scale of our model. We predict the elastic constants of the crystal and compare with experimental data. Our main findings are (1) crystal packing modifies the atomic fluctuations considerably and (2) thermal fluctuations are not the dominant contribution to crystallographic Debye-Waller factors. AVAILABILITY: The programs developed for this work are available as supplementary material at Bioinformatics Online.     

The packing density in proteins: standard radii and volumes.

The sizes of atomic groups are a fundamental aspect of protein structure. They are usually expressed in terms of standard sets of radii for atomic groups and of volumes for both these groups and whole residues. Atomic groups, which subsume a heavy-atom and its covalently attached hydrogen atoms into one moiety, are used because the positions of hydrogen atoms in protein structures are generally not known. We have calculated new values for the radii of atomic groups and for the volumes of atomic groups. These values should prove useful in the analysis of protein packing, protein recognition and ligand design. Our radii for atomic groups were derived from intermolecular distance calculations on a large number (approximately 30,000) of crystal structures of small organic compounds that contain the same atomic groups to those found in proteins. Our radii show significant differences to previously reported values. We also use this new radii set to determine the packing efficiency in different regions of the protein interior. This analysis shows that, if the surface water molecules are included in the calculations, the overall packing efficiency throughout the protein interior is high and fairly uniform. However, if the water structure is removed, the packing efficiency in peripheral regions of the protein interior is underestimated, by approximately 3.5 %.     

Preferential orientation of an immunoglobulin in a glycolipid monolayer controlled by the disintegration kinetics of proteo-lipidic vesicles spread at an air-buffer interface

The insertion of immunoglobulin (IgG) in a glycolipid monolayer was achieved by using the ability of new proteo-glycolipid vesicles to disintegrate into a mixed IgG-glycolipid interfacial film after spreading at an air-buffer interface. The interfacial disintegration kinetics was shown to be directly dependent on the initial vesicle surface density and on the buffer ionic strength. The presence of the immunoglobulin in the glycolipid film was displayed by an increase of the lateral compressibility (Cs) during monolayer compression. Cs magnitude modifications, due to the antibody effect on the monolayer packing, decreases as the spread vesicle density increases. At interfacial saturation, the lateral compressibility profile becomes similar to that of a control monolayer without antibody. However, the careful analysis of the mixed monolayer after transfer by Langmuir-Blodgett technique (ATR-FTIR characterisation, enzyme immunoassociation) clearly demonstrated that the antibody was still present in such conditions and was not completely squeezed out from the interface as compressibility changes could have meant. At nonsaturating vesicle surface density, IgG molecules initially lying in the lipid matrix with the Y-shape plane parallel to the interface move to a standing-up position during the compression, leading to lateral compressibility modifications. For a saturating vesicle surface density, the glycolipid molecules force the IgG molecules to directly adopt a more vertical position in the interfacial film and, consequently, no lateral compressibility modification was recorded during the compression.     

Changes in crystallographic structure and thermostability of a Cu,Zn superoxide dismutase mutant resulting from the removal of a buried cysteine

In principle, protein thermostability depends on efficient interior packing of apolar residues and on avoidance of irreversible denaturation in the unfolded state. To study these effects, the single free cysteine in the highly stable enzyme bovine Cu,Zn superoxide dismutase was mutated to alanine (Cys6----Ala), and the recombinant protein was expressed in yeast, purified, characterized for reversible and irreversible denaturation, crystallized isomorphously to the wild-type enzyme, and used to determine the atomic structure. Removal of the chemically reactive thiol significantly decreased the rate of irreversible denaturation (as monitored by thermal inactivation at 70 degrees C), but the observed energetic cost (delta delta G of 0.7-1.3 kcal/mol as determined by differential scanning calorimetry) was much less than predicted from either the change in hydrophobicity or packing due to removal of the interior sulfur atom. X-ray diffraction data were collected to 2.1-A resolution using an area detector, and the atomic model for the mutant enzyme was determined by fitting to electron density difference maps, followed by reciprocal space refinement both with stereochemical restraints using PROLSQ and with molecular dynamics using X-PLOR. The refined 2.1-A resolution crystallographic structure suggests that small concerted and compensating shifts (less than 0.5 A) of the surrounding side chains and of the adjacent N- and C-terminal beta-strands significantly reduced the energetic cost of the interior mutation by improving packing and stereochemistry in the mutant enzyme. Taken together, these results differentiate between the effects of reversible and irreversible processes as they impact the design of thermostable proteins and suggest that relatively subtle concerted shifts can significantly reduce the energetic cost of evolutionary variation in internal residues of proteins with Greek key beta-barrel folds.     

Volumetric and spectroscopic characterizations of glucose-hexokinase association

The binding of D-glucose to hexokinase PII at 25 degrees C and pH 8.7 has been investigated by a combination of ultrasonic velocimetry, high precision densimetry, and fluorescence spectroscopy. The binding of glucose to the enzyme results in significant dehydration of the two interacting molecules, while the intrinsic coefficient of adiabatic compressibility of hexokinase slightly decreases. Glucose-hexokinase association is an entropy-driven process. The favorable change in entropy results from compensation between two large contributions. The binding-induced increase in hydrational entropy slightly prevails over the decrease in the configurational entropy of the enzyme. Taken together, our results emphasize the crucial role of water in modulating the energetics of protein recognition.