Prediction of nonlocal elasticity parameters using high-throughput molecular dynamics simulations and machine learning

Structure-based virtual screening relies on classical scoring functions that often fail to reliably discriminate binders from nonbinders. In this work, we present a high-throughput protein-ligand complex molecular dynamics (MD) simulation that uses the output from AutoDock Vina to improve docking results in distinguishing active from decoy ligands in a directory of useful decoy-enhanced (DUD-E) dataset. MD trajectories are processed by evaluating ligand-binding stability using root-mean-square deviations. We select 56 protein targets (of 7 different protein classes) and 560 ligands (280 actives, 280 decoys) and show 22% improvement in ROC AUC (area under the curve, receiver operating characteristics curve), from an initial value of 0.68 (AutoDock Vina) to a final value of 0.83. The MD simulation demonstrates a robust performance across all seven different protein classes. In addition, some predicted ligand-binding modes are moderately refined during MD simulations. These results systematically validate the reliability of a physics-based approach to evaluate protein-ligand binding interactions.

Prediction of material behavior using machine learning (ML) requires consistent, accurate, and, representative large data for training. However, such consistent and reliable experimental datasets are not always available for materials. To address this challenge, we synergistically integrate ML with high-throughput reactive molecular dynamics (MD) simulations to elucidate the constitutive relationship of calcium–silicate–hydrate (C–S–H) gel—the primary binding phase in concrete formed via the hydration of ordinary portland cement. Specifically, a highly consistent dataset on the nine elastic constants of more than 300 compositions of C–S–H gel is developed using high-throughput reactive simulations. From a comparative analysis of various ML algorithms including neural networks (NN) and Gaussian process (GP), we observe that NN provides excellent predictions. To interpret the predicted results from NN, we employ SHapley Additive exPlanations (SHAP), which reveals that the influence of silicate network on all the elastic constants of C–S–H is significantly higher than that of water and CaO content. Additionally, the water content is found to have a more prominent influence on the shear components than the normal components along the direction of the interlayer spaces within C–S–H. This result suggests that the in-plane elastic response is controlled by water molecules whereas the transverse response is mainly governed by the silicate network. Overall, by seamlessly integrating MD simulations with ML, this paper can be used as a starting point toward accelerated optimization of C–S–H nanostructures to design efficient cementitious binders with targeted properties.

Machine learning (ML) force fields are revolutionizing molecular dynamics (MD) simulations as they bypass the computational cost associated with ab initio methods but do not sacrifice accuracy in the process. In this work, the GPyTorch library is used to create Gaussian process regression (GPR) models that are interfaced with the next-generation ML force field FFLUX. These models predict atomic properties of different molecular configurations that appear in a progressing MD simulation. An improved kernel function is utilized to correctly capture the periodicity of the input descriptors. The first FFLUX molecular simulations of ammonia, methanol, and malondialdehyde with the updated kernel are performed. Geometry optimizations with the GPR models result in highly accurate final structures with a maximum root-mean-squared deviation of 0.064 Å and sub-kJ mol-1 total energy predictions. Additionally, the models are tested in 298 K MD simulations with FFLUX to benchmark for robustness. The resulting energy and force predictions throughout the simulation are in excellent agreement with ab initio data for ammonia and methanol but decrease in quality for malondialdehyde due to the increased system complexity. GPR model improvements are discussed, which will ensure the future scalability to larger systems.

Calcium-dependent lectin I from Pseudomonas aeruginosa (PA-IL) binds specifically to oligosaccharides presenting an alpha-galactose residue at their nonreducing end, such as the disaccharides alphaGal1-2betaGalOMe, alphaGal1-3betaGalOMe, and alphaGal1-4betaGalOMe. This provides a unique model for studying the effect of the glycosidic linkage of the ligands on structure and thermodynamics of the complexes by means of experimental and theoretical tools. The structural features of PA-IL in complex with the three disaccharides were established by docking and molecular dynamics simulations and compared with those observed in available crystal structures, including PA-IL.alphaGal1-2betaGalOMe complex, which was solved at 2.4 A resolution and reported herein. The role of a structural bridge water molecule in the binding site of PA-IL was also elucidated through molecular dynamics simulations and free energy calculations. This water molecule establishes three very stable hydrogen bonds with O6 of nonreducing galactose, oxygen from Pro-51 main chain, and nitrogen from Gln-53 main chain of the lectin binding site. Binding free energies for PA-IL in complex with the three disaccharides were investigated, and the results were compared with the experimental data determined by titration microcalorimetry. When the bridge water molecule was included in the free energy calculations, the simulations predicted the correct binding affinity trends with the 1-2-linked disaccharide presenting three times stronger affinity ligand than the other two. These results highlight the role of the water molecule in the binding site of PA-IL and indicate that it should be taken into account when designing glycoderivatives active against P. aeruginosa adhesion.

Although continuum theories have been proven quite robust to describe confined fluid flow at molecular length scales, molecular dynamics (MD) simulations reveal mechanistic insights into the interfacial dissipation processes. Most MD simulations of confined fluids have used setups in which the lateral box size is not much larger than the gap height, thus breaking thin-film assumptions usually employed in continuum simulations. Here we explicitly probe the long-wavelength hydrodynamic correlations in confined simple fluids with MD and compare to gap-averaged continuum theories as typically applied in, e.g., lubrication. Relaxation times obtained from equilibrium fluctuations interpolate between the theoretical limits from bulk hydrodynamics and continuum formulations with increasing wavelength. We show how to exploit this characteristic transition to measure viscosity and slip length in confined systems simultaneously from equilibrium MD simulations. Moreover, the gap-averaged theory describes a geometry-induced dispersion relation that leads to overdamped sound relaxation at large wavelengths, which is confirmed by our MD simulations. Our results add to the understanding of transport processes under strong confinement and might be of technological relevance for the design of nanofluidic devices. Published by the American Physical Society 2024

In present study, vibration analysis of functionally graded (FG) porous microbeam under thermal effects is investigated considering modified strain gradient theory (MSGT) of elasticity. As the main novelty of this study, in order to effectively acquire the effects of size of size-dependent porous structures, MSGT is employed (considers three material length scale parameters) rather than modified couple stress theory (MCST) (considers one material length scale parameter). The nonlinear governing equation of microbeam based on MSGT is derived from Hamilton’s principle and to determine the natural frequency of the system under simply supports, Navier solution is employed. Numerical results presented for different beam models and theories and are compared with those available in previous studies. The obtained results show that frequencies from MSGT are higher than those from MCST and also classical theory (CT), especially when the microbeam thickness is comparable to the microbeam length scale parameter. Also, increases in the temperature of microbeam working environment, lead to decrease of microbeam natural frequency. Finally, parametric study is performed on natural frequency of FG microbeams to indicate the effects of: temperature changes, length scale parameter, slender ratio, and gradient index and porosity volume fraction.

Prediction and screening of glass properties based on high-throughput molecular dynamics simulations and machine learning

The term molecular mechanics (MM) refers to the use of simple potential-energy functions (e.g., harmonic oscillator or Coulombic potentials) to model molecular systems. Molecular mechanics approaches are widely applied in molecular structure refinement, molecular dynamics (MD) simulations, Monte Carlo (MC) simulations, and ligand-docking simulations. Typically, molecular mechanics models consist of spherical atoms connected by springs which represent bonds. Internal forces experienced in the model structure are described using simple mathematical functions. For example, Hooke’s law is commonly used to describe bonded interactions, and the nonbonded atoms might be treated as inelastic hard spheres or may interact according to a Lennard-Jones potential. Using these simple models, a molecular dynamics simulation numerically solves Newton’s equations of motion, thus allowing structural fluctuations to be observed with respect to time. Dynamic simulation methods are widely used to obtain information on the time evolution of conformations of proteins and other biological macromolecules1–4 and also kinetic and thermodynamic information. Simulations can provide fine detail concerning the motions of individual particles as a function of time. They can be utilized to quantify the properties of a system at a precision and on a time scale that is otherwise inaccessible, and simulation is, therefore, a valuable tool in extending our understanding of model systems. Theoretical consideration of a system additionally allows one to investigate the specific contributions to a property through “computational alchemy”,5 that is, modifying the simulation in a way that is nonphysical but nonetheless allows a model’s characteristics to be probed. One particular example is the artificial conversion of the energy function from that representing one system to that of another during a simulation. This is an important technique in free-energy calculations.6 Thus, molecular dynamics simulations, along with a range of complementary computational approaches, have become valuable tools for investigating the basis of protein structure and function. This review offers an outline of the origin of molecular dynamics simulation for protein systems and how it has developed into a robust and trusted tool. This review then covers more recent advances in theory and an illustrative selection of practical studies in which it played a central role. The range of studies in which MD has played a considerable or pivotal role is immense, and this review cannot do justice to them; MD simulations of biomedical importance were recently reviewed.4 Particular emphasis will be placed on the study of dynamic aspects of protein recognition, an area where molecular dynamics has scope to provide broad and far-ranging insights. This review concludes with a brief discussion of the future potential offered to advancement of the biological and biochemical sciences and the remaining issues that must be overcome to allow the full extent of this potential to be realized. 1.1. Historical Background MD methods were originally conceived within the theoretical physics community during the 1950s. In 1957, Alder and Wainwright7 performed the earliest MD simulation using the so-called hard-sphere model, in which the atoms interacted only through perfect collisions. Rahman8 subsequently applied a smooth, continuous potential to mimic real atomic interactions. During the 1970s, as computers became more widespread, MD simulations were developed for more complex systems, culminating in 1976 with the first simulation of a protein9,10 using an empirical energy function constructed using physics-based first-principles assumptions. MD simulations are now widely and routinely applied and especially popular in the fields of materials science11,12 and biophysics. As will be discussed later in this review, a variety of experimental conditions may be simulated with modern theories and algorithms. The initial simulations only considered single molecules in vacuo. Over time, more realistic or at least biologically relevant simulations could be performed. This trend is continuing today. The initial protein MD simulation, of the small bovine pancreatic trypsin inhibitor (BPTI), covered only 9.2 ps of simulation time. Modern simulations routinely have so-called equilibration periods much longer than that, and production simulations of tens of nanoseconds are routine, with the first microsecond MD simulation being reported in 1998.13 In addition, the original BPTI simulation included only about 500 atoms rather than the 104-106 atoms that are common today. While much of this advancement results from an immense increase in availability of computing power, major theoretical and methodological developments also contribute significantly. The number of publications regarding MD theory and application of MD to biological systems is growing at an extraordinary pace. A single review cannot do justice to the recent applications of MD. Using data from ISI Web of Science, the authors estimate that during 2005 at least 800 articles will be published that discuss molecular dynamics and proteins. The historical counts are shown in Figure 1. Open in a separate window Figure 1 Articles matching ISI Web of Science query “TS=(protein) AND TS=(molecular dynamics)”.

Both machine learning (ML) and molecular dynamics (MD) techniques have benefited substantially from the growing accessibility to computational resources, enabling more sophisticated ML models and more extensive MD simulations. Over recent years, the use of ML and MD methods has become increasingly prevalent in computer-aided drug discovery. To exploit the synergies between these methodologies, researchers have started to explore approaches that combine ML techniques and MD simulations. We start this chapter with a brief introduction of the basics of MD simulations as well as ML applications in drug discovery. Next, three broad ways of marrying MD with ML are discussed: (1) use of ML to parameterise systems for MD simulations with increased accuracy, (2) use of ML to improve sampling efficiency during MD simulations, and (3) use of ML to post-process MD trajectories in order to gain quantitative insights into molecular processes. We conclude the chapter with our perspective on areas for further development of this emerging interdisciplinary field.

In this study, the buckling behavior of functionalized single- and double-walled boron-nitride nanotubes (SWBNNTs and DWBNNTs) with a monomer of chitosan using molecular dynamics (MD) simulations is explored. The effect of chemical adsorption of chitosan molecule on the critical buckling force and strain is investigated. The results show that the critical buckling force considerably increases as the chitosan is attached to the side wall of boron-nitride nanotube which is more considerable for larger radii of nanotube. Moreover, increasing the number of walls reduces the sensitivity of boron-nitride nanotube to the functionalization compared with similar SWBNNTs. Further, it is shown that critical buckling of functionalized BNNTs increases by rising the weight percentage of chitosan. Considering the critical strain, it is observed that functionalization reduces the critical strain of functionalized BNNTs which is more pronounced in the case of SWBNNTs with bigger radii. Moreover, the buckling mode shape of functionalized BNNTs is presented.

The advent of Severe Acute Respiratory Syndrome coronavirus 2 (SARS-CoV-2), the etiological agent of the Coronavirus Disease 2019 (COVID-19) pandemic, has impacted physical and mental health worldwide. The lack of effective antiviral drugs necessitates a robust therapeutic approach to develop anti-SARS-CoV-2 drugs. Various investigations have recognized ACE2 as the primary receptor of SARS-CoV-2, and this amalgamation of ACE2 with the spike protein of the coronavirus is paramount for viral entry into the host cells and inducing infection. Consequently, restricting the virus's accessibility to ACE2 offers an alternative therapeutic approach to averting this illness. The study aimed to identify potent inhibitors with enhanced affinity for the ACE2 protein and validate their stability and efficacy against established inhibitors via molecular docking, machine learning, and MD simulations. 202 ACE2 inhibitors (PDB ID and 6LZG), comprising repurposed antiviral compounds and specific ACE2 inhibitors, were selected for molecular docking. The two most effective compounds obtained from docking were further analyzed using machine learning to identify potential compounds with enhanced ACE2-binding affinity. To refine the dataset, molecular decoys were generated through the Database of Useful Decoys: Enhanced (DUD-E) server, and Singular Value Decomposition (SVD) was applied for data preprocessing. The Tree-based Pipeline Optimization Tool (TPOT) was then utilized to optimize the machine learning pipeline. The most promising ML-predicted compounds were re-evaluated through docking and subjected to Molecular Dynamics (MD) simulations to evaluate their structural stability and interactions with ACE2. Finally, these compounds were evaluated against the top two pre-established inhibitors using various computational tools. The two best pre-established inhibitors were identified as Birinapant and Elbasvir, while the best machine-learning-predicted compounds were PubChem ID: 23658468 and PubChem ID: 117637105. Pharmacophore studies were conducted on the most effective machine-learning-predicted compounds, followed by a comparative ADME/T analysis between the best ML-screened and pre-established inhibitors. The results indicated that the top ML compound (PubChem ID: 23658468) demonstrated favorable BBB permeability and a high HIA index, highlighting its potential for therapeutic applications. The ML-screened ligand demonstrated structural stability with an RMSD (0.24 nm) and greater global stability (Rg: 2.08 nm) than Birinapant. Hydrogen bonding interactions further validated their strong binding affinity. MM/PBSA analysis confirmed the ML-screened compound's stronger binding affinity, with a binding free energy of - 132.90 kcal/mol, indicating enhanced stability in complex formation. The results emphasize the efficacy of integrating molecular docking, machine learning, and molecular dynamics simulations in facilitating the rapid identification of novel inhibitors. PubChem ID: 23658468 demonstrates robust binding affinity to ACE2 and favorable pharmacokinetic properties, establishing it as a promising candidate for further investigation.

Selectivity and transport in aquaporins from molecular simulation studies.

BackgroundArtificial sweeteners (AS) have been widely utilized in the food, beverage, and pharmaceutical industries for decades. While numerous publications have suggested a potential link between AS and diseases, particularly cancer, controversy still surrounds this issue. This study aims to investigate the association between AS consumption and cancer risk.MethodsTargets associated with commonly used AS were screened and validated using databases such as CTD, STITCH, Super-PRED, Swiss Target Prediction, SEA, PharmMapper, and GalaxySagittarius. Cancer-related targets were sourced from GeneCards, OMIM, and TTD databases. AS-cancer targets were identified through the intersection of these datasets. A network visualization (‘AS-targets-cancer’) was constructed using Cytoscape 3.9.0. Protein–protein interaction analysis was conducted using the STRING database to identify significant AS-cancer targets. GO and KEGG enrichment analyses were performed using the DAVID database. Core targets were identified from significant targets and genes involved in the ‘Pathways in cancer’ (map05200). Molecular docking and dynamics simulations were employed to verify interactions between AS and target proteins. Pan-cancer and univariate Cox regression analyses of core targets across 33 cancer types were conducted using GEPIA 2 and SangerBox, respectively. Gene chip datasets (GSE53757 for KIRC, GSE21354 for LGG, GSE42568 for BRCA, and GSE46602 for PRAD) were retrieved from the GEO database, while transcriptome and overall survival data were obtained from TCGA. Data normalization and identification of differentially expressed genes (DEGs) were performed on these datasets using R (version 4.3.2). Gene Set Enrichment Analysis (GSEA) was employed to identify critical pathways in the gene expression profiles between normal and cancer groups. A cancer risk prognostic model was constructed for key targets to further elucidate their significance in cancer initiation and progression. Finally, the HPA database was utilized to investigate variations in the expression of key AS-cancer target proteins across KIRC, LGG, BRCA, PRAD, and normal tissues.ResultsSeven commonly used AS (Aspartame, Acesulfame, Sucralose, NHDC, Cyclamate, Neotame, and Saccharin) were selected for study. A total of 368 AS-cancer intersection targets were identified, with 48 notable AS-cancer targets, including TP53, EGFR, SRC, PIK3R1, and EP300, retrieved. GO biological process analysis indicated that these targets are involved in the regulation of apoptosis, gene expression, and cell proliferation. Thirty-five core targets were identified from the intersection of the 48 significant AS-cancer targets and genes in the 'Pathways in cancer' (map05200). KEGG enrichment analysis of these core targets revealed associations with several cancer types and the PI3K-Akt signaling pathway. Molecular docking and dynamics simulations confirmed interactions between AS and these core targets. HSP90AA1 was found to be highly expressed across the 33 cancer types, while EGF showed the opposite trend. Univariate Cox regression analysis demonstrated strong associations of core targets with KIRC, LGG, BRCA, and PRAD. DEGs of AS-cancer core targets across these four cancers were analyzed. GSEA revealed upregulated and downregulated pathways enriched in KIRC, LGG, BRCA, and PRAD. Cancer risk prognostic models were constructed to elucidate the significant roles of key targets in cancer initiation and progression. Finally, the HPA database confirmed the crucial function of these targets in KIRC, LGG, BRCA, and PRAD.ConclusionThis study integrated data mining, machine learning, network toxicology, molecular docking, molecular dynamics simulations, and clinical sample analysis to demonstrate that AS increases the risk of kidney cancer, low-grade glioma, breast cancer, and prostate cancer through multiple targets and signaling pathways. This paper provides a valuable reference for the safety assessment and cancer risk evaluation of food additives. It urges food safety regulatory agencies to strengthen oversight and encourages the public to reduce consumption of foods and beverages containing artificial sweeteners and other additives.

Based on the Euler–Bernoulli beam model and the modified strain gradient theory, the size-dependent forced vibration of sandwich microbeams with a functionally graded (FG) core is presented. The equation of motion and the corresponding classical and nonclassical boundary conditions are derived using the Hamilton’s principle. An exact solution of the governing equation is developed for sandwich beams with various boundary conditions and subjected to an arbitrarily distributed harmonic transverse load. Finally, parametric studies are presented to investigate the effects of geometric ratios, length scale parameters, power index, boundary conditions, layup, and thickness of the FG layer on the frequency response of clamped and simply supported microbeams. Numerical results show that in the case of clamped microbeams, the essential and natural size-dependent boundary conditions have a significant effect on the resonance frequency and transverse deflection of microbeams. Also, it is seen that an optimal layup (without change in total volume of each material) can significantly improve the frequency characteristics of sandwich microbeams.

Introduction: Milk kefir, a fermenting milk made with kefir grains, has shown potential in promoting apoptosis, regulating the cell cycle, and reducing tumor growth in breast cancer cells. This study aimed to investigate the stability and potential of milk kefir metabolites as inhibitors of breast cancer growth by interacting with the estrogen receptor alpha (ER-α), a key protein involved in breast cancer cell proliferation. We used computational methods, specifically molecular docking simulations with AutoDock and molecular dynamics (MD) simulations with Gromacs, to analyze how these metabolites bind to ER-α. Methods: A combination of molecular docking and MD simulations was used to explore how metabolites derived from milk kefir interact with ER-α, a crucial target in breast cancer therapy. The methodology included multiple stages: preparation of target proteins, preparation and screening of the metabolites, geometry optimization, molecular docking, and MD simulations. Results: The molecular docking simulations of 43 metabolites revealed three promising candidates: 2-Methyl (S35), benzeneethanol (S42), and 2,6-dimethyl-4-heptanone (S54), with binding affinities (ΔG) of -5.08, -5.06, and -4.90 kcal/mol, respectively. MD simulations further showed that the selected metabolites stabilized the ER-α-metabolite complex, with the 2,6-dimethyl-4-heptanone (S54) metabolite demonstrating the most negative total MM-GBSA energy value (ΔG = -22.98 kcal/mol), indicating a strong and stable binding interaction. Conclusion: 2,6-Dimethyl-4-heptanone, a metabolite from milk kefir, showed promising potential as a candidate for further development as a breast cancer treatment, offering a novel alternative to conventional therapies.