Journal of Chemical Theory and Computation最新文献

A significant challenge in applying machine learning to computational chemistry, particularly considering the growing complexity of contemporary machine learning models, is the scarcity of available experimental data. To address this issue, we introduce an approach that derives molecular features from an intricate neural network-based model and applies them to a simpler conventional machine learning model that is robust to overfitting. This method can be applied to predict various properties of a liquid system, including viscosity or surface tension, based on molecular features drawn from the ab initio calculated free energy of solvation. Furthermore, we propose a modified kernel model that includes Arrhenius temperature dependence to incorporate theoretical principles and diminish extreme nonlinearity in the model. The modified kernel model demonstrated significant improvements in certain scenarios and possible extensions to various theoretical concepts of molecular systems.

A plethora of key biological events occur at the cellular membrane where the large spatiotemporal scales necessitate dimensionality reduction or coarse-graining approaches over conventional all-atom molecular dynamics simulation. Constructing coarse-grained descriptions of membranes systematically from statistical mechanical principles has largely remained challenging due to the necessity of capturing amphipathic self-assembling behavior in coarse-grained models. We show that bottom-up coarse-grained lipid models can possess metastable morphological behavior and that this potential metastability has ramifications for accurate development and training. We in turn develop a training algorithm which evades metastability issues by linking model training to self-assembling behavior, and demonstrate its robustness via construction of solvent-free coarse-grained models of various phospholipid membranes, including lipid species such as phosphatidylcholines, phosphatidylserines, sphingolipids, and cholesterol. The resulting coarse-grained lipid models are orders of magnitude faster than their atomistic counterparts while retaining structural fidelity and constitute a promising direction for the development of coarse-grained models of realistic cell membranes.

A multiscale force field (FF) is developed for an aqueous solution of trivalent actinide cations An³⁺ (An = U, Np, Pu, Am, Cm, Bk, and Cf) by using a 12-6-4 Lennard-Jones type potential considering ion-induced dipole interaction. Potential parameters are rigorously and automatically optimized by the meta-multilinear interpolation parametrization (meta-MIP) algorithm via matching the experimental properties, including ion-oxygen distance (IOD) and coordination number (CN) in the first solvation shell and hydration free energy (HFE). The water solvent models incorporate an especially developed polar coarse-grained (CG) water scheme named PW32 and three widely used all-atom (AA) level SPC/E, TIP3P, and TIP4P water schemes. Each PW32 is modeled as two bonded beads to represent three neighboring water molecules, the simulation efficiency of which is 1 to 2 orders of magnitude higher than that of AA waters. The newly developed FF shows high accuracy and transferability in reproducing the IOD, CN, and HFE of An³⁺. The molecular structure and water exchange dynamics of the first An³⁺ hydration shell and the ionic (van der Waals) radii are reinvestigated in this work. Moreover, the new FF can readily be transferred to other popular FFs, as it has practicably predicted the permeability of An³⁺ in a graphene oxide filter within the framework of optimized potentials for liquid simulations (OPLS)-AA FF. It holds promise for applications in exploring actinide aqueous solutions with multiscale computational overhead.

Biomolecular condensates are essential in various cellular processes, and their misregulation has been demonstrated to underlie disease. Small molecules that modulate condensate stability and material properties offer promising therapeutic approaches, but mechanistic insights into their interactions with condensates remain largely lacking. We employ a multiscale approach to enable long-time, equilibrated all-atom simulations of various condensate-ligand systems. Systematic characterization of the ligand binding poses reveals that condensates can form diverse and heterogeneous chemical environments with one or multiple chains to bind small molecules. Unlike traditional protein-ligand interactions, these chemical environments are dominated by nonspecific hydrophobic interactions. Nevertheless, the chemical environments feature unique amino acid compositions and physicochemical properties that favor certain small molecules over others, resulting in varied ligand partitioning coefficients within condensates. Notably, different condensates share similar sets of chemical environments but at different populations. This population shift drives ligand selectivity toward specific condensates. Our approach can enhance the interpretation of experimental screening data and may assist in the rational design of small molecules targeting specific condensates.

Due to their efficient molecular design, nonfullerene acceptors (NFAs) have significantly advanced organic photovoltaics (OPVs). However, the lack of models to screen and evaluate candidate NFAs based on the resulting device performance has impeded the rapid development of high-performance molecules. This work introduces a computational framework utilizing a kinetic Monte Carlo (kMC) model to derive device parameters from molecular properties computed through first principles. By analyzing the quantum chemical properties of diverse dimeric conformers, we estimate the relative probabilities of microscopic processes such as charge separation, recombination, and transport along with charge transfer state formation in the active layer of OPVs. These probabilities set up a random walk of charge carriers in a grid with disordered molecular sites, allowing us to track their average behavior and calculate key device parameters. Our model consistently predicts measured device parameters, including the short-circuit current and open-circuit voltage, for OPVs with diverse NFAs with high accuracy. Additionally, we applied the model to evaluate donor-acceptor combinations of known compounds and newly designed NFA molecules, identifying high-performing pairs. This model offers a computationally efficient approach for designing novel NFA molecules and optimizing the OPV performance.

Compared to ground-state electronic structure optimizations, accurate simulations of molecular real-time electron dynamics are usually much more difficult to perform. To simulate electron dynamics, the time-dependent density matrix renormalization group (TDDMRG) has been shown to offer an attractive compromise between accuracy and cost. However, many simulation parameters significantly affect the quality and efficiency of a TDDMRG simulation. So far, it is unclear whether common wisdom from ground-state DMRG carries over to the TDDMRG, and a guideline on how to choose these parameters is missing. Here, in order to establish such a guideline, we investigate the convergence behavior of the main TDDMRG simulation parameters, such as time integrator, the choice of orbitals, and the choice of matrix-product-state representation for complex-valued nonsinglet states. In addition, we propose a method to select orbitals that are tailored to optimize the dynamics. Lastly, we showcase the TDDMRG by applying it to charge migration ionization dynamics in furfural, where we reveal a rapid conversion from an ionized state with a σ character to one with a π character within less than a femtosecond.

In this work the 4-point polarizable SWM4 Drude water model is reparametrized. Multiple models were developed using different strategies toward reproduction of specific target data. Results indicate that no individual model can reproduce all the selected target data in the context of the present form of the potential energy function. The changes considered in the new models include, 1) variations in the gas phase dipole moment, 2) variations in the molecular polarizability, 3) variations of the distance between the oxygen and the M site, 4) variation of the oxygen Lennard-Jones (LJ) parameters, 5) introduction of a LJ potential to the hydrogen atoms, and 6) variations of the H-O-H angle. Detailed analysis is presented for 3 new water models from which a final model, SWM4-HLJ, is selected as the future default model for the Drude polarizable force field. The model maintains the gas phase dipole moment as the experimental value while the remaining listed terms were adjusted including a larger H-O-H angle (108.12). Compared to its predecessor, SWM4-NDP, the self-diffusion coefficient, water dimer properties, and water cluster energies are greatly improved. The temperature dependence of the density of the new model also performs better. Overall, the new SWM4-HLJ water model is a general improvement and a good balance between microscopic and bulk properties is achieved.

Attaining a complete thermodynamic and kinetic characterization for processes involving multiple interconnected rare-event transitions remains a central challenge in molecular biophysics. This challenge is amplified when the process must be understood under a range of reaction conditions. Herein, we present a novel condition-responsive kinetic modeling framework that can combine the strengths of bottom-up rate quantification from multiscale simulations with top-down solution refinement using both equilibrium and nonequilibrium experimental data. Although this framework can be applied to any process, we demonstrate its use for electrochemically driven transport through channels and transporters via the development of electrochemically responsive rates. Using the Cl^-/H⁺ antiporter ClC-ec1 as a model system, we show how optimal and predictive kinetic solutions can be obtained when the solution space is grounded by thermodynamic constraints, seeded through multiscale rate quantification, and further refined with experimental data, such as electrophysiology assays. Turning to the Shaker K⁺ channel, we demonstrate that optimal solutions and biophysical insights can also be obtained with sufficient experimental data. This multi-pathway method also proves capable of identifying single-pathway dominant channel mechanisms but reveals that competing and off-pathway flux is still essential to replicate experimental findings and to describe concentration-dependent channel rectification.

We investigate the relaxation dynamics of crystalline stems in relation to the molecular topology of the crystalline/amorphous interface, employing coarse-grained molecular dynamics. To efficiently generate model semicrystalline systems of linear polyethylene with a realistic interphase morphology, we simplified the Monte Carlo method by introducing molecular dynamics for faster relaxation. The structural properties of the generated systems are validated against experimental measurements, theoretical predictions, and existing simulation data. The models suggest that the probability distribution of loop-entry sites on the lamellar surface can be described by a power law in terms of the distance between the entry sites. By considering realistic interphase morphology, we are able to improve the prediction of the overall activation energy for the relaxation of crystalline stems, aligning it closely with experimental measurements. The largest model predicts that crystalline stems connected via large loops, i.e., those that exceed the entanglement length, and long tails are associated with increased activation energy; whereas stems connected to shorter tails show the lowest activation energy. These predictions can guide the future development of tougher semicrystalline polymers by providing insights into how amorphous chain morphology contributes to the activation energy and the relaxation dynamics of crystalline chains.

Radical pairs (electron-hole pairs, polaron pairs) are transient reaction intermediates that are found and exploited in all areas of science, from the hard realm of physics in the form of organic semiconductors, spintronics, quantum computing, and solar cells to the soft domain of chemistry and biology under the guise of chemical reactions in solution, biomimetic systems, and quantum biology. Quantitative analysis of radical pair phenomena has historically been successful by a few select groups. With this in mind, we present an intuitive open-source framework in the Python programming language that provides classical, semiclassical, and quantum simulation methodologies. A radical pair kinetic rate equation solver, Monte Carlo-based spin dephasing rate estimations, and molecule database functionalities are implemented. We introduce the kine-quantum method, a new approach that amalgamates classical rate equations, semiclassical, and quantum techniques. This method resolves the prohibitively large memory requirement issues of quantum approaches while achieving higher accuracy, and it also offers wavelength-resolved simulations, producing time- and wavelength-resolved magnetic field effect simulations. Model examples illustrate the versatility and ease of use of the software, including the new approach applied to the magnetosensitive absorption and fluorescence of flavin adenine dinucleotide photochemistry, spin-spin interaction estimation from molecular dynamics simulations on radical pairs inside reverse micelles, radical pair anisotropy inside proteins, and triplet exciton pairs in anthracene crystals. The intuitive interface also allows this software to be used as a teaching or learning aid for those interested in the field of spin chemistry. Furthermore, the software aims to be modular and extensible, with the aim to standardize how spin dynamics simulations are performed.