Molecular dynamics simulation study of association in trifluoroethanol/water mixtures

Molecular Aggregation Behavior and Microscopic Heterogeneity in Binary Osmolyte-Water Solutions

Osmolytes, small organic compounds, play a key role in modulating the protein stability in aqueous solutions, but the operating mechanism of the osmolyte remains inconclusive. Here, we attempt to clarify the mode of osmolyte action by quantitatively estimating the microheterogeneity of osmolyte-water mixtures with the aid of molecular dynamics simulation, graph theoretical analysis, and spatial distribution measurement in the four osmolyte solutions of trimethylamine-N-oxide (TMAO), tetramethylurea (TMU), dimethyl sulfoxide, and urea. TMAO, acting as a protecting osmolyte, tends to remain isolated with no formation of osmolyte aggregates while preferentially interacting with water, but there is a strong aggregation propensity in the denaturant TMU solution, characterized by favored hydrophobic interactions between TMU molecules. Taken together, the mechanism of osmolyte action on protein stability is proposed as a comprehensive one that encompasses the direct interactions between osmolytes and proteins and indirect interactions through the regulation of water properties in the osmolyte-water mixtures.

GaMD simulations as an alternative in the TFE-water mixture description

CONTEXT

2,2,2-Trifluoroethanol has been widely used to study the structure and dynamic properties of intrinsically disordered proteins. Experimentally, it is known that TFE-water mixtures stabilize secondary structures of IDPs, and therefore, it allows the studying of conformational ensembles of these proteins. In the last decades, molecular dynamic simulations have helped study the IDPs' conformational ensemble. Unfortunately, conventional MD requires very long simulation times to describe the properties of IDPs. Therefore, a variety of accelerated sampling techniques have been developed and employed. The TFE-water mixture arrangement description through MD has faced substantial difficulties since emulating the TFE nanocrowding at certain TFE:H[Formula: see text]O ratios (around 15-40% of TFE). In this work, we determine the most suitable conditions that reproduce experimentally reported properties of TFE-water mixtures. We compared the employment of conventional MD and GaMD simulations and various water parameters. Our results show that the combination of parameters that better reproduce the experimental information is the combination of the TIP4PD water model and GaMD simulations. Therefore, these conditions help accurately describe the structural ensemble of IDPs in TFE-water mixtures.

METHODS

Conventional MD and GaMD simulations were performed under AMBER 18 software. The TFE and water molecules were described using GAFF2 and a variety of water models, such as TIP3P, TIP4P2005, TIP4PD, and TIP5P, respectively. The systems were simulated a 100 ns at 298 K.

Aggregation in Aqueous Solutions of 2-(Tetrafluoro(trifluoromethyl)-λ6-sulfanyl-ethan-1-ol (CF3SF4-ethanol)): A Comparison with Aqueous Trifluoroethanol and Hexafluoroisopropanol Using Molecular Dynamics Simulations and Dynamic Light Scattering Experiments

2-Tetrafluoro(trifluoromethyl)-λ6-sulfanylethan-1-ol (CF3SF4-ethanol) combines the polar hydrophobicity of tetrafluoro(trifluoromethyl)-λ6-sulfanyl (CF3SF4) group with the polarity of simple alcohols. The properties of aqueous solutions of the well-known fluorinated alcohols 2,2,2-trifluoroethanol (TFE) and 1,1,1,3,3,3-hexafluoroisopropanol (HFIP) were compared with those of aqueous solutions of the novel CF3SF4-ethanol. Those properties were computed using all atom molecular dynamics simulations with OPLS-compatible parameters. DFT ab initio calculations were used to accurately describe the nonsymmetrical, hypervalent sulfur in CF3SF4-ethanol. Although the molecular and conformational characteristics of CF3SF4-ethanol are like those of both TFE and HFIP, the greater hydrophobicity and lower polarity of CF3SF4-ethanol resulted in solution phase aggregation at a much lower concentration. The properties computed for TFE and HFIP in this work were consistent with published computational and experimental studies. CF3SF4-ethanol is predicted to be environmentally benign and hence an excellent green solvent candidate while possessing many of the same properties as TFE or HFIP.

Anionic Dimers of Fluorinated Alcohols

Negative-ion mode electrospray ionization of solutions of ethanol (R^F0OH), 2-fluoroethanol (R^F1OH), 2,2-difluoroethanol (R^F2OH), and/or 2,2,2-trifluoroethanol (R^F3OH) produces anionic dimers of the types (R^FnO)₂H^- and (R^FnO)(R^Fn+1O)H^-. The exchange reactions of these anionic dimers with the neutral alcohols are examined in a quadrupole-ion trap to extract kinetic data, from which the reaction Gibbs energies are obtained. In all cases, the formation of anionic dimers containing the more highly fluorinated alcohols is favored. Quantum chemical calculations confirm this trend and, besides affording structural data, also determine the dissociation energies of the anionic dimers. These dissociation energies are much higher than those of the corresponding neutral dimers and increase further for the more highly fluorinated alcohols due to the stronger hydrogen-bond donor ability of the latter. The present results on the interaction of individual alkoxide anions and neutral alcohol molecules contribute to a better understanding of the association of the fluorinated alcohols in solution.

Computational and Experimental Evaluation of the Stability of a GLP-1-like Peptide in Ethanol–Water Mixtures

Aggregation resulting from the self-association of peptide molecules remains a major challenge during preformulation. Whereas certain organic solvents are known to promote aggregation, ethanol (EtOH) is capable of disrupting interactions between peptide molecules. It is unclear whether it is beneficial or counterproductive to include EtOH in formulations of short peptides. Here, we employed molecular dynamics simulations using the DAFT protocol and MARTINI force field to predict the formation of self-associated dimers and to estimate the stability of a GLP-1-like peptide (G48) in 0–80% aqueous EtOH solutions. Both simulation and experimental data reveal that EtOH leads to a remarkable increase in the conformational stability of the peptide when stored over 15 days at 27 °C. In the absence of EtOH, dimerisation and subsequent loss in conformational stability (α-helix → random coil) were observed. EtOH improved conformational stability by reducing peptide–peptide interactions. The data suggest that a more nuanced approach may be applied in formulation decision making and, if the native state of the peptide is an α-helix organic solvent, such as EtOH, may enhance stability and improve prospects of long-term storage.

Microheterogeneity and CO2 Switchability of N, N-Dimethylcyclohexylamine-Water Binary Mixtures

Binary mixtures of water and organic solvents are described as the aqueous solutions of organic solvents, which are usually spatially heterogeneous on the scale of a few molecular sizes but homogeneous on longer length scales, that is, microheterogeneity. For the water-organic solvent binary mixtures with microheterogeneity, most organic solvents are miscible with water at any ratio. Interestingly, some slightly water-miscible organic solvents can also be used to prepare binary mixtures with microheterogeneity. In this study, N, N-dimethylcyclohexylamine (DMCHA) was used to prepare binary mixtures with microheterogeneity and CO₂ switchability. With the help of conductivity, Fourier-transform infrared spectroscopy, ultraviolet-visible spectroscopy, and dynamic light scattering measurements, we found that water molecules are hydrogen-bonded together to form clusters over the water content range of 9 to 27 wt %, exhibiting microheterogeneity in the binary mixture. The size of the water clusters increases slightly with increasing water content. What is more, the DMCHA-water mixtures can be reversibly split into two phases by alternate bubbling of CO₂ and N₂, exhibiting excellent CO₂ switchability. The binary mixtures can be used as reaction media for the synthesis of CaCO₃ nanoparticles. Binary mixtures with microheterogeneity can also be formed under high salinity or high temperature conditions or be prepared using other slightly water-miscible organic solvents, opening up more interesting possibilities for binary mixtures with microheterogeneity.

Examination of Trifluoroethanol Interactions with Trp-Cage in Trifluoroethanol-Water at 298 K through Molecular Dynamics Simulations and Intermolecular Nuclear Overhauser Effects

Molecular dynamics simulations of the protein model Trp-cage in 42% trifluoroethanol (TFE)-water at 298 K have been carried out with the goal of exploring peptide hydrogen-solvent fluorine nuclear spin cross-relaxation. The TFE5 model of TFE developed in a previous work was used with the TIP5P-Ew model of water. System densities and component translational diffusion coefficients predicted by the simulations were within 20% of the experimental values. Consideration of the calculated relative amounts of TFE and water surrounding the hydrogens of Trp-cage indicated that the composition of the solvent mixture beyond ∼1.5 nm from the van der Waals surface of the peptide is close to the composition of the bulk solvent, but as observed by others, TFE accumulates preferentially near the peptide surface. In the simulations, both TFE and water molecules make contacts with the peptide surface; water molecules predominate in contacts with the peptide backbone atoms and TFE molecules generally preferentially interact with side chains. Translational diffusion of solvent molecules appears to be slowed near the surface of the peptide. Depending on the location in the structure, TFE molecules form complexes with the peptide that may persist for up to ∼7 ns. Many of the peptide spin-solvent fluorine cross-relaxation parameters (Σ_HF) for which experimental values are available are reasonably well-predicted from the simulations. However, the calculated Σ_HF values were too small for some hydrogens of the 6Trp indole ring and the amino acid hydrogens near this residue in the native structure, whereas Σ_HF values for hydrogens on the side chains of 1Asn, 4Ile, and 7Leu are too large. In 42% TFE-water, persistent conformations of Trp-cage are found, which differ from the conformation found in water by the orientation of the 3Tyr ring.

Experimental validation of the role of trifluoroethanol as a nanocrowder

Trifluoroethanol (TFE) is commonly used to induce protein secondary structure, especially α-helix formation. Due to its amphiphilic nature, however, TFE can also self-associate to form micellelike, nanometer-sized clusters. Herein, we hypothesize that such clusters can act as nanocrowders to increase protein folding rates via the excluded volume effect. To test this hypothesis, we measure the conformational relaxation kinetics of an intrinsically disordered protein, the phosphorylated kinase inducible domain (pKID), which forms a helix–turn–helix in TFE solutions. We find that the conformational relaxation rate of pKID displays a rather complex dependence on TFE percentage (v/v): while it first decreases between 0 and 5%, between 5 and 15% the rate increases and then remains relatively unchanged between 15 and 30% and finally decreases again at higher percentages (i.e., 50%). This trend coincides with the fact that TFE clustering is maximized in the range of 15–30%, thus providing validation of our hypothesis. Another line of supporting evidence comes from the observation that the relaxation rate of a monomeric helical peptide, which due to its predominantly local interactions in the folded state is less affected by crowding, does not show a similar TFE dependence.

Structural and functional characterization of a multifunctional alanine-rich peptide analogue from Pleuronectes americanus

Recently, defense peptides that are able to act against several targets have been characterized. The present work focuses on structural and functional evaluation of the peptide analogue Pa-MAP, previously isolated as an antifreeze peptide from Pleuronectes americanus. Pa-MAP showed activities against different targets such as tumoral cells in culture (CACO-2, MCF-7 and HCT-116), bacteria (Escherichia coli ATCC 8739 and Staphylococcus aureus ATCC 25923), viruses (HSV-1 and HSV-2) and fungi (Candida parapsilosis ATCC 22019, Trichophyton mentagrophytes (28d&E) and T. rubrum (327)). This peptide did not show toxicity against mammalian cells such as erythrocytes, Vero and RAW 264.7 cells. Molecular mechanism of action was related to hydrophobic residues, since only the terminal amino group is charged at pH 7 as confirmed by potentiometric titration. In order to shed some light on its structure-function relations, in vitro and in silico assays were carried out using circular dichroism and molecular dynamics. Furthermore, Pa-MAP showed partial unfolding of the peptide changes in a wide pH (3 to 11) and temperature (25 to 95°C) ranges, although it might not reach complete unfolding at 95°C, suggesting a high conformational stability. This peptide also showed a conformational transition with a partial α-helical fold in water and a full α-helical core in SDS and TFE environments. These results were corroborated by spectral data measured at 222 nm and by 50 ns dynamic simulation. In conclusion, data reported here show that Pa-MAP is a potential candidate for drug design against pathogenic microorganisms due to its structural stability and wide activity against a range of targets.

The R.E.D. tools: advances in RESP and ESP charge derivation and force field library building

Deriving atomic charges and building a force field library for a new molecule are key steps when developing a force field required for conducting structural and energy-based analysis using molecular mechanics. Derivation of popular RESP charges for a set of residues is a complex and error prone procedure because it depends on numerous input parameters. To overcome these problems, the R.E.D. Tools (RESP and ESP charge Derive, ) have been developed to perform charge derivation in an automatic and straightforward way. The R.E.D. program handles chemical elements up to bromine in the periodic table. It interfaces different quantum mechanical programs employed for geometry optimization and computing molecular electrostatic potential(s), and performs charge fitting using the RESP program. By defining tight optimization criteria and by controlling the molecular orientation of each optimized geometry, charge values are reproduced at any computer platform with an accuracy of 0.0001 e. The charges can be fitted using multiple conformations, making them suitable for molecular dynamics simulations. R.E.D. allows also for defining charge constraints during multiple molecule charge fitting, which are used to derive charges for molecular fragments. Finally, R.E.D. incorporates charges into a force field library, readily usable in molecular dynamics computer packages. For complex cases, such as a set of homologous molecules belonging to a common family, an entire force field topology database is generated. Currently, the atomic charges and force field libraries have been developed for more than fifty model systems and stored in the RESP ESP charge DDataBase. Selected results related to non-polarizable charge models are presented and discussed.