Protein-solvent Interface Research Articles

Temporal and spatial resolution of distal protein motions that activate hydrogen tunneling in soybean lipoxygenase

The enzyme soybean lipoxygenase (SLO) provides a prototype for deep tunneling mechanisms in hydrogen transfer catalysis. This work combines room temperature X-ray studies with extended hydrogen-deuterium exchange experiments to define a catalytically-linked, radiating cone of aliphatic side chains that connects an active site iron center of SLO to the protein-solvent interface. Employing eight variants of SLO that have been appended with a fluorescent probe at the identified surface loop, nanosecond fluorescence Stokes shifts have been measured. We report a remarkable identity of the energies of activation (Ea) for the Stokes shifts decay rates and the millisecond C-H bond cleavage step that is restricted to side chain mutants within an identified thermal network. These findings implicate a direct coupling of distal protein motions surrounding the exposed fluorescent probe to active site motions controlling catalysis. While the role of dynamics in enzyme function has been predominantly attributed to a distributed protein conformational landscape, the presented data implicate a thermally initiated, cooperative protein reorganization that occurs on a timescale faster than nanosecond and represents the enthalpic barrier to the reaction of SLO.

Molecular-dynamics simulation methods for macromolecular crystallography.

It is investigated whether molecular-dynamics (MD) simulations can be used to enhance macromolecular crystallography (MX) studies. Historically, protein crystal structures have been described using a single set of atomic coordinates. Because conformational variation is important for protein function, researchers now often build models that contain multiple structures. Methods for building such models can fail, however, in regions where the crystallographic density is difficult to interpret, for example at the protein-solvent interface. To address this limitation, a set of MD-MX methods that combine MD simulations of protein crystals with conventional modeling and refinement tools have been developed. In an application to a cyclic adenosine monophosphate-dependent protein kinase at room temperature, the procedure improved the interpretation of ambiguous density, yielding an alternative water model and a revised protein model including multiple conformations. The revised model provides mechanistic insights into the catalytic and regulatory interactions of the enzyme. The same methods may be used in other MX studies to seek mechanistic insights.

Temperature-dependent dynamics at protein-solvent interfaces.

We perform differential scanning calorimetry, broadband dielectric spectroscopy (BDS), and nuclear magnetic resonance (NMR) studies to understand the molecular dynamics in mixtures of ethylene glycol with elastin or lysozyme over broad temperature ranges. To focus on the protein-solvent interface, we use mixtures with about equal numbers of amino acids and solvent molecules. The elastin and lysozyme mixtures show similar glass transition steps, which extend over a broad temperature range of 157-185K. The BDS and NMR studies yield fully consistent results for the fastest process P1, which is caused by the structural relaxation of ethylene glycol between the protein molecules and follows an Arrhenius law with an activation energy of Ea = 0.63eV. It involves quasi-isotropic reorientation and is very similar in the elastin and lysozyme matrices but different from the α and β relaxations of bulk ethylene glycol. Two slower BDS processes, viz., P2 and P3, have protein-dependent time scales, but they exhibit a similar Arrhenius-like temperature dependence with an activation energy of Ea ∼ 0.81eV. However, P2 and P3 do not have a clear NMR signature. In particular, the NMR results for the lysozyme mixture reveal that the protein backbone does not show isotropic α-like motion on the P2 and P3 time scales but only restricted β-like reorientation. The different activation energies of the P1 and P2/P3 processes do not support an intimate coupling of protein and ethylene glycol dynamics. The present results are compared with previous findings for mixtures of proteins with water or glycerol, implying qualitatively different dynamical couplings at various protein-solvent interfaces.

Design, synthesis, and antiviral evaluation of novel piperidine-substituted arylpyrimidines as HIV-1 NNRTIs by exploring the hydrophobic channel of NNIBP

Herein, alkenylpiperidine and alkynylpiperidine moieties were introduced into the left wing of DAPYs (diarylpyrimidines) to explore the new site of the NNIBP (non-nucleoside inhibitor binding pocket) protein-solvent interface region via the structure-based drug design strategy. All the synthesized compounds displayed nanomolar to submicromolar activity against WT (wild-type) HIV-1. Among all, compound FT1 (EC50 = 19 nM) was found to be the most active molecule, which is better than NVP (EC50 = 0.10 μM). In addition, most of the compounds displayed micromolar activity against K103N and E138K mutant strains, while FT1 (EC50(K103N) = 50 nM, EC50(E138K) = 0.19 µM) still has the most effective activity. The molecular dynamics simulation studies revealed that the presence of pyridine moiety of FT1 was essential and played a significant role in its binding with RT (reverse transcriptase).

Similarities and Differences among Protein Dynamics Studied by Variable Temperature Nuclear Magnetic Resonance Relaxation.

Understanding and describing the dynamics of proteins is one of the major challenges in biology. Here, we use multifield variable-temperature NMR longitudinal relaxation (R1) measurements to determine the hierarchical activation energies of motions of four different proteins: two small globular proteins (GB1 and the SH3 domain of α-spectrin), an intrinsically disordered protein (the C-terminus of the nucleoprotein of the Sendai virus, Sendai Ntail), and an outer membrane protein (OmpG). The activation energies map the motions occurring in the side chains, in the backbone, and in the hydration shells of the proteins. We were able to identify similarities and differences in the average motions of the proteins. We find that the NMR relaxation properties of the four proteins do share similar features. The data characterizing average backbone motions are found to be very similar, the same for methyl group rotations, and similar activation energies are measured. The main observed difference occurs for the intrinsically disordered Sendai Ntail, where we observe much lower energy of activation for motions of protons associated with the protein-solvent interface as compared to the others. We also observe variability between the proteins regarding side chain 15N relaxation of lysine residues, with a higher activation energy observed in OmpG. This hints at strong interactions with negatively charged lipids in the bilayer and provides a possible mechanistic clue for the "positive-inside" rule for helical membrane proteins. Overall, these observations refine the understanding of the similarities and differences between hierarchical dynamics in proteins.

Identification of Thermal Conduits That Link the Protein-Water Interface to the Active Site Loop and Catalytic Base in Enolase.

We report here on the salient role of protein mobility in accessing conformational landscapes that enable efficient enzyme catalysis. We are focused on yeast enolase, a highly conserved lyase with a TIM barrel domain and catalytic loop, as part of a larger study of the relationship of site selective protein motions to chemical reactivity within superfamilies. Enthalpically hindered variants were developed by replacement of a conserved hydrophobic side chain (Leu 343) with smaller side chains. Leu343 is proximal to the active site base in enolase, and comparative pH rate profiles for the valine and alanine variants indicate a role for side chain hydrophobicity in tuning the pKa of the catalytic base. However, the magnitude of a substrate deuterium isotope effect is almost identical for wild-type (WT) and Leu343Ala, supporting an unchanged rate-determining proton abstraction step. The introduced hydrophobic side chains at position 343 lead to a discontinuous break in both activity and activation energy as a function of side chain volume. Hydrogen-deuterium exchange mass spectrometry (HDX-MS) experiments were performed as a function of time and temperature for WT and Leu343Ala, and provide a spatially resolved map of changes in protein flexibility following mutation. Impacts on protein flexibility are localized to specific networks that arise at the protein-solvent interface and terminate in a loop that has been shown by X-ray crystallography to close over the active site. These interrelated effects are discussed in the context of long-range, solvent-accessible and thermally activated networks that play key roles in tuning the precise distances and interactions among reactants.

Biophysical Characterization of a Disabled Double Mutant of Soybean Lipoxygenase: The "Undoing" of Precise Substrate Positioning Relative to Metal Cofactor and an Identified Dynamical Network.

Soybean lipoxygenase (SLO) has served as a prototype for understanding the molecular origin of enzymatic rate accelerations. The double mutant (DM) L546A/L754A is considered a dramatic outlier, due to the unprecedented size and near temperature-independence of its primary kinetic isotope effect, low catalytic efficiency, and elevated enthalpy of activation. To uncover the physical basis of these features, we herein apply three structural probes: hydrogen-deuterium exchange mass spectrometry, room-temperature X-ray crystallography and EPR spectroscopy on four SLO variants (wild-type (WT) enzyme, DM, and the two parental single mutants, L546A and L754A). DM is found to incorporate features of each parent, with the perturbation at position 546 predominantly influencing thermally activated motions that connect the active site to a protein-solvent interface, while mutation at position 754 disrupts the ligand field and solvation near the cofactor iron. However, the expanded active site in DM leads to more active site water molecules and their associated hydrogen bond network, and the individual features from L546A and L754A alone cannot explain the aggregate kinetic properties for DM. Using recently published QM/MM-derived ground-state SLO-substrate complexes for WT and DM, together with the thorough structural analyses presented herein, we propose that the impairment of DM is the combined result of a repositioning of the reactive carbon of linoleic acid substrate with regard to both the iron cofactor and a catalytically linked dynamic region of protein.

Paying the Price of Desolvation in Solvent-Exposed Protein Pockets: Impact of Distal Solubilizing Groups on Affinity and Binding Thermodynamics in a Series of Thermolysin Inhibitors.

In lead optimization, open, solvent-exposed protein pockets are often disregarded as prospective binding sites. Because of bulk-solvent proximity, researchers are instead enticed to attach charged polar groups at inhibitor scaffolds to improve solubility and pharmacokinetic properties. It is rarely considered that solvent effects from water reorganization in the first hydration shell of protein-ligand complexes can have a significant impact on binding. We investigate the thermodynamic fingerprint of thermolysin inhibitors featuring terminal charged ammonium groups that are gradually pulled from a distal, solvent-exposed position into the flat, bowl-shaped S2' pocket. Even for the most remote attachment, costs for partial desolvation of the polar group next to the protein-solvent interface are difficult to compensate by interactions with the protein or surrounding water molecules. Through direct comparison with hydrophobic analogues, a significant 180-fold affinity loss was recorded, which questions popular strategies to attach polar ligand-solubilizing groups at the exposed terminus of substituents accommodated in flat open pockets.

The heat released during catalytic turnover enhances the diffusion of an enzyme.

Recent studies have shown that the diffusivity of enzymes increases in a substrate-dependent manner during catalysis. Although this observation has been reported and characterized for several different systems, the precise origin of this phenomenon is unknown. Calorimetric methods are often used to determine enthalpies from enzyme-catalysed reactions and can therefore provide important insight into their reaction mechanisms. The ensemble averages involved in traditional bulk calorimetry cannot probe the transient effects that the energy exchanged in a reaction may have on the catalyst. Here we obtain single-molecule fluorescence correlation spectroscopy data and analyse them within the framework of a stochastic theory to demonstrate a mechanistic link between the enhanced diffusion of a single enzyme molecule and the heat released in the reaction. We propose that the heat released during catalysis generates an asymmetric pressure wave that results in a differential stress at the protein-solvent interface that transiently displaces the centre-of-mass of the enzyme (chemoacoustic effect). This novel perspective on how enzymes respond to the energy released during catalysis suggests a possible effect of the heat of reaction on the structural integrity and internal degrees of freedom of the enzyme.

The R-factor gap in macromolecular crystallography: an untapped potential for insights on accurate structures.

In macromolecular crystallography, the agreement between observed and predicted structure factors (Rcryst and Rfree) is seldom better than 20%. This is much larger than the estimate of experimental error (Rmerge). The difference between Rcryst and Rmerge is the R-factor gap. There is no such gap in small-molecule crystallography, for which calculated structure factors are generally considered more accurate than the experimental measurements. Perhaps the true noise level of macromolecular data is higher than expected? Or is the gap caused by inaccurate phases that trap refined models in local minima? By generating simulated diffraction patterns using the program MLFSOM, and including every conceivable source of experimental error, we show that neither is the case. Processing our simulated data yielded values that were indistinguishable from those of real data for all crystallographic statistics except the final Rcryst and Rfree. These values decreased to 3.8% and 5.5% for simulated data, suggesting that the reason for high R-factors in macromolecular crystallography is neither experimental error nor phase bias, but rather an underlying inadequacy in the models used to explain our observations. The present inability to accurately represent the entire macromolecule with both its flexibility and its protein-solvent interface may be improved by synergies between small-angle X-ray scattering, computational chemistry and crystallography. The exciting implication of our finding is that macromolecular data contain substantial hidden and untapped potential to resolve ambiguities in the true nature of the nanoscale, a task that the second century of crystallography promises to fulfill.DatabaseCoordinates and structure factors for the real data have been submitted to the Protein Data Bank under accession 4tws.

Influence of Site-Dependent Pigment–Protein Interactions on Excitation Energy Transfer in Photosynthetic Light Harvesting

A site-dependent spectral density system-bath model of the Fenna-Matthews-Olsen (FMO) pigment-protein complex is developed using results from ground-state molecular mechanics simulations together with a partial charge difference model for how the long-range contributions to the chromophore excitation energies fluctuate with environmental configuration. A discussion of how best to consistently process the chromophore excitation energy fluctuation correlation functions calculated in these classical simulations to obtain reliable site-dependent spectral densities is presented. The calculations reveal that chromophores that are close to the protein-water interface can experience strongly dissipative environmental interactions characterized by reorganization energies that can be as much as 2-3 times those of chromophores that are buried deep in the hydrophobic protein scaffolding. Using a linearized density matrix quantum propagation method, we demonstrate that the inhomogeneous system-bath model obtained from our site-dependent spectral density calculations gives results consistent with experimental dissipation and dephasing rates. Moreover, we show that this model can simultaneously enhance the energy-transfer rate and extend the decoherence time. Finally, we explore the influence of initially exciting different chromophores and mutating local environments on energy transfer through the network. These studies suggest that different pathways, selected by varying initial photoexcitation, can exhibit significantly different relaxation times depending on whether the energy-transfer path involves chromophores at the protein-solvent interface or if all chromophores in the pathway are buried in the protein.

Relationship between Nanomechanical Properties of Proteins and Secondary Structure

There is tremendous interest in understanding the nanomechanical properties of proteins and how secondary structure influences these properties. Recent studies have analyzed protein rigidity using single molecule techniques such as AFM and MD simulations. These have focused on model systems such as single domain fibrils or amyloids. Here, we investigate these properties in native proteins chosen to represent particular protein classes: α (myoglobin and bovine serum albumin (BSA)), α+β (lysozyme), and β (green fluorescent protein (GFP)). We used Brillouin spectroscopy to provide information on the mechanical modulus on an ∼100nm length scale, accounting for contributions from the protein, solvent, protein-protein interface, and protein-solvent interface; as well as Raman and neutron scattering to analyze the low-frequency vibrational spectra, the boson peak which might provide estimates of rigidity of a single protein.The results clearly demonstrate that the α-helix proteins are softer than the proteins containing the β-motifs, both in dry and hydrated states (h∼0.4). GFP with its β-barrel structure is the most rigid protein studied. The Young's modulus of all the proteins decreases with increasing temperature. Dry proteins are rigid at low temperature, with Young's Modulus approaching ∼7GPa (myoglobin) to ∼13GPa (GFP). Hydrated proteins exhibit an increase in rigidity at low temperatures due to the freezing of hydrogen bond fluctuations between the protein surface and hydration water molecules. At higher temperatures, the hydrated proteins soften significantly, showing high conformational flexibility and experiencing larger amplitude dynamics strongly coupled to the solvent. Our results also suggest that this plasticizing effect by water is weaker for larger proteins, whose hydrophobic core is less affected by the presence of solvent molecules. Finally, we speculate about the physical origin of the boson peak in terms of collective vibrations of specific secondary structural units.

Hypothesis: The unfolding power of protein dielectricity.

A hypothesis is proposed on a potential role of protein dielectricity as an unfolding factor in protein-protein interactions. It is suggested that large protein complexes and aggregation seeds can unfold target proteins by virtue of their effect on the dielectric properties of water at the protein-solvent interface. Here, similar to the effect of membrane surfaces, protein surface can cause decrease in the local dielectric constant of solvent and thereby induce structural changes in a target protein approaching this surface. Some potential implementations of this hypothetical mechanism are also discussed.

Interaction of organic solvents with protein structures at protein-solvent interface

The effect of non-denaturing concentrations of three different organic solvents, formamide, acetone and isopropanol, on the structure of haloalkane dehalogenases DhaA, LinB, and DbjA at the protein-solvent interface was studied using molecular dynamics simulations. Analysis of B-factors revealed that the presence of a given organic solvent mainly affects the dynamical behavior of the specificity-determining cap domain, with the exception of DbjA in acetone. Orientation of organic solvent molecules on the protein surface during the simulations was clearly dependent on their interaction with hydrophobic or hydrophilic surface patches, and the simulations suggest that the behavior of studied organic solvents in the vicinity of hyrophobic patches on the surface is similar to the air/water interface. DbjA was the only dimeric enzyme among studied haloalkane dehalogenases and provided an opportunity to explore effects of organic solvents on the quaternary structure. Penetration and trapping of organic solvents in the network of interactions between both monomers depends on the physico-chemical properties of the organic solvents. Consequently, both monomers of this enzyme oscillate differently in different organic solvents. With the exception of LinB in acetone, the structures of studied enzymes were stabilized in water-miscible organic solvents.

Anisotropy of the Coulomb Interaction between Folded Proteins: Consequences for Mesoscopic Aggregation of Lysozyme

Toward quantitative description of protein aggregation, we develop a computationally efficient method to evaluate the potential of mean force between two folded protein molecules that allows for complete sampling of their mutual orientation. Our model is valid at moderate ionic strengths and accounts for the actual charge distribution on the surface of the molecules, the dielectric discontinuity at the protein-solvent interface, and the possibility of protonation or deprotonation of surface residues induced by the electric field due to the other protein molecule. We apply the model to the protein lysozyme, whose solutions exhibit both mesoscopic clusters of protein-rich liquid and liquid-liquid separation; the former requires that protein form complexes with typical lifetimes of approximately milliseconds. We find the electrostatic repulsion is typically lower than the prediction of the Derjaguin-Landau-Verwey-Overbeek theory. The Coulomb interaction in the lowest-energy docking configuration is nonrepulsive, despite the high positive charge on the molecules. Typical docking configurations barely involve protonation or deprotonation of surface residues. The obtained potential of mean force between folded lysozyme molecules is consistent with the location of the liquid-liquid coexistence, but produces dimers that are too short-lived for clusters to exist, suggesting lysozyme undergoes conformational changes during cluster formation.

Effect of the ordered interfacial water layer in protein complex formation: A nonlocal electrostatic approach

Using a nonlocal electrostatic approach that incorporates the short-range structure of the contacting media, we evaluated the electrostatic contribution to the energy of the complex formation of two model proteins. In this study, we have demonstrated that the existence of an ordered interfacial water layer at the protein-solvent interface reduces the charging energy of the proteins in the aqueous solvent, and consequently increases the electrostatic contribution to the protein binding (change in free energy upon the complex formation of two proteins). This is in contrast with the finding of the continuum electrostatic model, which suggests that electrostatic interactions are not strong enough to compensate for the unfavorable desolvation effects.

Discretization of the induced-charge boundary integral equation.

Boundary-element methods (BEMs) for solving integral equations numerically have been used in many fields to compute the induced charges at dielectric boundaries. In this paper, we consider a more accurate implementation of BEM in the context of ions in aqueous solution near proteins, but our results are applicable more generally. The ions that modulate protein function are often within a few angstroms of the protein, which leads to the significant accumulation of polarization charge at the protein-solvent interface. Computing the induced charge accurately and quickly poses a numerical challenge in solving a popular integral equation using BEM. In particular, the accuracy of simulations can depend strongly on seemingly minor details of how the entries of the BEM matrix are calculated. We demonstrate that when the dielectric interface is discretized into flat tiles, the qualocation method of Tausch [IEEE Trans Comput.-Comput.-Aided Des. 20, 1398 (2001)] to compute the BEM matrix elements is always more accurate than the traditional centroid-collocation method. Qualocation is not more expensive to implement than collocation and can save significant computational time by reducing the number of boundary elements needed to discretize the dielectric interfaces.

Thermostabilization Mechanism of Bovine Serum Albumin by Trehalose

Thermal denaturation of bovine serum albumin (BSA) is analyzed from differential scanning calorimetry (DSC) and Raman spectroscopy investigations. DSC curves exhibit a marked dependence on protein concentration. BSA thermal denaturation becomes broader and bimodal, and the temperature of denaturation increases with increasing protein concentration. Raman scattering investigations simultaneously carried out in the low-frequency range (10-350 cm(-1)) and in the amide I band region (1500-1800 cm(-1)) indicate that the denaturation process is described as a biphasic process independent of protein concentration. The dependence of the protein stability upon the protein concentration can be interpreted from the coupling of protein and solvent dynamics. The confrontation of previous results obtained from Raman investigations on lysozyme (LYS) and the present study of BSA brings out significant information on protein dynamics and the coupling of protein and hydration-water dynamics in relation with the solvent accessible surface area. Contrary to LYS, the modification of the dynamics of hydration water by the protein is clearly observed on BSA. The influence of trehalose on the protein dynamics was analyzed. We found that trehalose reduces the dynamic fluctuations of polar side chains at the protein-solvent interface. The mechanism of thermostabilization by trehalose is related to the reduction of the exposure of hydrophobic groups of BSA to the water molecules, and to a strengthening of intermolecular O-H interactions in the hydrogen-bond network of water, leading to the stabilization of the tertiary structure.

Structure and dynamics of the anti-AMCV scFv(F8): Effects of selected mutations on the antigen combining site

The recombinant antibody fragment scFv(F8), which recognizes the coat protein of the plant virus AMCV, is characterized by peculiar high in vitro stability and functional folding even in reducing environments, making it fit for designing stable antibodies with desired properties. Mutagenesis and functional analysis evidenced two residues, at positions 47 and 58 of the V(H) chain, playing a crucial role in the antigen binding recognition. Here, we used a computational procedure to assess the effects of these mutations on the stability, structure and dynamics of the antigen-binding site. Structural models of the wild type scFv(F8) and of its H47 and H58 mutants were built by homology modelling and assessed by multiple 15.5ns of molecular dynamics simulations. Computational results indicate that the 47H substitution strongly affects the CDR-H(2) conformation, destabilizes the V(H)/V(L) interface and confers high conformational flexibility to the antigen-binding site, leading the mutant to functional loss. The mutation at position H58 strenghtens the binding site, bestowing a high antigen specificity on the mutant. The essential dynamics and the analysis of the protein-solvent interface further corroborate the correspondence between the extent of the structurally-determined flexibility of the binding site with the different functional behaviours proved by the wild-type and its mutants. These results may have useful implications for structure-based design of antibody combining site.

A Billion-fold Range in Acidity for the Solvent-Exposed Amides of Pyrococcus furiosus Rubredoxin

The exchange rates of the static solvent-accessible amide hydrogens of Pyrococcus furiosus rubredoxin range from near the diffusion-limited rate to a billion-fold slower for the non-hydrogen-bonded Val 38 (eubacterial numbering). Hydrogen exchange directly monitors the kinetic acidity of the peptide nitrogen. Electrostatic solvation free energies were calculated by Poisson-Boltzmann methods for the individual peptide anions that form during the hydroxide-catalyzed exchange reaction to examine how well the predicted thermodynamic acidities match the experimentally determined kinetic acidities. With the exception of the Ile 12 amide, the differential exchange rate constant for each solvent-exposed amide proton that is not hydrogen bonded to a backbone carbonyl can be predicted within a factor of 6 (10 (0.78)) root-mean-square deviation (rmsd) using the CHARMM22 electrostatic parameter set and an internal dielectric value of 3. Under equivalent conditions, the PARSE parameter set yields a larger rmsd value of 1.28 pH units, while the AMBER parm99 parameter set resulted in a considerably poorer correlation. Either increasing the internal dielectric value to 4 or reducing it to a value of 2 significantly degrades the quality of the prediction. Assigning the excess charge of the peptide anion equally between the peptide nitrogen and the carbonyl oxygen also reduces the correlation to the experimental data. These continuum electrostatic calculations were further analyzed to characterize the specific structural elements that appear to be responsible for the wide range of peptide acidities observed for these solvent-exposed amides. The striking heterogeneity in the potential at sites along the protein-solvent interface should prove germane to the ongoing challenge of quantifying the contribution that electrostatic interactions make to the catalytic acceleration achieved by enzymes.