Biomolecular-Solvent Stereodynamic Coupling Probed by Deuteration

Abstract

Thermodynamic interpretation of experiments with isotopically perturbed solvent supports the view that solvent stereodynamics is directly relevant to thermodynamic stability of biomolecules. According with the current understanding of the structure of the aqueous solvent, in any stereodynamic configuration of the latter, connectivity pathways are identifiable for their topologie and order properties. Perturbing the solvent by isotopie substitution or, e.g., by addition of co-solvents, can therefore be viewed as reinforcing or otherwise perturbing these topologie structures.

This microscopic model readily visualizes thermodynamic interpretation. In conclusion, the topologie stereodynamic structures of connectivity pathways in the solvent, as modified by interaction with solutes, acquire a specific thermodynamic and biological significance, and the problem of thermodynamic and functional stability of biomolecules is seen in its full pertinent phase space.     

Applications of pressure perturbation calorimetry in biophysical studies

Pressure perturbation calorimetry (PPC) is a relatively new and efficient technique, to study the volumetric properties of biomolecules in solution. In PPC, the coefficient of thermal expansion of the partial volume of the biomolecule is deduced from the heat consumed or produced after small isothermal pressure jumps (typically ± 5 bar). This strongly depends on the interaction of the biomolecule with the solvent or cosolvent as well as on its packing and internal dynamic properties. This technique, complemented by ultrasound velocity and densitometry, provides valuable insight into the basic thermodynamic properties of solvation and volume effects accompanying phase transitions and interactions of biomolecular systems. Here we review data on protein folding, ligand binding processes, and phospholipid phase transitions, together with discussion of interpretation and further significant applications.     

Mapping the structural core of human cerebral cortex

Structurally segregated and functionally specialized regions of the human cerebral cortex are interconnected by a dense network of cortico-cortical axonal pathways. By using diffusion spectrum imaging, we noninvasively mapped these pathways within and across cortical hemispheres in individual human participants. An analysis of the resulting large-scale structural brain networks reveals a structural core within posterior medial and parietal cerebral cortex, as well as several distinct temporal and frontal modules. Brain regions within the structural core share high degree, strength, and betweenness centrality, and they constitute connector hubs that link all major structural modules. The structural core contains brain regions that form the posterior components of the human default network. Looking both within and outside of core regions, we observed a substantial correspondence between structural connectivity and resting-state functional connectivity measured in the same participants. The spatial and topological centrality of the core within cortex suggests an important role in functional integration.     

Multiscale simulations of protein folding: application to formation of secondary structures

A multiscale simulation method of protein folding is proposed, using atomic representation of protein and solvent, combing genetic algorithms to determine the key protein structures from a global view, with molecular dynamic simulations to reveal the local folding pathways, thus providing an integrated landscape of protein folding. The method is found to be superior to previously investigated global search algorithms or dynamic simulations alone. For secondary structure formation of a selected peptide, RN24, the structures and dynamics produced by this method agree well with corresponding experimental results. Three most populated conformations are observed, including hairpin, β-sheet and α-helix. The energetic barriers separating these three structures are comparable to the kinetic energy of the atoms of the peptide, implying that the transition between these states can be easily triggered by kinetic perturbations, mainly through electrostatic interactions between charged atoms. Transitions between α-helix and β-sheet should jump over at least two energy barriers and may stay in the energetic trap of hairpin. It is proposed that the structure of proteins should be jointly governed by thermodynamic and dynamic factors; free energy is not the exclusive dominant for stability of proteins.     

Biomolecule–nanoparticle interactions: Elucidation of the thermodynamics by isothermal titration calorimetry

BackgroundNanomaterials (NMs) are often exposed to a broad range of biomolecules of different abundances. Biomolecule sorption driven by various interfacial forces determines the surface structure and composition of NMs, subsequently governs their functionality and the reactivity of the adsorbed biomolecules. Isothermal titration calorimetry (ITC) is a nondestructive technique that quantifies thermodynamic parameters through in-situ measurement of the heat absorption or release associated with an interaction.Scope of reviewThis review highlights the recent applications of ITC in understanding the thermodynamics of interactions between various nanoparticles (NPs) and biomolecules. Different aspects of a typical ITC experiment that are crucial for obtaining accurate and meaningful data, as well as the strengths, weaknesses, and challenges of ITC applications to NP research were discussed.Major conclusionsITC reveals the driving forces behind biomolecule–NP interactions and the effects of the physicochemical properties of both NPs and biomolecules by quantifying the crucial thermodynamics parameters (e.g., binding stoichiometry, ΔH, ΔS, and ΔG). Complimentary techniques would strengthen the interpretation of ITC results for a more holistic understanding of biomolecule–NP interactions.General significanceThe thermodynamic information revealed by ITC and its complimentary characterizations is important for understanding biomolecule–NP interactions that are fundamental to the biomedical and environmental applications of NMs and their toxicological effects. This article is part of a Special Issue entitled Microcalorimetry in the BioSciences — Principles and Applications, edited by Fadi Bou-Abdallah.     

Fifty years of solvent denaturation

The author has worked in the area of solvent denaturation and stabilization, off and on, for approximately 50 years. This paper is a personal view of the progress which has been made since 1950. The topic is limited to the development of thermodynamic molecular models for the interpretation of the unfolding and stabilization of protein structures. The story starts with the work in Kauzmann's laboratory shortly after World War II and proceeds through models for multisite binding, the linear denaturation curve, the special considerations required for understanding weak solvent exchange and a new model for the balance of solvent contact interactions and excluded volume.     

Thermodynamics of protein folding: a microscopic view

Statistical thermodynamics provides a powerful theoretical framework for analyzing, understanding and predicting the conformational properties of biomolecules. The central quantity is the potential of mean force or effective energy as a function of conformation, which consists of the intramolecular energy and the solvation free energy. The intramolecular energy can be reasonably described by molecular mechanics-type functions. While the solvation free energy is more difficult to model, useful results can be obtained with simple approximations. Such functions have been used to estimate the intramolecular energy contribution to protein stability and obtain insights into the origin of thermodynamic functions of protein folding, such as the heat capacity. With reasonable decompositions of the various energy terms, one can obtain meaningful values for the contribution of one type of interaction or one chemical group to stability. Future developments will allow the thermodynamic characterization of ever more complex biological processes.     

Interpretation of Solution X-Ray Scattering by Explicit-Solvent Molecular Dynamics

Small- and wide-angle x-ray scattering (SWAXS) and molecular dynamics (MD) simulations are complementary approaches that probe conformational transitions of biomolecules in solution, even in a time-resolved manner. However, the structural interpretation of the scattering signals is challenging, while MD simulations frequently suffer from incomplete sampling or from a force-field bias. To combine the advantages of both techniques, we present a method that incorporates solution scattering data as a differentiable energetic restraint into explicit-solvent MD simulations, termed SWAXS-driven MD, with the aim to direct the simulation into conformations satisfying the experimental data. Because the calculations fully rely on explicit solvent, no fitting parameters associated with the solvation layer or excluded solvent are required, and the calculations remain valid at wide angles. The complementarity of SWAXS and MD is illustrated using three biological examples, namely a periplasmic binding protein, aspartate carbamoyltransferase, and a nuclear exportin. The examples suggest that SWAXS-driven MD is capable of refining structures against SWAXS data without foreknowledge of possible reaction paths. In turn, the SWAXS data accelerates conformational transitions in MD simulations and reduces the force-field bias.     

Thermodynamic databases for proteins and protein-nucleic acid interactions

Thermodynamic data regarding proteins and their interactions are important for understanding the mechanisms of protein folding, protein stability, and molecular recognition. Although there are several structural databases available for proteins and their complexes with other molecules, databases for experimental thermodynamic data on protein stability and interactions are rather scarce. Thus, we have developed two electronically accessible thermodynamic databases. ProTherm, Thermodynamic Database for Proteins and Mutants, contains numerical data of several thermodynamic parameters of protein stability, experimental methods and conditions, along with structural, functional, and literature information. ProNIT, Thermodynamic Database for Protein-Nucleic Acid Interactions, contains thermodynamic data for protein-nucleic acid binding, experimental conditions, structural information of proteins, nucleic acids and the complex, and literature information. These data have been incorporated into 3DinSight, an integrated database for structure, function, and properties of biomolecules. A WWW interface allows users to search for data based on various conditions, with different display and sorting options, and to visualize molecular structures and their interactions. These thermodynamic databases, together with structural databases, help researchers gain insight into the relationship among structure, function, and thermodynamics of proteins and their interactions, and will become useful resources for studying proteins in the postgenomic era.     

Alternatively folded states of an immunoglobulin

Well-defined, non-native protein structures of low stability have been increasingly observed as intermediates in protein folding or as equilibrium structures populated under specific solvent conditions. These intermediate structures, frequently referred to as molten globule states, are characterized by the presence of secondary structure, a lack of significant tertiary contacts, increased hydrophobicity and partial specific volume as compared to native structures, and low cooperativity in thermal unfolding. The present study demonstrates that under acidic conditions (pH less than 3) the antibody MAK33 can assume a folded stable conformation. This A-state is characterized by a high degree of secondary structure, increased hydrophobicity, a native-like maximum wavelength of fluorescence emission, and a tendency toward slow aggregation. A prominent feature of this low-pH conformation is the stability against denaturant and thermal unfolding that is manifested in highly cooperative reversible phase transitions indicative of the existence of well-defined tertiary contacts. These thermodynamic results are corroborated by the kinetics of folding from the completely unfolded chain to the alternatively folded state at pH 2. The given data suggest that MAK33 at pH 2 adopts a cooperative structure that differs from the native immunoglobulin fold at pH 7. This alternatively folded state exhibits certain characteristics of the molten globule but differs distinctly from it by its extraordinary structural stability that is characteristic for native protein structures.     

Structure determination of a DNA octamer in solution by NMR spectroscopy. Effect of fast local motions.

NMR structures of biomolecules are primarily based on nuclear Overhauser effects (NOEs) between protons. For the interpretation of NOEs in terms of distances, usually the assumption of a single rotational correlation time corresponding to a rigid molecule approximation is made. Here we investigate the effect of fast internal motions of the interproton vectors in the context of the relaxation matrix approach for structure determination of biomolecules. From molecular dynamics simulations generalized order parameters were calculated for the DNA octamer d(GCGTTCGC).d(CGCAACGC), and these were used in the calculation of NOE intensities. The magnitudes of the order parameters showed some variation for the different types of interproton vectors. The lowest values were observed for the interresidue base H6/H8-H2" proton vectors (S2 = 0.60), while the cytosine H5-H6 interproton vectors were among the most motionally restricted (S2 = 0.92). Inclusion of the motion of the interproton vectors resulted in a much better agreement between theoretically calculated NOE spectra and the experimental spectra measured by 2D NOE spectroscopy. The interproton distances changed only slightly, with a maximum of 10%; nevertheless, the changes were significant and resulted in constraints that were better satisfied. The structure of the DNA octamer was determined by using restrained molecular dynamics simulations with H2O as a solvent, with and without the inclusion of local internal motions. Starting from A- or B-DNA, the structures showed a high local convergence (0.86 A), while the global convergence for the octamer was ca. 2.6 A.     

Regulatory Networks and Complex Interactions between the Insulin and Angiotensin II Signalling Systems: Models and Implications for Hypertension and Diabetes

The cross-talk between insulin and angiotensin II signalling pathways plays a significant role in the co-occurrence of diabetes and hypertension. We developed a mathematical model of the system of interactions among the biomolecules that are involved in the cross-talk between the insulin and angiotensin II signalling pathways. We have identified several feedback structures that regulate the dynamic behavior of the individual signalling pathways and their interactions. Different scenarios are simulated and dominant steady-state, dynamic and stability characteristics are revealed. The proposed mechanistic model describes how angiotensin II inhibits the actions of insulin and impairs the insulin-mediated vasodilation. The model also predicts that poor glycaemic control induced by diabetes contributes to hypertension by activating the renin angiotensin aystem.     

Hydration changes upon DNA folding studied by osmotic stress experiments

The thermal stability of nucleic acid structures is perturbed under the conditions that mimic the intracellular environment, typically rich in inert components and under osmotic stress. We now describe the thermodynamic stability of DNA oligonucleotide structures in the presence of high background concentrations of neutral cosolutes. Small cosolutes destabilize the basepair structures, and the DNA structures consisting of the same nearest-neighbor composition show similar thermodynamic parameters in the presence of various types of cosolutes. The osmotic stress experiments reveal that water binding to flexible loops, unstable mismatches, and an abasic site upon DNA folding are almost negligible, whereas the binding to stable mismatch pairs is significant. The studies using the basepair-mimic nucleosides and the peptide nucleic acid suggest that the sugar-phosphate backbone and the integrity of the basepair conformation make important contributions to the binding of water molecules to the DNA bases and helical grooves. The study of the DNA hydration provides the basis for understanding and predicting nucleic acid structures in nonaqueous solvent systems.     

Stability of symmetric idiotypic networks—a critique of Hoffmann's analysis

Hoffmann (1982) analysed a very simple model of suppressive idiotypic immune networks and showed that idiotypic interactions are stabilizing. He concluded that immune networks provide a counterexample to the general analysis of large dynamic systems (Gardner and Ashby, 1970; May, 1972). The latter is often verbalized as: an increase in size and/or connectivity decreases the system stability. We here analyse this apparent contradiction by extending the Hoffmann model (with a decay term), and comparing it to an ecological model that was used as a paradigm in the general analysis. Our analysis confirms that the neighbourhood stability of such idiotypic networks increases with connectivity and/or size. However, the contradiction is one of interpretation, and is not due to exceptional properties of immune networks. The contradiction is caused by the awkward normalization used in the general analysis.     

Kinetic stability of membrane proteins

Although membrane proteins constitute an important class of biomolecules involved in key cellular processes, study of the thermodynamic and kinetic stability of their structures is far behind that of soluble proteins. It is known that many membrane proteins become unstable when removed by detergent extraction from the lipid environment. In addition, most of them undergo irreversible denaturation, even under mild experimental conditions. This process was found to be associated with partial unfolding of the polypeptide chain exposing hydrophobic regions to water, and it was proposed that the formation of kinetically trapped conformations could be involved. In this review, we will describe some of the efforts toward understanding the irreversible inactivation of membrane proteins. Furthermore, its modulation by phospholipids, ligands, and temperature will be herein discussed.     

Solvophobic forces and molecular surface area changes in drug-biomolecule associations as with actinomycin-deoxyguanosine in a wide range of methanol/water mixtures

A method is given to predict the unitary free energies of complexation between drug-like and nucleoside-like biomolecules in a range of mixed solvent compositions. A stability maximum for the actinomycin (A)-deoxyguanosine (D) complex at 8% MeOH (v/v) in methanol/water mixtures is correctly predicted by the theory in agreement with existing experimental data. The molecular surface areas of A and D exposed to the solvent are found to diminish by 36.4 A(2) upon association. The 'microthermodynamic differential surface tension' of the solvophobic theory obtained for nucleoside-like and organic molecules in contact with MeOH/H2O can be used to predict the solvent effect free energies in other such molecular or biopolymeric associations in solution.     

Kinetically controlled nanostructure formation in self-assembled globular protein-polymer diblock copolymers

Aqueous processing of globular protein-polymer diblock copolymers into solid-state materials and subsequent solvent annealing enables kinetic and thermodynamic control of nanostructure formation to produce block copolymer morphologies that maintain a high degree of protein fold and function. When model diblock copolymers composed of mCherry-b-poly(N-isopropylacrylamide) are used, orthogonal control over solubility of the protein block through changes in pH and the polymer block through changes in temperature is demonstrated during casting and solvent annealing. Hexagonal cylinders, perforated lamellae, lamellae, or hexagonal and disordered micellar phases are observed, depending on the coil fraction of the block copolymer and the kinetic pathway used for self-assembly. Good solvents for the polymer block produce ordered structures reminiscent of coil-coil diblock copolymers, while an unfavorable solvent results in kinetically trapped micellar structures. Decreasing solvent quality for the protein improves long-range ordering, suggesting that the strength of protein interactions influences nanostructure formation. Subsequent solvent annealing results in evolution of the nanostructures, with the best ordering and the highest protein function observed when annealing in a good solvent for both blocks. While protein secondary structure was found to be almost entirely preserved for all processing pathways, UV-vis spectroscopy of solid-state films indicates that using a good solvent for the protein block enables up to 70% of the protein to be retained in its functional form.     

Order through disorder: hyper-mobile C-terminal residues stabilize the folded state of a helical peptide. a molecular dynamics study

Conventional wisdom has it that the presence of disordered regions in the three-dimensional structures of polypeptides not only does not contribute significantly to the thermodynamic stability of their folded state, but, on the contrary, that the presence of disorder leads to a decrease of the corresponding proteins' stability. We have performed extensive 3.4 μs long folding simulations (in explicit solvent and with full electrostatics) of an undecamer peptide of experimentally known helical structure, both with and without its disordered (four residue long) C-terminal tail. Our simulations clearly indicate that the presence of the apparently disordered (in structural terms) C-terminal tail, increases the thermodynamic stability of the peptide's folded (helical) state. These results show that at least for the case of relatively short peptides, the interplay between thermodynamic stability and the apparent structural stability can be rather subtle, with even disordered regions contributing significantly to the stability of the folded state. Our results have clear implications for the understanding of peptide energetics and the design of foldable peptides.