Search

The presence of antibiotics in soils is increasing drastically in last decades due to the intensive farming industry and excessive human consumption. Clay minerals are one of the soil components with great adsorption capacity for organic pollutants. The study of interactions between antibiotics and mineral surfaces will give us scientific knowledge of these pollutants through soils. In this work, we study the adsorption of the antibiotic ciprofloxacin in the clay mineral fraction of soils from the Argentinian zone of Santa Rosa (Corrientes), in a collaborative research of experiments and atomistic modelling calculations of the intercalation of ciprofloxacin in the interlayer space of montmorillonite. Adsorption and desorption isotherms were performed and compared with different isotherm models. Additionally, enthalpy, entropy, and free energy were determined from equilibrium constants at a function of temperature. All these experiments and calculations lead to the conclusions that two adsorption types of ciprofloxacin are found on clay minerals: one weakly sorbed that is released during the desorption experiments, and other one strongly joined that remains in the soil. [ABSTRACT FROM AUTHOR]

Machine learning potentials, particularly the deep potential (DP) model, have revolutionized molecular dynamics (MD) simulations, striking a balance between accuracy and computational efficiency. To facilitate the DP model's integration with the popular MD engine OpenMM, we have developed a versatile OpenMM plugin. This plugin supports a range of applications, from conventional MD simulations to alchemical free energy calculations and hybrid DP/MM simulations. Our extensive validation tests encompassed energy conservation in microcanonical ensemble simulations, fidelity in canonical ensemble generation, and the evaluation of the structural, transport, and thermodynamic properties of bulk water. The introduction of this plugin is expected to significantly expand the application scope of DP models within the MD simulation community, representing a major advancement in the field. [ABSTRACT FROM AUTHOR]

Thermal conductivities of polybutadiene crosslinked with sulfur as a function of the heat flux autocorrelation function by using an equilibrium molecular dynamic (EMD) simulation were investigated. The Green–Kubo method was used to calculate thermal conductivities. All simulations were performed by applying the LAMMPS software (version 3 Mar 2020) package. The united-atom force field (OPLS-UA) from the Moltemplate software (version 2.20.3) was applied in the simulations. The influence of uniform and random distributions of sulfur in polybutadiene on the final value of thermal conductivities was studied by polymeric model structures with similar and variable degrees of crosslinking. The results showed that for identical degrees of crosslinking, the distribution of crosslinkers in the polymeric model structures significantly influenced the final value of thermal conductivity. Moreover, the influence of the crosslinking degree on the final value of thermal conductivity was studied by considering polymeric model structures with different degrees of crosslinking. The results demonstrate that by having a random distribution of sulfur, the thermal conductivity will be enhanced. However, by increasing the degree of crosslinking to the higher percentage in random crosslinked model structures, the value of thermal conductivity drops significantly due to possible higher crystallization of the model structures, which decrease the degree of freedom for phonon contributions. [ABSTRACT FROM AUTHOR]

This paper reviews molecular dynamics (MD) concepts on heat transfer analysis of supercritical CO2, and highlights the major parameters that can affect the accuracy of respective thermal coefficients. Subsequently, the prime aspects of construction, transfer identification, and thermal performance are organized according to their challenges and prospective solutions associated with the mutability of supercritical CO2 properties. Likewise, the characteristics of bound force field schemes and thermal relaxation approaches are discussed on a case-by-case basis. Both convective and diffusive states of trans- and supercritical CO2 are debated, given their magnitude effects on molecular interactions. Following the scarcity of literature on similar enquiries, this paper recommended a future series of studies on molecular dynamics models in a large region of supercriticality and phase-interactions for coupled heat and mass transfer systems. This review recognizes that the foremost undertaking is to ascertain the thermo-hydraulic identity of supercritical CO2 for process feasibility of developed technology. [ABSTRACT FROM AUTHOR]

Boron containing compounds (BCCs) aroused increasing interest in the scientific community due to their wide application as drugs in various fields. In order to design new compounds hopefully endowed with pharmacological activity and also investigate their conformational behavior, the support of computational studies is crucial. Nevertheless, the suitable molecular mechanics parameterization and the force fields needed to perform these simulations are not completely available for this class of molecules. In this paper, Amber force field parameters for phenyl-, benzyl-, benzylamino-, and methylamino-boronates, a group of boron-containing compounds involved in different branches of the medicinal chemistry, were created. The robustness of the obtained data was confirmed through molecular dynamics simulations on ligand/β-lactamases covalent complexes. The ligand torsional angles, populated over the trajectory frames, were confirmed by values found in the ligand geometries, located through optimizations at the DFT/B3LYP/6-31g(d) level, using water as a solvent. In summary, this study successfully provided a library of parameters, opening the possibility to perform molecular dynamics simulations of this class of boron-containing compounds. [ABSTRACT FROM AUTHOR]

Host–guest binding, despite the relatively simple structural and chemical features of individual components, still poses a challenge in computational modelling. The extreme underperformance of standard end-point methods in host–guest binding makes them practically useless. In the current work, we explore a potentially promising modification of the three-trajectory realization. The alteration couples the binding-induced structural reorganization into free energy estimation and suffers from dramatic fluctuations in internal energies in protein–ligand situations. Fortunately, the relatively small size of host–guest systems minimizes the magnitude of internal fluctuations and makes the three-trajectory realization practically suitable. Due to the incorporation of intra-molecular interactions in free energy estimation, a strong dependence on the force field parameters could be incurred. Thus, a term-specific investigation of transferable GAFF derivatives is presented, and noticeable differences in many aspects are identified between commonly applied GAFF and GAFF2. These force-field differences lead to different dynamic behaviors of the macrocyclic host, which ultimately would influence the end-point sampling and binding thermodynamics. Therefore, the three-trajectory end-point free energy calculations are performed with both GAFF versions. Additionally, due to the noticeable differences between host dynamics under GAFF and GAFF2, we add additional benchmarks of the single-trajectory end-point calculations. When only the ranks of binding affinities are pursued, the three-trajectory realization performs very well, comparable to and even better than the regressed PBSA_E scoring function and the dielectric constant-variable regime. With the GAFF parameter set, the TIP3P water in explicit solvent sampling and either PB or GB implicit solvent model in free energy estimation, the predictive power of the three-trajectory realization in ranking calculations surpasses all existing end-point methods on this dataset. We further combine the three-trajectory realization with another promising modified end-point regime of varying the interior dielectric constant. The combined regime does not incur sizable improvements for ranks and deviations from experiment exhibit non-monotonic variations. [ABSTRACT FROM AUTHOR]

Hydration free energies of small molecules are commonly used as benchmarks for solvation models. However, errors in predicting hydration free energies are partially due to the force fields used and not just the solvation model. To address this, we have used the 3D reference interaction site model (3D-RISM) of molecular solvation and existing benchmark explicit solvent calculations with a simple element count correction (ECC) to identify problems with the non-bond parameters in the general AMBER force field (GAFF). 3D-RISM was used to calculate hydration free energies of all 642 molecules in the FreeSolv database, and a partial molar volume correction (PMVC), ECC, and their combination (PMVECC) were applied to the results. The PMVECC produced a mean unsigned error of 1.01 ± 0.04 kcal / mol and root mean squared error of 1.44 ± 0.07 kcal / mol , better than the benchmark explicit solvent calculations from FreeSolv, and required less than 15 s of computing time per molecule on a single CPU core. Importantly, parameters for PMVECC showed systematic errors for molecules containing Cl, Br, I, and P. Applying ECC to the explicit solvent hydration free energies found the same systematic errors. The results strongly suggest that some small adjustments to the Lennard–Jones parameters for GAFF will lead to improved hydration free energy calculations for all solvent models. [ABSTRACT FROM AUTHOR]

Charge scaling as an effective solution to the experiment–computation disagreement in molecular modelling of ionic liquids (ILs) could bring the computational results close to the experimental reference for various thermodynamic properties. According to the large-scale benchmark calculations of mass density, solvation, and water-ILs transfer-free energies in our series of papers, the charge-scaling factor of 0.8 serves as a near-optimal option generally applicable to most ILs, although a system-dependent parameter adjustment could be attempted for further improved performance. However, there are situations in which such a charge-scaling treatment would fail. Namely, charge scaling cannot really affect the simulation outcome, or minimally perturbs the results that are still far from the experimental value. In such situations, the vdW radius as an additional adjustable parameter is commonly tuned to minimize the experiment–calculation deviation. In the current work, considering two ILs from the quinuclidinium family, we investigate the impacts of this vdW-scaling treatment on the mass density and the solvation/partition thermodynamics in a fashion similar to our previous charge-scaling works, i.e., scanning the vdW-scaling factor and computing physical properties under these parameter sets. It is observed that the mass density exhibits a linear response to the vdW-scaling factor with slopes close to −1.8 g/mL. By further investigating a set of physiochemically relevant temperatures between 288 K and 348 K, we confirm the robustness of the vdW-scaling treatment in the estimation of bulk properties. The best vdW-scaling parameter for mass density would worsen the computation of solvation/partition thermodynamics, and a marginal decrease in the vdW-scaling factor is considered as an intermediate option balancing the reproductions of bulk properties and solvation thermodynamics. These observations could be understood in a way similar to the charge-scaling situation. i.e., overfitting some properties (e.g., mass density) would degrade the accuracy of the other properties (e.g., solvation free energies). Following this principle, the general guideline for applying this vdW-tuning protocol is by using values between the density-derived choice and the solvation/partition-derived solution. The charge and current vdW scaling treatments cover commonly encountered ILs, completing the protocol for accurate modelling of ILs with fixed-charge force fields. [ABSTRACT FROM AUTHOR]

"Plant-type" ferredoxins (Fds) in the thylakoid membranes of plants, algae, and cyanobacteria possess a single [2Fe-2S] cluster in active sites and mediate light-induced electron transfer from Photosystem I reaction centers to various Fd-dependent enzymes. Structural knowledge of plant-type Fds is relatively limited to static structures, and the detailed behavior of oxidized and reduced Fds has not been fully elucidated. It is important that the investigations of the effects of active-center reduction on the structures and dynamics for elucidating electron-transfer mechanisms. In this study, model systems of oxidized and reduced Fds were constructed from the high-resolution crystal structure of Chlamydomonas reinhardtii Fd1, and three 200 ns molecular dynamics simulations were performed for each system. The force field parameters of the oxidized and reduced active centers were independently obtained using quantum chemical calculations. There were no substantial differences in the global conformations of the oxidized and reduced forms. In contrast, active-center reduction affected the hydrogen-bond network and compactness of the surrounding residues, leading to the increased flexibility of the side chain of Phe61, which is essential for the interaction between Fd and the target protein. These computational results will provide insight into the electron-transfer mechanisms in the Fds. [ABSTRACT FROM AUTHOR]

Phosphorylation plays a key role in plant biology, such as the accumulation of plant cells to form the observed proteome. Statistical analysis found that many phosphorylation sites are located in disordered regions. However, current force fields are mainly trained for structural proteins, which might not have the capacity to perfectly capture the dynamic conformation of the phosphorylated proteins. Therefore, we evaluated the performance of ff03CMAP, a balanced force field between structural and disordered proteins, for the sampling of the phosphorylated proteins. The test results of 11 different phosphorylated systems, including dipeptides, disordered proteins, folded proteins, and their complex, indicate that the ff03CMAP force field can better sample the conformations of phosphorylation sites for disordered proteins and disordered regions than ff03. For the solvent model, the results strongly suggest that the ff03CMAP force field with the TIP4PD water model is the best combination for the conformer sampling. Additional tests of CHARMM36m and FB18 force fields on two phosphorylated systems suggest that the overall performance of ff03CMAP is similar to that of FB18 and better than that of CHARMM36m. These results can help other researchers to choose suitable force field and solvent models to investigate the dynamic properties of phosphorylation proteins. [ABSTRACT FROM AUTHOR]

The 3D-reference interaction site model (3D-RISM) molecular solvation theory in combination with the Kovalenko–Hirata (KH) closure is extended to seven heterocyclic liquids to understand their liquid states and to test the performance of the theory in solvation free energy (SFE) calculations of solutes in select solvents. The computed solvent site distribution profiles were compared with the all-atom molecular dynamics (MD) simulations, showing comparable performances. The computational results were compared against the structural parameters for liquids, whenever available, as well as against the experimental SFEs. The liquids are found to have local ordered structures held together via weak interactions in both the RISM and MD simulations. The 3D-RISM-KH computed SFEs are in good agreement with the benchmark values for the tetrahydrothiophene-S,S-dioxide, and showed comparatively larger deviations in the case of the SFEs in the tetrahydrofuran continuum. [ABSTRACT FROM AUTHOR]

Metal–organic frameworks (MOFs) are a highly versatile group of porous materials suitable for a broad range of applications, which often crucially depend on the MOFs' heat transport properties. Nevertheless, detailed relationships between the chemical structure of MOFs and their thermal conductivities are still largely missing. To lay the foundations for developing such relationships, we performed non-equilibrium molecular dynamics simulations to analyze heat transport in a selected set of materials. In particular, we focus on the impact of organic linkers, the inorganic nodes and the interfaces between them. To obtain reliable data, great care was taken to generate and thoroughly benchmark system-specific force fields building on ab-initio-based reference data. To systematically separate the different factors arising from the complex structures of MOF, we also studied a series of suitably designed model systems. Notably, besides the expected trend that longer linkers lead to a reduction in thermal conductivity due to an increase in porosity, they also cause an increase in the interface resistance between the different building blocks of the MOFs. This is relevant insofar as the interface resistance dominates the total thermal resistance of the MOF. Employing suitably designed model systems, it can be shown that this dominance of the interface resistance is not the consequence of the specific, potentially weak, chemical interactions between nodes and linkers. Rather, it is inherent to the framework structures of the MOFs. These findings improve our understanding of heat transport in MOFs and will help in tailoring the thermal conductivities of MOFs for specific applications. [ABSTRACT FROM AUTHOR]

Molecular mechanics force field calculations have historically shown significant limitations in modeling the energetic and conformational interconversions of highly substituted furanose rings. This is primarily due to the gauche effect that is not easily captured using pairwise energy potentials. In this study, we present a refinement to the set of torsional parameters in the General Amber Force Field (gaff) used to calculate the potential energy of mono, di-, and gem-fluorinated nucleosides. The parameters were optimized to reproduce the pseudorotation phase angle and relative energies of a diverse set of mono- and difluoro substituted furanose ring systems using quantum mechanics umbrella sampling techniques available in the IpolQ engine in the Amber suite of programs. The parameters were developed to be internally consistent with the gaff force field and the TIP3P water model. The new set of angle and dihedral parameters and partial charges were validated by comparing the calculated phase angle probability to those obtained from experimental nuclear magnetic resonance experiments. [ABSTRACT FROM AUTHOR]

Self-assembled cyclic peptide nanotubes with alternating D- and L-amino acid residues in the sequence of each subunit have attracted a great deal of attention due to their potential for new nanotechnology and biomedical applications, mainly in the field of antimicrobial peptides. Molecular dynamics simulations can be used to characterize these systems with atomic resolution at different time scales, providing information that is difficult to obtain via wet lab experiments. However, the performance of classical force fields typically employed in the simulation of biomolecules has not yet been extensively tested with this kind of highly constrained peptide. Four different classical force fields (AMBER, CHARMM, OPLS, and GROMOS), using a nanotube formed by eight D,L-α-cyclic peptides inserted into a lipid bilayer as a model system, were employed here to fill this gap. Significant differences in the pseudo-cylindrical cavities formed by the nanotubes were observed, the most important being the diameter of the nanopores, the number and location of confined water molecules, and the density distribution of the solvent molecules. Furthermore, several modifications were performed on GROMOS54a7, aiming to explore acceleration strategies of the MD simulations. The hydrogen mass repartitioning (HMR) and hydrogen isotope exchange (HIE) methods were tested to slow down the fastest degrees of freedom. These approaches allowed a significant increase in the time step employed in the equation of the motion integration algorithm, from 2 fs up to 5–7 fs, with no serious changes in the structural and dynamical properties of the nanopores. Subtle differences with respect to the simulations with the unmodified force fields were observed in the concerted movements of the cyclic peptides, as well as in the lifetime of several H-bonds. All together, these results are expected to contribute to better understanding of the behavior of self-assembled cyclic peptide nanotubes, as well as to support the methods tested to speed up general MD simulations; additionally, they do provide a number of quantitative descriptors that are expected to be used as a reference to design new experiments intended to validate and complement computational studies of antimicrobial cyclic peptides. [ABSTRACT FROM AUTHOR]

We investigated the effect of fluorinated molecules on dipalmitoylphosphatidylcholine (DPPC) bilayers by force-field molecular dynamics simulations. In the first step, we developed all-atom force-field parameters for additive molecules in membranes to enable an accurate description of those systems. On the basis of this force field, we performed extensive simulations of various bilayer systems containing different additives. The additive molecules were chosen to be of different size and shape, and they included small molecules such as perfluorinated alcohols, but also more complex molecules. From these simulations, we investigated the structural and dynamic effects of the additives on the membrane properties, as well as the behavior of the additive molecules themselves. Our results are in good agreement with other theoretical and experimental studies, and they contribute to a microscopic understanding of interactions, which might be used to specifically tune membrane properties by additives in the future. [ABSTRACT FROM AUTHOR]

To study the microscopic thermal aging mechanism of insulating paper cellulose through molecular dynamics simulation, it is important to select suitable DP (Degree of Polymerization) and force field for the cellulose model to shorten the simulation time and obtain correct and objective simulation results. Here, the variation of the mechanical properties and solubility parameters of models with different polymerization degrees and force fields were analyzed. Numerous cellulose models with different polymerization degrees were constructed to determine the relative optimal force field from the perspectives of the similarity of the density of cellulose models in equilibrium to the actual cellulose density, and the volatility and repeatability of the mechanical properties of the models through the selection of a stable polymerization degree using the two force fields. The results showed that when the polymerization degree was more than or equal to 10, the mechanical properties and solubility of cellulose models with the COMPASS (Condensed-phase Optimized Molecular Potential for Atomistic Simulation Studies) and PCFF (Polymer Consistent Force Field) force fields were in steady states. The steady-state density of the cellulose model using the COMPASS force field was closer to the actual density of cellulose. Thus, the COMPASS force field is favorable for molecular dynamics simulation of amorphous cellulose. [ABSTRACT FROM AUTHOR]

This review discusses the use of molecular modeling tools, together with existing experimental findings, to provide a complete atomic-level description of enzyme dynamics and function. We focus on functionally relevant conformational dynamics of enzymes and the protonation states of substrates. The conformational fluctuations of enzymes usually play a crucial role in substrate recognition and catalysis. Protein dynamics can be altered by a tiny change in a molecular system such as different protonation states of various intermediates or by a significant perturbation such as a ligand association. Here we review recent advances in applying atomistic molecular dynamics (MD) simulations to investigate allosteric and network regulation of tryptophan synthase (TRPS) and protonation states of its intermediates and catalysis. In addition, we review studies using quantum mechanics/molecular mechanics (QM/MM) methods to investigate the protonation states of catalytic residues of β-Ketoacyl ACP synthase I (KasA). We also discuss modeling of large-scale protein motions for HIV-1 protease with coarse-grained Brownian dynamics (BD) simulations. [ABSTRACT FROM AUTHOR]

The quantum quasi-docking procedure is used to compare the docking accuracies of two quantum-chemical semiempirical methods, namely, PM6-D3H4X and PM7. Quantum quasi-docking is an approximation to quantum docking. In quantum docking, it is necessary to search directly for the global minimum of the energy of the protein-ligand complex calculated by the quantum-chemical method. In quantum quasi-docking, firstly, we look for a wide spectrum of low-energy minima, calculated using the MMFF94 force field, and secondly, we recalculate the energies of all these minima using the quantum-chemical method, and among these recalculated energies we determine the lowest energy and the corresponding ligand position. Both PM6-D3H4X and PM7 are novel methods that describe well-dispersion interactions, hydrogen and halogen bonds. The PM6-D3H4X and PM7 methods are used with the COSMO implicit solvent model as it is implemented in the MOPAC program. The comparison is made for 25 high quality protein-ligand complexes. Firstly, the docking positioning accuracies have been compared, and we demonstrated that PM7+COSMO provides better positioning accuracy than PM6-D3H4X. Secondly, we found that PM7+COSMO demonstrates a much higher correlation between the calculated and measured protein–ligand binding enthalpies than PM6-D3H4X. For future quantum docking PM7+COSMO is preferable, but the COSMO model must be improved. [ABSTRACT FROM AUTHOR]

Thermal conductivities of polybutadiene rubbers crosslinked by 2.4 and 2.8 phr of sulfur have been found to be functions of temperature via molecular dynamics (MD) simulations using the Green–Kubo method. From an analysis of the heat flux autocorrelation functions, it has been revealed that the dominant means of heat transport in rubbers is governed by deformations of polymeric chains. Thermal conductivities of rubber samples vulcanized by 2.4 and 2.8 phr of sulfur have been measured by the heat flow meter method between 0 ∘ C and 60 ∘ C at atmospheric pressure. The temperature dependencies of the thermal conductivities of rubbers and their glass transition temperatures derived from MD simulations are in good agreement with the literature and experimental data. Details are discussed in the paper. [ABSTRACT FROM AUTHOR]

The development of novel materials has challenges besides their synthesis. Materials such as novel MXenes are difficult to probe experimentally due to their reduced size and low stability under ambient conditions. Quantum mechanics and molecular dynamics simulations have been valuable options for material properties determination. However, computational materials scientists may still have difficulty finding specific force field models for their simulations. Force fields are usually hard to parametrize, and their parameters' determination is computationally expensive. We show the Lennard-Jones (2-body interactions) combined with the Axilrod-Teller (3-body interactions) parametrization process' applicability for metals and new classes of materials (MXenes). Because this parametrization process is simple and computationally inexpensive, it allows users to predict materials' behaviors under close-to-ambient conditions in molecular dynamics, independent of pre-existing potential files. Using the process described in this work, we have made the T i 2 C parameters set available for the first time in a peer-reviewed work. [ABSTRACT FROM AUTHOR]

While the construction of a dependable force field for performing classical molecular dynamics (MD) simulation is crucial for elucidating the structure and function of biomolecular systems, the attempts to do this for glycans are relatively sparse compared to those for proteins and nucleic acids. Currently, the use of GLYCAM06 force field is the most popular, but there have been a number of concerns about its accuracy in the systematic description of structural changes. In the present work, we focus on the improvement of the GLYCAM06 force field for β -d-glucose, a simple and the most abundant monosaccharide molecule, with the aid of machine learning techniques implemented with the TensorFlow library. Following the pre-sampling over a wide range of configuration space generated by MD simulation, the atomic charge and dihedral angle parameters in the GLYCAM06 force field were re-optimized to accurately reproduce the relative energies of β -d-glucose obtained by the density functional theory (DFT) calculations according to the structural changes. The validation for the newly proposed force-field parameters was then carried out by verifying that the relative energy errors compared to the DFT value were significantly reduced and that some inconsistencies with experimental (e.g., NMR) results observed in the GLYCAM06 force field were resolved relevantly. [ABSTRACT FROM AUTHOR]

The spreading of antibiotic-resistant bacteria strains is one of the most serious problem in medicine to struggle nowadays. This triggered the development of alternative antimicrobial agents in recent years. One of such group is Gemini surfactants which are massively synthesised in various structural configurations to obtain the most effective antibacterial properties. Unfortunately, the comparison of antimicrobial effectiveness among different types of Gemini agents is unfeasible since various protocols for the determination of Minimum Inhibitory Concentration are used. In this work, we proposed alternative, computational, approach for such comparison. We designed a comprehensive database of 250 Gemini surfactants. Description of structure parameters, for instance spacer type and length, are included in the database. We parametrised modelled molecules to obtain force fields for the entire Gemini database. This was used to conduct in silico studies using the molecular dynamics to investigate the incorporation of these agents into model E. coli inner membrane system. We evaluated the effect of Gemini surfactants on structural, stress and mechanical parameters of the membrane after the agent incorporation. This enabled us to select four most likely membrane properties that could correspond to Gemini's antimicrobial effect. Based on our results we selected several types of Gemini spacers which could demonstrate a particularly strong effect on the bacterial membranes. [ABSTRACT FROM AUTHOR]

Organophosphate hydrolases are promising as potential biotherapeutic agents to treat poisoning with pesticides or nerve gases. However, these enzymes often need to be further engineered in order to become useful in practice. One example of such enhancement is the alteration of enantioselectivity of diisopropyl fluorophosphatase (DFPase). Molecular modeling techniques offer a unique opportunity to address this task rationally by providing a physical description of the substrate-binding process. However, DFPase is a metalloenzyme, and correct modeling of metal cations is a challenging task generally coming with a tradeoff between simulation speed and accuracy. Here, we probe several molecular mechanical parameter combinations for their ability to empower long simulations needed to achieve a quantitative description of substrate binding. We demonstrate that a combination of the Amber19sb force field with the recently developed 12-6 Ca2+ models allows us to both correctly model DFPase and obtain new insights into the DFP binding process. [ABSTRACT FROM AUTHOR]

Interactions of bioorganic moieties with clay minerals have attracted attention not only from the perspective of novel bioclay materials but also because they play a crucial role in our understanding of physical and chemical processes in soils. The aim of the present article is to explore the interactions responsible for the formation of a phosphatidylcholine-kaolinite bioclay by employing a series of classical molecular dynamic simulations. Detailed analysis of the structure and energies of the resulting bioclays reveals that the phosphatidylcholine molecules bind to the kaolinite surface either via their zwitterionic heads or hydrophobic aliphatic tails, depending on the kaolinite surface characteristics and the density of organic coating. The phosphatidylcholine molecules have a tendency to form irregular layers with a preferred parallel orientation of molecules with respect to the kaolinite surface. The tails exhibit varying degrees of flexibility and disorder depending on their distance from the surface and the density of surface coating. Significant differences in the binding can be spotted with respect to the two types of kaolinite basal surfaces, i.e., the hydrophobic siloxane surface, which possesses a considerable dispersion character, and the hydrophilic alumina surface, polarized by the surface hydroxyl groups. [ABSTRACT FROM AUTHOR]

Coarse-grained (CG) force fields have become promising tools for studies of protein behavior, but the balance of speed and accuracy is still a challenge in the research of protein coarse graining methodology. In this work, 20 CG beads have been designed based on the structures of amino acid residues, with which an amino acid can be represented by one or two beads, and a CG solvent model with five water molecules was adopted to ensure the consistence with the protein CG beads. The internal interactions in protein were classified according to the types of the interacting CG beads, and adequate potential functions were chosen and systematically parameterized to fit the energy distributions. The proposed CG force field has been tested on eight proteins, and each protein was simulated for 1000 ns. Even without any extra structure knowledge of the simulated proteins, the Ca root mean square deviations (RMSDs) with respect to their experimental structures are close to those of relatively short time all atom molecular dynamics simulations. However, our coarse grained force field will require further refinement to improve agreement with and persistence of native-like structures. In addition, the root mean square fluctuations (RMSFs) relative to the average structures derived from the simulations show that the conformational fluctuations of the proteins can be sampled. [ABSTRACT FROM AUTHOR]

Structural characterization of poly-L-lactic acid (PLLA) and poly(glycolic acid) (PGA) oligomers containing three units was carried out with an atomistic approach. Oligomer structures were first optimized through quantum chemical calculations, using density functional theory (DFT); rotational barriers concerning dihedral angles along the chain were then investigated. Diffusion coefficients of L-lactic acid and glycolic acid in pure water were estimated through molecular dynamic (MD) simulations. Monomer structures were obtained with quantum chemical computation in implicit water using DFT method; atomic charges were fitted with Restrained Electrostatic Potentials (RESP) formalism, starting from electrostatic potentials calculated with quantum chemistry. MD simulations were carried out in explicit water, in order to take into account solvent presence. [ABSTRACT FROM AUTHOR]

G-quadruplexes have drawn widespread attention for serving as a potential anti-cancer target and their application in material science. Molecular dynamics (MD) simulation is the key theoretical tool in the study of GQ's structure-function relationship. In this article, we systematically benchmarked the five force fields of parmbsc0, parmbsc1, OL15, AMOEBA, and Drude2017 on the MD simulation of G-quadruplex from four aspects: structural stability, central ion channel stability, description of Hoogsteen hydrogen bond network, and description of the main chain dihedral angle. The results show that the overall performance of the Drude force field is the best. Although there may be a certain over-polarization effect, it is still the best choice for the MD simulation of G-quadruplexes. [ABSTRACT FROM AUTHOR]

For the first time, the thermal conductivities of vulcanized polybutadiene and polyisoprene have been investigated according to their degree of crosslinking. The C-C and C-S-S-C crosslink bridges, which can be obtained via vulcanization processes using peroxides and sulfur, respectively, are considered. The temperature dependence of the thermal conductivity of soft rubber derived from molecular dynamics (MD) simulations is in very good agreement with the experimental results. The contributions of bonded and non-bonded interactions in the MD simulations and their influence on the thermal conductivities of polyisoprene and polybutadiene are presented. The details are discussed in this paper. [ABSTRACT FROM AUTHOR]

Glycosaminoglycans (GAGs) are the linear carbohydrate components of proteoglycans (PGs) and are key mediators in the bioactivity of PGs in animal tissue. GAGs are heterogeneous, conformationally complex, and polydisperse, containing up to 200 monosaccharide units. These complexities make studying GAG conformation a challenge for existing experimental and computational methods. We previously described an algorithm we developed that applies conformational parameters (i.e., all bond lengths, bond angles, and dihedral angles) from molecular dynamics (MD) simulations of nonsulfated chondroitin GAG 20-mers to construct 3-D atomic-resolution models of nonsulfated chondroitin GAGs of arbitrary length. In the current study, we applied our algorithm to other GAGs, including hyaluronan and nonsulfated forms of dermatan, keratan, and heparan and expanded our database of MD-generated GAG conformations. Here, we show that individual glycosidic linkages and monosaccharide rings in 10- and 20-mers of hyaluronan and nonsulfated dermatan, keratan, and heparan behave randomly and independently in MD simulation and, therefore, using a database of MD-generated 20-mer conformations, that our algorithm can construct conformational ensembles of 10- and 20-mers of various GAG types that accurately represent the backbone flexibility seen in MD simulations. Furthermore, our algorithm efficiently constructs conformational ensembles of GAG 200-mers that we would reasonably expect from MD simulations. [ABSTRACT FROM AUTHOR]

Six empirical force fields were tested for applicability to calculations for automated carbohydrate database filling. They were probed on eleven disaccharide molecules containing representative structural features from widespread classes of carbohydrates. The accuracy of each method was queried by predictions of nuclear Overhauser effects (NOEs) from conformational ensembles obtained from 50 to 100 ns molecular dynamics (MD) trajectories and their comparison to the published experimental data. Using various ranking schemes, it was concluded that explicit solvent MM3 MD yielded non-inferior NOE accuracy with newer GLYCAM-06, and ultimately PBE0-D3/def2-TZVP (Triple-Zeta Valence Polarized) Density Functional Theory (DFT) simulations. For seven of eleven molecules, at least one empirical force field with explicit solvent outperformed DFT in NOE prediction. The aggregate of characteristics (accuracy, speed, and compatibility) made MM3 dynamics with explicit solvent at 300 K the most favorable method for bulk generation of disaccharide conformation maps for massive database filling. [ABSTRACT FROM AUTHOR]

Intrinsically Disordered Peptides and Proteins (IDPs) in solution can span a broad range of conformations that often are hard to characterize by both experimental and computational methods. However, obtaining a significant representation of the conformational space is important to understand mechanisms underlying protein functions such as partner recognition. In this work, we investigated the behavior of the Sic1 Kinase-Inhibitor Domain (KID) in solution by Molecular Dynamics (MD) simulations. Our results point out that application of common descriptors of molecular shape such as Solvent Accessible Surface (SAS) area can lead to misleading outcomes. Instead, more appropriate molecular descriptors can be used to define 3D structures. In particular, we exploited Weighted Holistic Invariant Molecular (WHIM) descriptors to get a coarse-grained but accurate definition of the variegated Sic1 KID conformational ensemble. We found that Sic1 is able to form a variable amount of folded structures even in absence of partners. Among them, there were some conformations very close to the structure that Sic1 is supposed to assume in the binding with its physiological complexes. Therefore, our results support the hypothesis that this protein relies on the conformational selection mechanism to recognize the correct molecular partners. [ABSTRACT FROM AUTHOR]

This work deals with molecular dynamics analysis of properties of systems composed of carbon nanotubes and short telomeric DNA strands able to fold into i-motif structures at slightly acidic pH conditions. The studies are focused on possible application of such constructs as pH-controlled drug delivery and release systems. We study two different approaches. The first assumes that folding/unfolding property of these DNA strands might realize a gate closing/opening mechanism with carbon nanotube as a container for drug molecules. The second approach assumes that these DNA strands can modulate the drug intercalating property as a function of pH. As a model drug molecule we used doxorubicin. We found that the first approach is impossible to realize because doxorubicin is not effectively locked in the nanotube interior by DNA oligonuceotides. The second approach is more promising though direct drug release was not observed in unbiased molecular dynamics simulations. However, by applying detailed analysis of pair interaction energies, mobilities and potential of mean force we can show that doxorubicin can be released when the DNA strands fold into i-motifs. Carbon nanotube in that latter case acts mainly as a carrier for active phase which is composed of DNA fragments able to fold into noncanonical tetraplexes (i-motif). [ABSTRACT FROM AUTHOR]

The thermal conductivities of untreated polyisoprene and polybutadiene were calculated by molecular dynamics (MD) simulations using a Green-Kubo approach between −10 °C and 50 °C at atmospheric pressure. For comparison, the thermal conductivities of untreated polyisoprene with a molecular weight of 54,000 g/mol and untreated polybutadiene with a molecular weight of 45,000 g/mol were measured by the transient hot wire method in similar conditions. The simulation results of both polymers are in good agreement with the experimental data. We observed that the MD simulations slightly overestimate the thermal conductivity due to the chosen force field description. Details are discussed in the paper. [ABSTRACT FROM AUTHOR]

Boric acid, borate esters, and hydroxy derivatives are biologically active molecules. Thus, performing molecular dynamics simulations of these molecules is vital in terms of drug design, but it is difficult to find directly generated Amber parameters based on an ab initio method for these kinds of molecules in the literature. In this study, Amber parameters for such molecules containing boron were generated based on ab initio calculations using the paramfit program, which applies a combination of genetic and simplex algorithms, and the Visual Force Field Derivation Toolkit (VFFDT) program containing the Seminario method. The minimized structure, after obtaining novel parameters and using the sander program, was compared with the experimental crystallographic structures, and it was observed that the root-mean-square deviation (RMSD) value between the experimental structure and minimized structure agreed reasonably well. In addition, the molecule was heated, and the molecular dynamics simulation was successfully obtained with the novel parameters. [ABSTRACT FROM AUTHOR]

Glycosaminoglycans (GAGs) are linear, structurally diverse, conformationally complex carbohydrate polymers that may contain up to 200 monosaccharides. These characteristics present a challenge for studying GAG conformational thermodynamics at atomic resolution using existing experimental methods. Molecular dynamics (MD) simulations can overcome this challenge but are only feasible for short GAG polymers. To address this problem, we developed an algorithm that applies all conformational parameters contributing to GAG backbone flexibility (i.e., bond lengths, bond angles, and dihedral angles) from unbiased all-atom explicit-solvent MD simulations of short GAG polymers to rapidly construct models of GAGs of arbitrary length. The algorithm was used to generate non-sulfated chondroitin 10- and 20-mer ensembles which were compared to MD-generated ensembles for internal validation. End-to-end distance distributions in constructed and MD-generated ensembles have minimal differences, suggesting that our algorithm produces conformational ensembles that mimic the backbone flexibility seen in simulation. Non-sulfated chondroitin 100- and 200-mer ensembles were constructed within a day, demonstrating the efficiency of the algorithm and reduction in time and computational cost compared to simulation. [ABSTRACT FROM AUTHOR]

Our recently developed optimized potentials for liquid simulations-virtual site ionic liquid (OPLS-VSIL) force field has been shown to provide accurate bulk phase properties and local ion-ion interactions for a wide variety of imidazolium-based ionic liquids. The force field features a virtual site that offloads negative charge to inside the plane of the ring with careful attention given to hydrogen bonding interactions. In this study, the Diels-Alder reaction between cyclopentadiene and methyl acrylate was computationally investigated in the ionic liquid 1-butyl-3-methylimidazolium hexafluorophosphate, [BMIM][PF6], as a basis for the validation of the OPLS-VSIL to properly reproduce a reaction medium environment. Mixed ab initio quantum mechanics and molecular mechanics (QM/MM) calculations coupled to free energy perturbation and Monte Carlo sampling (FEP/MC) that utilized M06-2X/6-31G(d) and OPLS-VSIL gave activation free energy barriers of 14.9 and 16.0 kcal/mol for the endo-cis and exo-cis Diels-Alder reaction pathways, respectively (exptl. ΔH‡ of 14.6 kcal/mol). The endo selectivity trend was correctly predicted with a calculated 73% endo preference. The rate and selectivity enhancements present in the endo conformation were found to arise from preferential hydrogen bonding with the exposed C4 ring hydrogen on the BMIM cation. Weaker electronic stabilization of the exo transition state was predicted. For comparison, our earlier ±0.8 charge-scaled OPLS-2009IL force field also yielded a ΔG‡ of 14.9 kcal/mol for the favorable endo reaction pathway but did not adequately capture the highly organized solvent interactions present between the cation and Diels-Alder transition state. [ABSTRACT FROM AUTHOR]

The solenoid is a highly ordered structure observed in proteins, characterized by a set of symmetries. A group of enzymes—lyases containing solenoid fragments—was subjected to analysis with focus on their distribution of hydrophobicity/hydrophilicity, applying the fuzzy oil drop model. The model differentiates between a monocentric distribution hydrophobic core (spherical symmetry—mathematically modeled by a 3D Gaussian) and linear propagation of hydrophobicity (symmetry based on translation of structural units, i.e., chains—evident in amyloids). The linearly ordered solenoid carries information that affects the structure of the aqueous solvent in its neighborhood. Progressive disruption of its symmetry (via incorporation of asymmetrical fragments of varying size) appears to facilitate selective interaction with the intended substrate during enzymatic catalysis. [ABSTRACT FROM AUTHOR]

Protein structure is the result of the high synergy of all amino acids present in the protein. This synergy is the result of an overall strategy for adapting a specific protein structure. It is a compromise between two trends: The optimization of non-binding interactions and the directing of the folding process by an external force field, whose source is the water environment. The geometric parameters of the structural form of the polypeptide chain in the form of a local radius of curvature that is dependent on the orientation of adjacent peptide bond planes (result of the respective Phi and Psi rotation) allow for a comparative analysis of protein structures. Certain levels of their geometry are the criteria for comparison. In particular, they can be used to assess the differences between the structural form of biologically active proteins and their amyloid forms. On the other hand, the application of the fuzzy oil drop model allows the assessment of the role of amino acids in the construction of tertiary structure through their participation in the construction of a hydrophobic core. The combination of these two models—the geometric structure of the backbone and the determining of the participation in the construction of the tertiary structure that is applied for the comparative analysis of biologically active and amyloid forms—is presented. [ABSTRACT FROM AUTHOR]