How does temperature modulate the structural properties of aggregated melamine in aqueous solution—An answer from classical molecular dynamics simulation

Influence of salt and temperature on the self-assembly of cyclic peptides in water: a molecular dynamics study

It is found in the literature that cyclic peptides (CPs) are able to self-assemble in water to form cyclic peptide nanotubes (CPNTs) and are used extensively in the field of nanotechnology. Several factors influence the formation and stability of these nanotubes in water. However, an extensive study of the contribution of several important factors is still lacking. The purpose of this study is to explore the effect of temperature and salt (NaCl) on the association tendency of CPs. Furthermore, the self-association behavior of CPs in aqueous solutions at various temperatures is also thoroughly discussed. Cyclo-[(Asp-D-Leu-Lys-D-Leu)₂] is considered for this study and a series of classical molecular dynamics (MD) simulations at three different temperatures, viz. 280 K, 300 K, and 320 K, both in pure water and in NaCl solutions of different concentrations are carried out. The calculations of radial distribution functions, preferential interaction parameters, cluster formation and hydrogen bonding properties suggest a strong influence of NaCl concentration on the association propensity of CPs. Low NaCl concentration hinders CP association while high NaCl concentration facilitates the association of CPs. Besides this, the association of CPs is found to be enhanced at low temperature. Furthermore, the thermodynamics of CP association is predominantly found to be enthalpy driven in both the presence and absence of salt. No crossover between enthalpy and entropy in CP association is observed. In addition, the MM-GBSA method is used to investigate the binding free energies of the CP rings that self-assembled to form nanotube like structures at all three temperatures.

The chaotropic effect of ions on the self-aggregating propensity of Whitlock's molecular tweezers

Molecular tweezers feature the first class of artificial receptors to pique the interest of researchers and emerge as an effective therapeutic candidate. The exceptional structure and exquisite binding specificity of tweezers establish this overall class of receptors as a promising tool, with abundant applications. However, their inclination to self-aggregate by mutual π-π stacking interactions of their aromatic arms diminishes their efficacy as a therapeutic candidate. Therefore, following up on sporadic studies, since the discovery of the Hofmeister series, on the ability of ions to either solvate (salting-in) or induce aggregation (salting-out) of hydrophobic solutes, the notions of ion-specificity effects are utilized on tweezer moieties. The impacts of three different aluminum salts bearing anions Cl^-, ClO₄^- and SCN^- on the self-association propensity of Whitlock's caffeine-pincered molecular tweezers are investigated, with a specific emphasis placed on elucidating the varied behavior of the ions on the hydration ability of tweezers. The comparative investigation is conducted employing a series of all-atom molecular dynamics simulations of five tweezer molecules in pure water and three salt solutions, at two different concentrations each, maintaining a temperature of 300 K and a pressure of 1 atm, respectively. Radial distribution functions, coordination numbers, and SASA calculations display a steady reduction in the aggregation proclivity of the receptor molecules with an increase in salt concentration, as progressed along the Hofmeister series. Orientational preferences between the tweezer arms reveal a disruptive effect in the regular π-π stacking interactions, in the presence of high concentrations of ClO₄^- and SCN^- ions, while preferential interactions and tetrahedral order parameters unveil the underlying mechanism, by which the anions alter the solubility of the hydrophobic molecules. Overall, it is observed that SCN^- exhibits the highest salting-in effect, followed by ClO₄^-, with both anions inhibiting tweezer aggregation through different mechanisms. ClO₄^- ions impart an effect by moderately interacting with the solute molecules as well as modifying the water structure of the bulk solution promoting solvation, whereas, SCN^- ions engage entirely in interaction with specific tweezer sites. Cl^- being the most charge-dense of the three anionic species experiences stronger hydration and therefore, imparts a very negligible salting-in effect.

Molecular insights into the urea-choline-O-sulfate interactions in aqueous solution

Urea and choline-O-sulfate (COS) are both osmolytes, but have opposite effects on protein structure. Urea has been well-known for years to destabilize protein structure. Though COS has been revealed as an osmoprotective molecule against urea induced denaturation of proteins, the mechanism of this compensation is still unexplored. This study focuses on a theoretical investigation of the interdependent behavior of urea and COS in a mixture, to explore how urea becomes a weaker denaturing agent in the presence of COS. In this study, we have considered every possible interaction among the solute (urea and COS) and solvent (water) both at room temperature and high temperature, employing two different force field parameters i.e., CHARMM General Force Field parameters (CGenFF) and General AMBER Force Field (GAFF) parameters through classical molecular dynamics simulation studies. Different techniques have been used to analyze the average interactions between COS and urea as well as their solvation properties, which show that in the presence of COS, urea becomes a less effective denaturant than when alone. The water-water interaction shows that the mixed osmolyte solution of urea and COS strengthens the water hydrogen bonding network. The enhanced solvation of urea and COS in the urea-COS mixture and their mutual interactions, results in the exclusion of free urea as well as COS from the solution. This synergistic behavior of urea and COS could be the major reason behind COS counteracting urea's denaturation of proteins.

A computational approach on the stereoselective binding of peptides from aqueous medium with endo-functionalized molecular tubes

The need to obtain enantiomerically pure isomers of amino acids and peptides is often realized in the field of biology and in the pharmaceutical industry. Research is underway to devise simple methods for the chiral resolution of amino acids from their racemic mixtures. Inspired by this objective, in our present work, we have computationally shown the possibility of chiral separation of the enantiomeric pairs of two model peptides, namely, (D,L)-aspargine and (D,L)-phenylalanine, in the presence of water. For this purpose, we have used two synthetic supramolecular receptors named host-1a and host-1b, respectively. Molecular dynamics simulations and quantum chemical methods are employed to analyze the structural features and the energy aspects involved in the separation process. The information obtained at the molecular level helps us gain better insights into the key interactions that operate to produce such enantioselectivity. We have also investigated the dynamics and changes in the water structure in the vicinity of the host molecules, both in the presence and absence of the model peptides. The D- and L-isomers of the same peptide undergo complexation with a particular host molecule registering a difference of more than 1.5 kcal mol^-1 (obtained from PMF and MM-PBSA analyses) in their respective energies. This indicates that the chiral separation of the peptides with the help of these endo-fuctionalized molecular tube receptors may be energetically feasible. The connection between the peptide stereochemistry and its interaction with the endo-functionalized hosts would be instrumental in designing novel segregation techniques that can be further extended to separate larger model peptides or proteins.

Phase separation property of a hydrophobic deep eutectic solvent-water binary mixture: A molecular dynamics simulation study

Over the past decade, deep eutectic solvents (DESs) have earned applicability in numerous fields as non-flammable, non-volatile, and greener alternatives to conventional organic solvents. In a first of its kind, a hydrophobic DES composed of a 1:1 mixture of oleic acid and lidocaine was recently reported, possessing a lower critical solution temperature in water. The thermoreversible phase property of this DES-water system was utilized to sequester out dye molecules from their aqueous solutions. In this article, we explore the phase separation phenomena for this particular DES in its aqueous solution using an all-atom molecular dynamics simulation. A 50 wt. % solution of the DES in water was studied at three different temperatures (253, 293, and 313 K) to understand the various molecular interactions that dictate the phase segregation property of these systems. In this work, we have elaborated on the importance of hydrogen bonding interactions and the non-bonding interactions between the components and the competition between the two that leads to phase separation. Overall, we observe that the increase in unfavorable interaction between the DES components and water with increasing temperature determines the phase separation behavior. We have also studied the modification in the dynamical properties of water molecules close to the phase boundary. Such molecular insights would be beneficial for designing novel solvent systems that can be used as extraction-based media in industries.

Molecular Insight into the Effects of Enhanced Hydrophobicity on Amyloid-like Aggregation

Generally, hydrophobic amino acids provide hydrophobic interactions during peptide aggregation. However, besides the hydrophobic amino acids, some hydrophilic amino acids, such as glutamine, are also considered to be essential elements in many self-aggregating peptides. For example, huntingtin contains polyglutamine at its N-terminus and the yeast prion Sup35 protein has the GNNQQNY sequence, a peptide well-known for its ability for amyloid fibril formation. However, despite the frequent emergence of glutamine in self-assembling systems, the molecular mechanism of amyloid formation involving this unique amino acid has not been well documented. It is still not clear how this hydrophilic amino acid is responsible for the hydrophobic interaction in the self-association process. Therefore, in this study, we have carried out classical molecular dynamics simulations of the GNNQQNY peptide and its derivatives in pure water. We quantify the propensity for the formation of β-sheet conformation with an increasing glutamine number in the peptide sequence. In addition, we assess the importance of the hydrophobicity of the dimethanediyl group present in glutamine (as well as in glutamic acid) for the self-association of the peptides through nonpolar solvent medium simulations.

Investigation on the Mechanisms of Synchronous Interaction of K3Cit with Melamine and Uric Acid That Avoids the Formation of Large Clusters

Uric acid (UA) has an enormous competence to aggregate over melamine (Mel), producing large UA clusters that "drag" Mel toward them. Such a combination of donor-acceptor pairs provides a robust Mel-UA composite, thereby denoting a high complexity. Thus, a straightforward but pragmatic methodology might indeed require either destruction of the aggregation of UA or impediment of the hydrogen-bonded cluster of Mel and UA. Here, potassium citrate (K₃Cit) is used as a potent inhibitor for a significant decrease of large UA-Mel clusters. The underlying mechanisms of synchronous interactions between K₃Cit and the Mel-UA pair are examined by the classical molecular dynamics simulation coupled with the enhanced sampling method. K₃Cit binds to the Mel-UA pair profoundly to produce a Mel-UA-K₃Cit complex with favorable complexation energy (as indicated by the reckoning of pairwise ΔG_bind° employing the molecular mechanics Poisson-Boltzmann surface area (MM-PBSA) method). The strength of interaction follows the order UA-K₃Cit > Mel-K₃Cit > Mel-UA, thus clearly demonstrating the instability caused by upsetting the π-stacking of UA and hydrogen bonding of Mel-UA simultaneously. The comprehensive, strategically designed "direct approach" and "indirect approach" cluster structure analysis shows that K₃Cit reduces the direct approach Mel-UA cluster size significantly irrespective of ensemble variation. Furthermore, the estimation of potentials of mean force (PMFs) reveals that the (UA)_decamer-Mel interaction prevails over (UA)_tetramer-Mel. The dynamic property (dimer existence autocorrelation functions) proves the essence of dimerization between Mel and UA in the absence and presence of K₃Cit. Moreover, the calculation of the preferential interaction parameter provides the concentration at which Mel-K₃Cit and UA-K₃Cit interactions are predominant over the interaction of Mel and UA.

Switching of Self-Assembly to Solvent-Assisted Assembly of Molecular Motor: Unveiling the Mechanisms of Dynamic Control on Solvent Exchange

Natural biological molecular motors are capable of performing several biological functions, such as fuel production, mobility, transport, and many other dynamic features. Inspired by these biological motors, scientists effectively synthesized artificial molecular motors to mimic several biological functionalities. Several molecular systems, from sensitive materials to molecular motors, are essential for controlling dynamic processes in larger assemblies. In this work, we discuss the self-assembly of molecular motors in water and how this self-assembly switches to the solvent-assisted assembly as solvent changes to a water-THF (tetrahydrofuran) mixture. We present an elaborate description of the morphological changes of molecular motor assemblies from pure water to a water-THF mixture to pure THF. Under the influence of THF solvent, molecular motors form an assembled structure by taking a sufficient number of THF molecules in between themselves, resulting in an assembled molecular motor with a softened core. So, molecular motor assembly swells in the water-THF mixture, and in pure water, it shrinks. This solvent-assisted assembled structure has a specific shape. We have confirmed this assembly as a solvent-assisted assembly with the help of molecular dynamics simulation and quantum chemical analysis. Molecular motor-THF and THF-THF interactions are the main responsible interactions for solvent-assisted assembly over self-assembly. This work is a perfect example of conversion between self-assembly (shrinking) and solvent-assisted assembly (swelling) of molecular motors by adding THF into water or vice versa. A spectacular check on the shrinking and swelling by merely altering solvents illustrates so many intriguing possibilities for an alteration of dynamic processes at the nanoscale.

Computational Study of Encapsulation of Polyaromatic Hydrocarbons by Endo-Functionalized Receptors in Nonpolar Medium

Polycyclic aromatic hydrocarbons (PAHs) constitute a large group of organic pollutants produced from either natural or artificial sources during the incomplete combustion of fossil fuels or derived from various industrial processes (such as refinery processes of crude petroleum). They are seriously hazardous to human health, and removing them is of major importance. The complexation likeliness with and selective recognition of PAH guests by endo-functionalized molecular tube hosts (host-abu and host-abtu) in a nonpolar medium are investigated using classical molecular dynamics simulation and quantum calculation to probe the factors and the molecular mechanism involved in complexation processes. We examine the role of different guest molecules in the structural changes of hosts, a prelude to van der Waals interactions and binding free energy in the complexation process. These types of host-guest interactions depend on various factors. We find that (i) both the host molecules (host-abtu and host-abu) interact with the guest π-electron cloud almost equally and (ii) these interactions also depend on the molecular size of PAHs. The larger the nonpolar surface area of PAHs, the greater the interactions with the host, and the more extensive the π-electron cloud of the guest, the stronger the interactions. The linear PAHs interact more strongly than isomeric branched/curved PAHs, and the presence of heteroatoms on PAHs decreases the interactions with the host by creating repulsion between the lone pairs of heteroatoms and the π-electron cloud of the host. Noncovalent van der Waals interactions and N-H···π interactions dominate the high affinities of PAHs toward host-abu and host-abtu. The potential of mean force and molecular mechanics Poisson-Boltzmann surface area calculations reveal that all host-guest complexes are energetically stable.