Collection of Partition Coefficients in Hexadecyltrimethylammonium Bromide, Sodium Cholate, and Lithium Perfluorooctanesulfonate Micellar Solutions: Experimental Determination and Computational Predictions.

This study focuses on determining the partition coefficients (logP) of a diverse set of 63 molecules in three distinct micellar systems: hexadecyltrimethylammonium bromide (HTAB), sodium cholate (SC), and lithium perfluorooctanesulfonate (LPFOS). The experimental log p values were obtained through micellar electrokinetic chromatography (MEKC) experiments, conducted under controlled pH conditions. Then, Quantum Mechanics (QM) and machine learning approaches are proposed for the prediction of the partition coefficients in these three micellar systems. In the applied QM approach, the experimentally obtained partition coefficients were correlated with the calculated values for the case of the 15 solvent mixtures. Using Density Function Theory (DFT) with the B3LYP functional, we calculated the solvation free energies of 63 molecules in these 16 solvents. The combined data from the experimental partition coefficients in the three micellar formulations showed that the 1-propanol/water combination demonstrated the best agreement with the experimental partition coefficients for the SC and HTAB micelles. Moreover, we employed the SVM approach and k-means clustering based on the generation of the chemical descriptor space. The analysis revealed distinct partitioning patterns associated with specific characteristic features within each identified class. These results indicate the utility of the combined techniques when we want an efficient and quicker model for predicting partition coefficients in diverse micelles.

Unusual Aspects of Charge Regulation in Flexible Weak Polyelectrolytes.

This article reviews the state of the art of the studies on charge regulation (CR) effects in flexible weak polyelectrolytes (FWPE). The characteristic of FWPE is the strong coupling of ionization and conformational degrees of freedom. After introducing the necessary fundamental concepts, some unconventional aspects of the the physical chemistry of FWPE are discussed. These aspects are: (i) the extension of statistical mechanics techniques to include ionization equilibria and, in particular, the use of the recently proposed Site Binding-Rotational Isomeric State (SBRIS) model, which allows the calculation of ionization and conformational properties on the same foot; (ii) the recent progresses in the inclusion of proton equilibria in computer simulations; (iii) the possibility of mechanically induced CR in the stretching of FWPE; (iv) the non-trivial adsorption of FWPE on ionized surfaces with the same charge sign as the PE (the so-called "wrong side" of the isoelectric point); (v) the influence of macromolecular crowding on CR.

Molecular dynamics simulations of an α-synuclein NAC domain fragment with a ff14IDPSFF IDP-specific force field suggest ß-sheet intermediate states of fibrillation.

For the discovery of treatments against synucleinopathies, it is necessary to unravel and fully understand the mechanism of fibrillation of proteins involved. Among them, α-synuclein (αS) plays a key role in the development of these diseases through its aggregation into oligomers found in Lewy bodies. However, its structural disorder as an intrinsically disordered protein (IDP) makes its characterization by experimental techniques arduously difficult. Atomistic simulations aim to provide insights into this blank canvas and, fortunately, some studies have already suggested promising mechanisms. Still, it is urgent to consider the IDP features in simulations, so recently a lot of force fields designed to deal with IDPs have been developed. In this study, we have carried out a total of 12 µs simulations of an αS core fragment using a popular ff14SB AMBER force field and the ff14IDPSFF variation that includes a grid-based energy correction map (CMAP) method. The predicted chemical shifts from the simulations and those measured from the αS protein in the NMR solution indicate that ff14IDPSFF reproduces the experimental data more accurately. Moreover, structural analysis exhibits opposite trends between secondary structure propensities. The ff14SB force field preserves the α-helices found in the micelle-bound αS structure, which is used as an initial conformation, while ff14IDPSFF stands out with increased structural disorder and the formation of ß-sheets, which suggests that the IDP-specific force field can capture more suitable conformations representing the possible intermediate states of the fibrillation process.

Adsorption of flexible proteins in the 'wrong side' of the isoelectric point: Casein macropeptide as a model system.

We analyze the conditions of the adsorption of a flexible peptide onto a charged substrate in the 'wrong side' of the isoelectric point (WSIP), i.e. when surface and peptide charges have the same sign. As a model system, we focus on the casein macropeptide (CMP), both in the aglycosylated (aCMP) and fully glycosydated (gCMP) forms. We model the substrate as a uniformly charged plane while CMP is treated as a bead-and-spring model including electrostatic interactions, excluded volume effects and acid/base equilibria. Adsorption coverage, aminoacid charges and concentration profiles are computed by means of Monte Carlo simulations at fixed pH and salt concentration. We conclude that for different reasons the CMP can be adsorbed to both positively and negatively charged surfaces in the WSIP. For negatively charged surfaces, WSIP adsorption is due to the patchy distribution of charges: the peptide is attached to the surface by the positively charged end of the chain, while the repulsion of the surface for the negatively charged tail is screened by the small ions of the added salt. This effect increases with salt concentration. Conversely, a positively charged substrate induces strong charge regulation of the peptide: the acidic groups are deprotonated, and the peptide becomes negatively charged. This effect is stronger at low salt concentrations and it is more intense for gCMP than for aCMP, due to the presence of the additional sialic groups in gCMP.

Inhibition Ability of Natural Compounds on Receptor-Binding Domain of SARS-CoV2: An In Silico Approach.

The lack of medication to treat COVID-19 is still an obstacle that needs to be addressed by all possible scientific approaches. It is essential to design newer drugs with varied approaches. A receptor-binding domain (RBD) is a key part of SARS-CoV-2 virus, located on its surface, that allows it to dock to ACE2 receptors present on human cells, which is followed by admission of virus into cells, and thus infection is triggered. Specific receptor-binding domains on the spike protein play a pivotal role in binding to the receptor. In this regard, the in silico method plays an important role, as it is more rapid and cost effective than the trial and error methods using experimental studies. A combination of virtual screening, molecular docking, molecular simulations and machine learning techniques are applied on a library of natural compounds to identify ligands that show significant binding affinity at the hydrophobic pocket of the RBD. A list of ligands with high binding affinity was obtained using molecular docking and molecular dynamics (MD) simulations for protein-ligand complexes. Machine learning (ML) classification schemes have been applied to obtain features of ligands and important descriptors, which help in identification of better binding ligands. A plethora of descriptors were used for training the self-organizing map algorithm. The model brings out descriptors important for protein-ligand interactions.

Unveiling the Effect of Low pH on the SARS-CoV-2 Main Protease by Molecular Dynamics Simulations.

(1) Background: Main Protease (Mpro) is an attractive therapeutic target that acts in the replication and transcription of the SARS-CoV-2 coronavirus. Mpro is rich in residues exposed to protonation/deprotonation changes which could affect its enzymatic function. This work aimed to explore the effect of the protonation/deprotonation states of Mpro at different pHs using computational techniques. (2) Methods: The different distribution charges were obtained in all the evaluated pHs by the Semi-Grand Canonical Monte Carlo (SGCMC) method. A set of Molecular Dynamics (MD) simulations was performed to consider the different protonation/deprotonation during 250 ns, verifying the structural stability of Mpro at different pHs. (3) Results: The present findings demonstrate that active site residues and residues that allow Mpro dimerisation was not affected by pH changes. However, Mpro substrate-binding residues were altered at low pHs, allowing the increased pocket volume. Additionally, the results of the solvent distribution around Sγ, Hγ, Nδ1 and Hδ1 atoms of the catalytic residues Cys145 and His41 showed a low and high-water affinity at acidic pH, respectively. It which could be crucial in the catalytic mechanism of SARS-CoV-2 Mpro at low pHs. Moreover, we analysed the docking interactions of PF-00835231 from Pfizer in the preclinical phase, which shows excellent affinity with the Mpro at different pHs. (4) Conclusion: Overall, these findings indicate that SARS-CoV-2 Mpro is highly stable at acidic pH conditions, and this inhibitor could have a desirable function at this condition.

Unravelling Constant pH Molecular Dynamics in Oligopeptides with Explicit Solvation Model.

An accurate description of the protonation state of amino acids is essential to correctly simulate the conformational space and the mechanisms of action of proteins or other biochemical systems. The pH and the electrochemical environments are decisive factors to define the effective pKa of amino acids and, therefore, the protonation state. However, they are poorly considered in Molecular Dynamics (MD) simulations. To deal with this problem, constant pH Molecular Dynamics (cpHMD) methods have been developed in recent decades, demonstrating a great ability to consider the effective pKa of amino acids within complex structures. Nonetheless, there are very few studies that assess the effect of these approaches in the conformational sampling. In a previous work of our research group, we detected strengths and weaknesses of the discrete cpHMD method implemented in AMBER when simulating capped tripeptides in implicit solvent. Now, we progressed this assessment by including explicit solvation in these peptides. To analyze more in depth the scope of the reported limitations, we also carried out simulations of oligopeptides with distinct positions of the titratable amino acids. Our study showed that the explicit solvation model does not improve the previously noted weaknesses and, furthermore, the separation of the titratable amino acids in oligopeptides can minimize them, thus providing guidelines to improve the conformational sampling in the cpHMD simulations.

Combinatorial K-Means Clustering as a Machine Learning Tool Applied to Diabetes Mellitus Type 2.

A new original procedure based on k-means clustering is designed to find the most appropriate clinical variables able to efficiently separate into groups similar patients diagnosed with diabetes mellitus type 2 (DMT2) and underlying diseases (arterial hypertonia (AH), ischemic heart disease (CHD), diabetic polyneuropathy (DPNP), and diabetic microangiopathy (DMA)). Clustering is a machine learning tool for discovering structures in datasets. Clustering has been proven to be efficient for pattern recognition based on clinical records. The considered combinatorial k-means procedure explores all possible k-means clustering with a determined number of descriptors and groups. The predetermined conditions for the partitioning were as follows: every single group of patients included patients with DMT2 and one of the underlying diseases; each subgroup formed in such a way was subject to partitioning into three patterns (good health status, medium health status, and degenerated health status); optimal descriptors for each disease and groups. The selection of the best clustering is obtained through the parameter called global variance, defined as the sum of all variance values of all clinical variables of all the clusters. The best clinical parameters are found by minimizing this global variance. This methodology has to identify a set of variables that are assumed to separate each underlying disease efficiently in three different subgroups of patients. The hierarchical clustering obtained for these four underlying diseases could be used to build groups of patients with correlated clinical data. The proposed methodology gives surmised results from complex data based on a relationship with the health status of the group and draws a picture of the prediction rate of the ongoing health status.

Influence of macromolecular crowding on the charge regulation of intrinsically disordered proteins.

In this work we study the coupling between ionization and conformational properties of two IDPs, histatin-5 and ß-amyloid 42, in the presence of neutral and charged crowders. The latter is modeled to resemble bovine serum albumin (BSA). With this aim, semi-grand canonical Monte Carlo simulations are performed, so that the IDP charge is a dynamic property, undergoing protonation/deprotonation processes. Both ionization properties (global and specific amino acid charge and binding capacitance) and radius of gyration are analyzed in a large range of pH values and salt concentrations. Without crowder agents, the titration curve of histatin-5, a polycation, is salt-dependent while that of ß-amyloid 42, a polyampholyte, is almost unaffected. The salt concentration is found to be particularly relevant at pH values where the protein binding capacitance (directly linked with charge fluctuation) is larger. Upon addition of neutral crowders, charge regulation is observed in histatin-5, while for ß-amyloid 42 this effect is very small. The main mechanism for charge regulation is found to be the effective increase in the ionic strength due to the excluded volume. In the presence of charged crowders, a significant increase in the charge of both IDPs is observed in almost all the pH range. In this case, the IDP charge is altered not only by the increase in the effective ionic strength but also by its direct electrostatic interaction with the charged crowders.

Effect of pH on the Supramolecular Structure of Helicobacter pylori Urease by Molecular Dynamics Simulations.

The effect of pH on the supramolecular structure of Helicobacter pylori urease was studied by means of molecular dynamics simulations at seven different pHs. Appropriate urease charge distributions were calculated using a semi-grand canonical Monte Carlo (SGCMC) procedure that assigns each residue's charge state depending on the assigned individual pKa obtained by PROPKA. The effect of pH on protein stability has been analyzed through root-mean-square deviation (RMSD), radius of gyration (RG), solvent-accessible surface area (SASA), hydrogen bonds (HB) and salt bridges (SB). Urease catalyses the hydrolysis of urea in 12 active sites that are covered by mobile regions that act like flaps. The mobility of these flaps is increased at acidic pHs. However, extreme acidic conditions cause urease to have the least number of stabilizing interactions. This initiates the process of denaturalization, wherein the four (αß)3 subunits of the global structure ((αß)3)4 of urease start to separate.

Fuzzy partitioning of clinical data for DMT2 patients.

The present study represents an original approach to data interpretation of clinical data for patients with diagnosis diabetes mellitus type 2 (DMT2) using fuzzy clustering as a tool for intelligent data analysis. Fuzzy clustering is often used in classification and interpretation of medical data (including in medical diagnosis studies) but in this study it is applied with a different goal: to separate a group of 100 patients with DMT2 from a control group of healthy volunteers and, further, to reveal three different patterns of similarity between the patients. Each pattern is described by specific descriptors (variables), which ensure pattern interpretation by appearance of underling disease to DMT2.

Vibrational Analysis of Manganese(II) Oxalates Hydrates: An In Silico Statistical Approach.

The experimental and computational vibrational study for three different manganese(II) oxalates hydrates was explored. The elucidation of IR and Raman spectra were discussed based on their structural singularity; in the same way, they establish some interesting relations between them in the field of computational and statistical approaches. The density functional theory (DFT) computational approach was conducted for accurate prediction and interpretation of the intermolecular effects based on experimental and calculated IR and Raman spectra in the solid-state data in combination with multivariate statistical technique. The proposed computational scheme was also explored for the case of the isolated-molecule model. The goals of the study were to access the accuracy of the proposed procedure for solid-state calculations along with electron calculations for the isolated molecules and to reveal the similarities within the groups of objects by the cluster analysis (CA) techniques and two-way CA for the data. The presented simulation procedure should be very valuable for exploring and to classify other oxalate compounds.

On the Use of the Discrete Constant pH Molecular Dynamics to Describe the Conformational Space of Peptides.

Solvent pH is an important property that defines the protonation state of the amino acids and, therefore, modulates the interactions and the conformational space of the biochemical systems. Generally, this thermodynamic variable is poorly considered in Molecular Dynamics (MD) simulations. Fortunately, this lack has been overcome by means of the Constant pH Molecular Dynamics (CPHMD) methods in the recent decades. Several studies have reported promising results from these approaches that include pH in simulations but focus on the prediction of the effective pKa of the amino acids. In this work, we want to shed some light on the CPHMD method and its implementation in the AMBER suitcase from a conformational point of view. To achieve this goal, we performed CPHMD and conventional MD (CMD) simulations of six protonatable amino acids in a blocked tripeptide structure to compare the conformational sampling and energy distributions of both methods. The results reveal strengths and weaknesses of the CPHMD method in the implementation of AMBER18 version. The change of the protonation state according to the chemical environment is presumably an improvement in the accuracy of the simulations. However, the simulations of the deprotonated forms are not consistent, which is related to an inaccurate assignment of the partial charges of the backbone atoms in the CPHMD residues. Therefore, we recommend the CPHMD methods of AMBER program but pointing out the need to compare structural properties with experimental data to bring reliability to the conformational sampling of the simulations.

Role of Charge Regulation and Fluctuations in the Conformational and Mechanical Properties of Weak Flexible Polyelectrolytes.

This work addresses the role of charge regulation (CR) and the associated fluctuations in the conformational and mechanical properties of weak polyelectrolytes. Due to CR, changes in the pH-value modifies the average macromolecular charge and conformational equilibria. A second effect is that, for a given average charge per site, fluctuations can alter the intensity of the interactions by means of correlation between binding sites. We investigate both effects by means of Monte Carlo simulations at constant pH-value, so that the charge is a fluctuating quantity. Once the average charge per site is available, we turn off the fluctuations by assigning the same average charge to every site. A constant charge MC simulation is then performed. We make use of a model which accounts for the main fundamental aspects of a linear flexible polyelectrolyte that is, proton binding, angle internal rotation, bond stretching and bending. Steric excluded volume and differentiated treatment for short-range and long-range interactions are also included. This model can be regarded as a kind of "minimal" in the sense that it contains a minimum number of parameters but still preserving the atomistic detail. It is shown that, if fluctuations are activated, gauche state bond probabilities increase and the persistence length decreases, so that the polymer becomes more folded. Macromolecular stretching is also analyzed in presence of CR (the charge depends on the applied force) and without CR (the charge is fixed to the value at zero force). The analysis of the low force scaling behavior concludes that Pincus exponent becomes pH-dependent. Both, with and without CR, a transition from 1/2 at high pH-values (phantom chain) to 3/5 at low pH-values (Pincus regime) is observed. Finally, the intermediate force stretching regime is investigated. It is found that CR induces a moderate influence in the force-extension curves and persistence length (which in this force regime becomes force-dependent). It is thus concluded that the effect of CR on the stretching curves is mainly due to the changes in the average charge at zero force. It is also found that, for the cases studied, the effect of steric excluded volume is almost irrelevant compared to electrostatic interactions.

Calculating the Partition Coefficients of Organic Solvents in Octanol/Water and Octanol/Air.

Partition coefficients define how a solute is distributed between two immiscible phases at equilibrium. The experimental estimation of partition coefficients in a complex system can be an expensive, difficult, and time-consuming process. Here a computational strategy to predict the distributions of a set of solutes in two relevant phase equilibria is presented. The octanol/water and octanol/air partition coefficients are predicted for a group of polar solvents using density functional theory (DFT) calculations in combination with a solvation model based on density (SMD) and are in excellent agreement with experimental data. Thus, the use of quantum-chemical calculations to predict partition coefficients from free energies should be a valuable alternative for unknown solvents. The obtained results indicate that the SMD continuum model in conjunction with any of the three DFT functionals (B3LYP, M06-2X, and M11) agrees with the observed experimental values. The highest correlation to experimental data for the octanol/water partition coefficients was reached by the M11 functional; for the octanol/air partition coefficient, the M06-2X functional yielded the best performance. To the best of our knowledge, this is the first computational approach for the prediction of octanol/air partition coefficients by DFT calculations, which has remarkable accuracy and precision.

Macromolecular diffusion in crowded media beyond the hard-sphere model.

The effect of macromolecular crowding on diffusion beyond the hard-core sphere model is studied. A new coarse-grained model is presented, the Chain Entanglement Softened Potential (CESP) model, which takes into account the macromolecular flexibility and chain entanglement. The CESP model uses a shoulder-shaped interaction potential that is implemented in the Brownian Dynamics (BD) computations. The interaction potential contains only one parameter associated with the chain entanglement energetic cost (Ur). The hydrodynamic interactions are included in the BD computations via Tokuyama mean-field equations. The model is used to analyze the diffusion of a streptavidin protein among different sized dextran obstacles. For this system, Ur is obtained by fitting the streptavidin experimental long-time diffusion coefficient Dlongversus the macromolecular concentration for D50 (indicating their molecular weight in kg mol-1) dextran obstacles. The obtained Dlong values show better quantitative agreement with experiments than those obtained with hard-core spheres. Moreover, once parametrized, the CESP model is also able to quantitatively predict Dlong and the anomalous exponent (α) for streptavidin diffusion among D10, D400 and D700 dextran obstacles. Dlong, the short-time diffusion coefficient (Dshort) and α are obtained from the BD simulations by using a new empirical expression, able to describe the full temporal evolution of the diffusion coefficient.

Coupling of Charge Regulation and Conformational Equilibria in Linear Weak Polyelectrolytes: Treatment of Long-Range Interactions via Effective Short-Ranged and pH-Dependent Interaction Parameters.

The classical Rotational Isomeric State (RIS) model, originally proposed by Flory, has been used to rationalize a wide range of physicochemical properties of neutral polymers. However, many weak polyelectrolytes of interest are able to regulate their charge depending on the conformational state of the bonds. Recently, it has been shown that the RIS model can be coupled with the Site Binding (SB) model, for which the ionizable sites can adopt two states: protonated or deprotonated. The resulting combined scheme, the SBRIS model, allows for analyzing ionization and conformational equilibria on the same foot. In the present work, this approach is extended to include pH-dependent electrostatic Long-Range (LR) interactions, ubiquitous in weak polyelectrolytes at moderate and low ionic strengths. With this aim, the original LR interactions are taken into account by defining effective Short-Range (SR) and pH-dependent parameters, such as effective microscopic protonation constants and rotational bond energies. The new parameters are systematically calculated using variational methods. The machinery of statistical mechanics for SR interactions, including the powerful and fast transfer matrix methods, can then be applied. The resulting technique, which we will refer to as the Local Effective Interaction Parameters (LEIP) method, is illustrated with a minimal model of a flexible linear polyelectrolyte containing only one type of rotating bond. LEIP reproduces very well the pH dependence of the degree of protonation and bond probabilities obtained by semi-grand canonical Monte Carlo simulations, where LR interactions are explicitly taken into account. The reduction in the computational time in several orders of magnitude suggests that the LEIP technique could be useful in a range of areas involving linear weak polyelectrolytes, allowing direct fitting of the relevant physical parameters to the experimental quantities.

Pre-selection and assessment of green organic solvents by clustering chemometric tools.

The study presents the result of the application of chemometric tools for selection of physicochemical parameters of solvents for predicting missing variables - bioconcentration factors, water-octanol and octanol-air partitioning constants. EPI Suite software was successfully applied to predict missing values for solvents commonly considered as "green". Values for logBCF, logKOW and logKOA were modelled for 43 rather nonpolar solvents and 69 polar ones. Application of multivariate statistics was also proved to be useful in the assessment of the obtained modelling results. The presented approach can be one of the first steps and support tools in the assessment of chemicals in terms of their greenness.

Ionization and Conformational Equilibria of Citric Acid: Delocalized Proton Binding in Solution.

The microspeciation of citric acid is studied by analyzing NMR titration data. When the site binding (SB) model, which assumes fully localized proton binding to the carboxylic groups, is used to obtain microscopic energy parameters (dissociation constants, pair and triplet interaction energies between charged carboxylate groups), contradictory results are obtained. The resulting macroscopic constants are in very good agreement with the values reported in the literature using potentiometry. However, the found pair interaction energy between the terminal carboxylates and the triplet interaction energy are physically meaningless. To solve this apparent contradiction, we consider the possibility of delocalized proton binding, so that the proton can be exchanged at high velocity in the NMR time scale through short, strong, low-barrier (SSLB) hydrogen bonds. With this aim, ab initio MP2 calculations using the SMD polarizable continuum model for the solvent were performed and the fully roto-microspeciation elucidated. First, fully localized proton binding was assumed, and the resulting microstate probabilities are in reasonable agreement with those reported in previous works that use selective blocking of the carboxylic groups. They are, however, in clear disagreement with the microstate probabilities derived from the NMR titration data, which predict, within a very narrow confidence interval, a unique microspecies for the symmetric di-ionized form. Moreover, counterintuitively, the interaction between terminal charged groups is much larger than that between central and terminal groups. As a consequence, we have explored the possibility of delocalized proton binding by calculating the energy of intermediate proton positions between two carbolxylic groups. The results reveal that the exchange of the proton through the hydrogen bonds is in some cases produced without energetic barrier. This effect is specially relevant in the di-ionized form, with all the most stable conformations forming a SSLB, which together would constitute the only microstate detected by NMR. An alternative reaction scheme for the ionization process, based on proton delocalization, is proposed.

Insight into electric field-induced rupture mechanism of water-in-toluene emulsion films from a model system.

This paper presents a model, which we have designed to get insight into the development of electro-induced instability of a thin toluene emulsion film in contact with the saline aqueous phase. Molecular dynamics (MD) simulations demonstrate the role of charge accumulation in the toluene-film rupture induced by a DC electric field. Two ensembles-NVT and NPT-are used to determine the critical value of the external field at which the film ruptures, the charge distribution and capacitance of the thin film, number densities, and the film structure. The rupture mechanism as seen from this model is the following: in both NVT and NPT ensembles, condenser plates, where the charge density is maximal, are situated at the very border between the bulk aqueous (water) phase and the mixed layer. No ion penetration is observed within the toluene core, thus leaving all the distribution of charges within the mixed zone and the bulk phase that could be attributed to the formation of hydration shells. When the critical electric field is reached within a certain time after the field application, electric discharge occurs indicating the beginning of the rupturing process. The MD simulations indicate that the NPT ensemble predicts a value of the critical field that is closer to the experimental finding.