Exploring the structure and dynamics of proteins in soil organic matter.

Alongside inorganic materials, water, and air, soil organic matter (SOM) is one of the major components of soil and has tremendous influence on the environment given its vital role in the carbon cycle. Many soil dwelling organisms like plants, fungi and bacteria excrete proteins, whose interaction with SOM is poorly understood on an atomistic level. In this study, molecular dynamics simulations were used to investigate selected proteins in soil models of different complexity from simple co-solvent molecules to Leonardite humic acids (LHA). We analyzed the proteins in terms of their structural stability, the nature and strength of the interactions with their surroundings, as well as their aggregation behavior. Upon insertion of proteins in complex SOM models, their structural stability decreased, although no unfolding or disruption of secondary structure was observed. The interactions of proteins and SOM were primarily governed by electrostatic forces, often in form of hydrogen bonds. However, also weaker van der Waals forces made a significant contribution to the total interaction energies. Moreover, we showed that even though the molecular structure and size of SOM molecules varied, the functional groups of SOM ordered around the protein in a similar pattern. Finally, the number of aggregates formed by proteins and SOM molecules was shown to be primarily proportional to the size of the latter. Strikingly, for varying protein net charges no changes in the formation of aggregates with the strongly negatively charged LHA were observed.

Exploring the Cold-Adaptation Mechanism of Serine Hydroxymethyltransferase by Comparative Molecular Dynamics Simulations.

Cold-adapted enzymes feature a lower thermostability and higher catalytic activity compared to their warm-active homologues, which are considered as a consequence of increased flexibility of their molecular structures. The complexity of the (thermo)stability-flexibility-activity relationship makes it difficult to define the strategies and formulate a general theory for enzyme cold adaptation. Here, the psychrophilic serine hydroxymethyltransferase (pSHMT) from Psychromonas ingrahamii and its mesophilic counterpart, mSHMT from Escherichia coli, were subjected to µs-scale multiple-replica molecular dynamics (MD) simulations to explore the cold-adaptation mechanism of the dimeric SHMT. The comparative analyses of MD trajectories reveal that pSHMT exhibits larger structural fluctuations and inter-monomer positional movements, a higher global flexibility, and considerably enhanced local flexibility involving the surface loops and active sites. The largest-amplitude motion mode of pSHMT describes the trends of inter-monomer dissociation and enlargement of the active-site cavity, whereas that of mSHMT characterizes the opposite trends. Based on the comparison of the calculated structural parameters and constructed free energy landscapes (FELs) between the two enzymes, we discuss in-depth the physicochemical principles underlying the stability-flexibility-activity relationships and conclude that (i) pSHMT adopts the global-flexibility mechanism to adapt to the cold environment and, (ii) optimizing the protein-solvent interactions and loosening the inter-monomer association are the main strategies for pSHMT to enhance its flexibility.

Determinants of human glucokinase activation and implications for small molecule allosteric control.

Glucokinase (GK) is an enzyme that catalyzes the ATP-dependent phosphorylation of glucose to form glucose-6-phosphate, and it is a tightly regulated checkpoint in glucose homeostasis. GK is known to undergo substantial conformational changes upon glucose binding. The monomeric enzyme possesses a highly exotic kinetic activity profile with an unusual sigmoidal dependence on glucose concentration. In this interdisciplinary study, which draws on small angle X-ray scattering (SAXS) integrated with 250 ns of atomistic molecular dynamics (MD) simulations and experimental glucose binding thermodynamics, we reveal that the critical regulation of this glucose sensor is due to a solvent controlled "switch". We demonstrate that the "solvent switch" is driven by specific protein structural dynamics, which leads to an enzyme structure that has a much more favorable solvation energy than most of the protein ensemble. These findings uncover the physical workings of an agile and flexible protein scaffold, which derives its long-range allosteric control through specific regions with favorable solvation energy. The physiological framework presented herein provides insights that have direct implications for the design of small molecule GK activators as anti-diabetes therapeutics as well as for understanding how proteins can be designed to have built-in regulatory functions via solvation energy dynamics.

A Reappraisal of Sedimentation Nonideality Coefficients for the Analysis of Weak Interactions of Therapeutic Proteins.

The study of weak or colloidal interactions of therapeutic proteins in different formulations allows prediction and optimization of protein stability. Various biophysical techniques have been applied to determine the second osmotic virial coefficient B2 as it reflects on the macromolecular distance distribution that governs solution behavior at high concentration. In the present work, we exploit a direct link predicted by hydrodynamic theory between B2 and the nonideality of sedimentation, commonly measured in sedimentation velocity analytical ultracentrifugation through the nonideality coefficient of sedimentation, kS. Using sedimentation equilibrium analytical ultracentrifugation for independent measurement of B2, we have examined the dependence of kS on B2 for model proteins in different buffers. The data exhibit the expected linear relationship and highlight the impact of protein shape on the magnitude of the nonideality coefficient kS. Recently, measurements of kS have been considerably simplified allowing higher throughput and simultaneous polydispersity assessment at higher protein concentrations. Thus, sedimentation velocity may offer a useful approach to compare the impact of formulation conditions on weak interactions and simultaneously on higher-order structure of therapeutic proteins.

Effects of protein-crystal hydration and temperature on side-chain conformational heterogeneity in monoclinic lysozyme crystals.

The modulation of main-chain and side-chain conformational heterogeneity and solvent structure in monoclinic lysozyme crystals by dehydration (related to water activity) and temperature is examined. Decreasing the relative humidity (from 99 to 11%) and decreasing the temperature both lead to contraction of the unit cell, to an increased area of crystal contacts and to remodeling of primarily contact and solvent-exposed residues. Both lead to the depopulation of some minor side-chain conformers and to the generation of new conformations. Side-chain modifications and main-chain r.m.s.d.s associated with cooling from 298 to 100 K depend on relative humidity and are minimized at 85% relative humidity (r.h.). Dehydration from 99 to 93% r.h. and cooling from 298 to 100 K result in a comparable number of remodeled residues, with dehydration-induced remodeling somewhat more likely to arise from contact interactions. When scaled to equivalent temperatures based on unit-cell contraction, the evolution of side-chain order parameters with dehydration shows generally similar features to those observed on cooling to T = 100 K. These results illuminate the qualitative and quantitative similarities between structural perturbations induced by modest dehydration, which routinely occurs in samples prepared for 298 and 100 K data collection, and cryocooling. Differences between these perturbations in terms of energy landscapes and occupancies, and implications for variable-temperature crystallography between 180 and 298 K, are discussed. It is also noted that remodeling of a key lysozyme active-site residue by dehydration, which is associated with a radical decrease in the enzymatic activity of lysozyme powder, arises due to a steric clash with the residue of a symmetry mate.

Comparative thermal unfolding study of psychrophilic and mesophilic subtilisin-like serine proteases by molecular dynamics simulations.

Molecular dynamics (MD) simulations of a subtilisin-like serine protease VPR from the psychrophilic marine bacterium Vibrio sp. PA-44 and its mesophilic homologue, proteinase K (PRK), have been performed for 20 ns at four different temperatures (300, 373, 473, and 573 K). The comparative analyses of MD trajectories reveal that at almost all temperatures, VPR exhibits greater structural fluctuations/deviations, more unstable regular secondary structural elements, and higher global flexibility than PRK. Although these two proteases follow similar unfolding pathways at high temperatures, VPR initiates unfolding at a lower temperature and unfolds faster at the same high temperatures than PRK. These observations collectively indicate that VPR is less stable and more heat-labile than PRK. Analyses of the structural/geometrical properties reveal that, when compared to PRK, VPR has larger radius of gyration (Rg), less intramolecular contacts and hydrogen bonds (HBs), more protein-solvent HBs, and smaller burial of nonpolar area and larger exposure of polar area. These suggest that the increased flexibility of VPR would be most likely caused by its reduced intramolecular interactions and more favourable protein-solvent interactions arising from the larger exposure of the polar area, whereas the enhanced stability of PRK could be ascribed to its increased intramolecular interactions arising from the better optimized hydrophobicity. The factors responsible for the significant differences in local flexibility between these two proteases were also analyzed and ascertained. This study provides insights into molecular basis of thermostability of homologous serine proteases adapted to different temperatures.

Enhancing Structure Prediction and Design of Soluble and Membrane Proteins with Explicit Solvent-Protein Interactions.

Solvent molecules interact intimately with proteins and can profoundly regulate their structure and function. However, accurately and efficiently modeling protein solvation effects at the molecular level has been challenging. Here, we present a method that improves the atomic-level modeling of soluble and membrane protein structures and binding by efficiently predicting de novo protein-solvent molecule interactions. The method predicted with unprecedented accuracy buried water molecule positions, solvated protein conformations, and challenging mutational effects on protein binding. When applied to homology modeling, solvent-bound membrane protein structures, pockets, and cavities were recapitulated with near-atomic precision even from distant homologs. Blindly refined atomic-level structures of evolutionary distant G protein-coupled receptors imply strikingly different functional roles of buried solvent between receptor classes. The method should prove useful for refining low-resolution protein structures, accurately modeling drug-binding sites in structurally uncharacterized receptors, and designing solvent-mediated protein catalysis, recognition, ligand binding, and membrane protein signaling.

The role of solvent in protein folding and in aggregation.

We discuss features of the effect of solvent on protein folding andaggregation, highlighting the physics related to the particulate nature and the peculiar structure of the aqueous solvent, and the biological significance of interactions between solvent and proteins. To this purpose we use a generalized energy landscape of extended dimensionality. A closer look at the properties of solvent induced interactions and forces proves useful for understanding the physical grounds of `ad hoc' interactions and for devising realistic ways of accounting for solvent effects. The solvent has long been known to be a crucially important part of biological systems, and times appear mature for it to be adequately accounted for in the protein folding problem. Use of the extended dimensionality energy landscape helpseliciting the possibility of coupling among conformational changes and aggregation, such as proved by experimental data in the literature.