The amounts of thermal vibrations and static disorder in protein X-ray crystallographic B-factors

Differential Responses in the Core, Active Site and Peripheral Regions of Cytochrome c Peroxidase to Extreme Pressure and Temperature

In consideration of life in extreme environments, the effects of hydrostatic pressure on proteins at the atomic level have drawn substantial interest. Large deviations of temperature and pressure from ambient conditions can shift the free energy landscape of proteins to reveal otherwise lowly populated structural states and even promote unfolding. We report the crystal structure of the heme-containing peroxidase, cytochrome c peroxidase (CcP) at 1.5 and 3.0 kbar and make comparisons to structures determined at 1.0 bar and cryo-temperatures (100 K). Pressure produces anisotropic changes in CcP, but compressibility plateaus after 1.5 kbar. CcP responds to pressure with volume declines at the periphery of the protein where B-factors are relatively high but maintains nearly intransient core structure, hydrogen bonding interactions and active site channels. Changes in active-site solvation and heme ligation reveal pressure sensitivity to protein-ligand interactions and a potential docking site for the substrate peroxide. Compression at the surface affects neither alternate side-chain conformers nor B-factors. Thus, packing in the core, which resembles a crystalline solid, limits motion and protects the active site, whereas looser packing at the surface preserves side-chain dynamics. These data demonstrate that conformational dynamics and packing densities are not fully correlated in proteins and that encapsulation of cofactors by the polypeptide can provide a precisely structured environment resistant to change across a wide range of physical conditions.

Computational insight into stability-enhanced systems of anthocyanin with protein/peptide

Anthocyanins, which belong to the flavonoid group, are commonly found in the organs of plants native to South and Central America. However, these pigments are unstable under conditions of varying pH, heat, etc., which limits their potential applications. One method for preserving the stability of anthocyanins is through encapsulation using proteins or peptides. Nevertheless, the complex and diverse structure of these molecules, as well as the limitation of experimental technologies, have hindered a comprehensive understanding of the encapsulation processes and the mechanisms by which stability is enhanced. To address these challenges, computational methods, such as molecular docking and molecular dynamics simulation have been used to study the binding affinity and dynamics of interactions between proteins/peptides and anthocyanins. This review summarizes the mechanisms of interaction between these systems, based on computational approaches, and highlights the role of proteins and peptides in the stability enhancement of anthocyanins. It also discusses the current limitations of these methods and suggests possible solutions.

Cooperative mechanics of PR65 scaffold underlies the allosteric regulation of the phosphatase PP2A

PR65, a horseshoe-shaped scaffold composed of 15 HEAT (observed in Huntingtin, elongation factor 3, protein phosphatase 2A, and the yeast kinase TOR1) repeats, forms, together with catalytic and regulatory subunits, the heterotrimeric protein phosphatase PP2A. We examined the role of PR65 in enabling PP2A enzymatic activity with computations at various levels of complexity, including hybrid approaches that combine full-atomic and elastic network models. Our study points to the high flexibility of this scaffold allowing for end-to-end distance fluctuations of 40-50 Å between compact and extended conformations. Notably, the intrinsic dynamics of PR65 facilitates complexation with the catalytic subunit and is retained in the PP2A complex enabling PR65 to engage the two domains of the catalytic subunit and provide the mechanical framework for enzymatic activity, with support from the regulatory subunit. In particular, the intra-repeat coils at the C-terminal arm play an important role in allosterically mediating the collective dynamics of PP2A, pointing to target sites for modulating PR65 function.

Simulating diffraction photographs based on molecular dynamics trajectories of a protein crystal: a new option to examine structure-solving strategies in protein crystallography

A molecular dynamics (MD)-based pipeline has been designed and implemented to emulate the entire process of collecting diffraction photographs and calculating crystallographic structures of proteins from them. Using a structure of lysozyme solved in-house, a supercell comprising 125 (5 × 5 × 5) crystal unit cells containing a total of 1000 protein molecules and explicit interstitial solvent was constructed. For this system, two 300 ns MD trajectories at 298 and 250 K were recorded. A series of snapshots from these trajectories were then used to simulate a fully realistic set of diffraction photographs, which were further fed into the standard pipeline for structure determination. The resulting structures show very good agreement with the underlying MD model not only in terms of coordinates but also in terms of B factors; they are also consistent with the original experimental structure. The developed methodology should find a range of applications, such as optimizing refinement protocols to solve crystal structures and extracting dynamics information from diffraction data or diffuse scattering.

Protein Fluctuations in Response to Random External Forces

Elastic network models (ENMs) have been widely used in the last decades to investigate protein motions and dynamics. There the intrinsic fluctuations based on the isolated structures are obtained from the normal modes of these elastic networks, and they generally show good agreement with the B-factors extracted from X-ray crystallographic experiments, which are commonly considered to be indicators of protein flexibility. In this paper, we propose a new approach to analyze protein fluctuations and flexibility, which has a more appropriate physical basis. It is based on the application of random forces to the protein ENM to simulate the effects of collisions of solvent on a protein structure. For this purpose, we consider both the Cα-atom coarse-grained anisotropic network model (ANM) and an elastic network augmented with points included for the crystallized waters. We apply random forces to these protein networks everywhere, as well as only on the protein surface alone. Despite the randomness of the directions of the applied perturbations, the computed average displacements of the protein network show a remarkably good agreement with the experimental B-factors. In particular, for our set of 919 protein structures, we find that the highest correlation with the B-factors is obtained when applying forces to the external surface of the water-augmented ANM (an overall gain of 3% in the Pearson’s coefficient for the entire dataset, with improvements up to 30% for individual proteins), rather than when evaluating the fluctuations obtained from the normal modes of a standard Cα-atom coarse-grained ANM. It follows that protein fluctuations should be considered not just as the intrinsic fluctuations of the internal dynamics, but also equally well as responses to external solvent forces, or as a combination of both.

Low-Frequency Harmonic Perturbations Drive Protein Conformational Changes

Protein dynamics has been investigated since almost half a century, as it is believed to constitute the fundamental connection between structure and function. Elastic network models (ENMs) have been widely used to predict protein dynamics, flexibility and the biological mechanism, from which remarkable results have been found regarding the prediction of protein conformational changes. Starting from the knowledge of the reference structure only, these conformational changes have been usually predicted either by looking at the individual mode shapes of vibrations (i.e., by considering the free vibrations of the ENM) or by applying static perturbations to the protein network (i.e., by considering a linear response theory). In this paper, we put together the two previous approaches and evaluate the complete protein response under the application of dynamic perturbations. Harmonic forces with random directions are applied to the protein ENM, which are meant to simulate the single frequency-dependent components of the collisions of the surrounding particles, and the protein response is computed by solving the dynamic equations in the underdamped regime, where mass, viscous damping and elastic stiffness contributions are explicitly taken into account. The obtained motion is investigated both in the coordinate space and in the sub-space of principal components (PCs). The results show that the application of perturbations in the low-frequency range is able to drive the protein conformational change, leading to remarkably high values of direction similarity. Eventually, this suggests that protein conformational change might be triggered by external collisions and favored by the inherent low-frequency dynamics of the protein structure.