Kirkwood–Buff-Derived Alcohol Parameters for Aqueous Carbohydrates and Their Application to Preferential Interaction Coefficient Calculations of Proteins

Solvent accessible surface area-assessed molecular basis of osmolyte-induced protein stability

In solvent-modulated protein folding, under certain physiological conditions, an equilibrium exists between the unfolded and folded states of the protein without any need to break or make a covalent bond. In this process, interactions between various protein groups (peptides) and solvent molecules are known to play a major role in determining the directionality of the chemical reaction. However, an understanding of the mechanism of action of the co(solvent) by a generic theoretical underpinning is lacking. In this study, a generic solvation model is developed based on statistical mechanics and the thermodynamic transfer free energy model by considering the microenvironment polarity of the interacting co(solvent)-protein system. According to this model, polarity and the fractional solvent-accessible surface areas contribute to the interaction energies. The present model includes various orientations of participating interactant solvent surfaces of suitable areas. As model systems, besides the backbone we consider naturally occurring amino acid residues solvated in ten different osmolytes, small organic compounds known to modulate protein stability. The present model is able to predict the correct trend of the osmolyte-peptide interactions ranging from stabilizing to destabilizing not only for the backbone but also for side chains. Our model predicts Asn, Gln, Asp, Glu, Arg and Pro to be highly stable in most of the protecting osmolytes while Ala, Val, Ile, Leu, Thr, Met, Lys, Phe, Trp and Tyr are predicted to be moderately stable, and Ser, Cys and Histidine are predicted to be least stable. However, in denaturing solvents, both backbone and side chain models show similar stabilities in urea and guanidine. One of the important aspects of this model is that it is essentially parameter-free and consistent with the electrostatics of the interaction partners that make this model suitable for estimating any solute-solvent interaction energies.

Preferential Solvation Phenomena as Solute-Induced Structure-Making/Breaking Processes: Linking Thermodynamic Preferential Interaction Parameters to Fundamental Structure Making/Breaking Functions

In this work, we identify the explicit macroscopic-to-microscopic rigorous links between existing thermodynamic preferential interaction parameters Γ Q α Q β ( χ i ) and microstructural descriptors based on total correlation function integrals, leading to their unambiguous characterization in terms of fundamental structure making/breaking functions S α β . First, we provide the statistics mechanical framework to identify a universal molecular-based signature for the preferential solvation P S phenomenon involving solutes at infinite dilution in mixed-solvent environments and discuss its fundamental properties. Then, we characterize the S α β functions relevant to the P S process, identify the microscopic markers for the existing preferential interaction parameters Γ Q α Q β ( χ i ) in terms of the S α β functions, and test their compliance with a pair of essential microstructural constraints linked to the properties of the universal P S signature. Moreover, we illustrate the analysis by probing the behavior of a representative ternary system comprising the solubility of methane in aqueous 1,4-dioxane mixed-solvent environments under ambient conditions. Finally, we discuss some relevant issues surrounding the statistical mechanical (microstructural) interpretation of the thermodynamic (macroscopic) preferential interaction parameters, review some pitfalls in their evaluation from molecular simulation, and provide an outlook.

Not Always Sticky: Specificity of Protein Stabilization by Sugars Is Conferred by Protein-Water Hydrogen Bonds

Solutes added to buffered solutions directly impact protein folding. Protein stabilization by cosolutes or crowders has been shown to be largely driven by protein-cosolute volume exclusion complemented by chemical and soft interactions. By contrast to previous studies that indicate the invariably destabilizing role of soft protein-sugar attractions, we show here that soft interactions with sugar cosolutes are protein-specific and can be stabilizing or destabilizing. We experimentally follow the folding of two model miniproteins that are only marginally stable but in the presence of sugars and polyols fold into representative and distinct secondary structures: β-hairpin or α-helix. Our mean-field model reveals that while protein-sugar excluded volume interactions have a similar stabilizing effect on both proteins, the soft interactions add a destabilizing contribution to one miniprotein but further stabilize the other. Using molecular dynamics simulations, we link the soft protein-cosolute interactions to the weakening of direct protein-water hydrogen bonding due to the presence of sugars. Although these weakened hydrogen bonds destabilize both the native and denatured states of the two proteins, the resulting contribution to the folding free energy can be positive or negative depending on the amino acid sequence. This study indicates that the significant variation between proteins in their soft interactions with sugar determines the specific response of different proteins, even to the same sugar.

Combining Molecular Dynamics Simulations and Biophysical Characterization to Investigate Protein-Specific Excipient Effects on Reteplase during Freeze Drying

We performed molecular dynamics simulations of Reteplase in the presence of different excipients to study the stabilizing mechanisms and to identify the role of excipients during freeze drying. To simulate the freeze-drying process, we divided the process into five distinct steps: (i) protein-excipient formulations at room temperature, (ii) the ice-growth process, (iii)-(iv) the partially solvated and fully dried formulations, and (v) the reconstitution. Furthermore, coarse-grained (CG) simulations were employed to explore the protein-aggregation process in the presence of arginine. By using a coarse-grained representation, we could observe the collective behavior and interactions between protein molecules during the aggregation process. The CG simulations revealed that the presence of arginine prevented intermolecular interactions of the catalytic domain of Reteplase, thus reducing the aggregation propensity. This suggests that arginine played a stabilizing role by interacting with protein-specific regions. From the freeze-drying simulations, we could identify several protein-specific events: (i) collapse of the domain structure, (ii) recovery of the drying-induced damages during reconstitution, and (iii) stabilization of the local aggregation-prone region via direct interactions with excipients. Complementary to the simulations, we employed nanoDSF, size-exclusion chromatography, and CD spectroscopy to investigate the effect of the freeze-drying process on the protein structure and stability.

Multifaceted Computational Modeling in Glycoscience

Glycoscience assembles all the scientific disciplines involved in studying various molecules and macromolecules containing carbohydrates and complex glycans. Such an ensemble involves one of the most extensive sets of molecules in quantity and occurrence since they occur in all microorganisms and higher organisms. Once the compositions and sequences of these molecules are established, the determination of their three-dimensional structural and dynamical features is a step toward understanding the molecular basis underlying their properties and functions. The range of the relevant computational methods capable of addressing such issues is anchored by the specificity of stereoelectronic effects from quantum chemistry to mesoscale modeling throughout molecular dynamics and mechanics and coarse-grained and docking calculations. The Review leads the reader through the detailed presentations of the applications of computational modeling. The illustrations cover carbohydrate-carbohydrate interactions, glycolipids, and N- and O-linked glycans, emphasizing their role in SARS-CoV-2. The presentation continues with the structure of polysaccharides in solution and solid-state and lipopolysaccharides in membranes. The full range of protein-carbohydrate interactions is presented, as exemplified by carbohydrate-active enzymes, transporters, lectins, antibodies, and glycosaminoglycan binding proteins. A final section features a list of 150 tools and databases to help address the many issues of structural glycobioinformatics.

A Transfer Free Energy Based Implicit Solvent Model for Protein Simulations in Solvent Mixtures: Urea-Induced Denaturation as a Case Study

We developed a method for implicit solvent molecular dynamics simulations of proteins in solvent mixtures (model with implicit solvation thermodynamics, MIST). The MIST method introduces experimental group transfer free energies to the generalized Born formulation for generating molecular trajectories without the need for developing rigorous explicit-solvent force fields for multicomponent solutions. As a test case, we studied the urea-induced denaturation of the Trp-cage miniprotein in water. We demonstrate that our method allows efficient exploration of the conformational space of the protein in only a few hundreds of nanoseconds of all-atom unbiased simulations. Furthermore, selective implementation of the transfer free energies of specific peptide groups, backbone, and side chains enables us to decouple their specific energetic contributions to the conformational changes of the protein. The approach herein developed can readily be extended to the investigation of complex matrices as well as to the characterization of protein aggregation. The MIST method is implemented in Plumed (ver. 2.8) as a separate module called SASA.

The Study of Protein-Cyclitol Interactions

Investigation of interactions between the target protein molecule and ligand allows for an understanding of the nature of the molecular recognition, functions, and biological activity of protein-ligand complexation. In the present work, non-specific interactions between a model protein (Bovine Serum Albumin) and four cyclitols were investigated. D-sorbitol and adonitol represent the group of linear-structure cyclitols, while shikimic acid and D-(-)-quinic acid have cyclic-structure molecules. Various analytical methods, including chromatographic analysis (HPLC-MS/MS), electrophoretic analysis (SDS-PAGE), spectroscopic analysis (spectrofluorimetry, Fourier transform infrared spectroscopy, and Raman spectroscopy), and isothermal titration calorimetry (ITC), were applied for the description of protein-cyclitol interactions. Additionally, computational calculations were performed to predict the possible binding places. Kinetic studies allowed us to clarify interaction mechanisms that may take place during BSA and cyclitol interaction. The results allow us, among other things, to evaluate the impact of the cyclitol's structure on the character of its interactions with the protein.

Force Field Parameterization for the Description of the Interactions between Hydroxypropyl-β-Cyclodextrin and Proteins

Cyclodextrins are cyclic oligosaccharides, widely used as drug carriers, solubilizers, and excipients. Among cyclodextrins, the functionalized derivative known as hydroxypropyl-β-cyclodextrin (HPβCD) offers several advantages due to its unique structural features. Its optimal use in pharmaceutical and medical applications would benefit from a molecular-level understanding of its behavior, as can be offered by molecular dynamics simulations. Here, we propose a set of parameters for all-atom simulations of HPβCD, based on the ADD force field for sugars developed in our group, and compare it to the original CHARMM36 description. Using Kirkwood-Buff integrals of binary HPβCD-water mixtures as target experimental data, we show that the ADD-based description results in a considerably improved prediction of HPβCD self-association and interaction with water. We then use the new set of parameters to characterize the behavior of HPβCD toward the different amino acids. We observe pronounced interactions of HPβCD with both polar and nonpolar moieties, with a special preference for the aromatic rings of tyrosine, phenylalanine, and tryptophan. Interestingly, our simulations further highlight a preferential orientation of HPβCD's hydrophobic cavity toward the backbone atoms of amino acids, which, coupled with a favorable interaction of HPβCD with the peptide backbone, suggest a propensity for HPβCD to denature proteins.

Kirkwood-Buff-Derived Force Field for Peptides and Proteins: Philosophy and Development of KBFF20

A new classical nonpolarizable force field, KBFF20, for the simulation of peptides and proteins is presented. The force field relies heavily on the use of Kirkwood-Buff theory to provide a comparison of simulated and experimental Kirkwood-Buff integrals for solutes containing the functional groups common in proteins, thus ensuring intermolecular interactions that provide a good balance between the peptide-peptide, peptide-solvent, and solvent-solvent distributions observed in solution mixtures. In this way, it differs significantly from other biomolecular force fields. Further development and testing of the intermolecular potentials are presented here. Subsequently, rotational potentials for the ϕ/ψ and χ dihedral degrees of freedom are obtained by analysis of the Protein Data Bank, followed by small modifications to provide a reasonable balance between simulated and observed α and β percentages for small peptides. This, the first of two articles, describes in detail the philosophy and development behind KBFF20.

Three-Dimensional Structures of Carbohydrates and Where to Find Them

Analysis and systematization of accumulated data on carbohydrate structural diversity is a subject of great interest for structural glycobiology. Despite being a challenging task, development of computational methods for efficient treatment and management of spatial (3D) structural features of carbohydrates breaks new ground in modern glycoscience. This review is dedicated to approaches of chemo- and glyco-informatics towards 3D structural data generation, deposition and processing in regard to carbohydrates and their derivatives. Databases, molecular modeling and experimental data validation services, and structure visualization facilities developed for last five years are reviewed.

Cyclodextrins in Parenteral Formulations

Most drugs have very limited solubility in water and some can be extremely difficult to formulate as parenteral solutions where the dose should preferably be dissolved in couple of ml of aqueous media without use of organic solvents and surface active agents, or application of somewhat extreme techniques such as prodrug formation. Thus, pharmaceutical formulators are constantly looking for new, biologically acceptable, and low-cost armamentarium for parenteral formulation development. Cyclodextrins (CDs) are enabling pharmaceutical excipients that can temporarily camouflage undesirable physiochemical drug properties such as low aqueous solubility through formation of drug/CD inclusion complexes. CDs are cyclic oligosaccharides that have similar physiological and biological properties like linear saccharides of comparable molecular weight. Due to their very favorable toxicological and pharmacokinetic profiles their usage in parenteral drug formulations is frequently preferred over other solubilizing techniques. Here the physiochemical and biological properties of CDs are reviewed as well as their pharmacokinetics after intravenous administration. Their regulatory status is briefly reviewed and their tendency to self-assemble to form clusters or aggregates discussed. Finally, some examples are given of how CDs are applied in aqueous parenteral formulations, how their solubilizing effect has been enhanced and how their target concentration is determined.

ADD Force Field for Sugars and Polyols: Predicting the Additivity of Protein-Osmolyte Interaction

The protein–osmolyte interaction has been shown experimentally to follow an additive construct, where the individual osmolyte–backbone and osmolyte–side-chain interactions contribute to the overall conformational stability of proteins. Here, we computationally reconstruct this additive relation using molecular dynamics simulations, focusing on sugars and polyols, including sucrose and sorbitol, as model osmolytes. A new set of parameters (ADD) is developed for this purpose, using the individual Kirkwood–Buff integrals for sugar–backbone and sugar–side-chain interactions as target experimental data. We show that the ADD parameters can reproduce the additivity of protein–sugar interactions and correctly predict sucrose and sorbitol self-association and their interaction with water. The accurate description of the separate osmolyte–backbone and osmolyte–side-chain contributions also automatically translates into a good prediction of preferential exclusion from the surface of ribonuclease A and α-chymotrypsinogen A. The description of sugar polarity is improved compared to previous force fields, resulting in closer agreement with the experimental data and better compatibility with charged groups, such as the guanidinium moiety. The ADD parameters are developed in combination with the CHARMM36m force field for proteins, but good compatibility is also observed with the AMBER 99SB-ILDN and the OPLS-AA force fields. Overall, exploiting the additivity of protein–osmolyte interactions is a promising approach for the development of new force fields.

Properties of Aqueous Trehalose Mixtures: Glass Transition and Hydrogen Bonding

Trehalose is a naturally occurring disaccharide known to remarkably stabilize biomacromolecules in the biologically active state. The stabilizing effect is typically observed over a large concentration range and affects many macromolecules including proteins, lipids, and DNA. Of special interest is the transition from aqueous solution to the dense and highly concentrated glassy state of trehalose that has been implicated in bioadaptation of different organisms toward desiccation stress. Although several mechanisms have been suggested to link the structure of the low water content glass with its action as an exceptional stabilizer, studies are ongoing to resolve which are most pertinent. Specifically, the role that hydrogen bonding plays in the formation of the glass is not well resolved. Here we model aqueous trehalose mixtures over a wide concentration range, using molecular dynamics simulations with two available force fields. Both force fields indicate glass transition temperatures and osmotic pressures that are close to experimental values, particularly at high trehalose contents. We develop and employ a methodology that allows us to analyze the thermodynamics of hydrogen bonds in simulations at different water contents and temperatures. Remarkably, this analysis is able to link the liquid to glass transition with changes in hydrogen bond characteristics. Most notably, the onset of the glassy state can be quantitatively related to the transition from weakly to strongly correlated hydrogen bonds. Our findings should help resolve the properties of the glass and the mechanisms of its formation in the presence of added macromolecules.

Molecular Computations of Preferential Interaction Coefficients of IgG1 Monoclonal Antibodies with Sorbitol, Sucrose, and Trehalose and the Impact of These Excipients on Aggregation and Viscosity

Preferential interactions of formulation excipients govern their overall interactions with protein molecules, and molecular dynamics simulations allow for the examination of the interactions at the molecular level. We used molecular dynamics simulations to examine the interactions of sorbitol, sucrose, and trehalose with three different IgG1 antibodies to gain insight into how these excipients impact aggregation and viscosity. We found that sucrose and trehalose reduce aggregation more than sorbitol because of their larger size and their stronger interactions with high-spatial aggregation propensity residues compared to sorbitol. Two of the antibodies had high viscosity in sodium acetate buffer, and for these, we found that sucrose and trehalose tended to have opposite effects on viscosity. The data presented here provide further insight into the mechanisms of interactions of these three carbohydrate excipients with the antibody surface and thus their impact on excipient stabilization of antibody formulations.

Self-association of cyclodextrins and cyclodextrin complexes in aqueous solutions

Cyclodextrins (CDs) are oligosaccharides that self-assemble in aqueous solutions to form transient clusters, nanoparticles and small microparticles. The critical aggregation concentration (cac) of the natural αCD, βCD and γCD in pure aqueous solutions was estimated to be 25, 8 and 9 mg/ml, respectively. The cac of 2-hydroxypropyl-β-cyclodextrin (HPβCD), that consists of mixture of isomers, was estimated to be significantly higher or 118 mg/ml. Addition of chaotropic agents (i.e. that disrupts non-covalent bonds such as hydrogen bonds) to the aqueous media increases the cac. Formation of drug/CD complexes can increase or decrease the cac. Due to the transient nature of the CD clusters and nanoparticles they can be difficult to detect and their presence is frequently ignored. However, they have profound effect on the physiochemical properties of CDs and their pharmaceutical applications. For example, the values of stability constants of drug/CD complexes can be both concentration dependent and method dependent. Like in the case of micelles water-soluble polymers can enhance the solubilizing effect of CDs. Also, formation of drug/CD complex nanoparticles appears to increase the ability of CDs to enhance drug delivery through some mucosal membranes.