Olanrewaju Ayodeji Durojaye, Sm Faysal Bellah, Henrietta Onyinye Uzoeto, Nkwachukwu Oziamara Okoro, Samuel Cosmas, Judith Nnedimkpa Ajima, Amarachukwu Vivian Arazu, Somtochukwu Precious Ezechukwu, Chiemekam Samuel Ezechukwu, Arome Solomon Odiba
{"title":"Harnessing AI-driven reverse docking in drug discovery: a comprehensive review of opportunities, challenges, and emerging trends","authors":"Olanrewaju Ayodeji Durojaye, Sm Faysal Bellah, Henrietta Onyinye Uzoeto, Nkwachukwu Oziamara Okoro, Samuel Cosmas, Judith Nnedimkpa Ajima, Amarachukwu Vivian Arazu, Somtochukwu Precious Ezechukwu, Chiemekam Samuel Ezechukwu, Arome Solomon Odiba","doi":"10.1007/s00894-025-06480-y","DOIUrl":"10.1007/s00894-025-06480-y","url":null,"abstract":"<p>The integration of artificial intelligence (AI) with reverse docking methodologies is reshaping drug discovery by streamlining the identification of drug targets and therapeutic interactions. This approach is pivotal in drug repurposing, safety profiling, and predicting off-target effects. Reverse docking uniquely identifies potential binding sites across diverse protein targets, providing insights into drug efficacy and adverse outcomes. AI technologies, such as machine learning, deep learning, and reinforcement learning, enhance this workflow by optimizing target selection, virtual screening, and conformational sampling. Despite challenges like data limitations and algorithmic complexities, AI-driven reverse docking has shown promise in drug repurposing and precision medicine, as illustrated by successful case studies. This review highlights its transformative potential and future prospects, including the incorporation of multi-omics data and real-time discovery pipelines for personalized medicine.</p><p>The computational strategies discussed leverage reverse docking platforms integrated with AI frameworks. Machine learning and deep learning models were employed for target selection and interaction prediction, while reinforcement learning facilitated advanced sampling techniques. Virtual screening workflows incorporated AI-driven optimizations for docking simulations. These methodologies were implemented using widely recognized computational tools, including AI libraries and molecular docking software, ensuring robust and reproducible results. Challenges in data integration were addressed by employing high-throughput pipelines capable of processing multi-omics datasets, thus supporting comprehensive drug discovery initiatives.\u0000</p>","PeriodicalId":651,"journal":{"name":"Journal of Molecular Modeling","volume":"31 9","pages":""},"PeriodicalIF":2.5,"publicationDate":"2025-08-25","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"144893920","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":4,"RegionCategory":"化学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Jianhua Wang, Shuya Wang, Jiaxuan Zhou, Hang Li, Yanwu Yu
{"title":"The influence of mixed solvents on the morphology of ε-CL-20 crystal","authors":"Jianhua Wang, Shuya Wang, Jiaxuan Zhou, Hang Li, Yanwu Yu","doi":"10.1007/s00894-025-06476-8","DOIUrl":"10.1007/s00894-025-06476-8","url":null,"abstract":"<div><h3>Context</h3><p>CL-20 (especially the ε-crystal form) is a high-energy density explosive, but the preparation of centimeter-sized ε-CL-20 crystals is complex, which limits the study of its properties and mechanisms. Since crystal morphology can significantly affect the performance and sensitivity of energetic materials, this study combines MD simulation and experimental techniques to study and examine the differences in ε-CL-20 crystal habit as well as the interplay between the principal growing crystal planes of ε-CL-20 and the mixed solvent molecules. The findings reveal that according to the AE model, ε-CL-20 consists of (0 1 1), (1 1 0), (1 0 -1), (1 1 -1), (0 0 2), (0 2 1), and (1 0 1) seven independent surfaces under vacuum, and the (0 1 1) surface is the most important growth surface. During the crystallization process of ε-CL-20 in three mixed solvents, the (0 1 1) surface remains an important crystal surface due to its low modified attachment energy. Also, the impact of the solvent component on the morphology of the crystal is more significant compared to that of the non-solvent component. The size of the crystal cultured in the experiment can reach centimeter level and is consistent with the simulated morphology. The X-ray diffraction patterns show that the synthesized crystals were in the <i>ε</i> form, and the sections are (0 1 1) surface. This study provides a novel approach for cultivating large ε-CL-20 crystals, paving the way for further studies on its performance and mechanism.</p><h3>Methods</h3><p>The ε-CL-20 crystal morphologies were obtained within a vacuum and different binary mixed solvent environments under COMPASS force field utilizing the AE model via the MD method on the Materials Studio 2019 platform. The complete simulation was generated using the NTV system. The temperature control method was chosen as Anderson, whereas the temperature of the system was set as 298 K. The data were acquired after every 5000 steps and the step size was 1 fs, while the total time of the simulation was equal to 100 ps.</p></div>","PeriodicalId":651,"journal":{"name":"Journal of Molecular Modeling","volume":"31 9","pages":""},"PeriodicalIF":2.5,"publicationDate":"2025-08-25","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"144893921","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":4,"RegionCategory":"化学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Peng Zhang, MengJie Bo, YuQin Chu, Yang Zhu, Peng Ma
{"title":"Electrostatic field modulation of high-nitrogen heterocycles: a computational study on energy and stability","authors":"Peng Zhang, MengJie Bo, YuQin Chu, Yang Zhu, Peng Ma","doi":"10.1007/s00894-025-06482-w","DOIUrl":"10.1007/s00894-025-06482-w","url":null,"abstract":"<div><h3>Context</h3><p>In this study, the effect of external electric field (EEF) on three kinds of high nitrogen heterocyclic compounds with different ring structures (<b>PA-1 ~ 3</b>, oxadiazole, triazole, and tetrazole) was investigated by quantum chemical calculation. By using B3LYP/6-311G + (d, p) theory and Laplace bond order (LBO) analysis, the trigger bond was determined, and further EEF was applied along the <i>X</i>/<i>Y</i>/<i>Z</i> axis (0–0.02 a.u.) to study its effects on molecular structure, electron distribution, and DOS. The results show that EEF in different directions has a significant effect on charge transfer and intermolecular interaction, which provides a theoretical basis for regulating energy output and thermal stability. This work highlights the potential of computational strategies in the design of energetic materials with controllable properties, and builds a bridge between theoretical prediction and practical application in high-energy systems.</p><h3>Methods</h3><p>The Gaussian 16 software has been chosen for simulation and computation in this study. The B3LYP functional and 6-311G + (d, p) basis set has been utilized for the calculation and simulation of external electric fields, with the strength ranging from 0 to 0.020 a.u. and an increment gradient of 0.005 a.u.</p></div>","PeriodicalId":651,"journal":{"name":"Journal of Molecular Modeling","volume":"31 9","pages":""},"PeriodicalIF":2.5,"publicationDate":"2025-08-25","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"144894058","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":4,"RegionCategory":"化学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Yang Zhu, Peng Zhang, YuQin Chu, Wen Jiang, Peng Ma, CongMing Ma
{"title":"Quantum chemical calculations based on 4,5-di(1,2,4-oxadiazol-3-yl)-2H-1,2,3-triazole and 3,3′ -(2H-1,2,3-triazole-4,5-diyl)-bis(1,2,4-oxadiazol-5(4H)-one) derivatives: a DFT study","authors":"Yang Zhu, Peng Zhang, YuQin Chu, Wen Jiang, Peng Ma, CongMing Ma","doi":"10.1007/s00894-025-06473-x","DOIUrl":"10.1007/s00894-025-06473-x","url":null,"abstract":"<div><h3>Context</h3><p>1,2,3-Triazole has excellent thermal stability, and both 1,2,3-triazole and 1,2,4-oxazolium are good in structural modifiability. This study employed density functional theory (DFT) to design 24 target compounds (designated <b>PA-1</b> to <b>PA-24</b>) based on the characteristics of these two heterocyclic systems. Using 4,5-diallyl-2H-1,2,3-triazole-3,3′-(2H-1,2,3-triazole-4,5-diyl)-bis(1,2,4-oxadiazol-5(4H)-one) as the parent compound (labeled <b>A1</b>), we introduced different energetic substituents at the C3 position of the oxadiazole ring and the N2 position of the triazole ring. Alternatively, using 3,3′-(2H-1,2,3-triazole-4,5-diyl)-bis(1,2,4-oxadiazol-5(4H)-one) as the parent compound (labeled <b>A2</b>), we introduced different energetic substituents at the N2 position of the triazole ring. The selection process considered both macroscopic properties (heat of formation and detonation performance) and microscopic factors including frontier molecular orbitals, energy gaps, weak interactions, surface electrostatic potential, and bond parameters. This approach aimed to identify energetic materials with both detonation performance and safety characteristics while investigating how introducing different substituents at these two positions affects detonation performance and safety. The research results indicate the following: First, introducing an azide group at the C3 position of the oxadiazole in <b>A1</b> more effectively enhances the compound’s heat of formation. <b>PA-13</b> (△<i>H</i><sub>f,sold</sub> = 972.97 kJ/mol), <b>PA-14</b> (△<i>H</i><sub>f,sold</sub> = 1023.62 kJ/mol), and <b>PA-15</b> (△<i>H</i><sub>f,sold</sub> = 1069.24 kJ/mol), which all contain azide groups, rank among the top three in heat of formation among the 24 compounds. Second, <b>PA-6</b> (<i>D</i> = 9.31 km/s, <i>P</i> = 39.93 GPa, <i>Q</i> = 7.39 kJ/g) exhibits the best detonation performance: its detonation pressure exceeds those of RDX (<i>G</i> = 36GPa), its detonation velocity and detonation heat surpass that of HMX (<i>D</i> = 9 km/s, <i>Q</i> = 6.4 kJ/g), and it achieves the ideal state of zero oxygen balance, resulting in the highest energy utilization efficiency. Third, <b>PA-19</b> showed the least impact sensitivity and showed the potential of an insensitive explosive. Fourth, hydrogen bonds are common in the molecular structure of <b>PA-8</b> ~ <b>PA-12</b> derivatives of <b>A2</b>, which may contribute to the stability of the molecules. Fifth, trinitromethyl and dinitromethyl groups are more prone to fragmentation. The substituents introduced at the C3 position of the oxadiazole have almost no effect on the ring length of the oxadiazole, but introducing trinitromethyl at the N2 position of the triazole causes the triazole ring structure to expand, resulting in an increased ring length.</p><h3>Method</h3><p>All calculations in this paper are based on density functional theory and were performed using the Gaussian16 software. Initi","PeriodicalId":651,"journal":{"name":"Journal of Molecular Modeling","volume":"31 9","pages":""},"PeriodicalIF":2.5,"publicationDate":"2025-08-21","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"144880944","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":4,"RegionCategory":"化学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
{"title":"Sensing attributes of ethylbenzene and methyl ethyl ketone vapours using novel β-arsenic nitride nanosheets based on first-principles study","authors":"R. Chandiramouli, A. Varshini, V. Nagarajan","doi":"10.1007/s00894-025-06449-x","DOIUrl":"10.1007/s00894-025-06449-x","url":null,"abstract":"<div><h3>Context</h3><p>In this research study, we employed a novel group VA-VA two-dimensional material β-arsenic nitride (β-AsN) nanosheet to explore the adsorption behavior of ethylbenzene and methyl ethyl ketone using the density functional theory (DFT) method. Initially, the structural stability of the β-AsN is validated by formation energy and phonon band spectrum. With the influence of band structure and projected density of states (PDOS) spectrum, we investigated the electronic characteristics of β-AsN monolayer. The computed energy gap value of β-AsN is 3.427 eV, which shows its semiconducting nature, and it can be utilized for numerous applications, viz., chemical sensors, bio-sensors, and optoelectronic devices. Using the most significant factors, namely relative band gap variation, Mulliken charge transfer, and adsorption energy, the adsorption behavior of ethylbenzene and methyl ethyl ketone on β-AsN is studied. The adsorption energy range is observed to be (− 0.134 eV to − 0.820 eV), which confirms that a weak van der Waals force acts between the base material and pollutants. The overall outcomes claimed that the β-AsN can be efficiently utilized for detecting both ethylbenzene and methyl ethyl ketone molecules in the common air environment.</p><h3>Methods</h3><p>The electronic and structural properties of β-AsN monolayer are calculated using the Quantum ATK package. We used a hybrid generalized gradient approximation (GGA) level of theory and Becke-3-Lee–Yang–Parr (B3LYP) exchange–correlation functional during the calculation. Also, methyl ethyl ketone and ethylbenzene adsorption on β-AsN monolayer is analyzed using the DFT calculations by employing the empirical dispersion correction of Grimme (DFT-D3) owing to weak van der Waals interactions.</p></div>","PeriodicalId":651,"journal":{"name":"Journal of Molecular Modeling","volume":"31 9","pages":""},"PeriodicalIF":2.5,"publicationDate":"2025-08-20","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"144880845","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":4,"RegionCategory":"化学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
F. E. Jorge, M. J. S. Matos, A. A. Vieira, T. Cazati, B. B. Postacchini
{"title":"Analysis and interpretation of experimental UV–Vis absorption spectra of benzochalcogenadiazoles and quinoxaline heterocycles through TDDFT","authors":"F. E. Jorge, M. J. S. Matos, A. A. Vieira, T. Cazati, B. B. Postacchini","doi":"10.1007/s00894-025-06460-2","DOIUrl":"10.1007/s00894-025-06460-2","url":null,"abstract":"<div><h3>Context</h3><p>Luminescent liquid crystals, particularly those derived from heterocyclic quinoxalines and benzochalcogenadiazoles, have garnered interest for their combined photophysical and mesomorphic properties. These compounds feature central heterocycles (selenium, oxygen, sulfur, or quinoxaline) connected to phenyl groups via triple bonds, with terminal alkoxy chains of uniform length. Their photophysical properties arise from conjugated systems involving delocalized electrons, while mesomorphic behavior depends on chain size and functional groups. This study performed a comparative analysis of experimental UV–vis absorption spectra against TDDFT computational calculations to clarify the photophysical properties of quinoxaline, 2,1,3-benzoxadiazole, 2,1,3-benzothiadiazole, and 2,1,3-benzoselenadiazole derivatives. Insights into the effects of scalar relativistic corrections and long-range functional adjustments on absorption energies are provided.</p><h3>Method</h3><p>Time-dependent density functional theory with the polarizable continuum model, available in the Gaussian 09 program, was used to simulate UV–Vis absorption spectra. Calculations utilized the CAM-B3LYP, B3LYP, and M062X functionals, along with basis sets incorporating scalar relativistic corrections via the Douglas-Kroll-Hess transformation. The impact of various functionals and the size of the basis set on the excitation energies was evaluated, especially for selenium-containing compound. The emission spectra in toluene were calculated with the ORCA 5.0.3 code. Experimental UV–Vis absorption and emission spectra were recorded in toluene for comparison with theoretical predictions to determine the most accurate computational approach.\u0000</p></div>","PeriodicalId":651,"journal":{"name":"Journal of Molecular Modeling","volume":"31 9","pages":""},"PeriodicalIF":2.5,"publicationDate":"2025-08-20","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"144868872","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":4,"RegionCategory":"化学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Giulio Poli, Salvatore Galati, Marco Macchia, Tiziano Tuccinardi
{"title":"ConPath 2.0: an optimized consensus strategy for assessing the potential pathogenicity of hRPE65 missense variants","authors":"Giulio Poli, Salvatore Galati, Marco Macchia, Tiziano Tuccinardi","doi":"10.1007/s00894-025-06481-x","DOIUrl":"10.1007/s00894-025-06481-x","url":null,"abstract":"<div><h3>Context</h3><p><i>h</i>RPE65 is an essential enzyme in the retinoid visual cycle. Numerous missense variants of <i>h</i>RPE65 have been linked to retinal disorders, such as retinitis pigmentosa and Leber congenital amaurosis. Moreover, many <i>h</i>RPE65 missense mutations are currently classified as variants of uncertain significance (VUS) due to insufficient evidence for a definitive pathogenicity classification. Addressing this limitation is critical for enabling accurate diagnoses and identifying suitable candidates for gene therapy. For this reason, we developed a <i>h</i>RPE65-tailored computational strategy, based on a consensus of multiple in silico pathogenicity predictions, enabling a rapid and reliable evaluation of the potential pathogenic effect of over 200 <i>h</i>RPE65 VUS. The analysis provided valuable insights to support the reclassification of these variants and assist clinicians in assessing their suitability for gene therapy.</p><h3>Methods</h3><p>In this study, we optimized our recently developed ConPath approach, which combines variant pathogenicity predictions from 19 different computational tools. In particular, we expanded the pool of predictive tools combined in the approach to 27, incorporated 3D-based methods that employ structural information for their prediction, and we refined the statistical framework for selecting the most reliable tools within an extended pool of more than 70 different methods. The tools were assessed based on their ability to discriminate pathogenic from benign <i>h</i>RPE65 missense mutations using an updated and expanded dataset of known <i>h</i>RPE65 variants. The resulting enhanced strategy, ConPath 2.0, was applied to the 210 <i>h</i>RPE65 VUS reported in the ClinVar database.</p></div>","PeriodicalId":651,"journal":{"name":"Journal of Molecular Modeling","volume":"31 9","pages":""},"PeriodicalIF":2.5,"publicationDate":"2025-08-20","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"144880846","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":4,"RegionCategory":"化学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
{"title":"Diffusion mechanism of dissolved gases in transformer mineral oil under furfural effect","authors":"Jia Tao, Xiongying Duan, Peng Yu, Xubin Li, Minfu Liao, Jinjin Li","doi":"10.1007/s00894-025-06475-9","DOIUrl":"10.1007/s00894-025-06475-9","url":null,"abstract":"<div><h3>Context</h3><p>Dissolved gas analysis (DGA) technology is a crucial technique for evaluating the operational state of transformers by detecting characteristic gases dissolved in insulating oil. Notably, the aging of insulation paper generates both characteristic gases and furfural, which also dissolves into the oil. However, the influence of furfural on gas diffusion behavior remains unclear. This study systematically investigates the effects of different furfural concentrations (0%, 1%, 3%, and 5% by mass) on the diffusion behavior of characteristic gases in insulating oil using molecular dynamics (MD) simulations. By calculating diffusion coefficients, free volume, and interaction energies, the underlying microscopic mechanisms were revealed. The results indicate that increasing furfural content slightly raises the free volume of gas molecules, thereby facilitating gas diffusion. Furthermore, the binding energies between furfural and insulating oil, as well as between furfural and gas molecules, increase with furfural concentration—mainly due to van der Waals interactions between polar and non-polar molecules. This study provides theoretical insight into the diffusion behavior of characteristic gases in insulating oil containing furfural.</p><h3>Methods</h3><p>The MD simulations were conducted using Materials Studio 2023. The amorphous cell module was employed to construct insulating oil mixture models incorporating varying concentrations of furfural molecules, utilizing the COMPASS III force field. Simulation results yielded key parameters including the mean square displacement (MSD) of characteristic gas molecules and the fractional free volume (FFV) of the system. Intermolecular interaction energies were computed using a self-developed Perl script. The simulation results demonstrate that increasing furfural concentration in insulating oil significantly enhances the diffusion characteristics of dissolved gases at the molecular level.</p></div>","PeriodicalId":651,"journal":{"name":"Journal of Molecular Modeling","volume":"31 9","pages":""},"PeriodicalIF":2.5,"publicationDate":"2025-08-18","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"144861534","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":4,"RegionCategory":"化学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
{"title":"Influence of SiO2 nanoparticles grafted with surfactant on coal wetting properties and microscopic mechanism of action","authors":"Xiangjun Chen, Jiahui Xu, Peiqi Zuo, Yibo Li, Liyang Li, San Zhao","doi":"10.1007/s00894-025-06466-w","DOIUrl":"10.1007/s00894-025-06466-w","url":null,"abstract":"<div><h3>Context</h3><p>In coal mine spray dust suppression, conventional studies on individual surfactants or surfactant mixtures often fail to achieve optimal coal wettability enhancement. This study innovatively introduces SiO₂ nanoparticles and demonstrates through simulations that the nanoparticle-surfactant composite system substantially enhances coal wettability. Electrostatic potential analysis indicates that the modified nanoparticle molecules formed by surfactant-grafted nanoparticles are more likely to form stable structures with coal molecules, enabling the nanoparticles to be better fixed on the coal molecules and establish a stable adsorption configuration. When the surfactant-grafted composite structure adsorbs onto the coal surface, it significantly enhances water molecular diffusivity, facilitates increased hydrogen bond formation with the composite surface, and substantially strengthens the interaction energy between coal and surfactants. Moreover, nanoparticle incorporation leads to noticeable thickening of the coal wetting layer and a marked increase in water molecule density at the interface. The results conclusively demonstrate that this composite structure enhances coal wettability through three synergistic mechanisms: interface optimization, hydrogen bond reinforcement, and water molecule diffusion enhancement.</p><h3>Methods</h3><p>This study innovatively employed molecular dynamics simulation methods to construct two adsorption systems—“water-surfactant-coal” and “water-nanoparticle-grafted surfactant-coal”—based on the Wiser coal molecular model. Through comprehensive analysis of key parameters including adsorption configurations, interaction energies, relative concentration profiles, hydrogen bond counts, water molecule diffusion behaviors, and detailed examination of surface electrostatic potentials of various components, we conducted comparative studies on the wetting behaviors of sodium dodecyl sulfate (SDS), alcohol ethoxylate (AEO<sub>9</sub>), and their nanoparticle-composite systems on coal. This investigation successfully revealed the microscopic mechanism by which composite structures enhance coal wettability at the molecular level.</p></div>","PeriodicalId":651,"journal":{"name":"Journal of Molecular Modeling","volume":"31 9","pages":""},"PeriodicalIF":2.5,"publicationDate":"2025-08-18","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"144861533","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":4,"RegionCategory":"化学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Nicole Rodgers, Christophe Lalaurie, Thomas Christopher Richard McDonnell
{"title":"Simulation of the trimeric globular head of C1q reveals temperature-sensitive network: implications for inflammation","authors":"Nicole Rodgers, Christophe Lalaurie, Thomas Christopher Richard McDonnell","doi":"10.1007/s00894-025-06464-y","DOIUrl":"10.1007/s00894-025-06464-y","url":null,"abstract":"<div><h3>Context</h3><p>C1q is an important protein in immune processes, driving complement activation through the classical pathway. Further to this, alterations in C1q either through SNPs or through autoantibodies can lead to systemic lupus erythematosus. Beyond these functions, C1q can also bind to other inflammatory proteins such as C-reactive protein (CRP) via its globular domain, when CRP is in the pentameric form. These interactions require specific structures to facilitate binding. Using molecular dynamics simulations, it is possible to measure the movements of proteins over time, with increasing temperatures allowing them to explore most of their available conformational space. Here, we describe using an increasing temperature simulation of C1q to identify potential structures generated during states of increased energy such as inflammation. Increasing temperature yielded significantly more movement of the monomeric and trimeric protein forms. Monomer A drove most movement within the molecule regardless of temperature, within the monomer and trimer. Further to this, novel structures were generated at higher temperatures, with significant movement of the CRP binding site. The altered movement in the CRP binding amino acids was correlative with increased temperature driving a loss of correlation between the different amino acids involved. Increased temperature and energy in the system leads to an alteration of C1q’s structure, which may leave it unable to bind to CRP in solution. This could have implications for the activity of the C1q/CRP complex as well as both proteins individually.</p><h3>Methods</h3><p>Models were generated using PDB:1PK6 and prepared using Charmm-GUI’s online platform. Protein simulations were run using NAMD on the UCL HPC facility (ARC). Trajectories were combined and aligned for analysis and visualised using Visual Molecular Dynamics (VMD). Analysis was carried out using VMD, R Studio, and Excel to identify novel structures of C1q, areas of increased flexibility, and potential protein networks.</p></div>","PeriodicalId":651,"journal":{"name":"Journal of Molecular Modeling","volume":"31 9","pages":""},"PeriodicalIF":2.5,"publicationDate":"2025-08-13","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"https://link.springer.com/content/pdf/10.1007/s00894-025-06464-y.pdf","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"144832224","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":4,"RegionCategory":"化学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"OA","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}