Search

The phonon-related properties of crystalline polymers are highly relevant for various applications. Their simulation is, however, particularly challenging, as the systems that need to be modeled are often too extended to be treated by ab initio methods, while classical force fields are too inaccurate. Machine-learned potentials parametrized against material-specific ab initio data hold the promise of being extremely accurate and also highly efficient. Still, for their successful application, protocols for their parametrization need to be established to ensure an optimal performance, and the resulting potentials need to be thoroughly benchmarked. These tasks are tackled in the current manuscript, where we devise a protocol for parametrizing moment tensor potentials (MTPs) to describe the structural properties, phonon band structures, elastic constants, and forces in molecular dynamics simulations for three prototypical crystalline polymers: polyethylene (PE), polythiophene (PT), and poly-3-hexylthiophene (P3HT). For PE, the thermal conductivity and thermal expansion are also simulated and compared to experiments. A central element of the approach is to choose training data in view of the considered use case of the MTPs. This not only yields a massive speedup for complex calculations while essentially maintaining DFT accuracy, but also enables the reliable simulation of properties that, so far, have been entirely out of reach. [ABSTRACT FROM AUTHOR]

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]

Based on the TraPPE force field, previously unknown values of the parameters of the intermolecular interaction potential of trans-decalin were determined. Parametrization was carried out using experimental data on saturated vapor pressure and density at atmospheric pressure. The found parameters make it possible to adequately describe the boiling and condensation lines of trans-decalin and also predict the critical values of pressure, density, and temperature with satisfactory accuracy. Calculations of vapor-liquid phase equilibrium conditions for a binary CO2—trans-decalin mixture in supercritical conditions for CO2 were carried out. When quantitatively comparing the calculated values with experimental data, an underestimation of pressure at a temperature of 345.4 K by 30% is observed, which decreases to 5% for temperatures up to 525 K. [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 research focuses on obtaining the behavior of chemical compounds with respect to their molecular weight and optimization energy based on the variation in properties in organic carbon links. Here, behavioral analysis of compounds is used in the application of a metal organic framework to denote the high-grade compounds. The grade was selected based on the essential measure of optimization energy and molecular weight, and in turn, depicts the stability of material. Computation of the optimization energy and molecular weight of chemical compounds was performed with Avogadro software. Several force fields can be considered to compute optimized energy. Exclusively, three force fields, namely, the Universal Force Field (UFF), the General Amber Force Field (GAFF), and the Ghemical force field (Ghemical) were selected from Avogadro as these were more relevant to compounds considered in this research. The various chemical compounds examined in this work are Aluminum (Al), Boron (Br), Calcium (Ca), Chlorine (Cl), Indium (In), Potassium (K), Scandium (Sc), Silicon (Si), and Tungsten (W). Hence, molecular modeling of different compounds incorporated with three different force fields was evaluated in this work. In this study, we found that the In structure has more energy reduction, of 22.673 kJ mol−1 in UFF, when compared with the other two force fields. Thus, In has higher potential with more stability. [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]

The longevity of asphalt pavements is a key focus of road engineering, which closely relates to the self-healing ability of bitumen. Our work aims to establish a CGMD model and matched force field for bitumen and break through the limitations of the research scale to further explore the microscopic mechanism of bitumen self-healing. In this study, a CGMD mapping scheme containing 16 kinds of beads is proposed, and the non-bond potential energy function and bond potential energy function are calculated based on all-atom simulation to construct and validate a coarse-grained model for bitumen. On this basis, a micro-crack model with a width of 36.6nm is simulated, and the variation laws of potential energy, density, diffusion coefficient, relative concentration and temperature in the process of bitumen self-healing are analyzed with the cracking rate parameter proposed to characterize the degree of bitumen crack healing. The results show that the computational size of the coarse-grained simulation is much larger than that of the all-atom, which can explain the self-healing mechanism at the molecular level. In the self-healing process, non-bonded interactions dominate the molecular movement, and differences in the decreased rate of diffusion among the components indicate that saturates and aromatics play a major role in self-healing. Meanwhile, the variations in crack rates reveal that healing time is inversely proportional to temperature. The impact of increasing temperature on reducing healing time is most obvious when the temperature approaches the glass transition temperature (300 K). [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]

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]

Infrared (IR) spectroscopy is a valuable tool to obtain information about protein secondary structure. The far-infrared (FIR) spectrum is characterized by a complex combination of different molecular contributions which, for small molecules, may be interpreted with the help of quantum-mechanical (QM) calculations. Unfortunately, the high computational cost of QM calculations makes them inapplicable to larger molecules, such as proteins and peptides. In this work, we present a theoretical study on the secondary structure, molecular properties, and vibrational spectra of different peptides, using both a classical and a QM approach. Our results show that the amide I main peak value, and related quantities, such as dipole strength (DS) and transition dipole moment (TDM), depends on protein secondary structure; in particular, from QM calculations arises that α -rich molecular systems present lower intensities than β -rich ones. Furthermore, it is possible to decouple and identify the intensity of the different contributions of the inter- and intra-molecular motions which characterize the FIR spectrum, starting from the results obtained with QM calculations. [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]

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]

Molecular dynamics (MD) simulation is a particularly useful technique in food processing. Normally, food processing techniques can be optimized to favor the creation of higher-quality, safer, more functional, and more nutritionally valuable food products. Modeling food processes through the application of MD simulations, namely, the Groningen Machine for Chemical Simulations (GROMACS) software package, is helpful in achieving a better understanding of the structural changes occurring at the molecular level to the biomolecules present in food products during processing. MD simulations can be applied to define the optimal processing conditions required for a given food product to achieve a desired function or state. This review presents the development history of MD simulations, provides an in-depth explanation of the concept and mechanisms employed through the running of a GROMACS simulation, and outlines certain recent applications of GROMACS MD simulations in the food industry for the modeling of proteins in food products, including peanuts, hazelnuts, cow's milk, soybeans, egg whites, PSE chicken breast, and kiwifruit. [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]

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 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]

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]

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]