Search

Computational molecular design—the endeavor to design molecules, with various missions, aided by machine learning and molecular dynamics approaches—has been widely applied to create valuable new molecular entities, from small molecule therapeutics to protein biologics. In the small data regime, physics-based approaches model the interaction between the molecule being designed and proteins of key physiological functions, providing structural insights into the mechanism. When abundant data have been collected, a quantitative structure–activity relationship can be more directly constructed from experimental data, from which machine learning can distill key insights to guide the design of the next round of experiment design. Machine learning methodologies can also facilitate physical modeling, from improving the accuracy of force fields and extending them to unseen chemical spaces to more directly enhancing the sampling on the conformational spaces. We argue that these techniques are mature enough to be applied to not just extend the longevity of life but the beauty it manifests. In this Perspective, we review the current frontiers in the research and development of skincare products, as well as the statistical and physical toolbox applicable to addressing the challenges in this industry. Feasible interdisciplinary research projects are proposed to harness the power of machine learning tools to design innovative, effective, and inexpensive skincare products. [ABSTRACT FROM AUTHOR]

High affinity is crucial for the efficacy and specificity of antibody. Due to involving high-throughput screens, biological experiments for antibody affinity maturation are time-consuming and have a low success rate. Precise computational-assisted antibody design promises to accelerate this process, but there is still a lack of effective computational methods capable of pinpointing beneficial mutations within the complementarity-determining region (CDR) of antibodies. Moreover, random mutations often lead to challenges in antibody expression and immunogenicity. In this study, to enhance the affinity of a human antibody against avian influenza virus, a CDR library was constructed and evolutionary information was acquired through sequence alignment to restrict the mutation positions and types. Concurrently, a statistical potential methodology was developed based on amino acid interactions between antibodies and antigens to calculate potential affinity-enhanced antibodies, which were further subjected to molecular dynamics simulations. Subsequently, experimental validation confirmed that a point mutation enhancing 2.5-fold affinity was obtained from 10 designs, resulting in the antibody affinity of 2 nM. A predictive model for antibody–antigen interactions based on the binding interface was also developed, achieving an Area Under the Curve (AUC) of 0.83 and a precision of 0.89 on the test set. Lastly, a novel approach involving combinations of affinity-enhancing mutations and an iterative mutation optimization scheme similar to the Monte Carlo method were proposed. This study presents computational methods that rapidly and accurately enhance antibody affinity, addressing issues related to antibody expression and immunogenicity. [ABSTRACT FROM AUTHOR]

DNA-PKcs is a crucial protein target involved in DNA repair and response pathways, with its abnormal activity closely associated with the occurrence and progression of various cancers. In this study, we employed a deep learning-based screening and molecular dynamics (MD) simulation-based pipeline, identifying eight candidates for DNA-PKcs targets. Subsequent experiments revealed the effective inhibition of DNA-PKcs-mediated cell proliferation by three small molecules (5025-0002, M769-1095, and V008-1080). These molecules exhibited anticancer activity with IC50 (inhibitory concentration at 50%) values of 152.6 μM, 30.71 μM, and 74.84 μM, respectively. Notably, V008-1080 enhanced homology-directed repair (HDR) mediated by CRISPR/Cas9 while inhibiting non-homologous end joining (NHEJ) efficiency. Further investigations into the structure-activity relationships unveiled the binding sites and critical interactions between these small molecules and DNA-PKcs. This is the first application of DeepBindGCN_RG in a real drug screening task, and the successful discovery of a novel DNA-PKcs inhibitor demonstrates its efficiency as a core component in the screening pipeline. Moreover, this study provides important insights for exploring novel anticancer therapeutics and advancing the development of gene editing techniques by targeting DNA-PKcs. [ABSTRACT FROM AUTHOR]

The interactions of the polar chemical bonds such as C=O and N–H with an external electric field were investigated, and a linear relationship between the QM/MM interaction energies and the electric field along the chemical bond is established in the range of weak to intermediate electrical fields. The linear relationship indicates that the electrostatic interactions of a polar group with its surroundings can be described by a simple model of a dipole with constant moment under the action of an electric field. This relationship is employed to develop a general approach to generating an electrostatic energy-based charge (EEC) model for molecules containing single or multiple polar chemical bonds. Benchmark test studies of this model were carried out for (CH3)2–CO and N-methyl acetamide in explicit water, and the result shows that the EEC model gives more accurate electrostatic energies than those given by the widely used charge model based on fitting to the electrostatic potential (ESP) in direct comparison to the energies computed by the QM/MM method. The MD simulations of the electric field at the active site of ketosteroid isomerase based on EEC demonstrated that EEC gave a better representation of the electrostatic interaction in the hydrogen-bonding environment than the Amber14SB force field by comparison with experiment. The current study suggests that EEC should be better suited for molecular dynamics study of molecular systems with polar chemical bonds such as biomolecules than the widely used ESP or RESP (restrained ESP) charge models. [ABSTRACT FROM AUTHOR]

Serine hydroxymethyltransferase 2 (SHMT2) has garnered significant attention as a critical catalytic regulator of the serine/glycine pathway in the one-carbon metabolism of cancer cells. Despite its potential as an anti-cancer target, only a limited number of inhibitors have been identified so far. In this study, we employed seven different scoring functions and skeleton clustering to screen the ChemDiv database for 38 compounds, 20 of which originate from the same skeleton structure. The most significant residues from SHMT2 and chemical groups from the inhibitors were identified using ASGBIE (Alanine Scanning with Generalized Born model and Interaction Entropy), and the binding energy of each residue was quantitatively determined, revealing the essential features of the protein–inhibitor interaction. The two most important contributing residues are TYR105 and TYR106 of the B chain followed by LEU166 and ARG425 of the A chain. The findings will be greatly helpful in developing a thorough comprehension of the binding mechanisms involved in drug–SHMT2 interactions and offer valuable direction for designing more potent inhibitors. [ABSTRACT FROM AUTHOR]

A practical approach to include the polarization effect in a molecular force field is the fluctuating charge method in which atomic charges vary as the configuration of the molecular system changes. However, the use of the Coulomb formula to evaluate energy in a fluctuating charge method is theoretically inconsistent with the forces given by the fluctuating method. In this work, we propose a force-consistent method to correctly calculate electrostatic energies of molecular systems using a fluctuating charge model (Effective Polarizable Bond or EPB). In this protocol, the electrostatic energy is obtained by numerical interaction of the atomic forces along the MD trajectory, rather than using the default Coulomb formula in the EPB model. Test study on the benchmark Barnase-Barstar protein-protein interaction system demonstrates that although the total electrostatic energy of the system shows little deviation due to the averaging effect, specific residue-residue electrostatic interaction energy is affected and the level of the effect depends on the charges of the interacting residues with charged residues showing pronounced differences in calculated energies between using the current protocol and the standard Coulomb formula. It is recommended that the proposed numerical interaction method should be preferred in the calculation of electrostatic energy in fluctuating charge models used in molecular dynamics simulations. [ABSTRACT FROM AUTHOR]

Human telomerase exhibits significant activity in cancer cells relative to normal cells, which contributes to the immortal proliferation of cancer cells. To counter this, the stabilization of G-quadruplexes formed in the guanine-rich sequence of the cancer cell chromosome has emerged as a promising avenue for anti-cancer therapy. Berberine (BER), an alkaloid that is derived from traditional Chinese medicines, has shown potential for stabilizing G-quadruplexes. To investigate the atomic interactions between G-quadruplexes and BER and its derivatives, molecular dynamics simulations were conducted. Modeling the interactions between G-quadruplexes and ligands accurately is challenging due to the strong negative charge of nucleic acids. Thus, various force fields and charge models for the G-quadruplex and ligands were tested to obtain precise simulation results. The binding energies were calculated by a combination of molecular mechanics/generalized Born surface area and interaction entropy methods, and the calculated results correlated well with experimental results. B-factor and hydrogen bond analyses demonstrated that the G-quadruplex was more stable in the presence of ligands than in the absence of ligands. Calculation of the binding free energy showed that the BER derivatives bind to a G-quadruplex with higher affinity than that of BER. The breakdown of the binding free energy to per-nucleotide energies suggested that the first G-tetrad played a primary role in binding. Additionally, energy and geometric properties analyses indicated that van der Waals interactions were the most favorable interactions between the derivatives and the G-quadruplexes. Overall, these findings provide crucial atomic-level insights into the binding of G-quadruplexes and their inhibitors. [ABSTRACT FROM AUTHOR]