Short Molecular Dynamics Simulations Research Articles

Extracting binding energies and binding modes from biomolecular simulations of fragment binding to endothiapepsin.

Fragment-based drug discovery (FBDD) aims to discover a set of small binding fragments that may be subsequently linked together. Therefore, in-depth knowledge of the individual fragments' structural and energetic binding properties is essential. In addition to experimental techniques, the direct simulation of fragment binding by molecular dynamics (MD) simulations became popular to characterize fragment binding. However, former studies showed that long simulation times and high computational demands per fragment are needed, which limits applicability in FBDD. Here, we performed short, unbiased MD simulations of direct fragment binding to endothiapepsin, a well-characterized model system of pepsin-like aspartic proteases. To evaluate the strengths and limitations of short MD simulations for the structural and energetic characterization of fragment binding, we predicted the fragments' absolute free energies and binding poses based on the direct simulations of fragment binding and compared the predictions to experimental data. The predicted absolute free energies are in fair agreement with the experiment. Combining the MD data with binding mode predictions from molecular docking approaches helped to correctly identify the most promising fragments for further chemical optimization. Importantly, all computations and predictions were done within 5 days, suggesting that MD simulations may become a viable tool in FBDD projects.

PigmentHunter: A point-and-click application for automated chlorophyll-protein simulations.

Chlorophyll proteins (CPs) are the workhorses of biological photosynthesis, working together to absorb solar energy, transfer it to chemically active reaction centers, and control the charge-separation process that drives its storage as chemical energy. Yet predicting CP optical and electronic properties remains a serious challenge, driven by the computational difficulty of treating large, electronically coupled molecular pigments embedded in a dynamically structured protein environment. To address this challenge, we introduce here an analysis tool called PigmentHunter, which automates the process of preparing CP structures for molecular dynamics (MD), running short MD simulations on the nanoHUB.org science gateway, and then using electrostatic and steric analysis routines to predict optical absorption, fluorescence, and circular dichroism spectra within a Frenkel exciton model. Inter-pigment couplings are evaluated using point-dipole or transition-charge coupling models, while site energies can be estimated using both electrostatic and ring-deformation approaches. The package is built in a Jupyter Notebook environment, with a point-and-click interface that can be used either to manually prepare individual structures or to batch-process many structures at once. We illustrate PigmentHunter's capabilities with example simulations on spectral line shapes in the light harvesting 2 complex, site energies in the Fenna-Matthews-Olson protein, and ring deformation in photosystems I and II.

Using machine learning with atomistic surface and local water features to predict heterogeneous ice nucleation.

Heterogeneous ice nucleation (HIN) has applications in climate science, nanotechnology, and cryopreservation. Ice nucleation on the earth's surface or in the atmosphere usually occurs heterogeneously involving foreign substrates, known as ice nucleating particles (INPs). Experiments identify good INPs but lack sufficient microscopic resolution to answer the basic question: What makes a good INP? We employ molecular dynamics (MD) simulations in combination with machine learning (ML) to address this question. Often, the large amount of computational cost required to cross the nucleation barrier and observe HIN in MD simulations is a practical limitation. We use information obtained from short MD simulations of atomistic surface and water models to predict the likelihood of HIN. We consider 153 atomistic substrates with some surfaces differing in elemental composition and others only in terms of lattice parameters, surface morphology, or surface charges. A range of water features near the surface (local) are extracted from short MD simulations over a time interval (≤300 ns) where ice nucleation has not initiated. Three ML classification models, Random Forest (RF), support vector machine, and Gaussian process classification are considered, and the accuracies achieved by all three approaches lie within their statistical uncertainties. Including local water features is essential for accurate prediction. The accuracy of our best RF classification model obtained including both surface and local water features is 0.89 ± 0.05. A similar accuracy can be achieved including only local water features, suggesting that the important surface properties are largely captured by the local water features. Some important features identified by ML analysis are local icelike structures, water density and polarization profiles perpendicular to the surface, and the two-dimensional lattice match to ice. We expect that this work, with its strong focus on realistic surface models, will serve as a guide to the identification or design of substrates that can promote or discourage ice nucleation.

Molecular dynamics simulations as a guide for modulating small molecule aggregation

Small colloidally aggregating molecules (SCAMs) can be problematic for biological assays in drug discovery campaigns. However, the self-associating properties of SCAMs have potential applications in drug delivery and analytical biochemistry. Consequently, the ability to predict the aggregation propensity of a small organic molecule is of considerable interest. Chemoinformatics-based filters such as ChemAGG and Aggregator Advisor offer rapid assessment but are limited by the assay quality and structural diversity of their training set data. Complementary to these tools, we explore here the ability of molecular dynamics (MD) simulations as a physics-based method capable of predicting the aggregation propensity of diverse chemical structures. For a set of 32 molecules, using simulations of 100 ns in explicit solvent, we find a success rate of 97% (one molecule misclassified) as opposed to 75% by Aggregator Advisor and 72% by ChemAGG. These short timescale MD simulations are representative of longer microsecond trajectories and yield an informative spectrum of aggregation propensities across the set of solutes, capturing the dynamic behaviour of weakly aggregating compounds. Implicit solvent simulations using the generalized Born model were less successful in predicting aggregation propensity. MD simulations were also performed to explore structure-aggregation relationships for selected molecules, identifying chemical modifications that reversed the predicted behaviour of a given aggregator/non-aggregator compound. While lower throughput than rapid cheminformatics-based SCAM filters, MD-based prediction of aggregation has potential to be deployed on the scale of focused subsets of moderate size, and, depending on the target application, provide guidance on removing or optimizing a compound’s aggregation propensity.Graphical

More with less: Harvesting transition mechanism, free energy, and rates from short parallel molecular dynamics simulations

More with less: Harvesting transition mechanism, free energy, and rates from short parallel molecular dynamics simulations

Assessing the performance of MM/PBSA and MM/GBSA methods. 10. Prediction reliability of binding affinities and binding poses for RNA-ligand complexes.

Ribonucleic acid (RNA)-ligand interactions play a pivotal role in a wide spectrum of biological processes, ranging from protein biosynthesis to cellular reproduction. This recognition has prompted the broader acceptance of RNA as a viable candidate for drug targets. Delving into the atomic-scale understanding of RNA-ligand interactions holds paramount importance in unraveling intricate molecular mechanisms and further contributing to RNA-based drug discovery. Computational approaches, particularly molecular docking, offer an efficient way of predicting the interactions between RNA and small molecules. However, the accuracy and reliability of these predictions heavily depend on the performance of scoring functions (SFs). In contrast to the majority of SFs used in RNA-ligand docking, the end-point binding free energy calculation methods, such as molecular mechanics/generalized Born surface area (MM/GBSA) and molecular mechanics/Poisson Boltzmann surface area (MM/PBSA), stand as theoretically more rigorous approaches. Yet, the evaluation of their effectiveness in predicting both binding affinities and binding poses within RNA-ligand systems remains unexplored. This study first reported the performance of MM/PBSA and MM/GBSA with diverse solvation models, interior dielectric constants (εin) and force fields in the context of binding affinity prediction for 29 RNA-ligand complexes. MM/GBSA is based on short (5 ns) molecular dynamics (MD) simulations in an explicit solvent with the YIL force field; the GBGBn2 model with higher interior dielectric constant (εin = 12, 16 or 20) yields the best correlation (Rp = -0.513), which outperforms the best correlation (Rp = -0.317, rDock) offered by various docking programs. Then, the efficacy of MM/GBSA in identifying the near-native binding poses from the decoys was assessed based on 56 RNA-ligand complexes. However, it is evident that MM/GBSA has limitations in accurately predicting binding poses for RNA-ligand systems, particularly compared with notably proficient docking programs like rDock and PLANTS. The best top-1 success rate achieved by MM/GBSA rescoring is 39.3%, which falls below the best results given by docking programs (50%, PLNATS). This study represents the first evaluation of MM/PBSA and MM/GBSA for RNA-ligand systems and is expected to provide valuable insights into their successful application to RNA targets.

Critical Assessment of Self-Consistency Checks in the All-Atom Molecular Dynamics Simulation of Intrinsically Disordered Proteins.

All atom simulations can be used to quantify conformational properties of Intrinsically Disordered Proteins (IDP). However, simulations must satisfy convergence checks to ensure observables computed from simulation are reliable and reproducible. While absolute convergence is purely a theoretical concept requiring infinitely long simulation, a more practical, yet rigorous, approach is to impose Self Consistency Checks (SCCs) to gain confidence in the simulated data. Currently there is no study of SCCs in IDPs, unlike their folded counterparts. In this paper, we introduce different criteria for self-consistency checks for IDPs. Next, we impose these SCCs to critically assess the performance of different simulation protocols using the N terminal domain of HIV Integrase and the linker region of SARS-CoV-2 Nucleoprotein as two model IDPs. All simulation protocols begin with all-atom implicit solvent Monte Carlo (MC) simulation and subsequent clustering of MC generated conformations to create the representative structures of the IDPs. These representative structures serve as the initial structure for subsequent molecular dynamics (MD) runs with explicit solvent. We conclude that generating multiple short (∼3 μs) MD simulation trajectories─all starting from the most representative MC generated conformation─and merging them is the protocol of choice due to (i) its ability to satisfy multiple SCCs, (ii) consistently reproducing experimental data, and (iii) the efficiency of running independent trajectories in parallel by harnessing multiple cores available in modern GPU clusters. Running one long trajectory (greater than 20 μs) can also satisfy the first two criteria but is less desirable due to prohibitive computation time. These findings help resolve the challenge of identifying a usable starting configuration, provide an objective measure of SCC, and establish rigorous criteria to determine the minimum length (for one long simulation) or number of trajectories needed in all-atom simulation of IDPs.

Enhancing Conformational Sampling for Intrinsically Disordered and Ordered Proteins by Variational Autoencoder.

Intrinsically disordered proteins (IDPs) account for more than 50% of the human proteome and are closely associated with tumors, cardiovascular diseases, and neurodegeneration, which have no fixed three-dimensional structure under physiological conditions. Due to the characteristic of conformational diversity, conventional experimental methods of structural biology, such as NMR, X-ray diffraction, and CryoEM, are unable to capture conformational ensembles. Molecular dynamics (MD) simulation can sample the dynamic conformations at the atomic level, which has become an effective method for studying the structure and function of IDPs. However, the high computational cost prevents MD simulations from being widely used for IDPs conformational sampling. In recent years, significant progress has been made in artificial intelligence, which makes it possible to solve the conformational reconstruction problem of IDP with fewer computational resources. Here, based on short MD simulations of different IDPs systems, we use variational autoencoders (VAEs) to achieve the generative reconstruction of IDPs structures and include a wider range of sampled conformations from longer simulations. Compared with the generative autoencoder (AEs), VAEs add an inference layer between the encoder and decoder in the latent space, which can cover the conformational landscape of IDPs more comprehensively and achieve the effect of enhanced sampling. Through experimental verification, the Cα RMSD between VAE-generated and MD simulation sampling conformations in the 5 IDPs test systems was significantly lower than that of AE. The Spearman correlation coefficient on the structure was higher than that of AE. VAE can also achieve excellent performance regarding structured proteins. In summary, VAEs can be used to effectively sample protein structures.

Deciphering Selectivity Mechanism of BRD9 and TAF1(2) toward Inhibitors Based on Multiple Short Molecular Dynamics Simulations and MM-GBSA Calculations

BRD9 and TAF1(2) have been regarded as significant targets of drug design for clinically treating acute myeloid leukemia, malignancies, and inflammatory diseases. In this study, multiple short molecular dynamics simulations combined with the molecular mechanics generalized Born surface area method were employed to investigate the binding selectivity of three ligands, 67B, 67C, and 69G, to BRD9/TAF1(2) with IC50 values of 230/59 nM, 1400/46 nM, and 160/410 nM, respectively. The computed binding free energies from the MM-GBSA method displayed good correlations with that provided by the experimental data. The results indicate that the enthalpic contributions played a critical factor in the selectivity recognition of inhibitors toward BRD9 and TAF1(2), indicating that 67B and 67C could more favorably bind to TAF1(2) than BRD9, while 69G had better selectivity toward BRD9 over TAF1(2). In addition, the residue-based free energy decomposition approach was adopted to calculate the inhibitor-residue interaction spectrum, and the results determined the gatekeeper (Y106 in BRD9 and Y1589 in TAF1(2)) and lipophilic shelf (G43, F44, and F45 in BRD9 and W1526, P1527, and F1528 in TAF1(2)), which could be identified as hotspots for designing efficient selective inhibitors toward BRD9 and TAF1(2). This work is also expected to provide significant theoretical guidance and insightful molecular mechanisms for the rational designs of efficient selective inhibitors targeting BRD9 and TAF1(2).

Generation of unfolded outer membrane protein ensembles defined by hydrodynamic properties.

Outer membrane proteins (OMPs) must exist as an unfolded ensemble while interacting with a chaperone network in the periplasm of Gram-negative bacteria. Here, we developed a method to model unfolded OMP (uOMP) conformational ensembles using the experimental properties of two well-studied OMPs. The overall sizes and shapes of the unfolded ensembles in the absence of a denaturant were experimentally defined by measuring the sedimentation coefficient as a function of urea concentration. We used these data to model a full range of unfolded conformations by parameterizing a targeted coarse-grained simulation protocol. The ensemble members were further refined by short molecular dynamics simulations to reflect proper torsion angles. The final conformational ensembles have polymer properties different from unfolded soluble and intrinsically disordered proteins and reveal inherent differences in the unfolded states that necessitate further investigation. Building these uOMP ensembles advances the understanding of OMP biogenesis and provides essential information for interpreting structures of uOMP-chaperone complexes.

Relation in ligand dynamics and protein mutation using unsupervised machine learning.

Drug resistance mutation, which decreases ligand binding affinity to the target biomolecule, is a crucial issue all over the world. To overcome this issue, many approaches are proposed and tried. Free energy calculations are an effective and widely used method to estimate the ligand binding affinity. Relative free energy calculations can estimate the influence of the mutation on the ligand binding. However, the accuracy of the estimation should be limited partly because of insufficient sampling and relaxation of the mutants in the short molecular dynamics (MD) simulations. On the other hand, longer MD simulations make the calculation cost incredibly expensive. To improve the accuracy and computational cost, we have performed a machine-learning-assisted analysis of MD simulations [Yasuda, Commun. Biol. (2022)] of mutant protein-ligand systems. We used nilotinib-Abl-kinase systems, which were previously reported by other groups [Hauser, Commun. biol. (2018), Aldeghi, ACS Cent. Sci. (2019) ]. Our approach can compare the difference in the protein dynamics on ligand binding using unsupervised learning and extracts features of the dynamics of each mutant protein. We will show the relation of the feature to the ligand-binding affinity. Furthermore, interpretation of the feature enables us to understand the part of the ligand that shows different dynamics induced by the mutants. This would provide hints for the mechanics of the binding and improvements of the ligand.

Ensemble-based, high-throughput virtual screening of potential inhibitor targeting putative farnesol dehydrogenase of Metisa plana (Lepidoptera: Psychidae)

Metisa plana (Lepidoptera: Psychidae) bagworm is a leaf-eater caterpillar ubiquitously found as a damaging pest in oil palm plantations, specifically in Malaysia. Various strategies have been implemented, including the usage of chemical insecticides. However, the main challenges include the development of insecticide resistance and its detrimental effects on the environment and non-target organisms. Therefore, a biorational insecticide is introduced by targeting the juvenile hormone (JH) biosynthetic pathway, which is mainly present in the insect and vital for the insect's growth, diapause, metamorphosis, and adult reproduction. This study aimed to investigate the potential inhibitor for the rate-limiting enzyme involved in the JH pathway known as farnesol dehydrogenase. A 255 amino acids sequence encoded for the putative M. plana farnesol dehydrogenase (MpFolDH) open reading frame had been identified and isolated. The three-dimensional structure of MpFolDH was predicted to have seven β- sheets with α-helices at both sides, showing typical characteristics for classical short-chain dehydrogenase and associated with oxidoreductase activity. Then, the ensemble-based virtual screening was conducted based on the ZINC20 database, in which 43 768 compounds that fulfilled pesticide-likeness criteria were screened by site-specific molecular docking. After a short molecular dynamics simulation (5 ns) was conducted towards 102 compounds, only the top 10 compounds based on their most favourable binding energy were selected for a more extended simulation (100 ns). Based on the protein-ligand stability, protein compactness, residues rigidity, binding interaction, binding energy throughout the 100 ns simulation, and physicochemical analysis, ZINC000408743205 was selected as a potential inhibitor for this enzyme. Amino acids decomposition analysis indicates Ile18, Ala95, Val198 and Val202 were the critical contributor residues for MpFolDH-inhibitors(s) complex.

Sequence Control of the Self-Assembly of Elastin-Like Polypeptides into Hydrogels with Bespoke Viscoelastic and Structural Properties.

The biofabrication of structural proteins with controllable properties via amino acid sequence design is interesting for biomedicine and biotechnology, yet a complete framework that connects amino acid sequence to material properties is unavailable, despite great progress to establish design rules for synthesizing peptides and proteins with specific conformations (e.g., unfolded, helical, β-sheets, or β-turns) and intermolecular interactions (e.g., amphipathic peptides or hydrophobic domains). Molecular dynamics (MD) simulations can help in developing such a framework, but the lack of a standardized way of interpreting the outcome of these simulations hinders their predictive value for the design of de novo structural proteins. To address this, we developed a model that unambiguously classifies a library of de novo elastin-like polypeptides (ELPs) with varying numbers and locations of hydrophobic/hydrophilic and physical/chemical-cross-linking blocks according to their thermoresponsiveness at physiological temperature. Our approach does not require long simulation times or advanced sampling methods. Instead, we apply (un)supervised data analysis methods to a data set of molecular properties from relatively short MD simulations (150 ns). We also experimentally investigate hydrogels of those ELPs from the library predicted to be thermoresponsive, revealing several handles to tune their mechanical and structural properties: chain hydrophilicity/hydrophobicity or block distribution control the viscoelasticity and thermoresponsiveness, whereas ELP concentration defines the network permeability. Our findings provide an avenue to accelerate the design of de novo ELPs with bespoke phase behavior and material properties.

Binding Kinetics Toolkit for Analyzing Transient Molecular Conformations and Computing Free Energy Landscapes.

Understanding ligand binding kinetics and thermodynamics, which involves investigating the free, transient, and final complex conformations, is important in fundamental studies and applications for chemical and biomedical systems. Examining the important but transient ligand-protein-bound conformations, in addition to experimentally determined structures, also provides a more accurate estimation for drug efficacy and selectivity. Moreover, obtaining the entire picture of the free energy landscape during ligand binding/unbinding processes is critical in understanding binding mechanisms. Here, we present a Binding Kinetics Toolkit (BKiT) that includes several utilities to analyze trajectories and compute a free energy and kinetics profile. BKiT uses principal component space to generate approximated unbinding or conformational transition coordinates for accurately describing and easily visualizing the molecular motions. We implemented a new partitioning approach to assign indexes along the approximated coordinates that can be used as milestones or microstates. The program can generate input files to run many short classical molecular dynamics simulations and uses milestoning theory to construct the free energy profile and estimate binding residence time. We first validated the method with a host-guest system, aspirin unbinding from β-cyclodextrin, and then applied the protocol to pyrazolourea compounds and cyclin-dependent kinase 8 and cyclin C complexes, a kinase system of pharmacological interest. Overall, our approaches yielded good agreement with published results and suggest ligand design strategies. The computed unbinding free energy landscape also provides a more complete picture of ligand-receptor binding barriers and stable local minima for deepening our understanding of molecular recognition. BKiT is easy to use and has extensible features for future expansion of utilities for postanalysis and calculations.

Binding selectivity-dependent molecular mechanism of inhibitors towards CDK2 and CDK6 investigated by multiple short molecular dynamics and free energy landscapes

Understanding selectivity-dependent molecular mechanism of inhibitors towards CDK2 over CDK6 is prominent for improving drug design towards the CDK family. Multiple short molecular dynamics (MD) simulations combined with MM-GBSA approach are adopted to investigate molecular mechanism on binding selectivity of inhibitors X64, X3A, and 4 AU to CDK2 and CDK6. The RMSF analysis and calculations of molecular surface areas indicate that local structural and global flexibility of CDK6 are stronger than that of CDK2. Based on dynamics cross-correlation maps (DCCMs), motion modes of CDK2 and CDK6 produce difference due to associations of X64, X3A, and 4 AU. The calculated binding free energies (BFEs) demonstrate that the compensation between binding enthalpy and entropy of X64, X34, and 4 AU is a key force driving selectivity of inhibitors towards CDK2 over CDK6. This work provides valuable information for designing highly selective inhibitors towards CDK2 and CDK6 and further promotes identification of efficient anticancer drugs in the future.

Decoding drug resistant mechanism of V32I, I50V and I84V mutations of HIV-1 protease on amprenavir binding by using molecular dynamics simulations and MM-GBSA calculations

ABSTRACT Mutations V32I, I50V and I84V in the HIV-1 protease (PR) induce drug resistance towards drug amprenavir (APV). Multiple short molecular dynamics (MSMD) simulations and molecular mechanics generalized Born surface area (MM-GBSA) method were utilized to investigate drug-resistant mechanism of V32I, I50V and I84V towards APV. Dynamic information arising from MSMD simulations suggest that V32I, I50V and I84V highly affect structural flexibility, motion modes and conformational behaviours of two flaps in the PR. Binding free energies calculated by MM-GBSA method suggest that the decrease in binding enthalpy and the increase in binding entropy induced by mutations V32I, I50V and I84V are responsible for drug resistance of the mutated PRs on APV. The energetic contributions of separate residues on binding of APV to the PR show that V32I, I50V and I84V highly disturb the interactions of two flaps with APV and mostly drive the decrease in binding ability of APV to the PR. Thus, the conformational changes of two flaps in the PR caused by V32I, I50V and I84V play key roles in drug resistance of three mutated PR towards APV. This study can provide useful dynamics information for the design of potent inhibitors relieving drug resistance.

CHARMM‐GUIhigh‐throughput simulator for efficient evaluation of protein–ligand interactions with different force fields

AbstractMolecular docking is one of the most popular computational tools for the hit discovery step in drug design. However, there is ample room for improvement of docking's ability to identify correct binding modes and discriminate active from decoy compounds. Molecular dynamics (MD) simulations of protein–ligand docking structures have been shown to be effective in improving docking results. Here, we present CHARMM‐GUI high‐throughput simulator (HTS) that prepares MD simulation systems and inputs for multiple protein–ligand complex structures in a high‐throughput manner. HTS supports commonly used MD programs (NAMD, GROMACS, AMBER, OpenMM, GENESIS, Desmond, LAMMPS, and Tinker) along with various force field combinations for protein and ligand, including CHARMM36m, Amber (ff19SB/ff14SB), OPLS‐AA/M, CGenFF, GAFF2, and OpenFF. Validation tests using Miller and the directory of useful decoys‐enhanced (DUD‐E) datasets demonstrate that short MD simulations using HTS‐generated systems and simple ligand RMSD calculations consistently outperform docking results. Specifically, MD simulations can better identify correct ligand‐binding modes among top 10 binding poses as compared to docking scores. In addition, MD simulations can better discriminate active from decoy compounds in the DUD‐E dataset than docking scores for both soluble and membrane proteins. We expect that HTS can be a useful tool to facilitate the hit discovery process in drug design by improving docking results.

Protein Model Quality Estimation Using Molecular Dynamics Simulation.

The estimation of protein model quality remains a challenging task and is important for protein structural model utilization. In the last decade, existing methods that rely on machine learning to deep learning have been developed and shown progressive improvement. Despite utilizing more sophisticated techniques and introducing new features, none of these methods employ explicit protein structure stability information. Hypothetically, protein model quality might be indicated by its structural stability in an in silico system disclosed by the structural difference from its initial structure. One of the possible methods to exploit such information is by implementing molecular dynamics simulations that have shown successful applications in many research fields. We present a novel approach by introducing explicit protein structure stability information using molecular dynamics simulation. Despite using only simple features, small data with no training process required, and a short molecular dynamics simulation time, our method shows comparable performance to the state-of-the-art deep learning-based method.

Artificial intelligence guided conformational mining of intrinsically disordered proteins

Artificial intelligence recently achieved the breakthrough of predicting the three-dimensional structures of proteins. The next frontier is presented by intrinsically disordered proteins (IDPs), which, representing 30% to 50% of proteomes, readily access vast conformational space. Molecular dynamics (MD) simulations are promising in sampling IDP conformations, but only at extremely high computational cost. Here, we developed generative autoencoders that learn from short MD simulations and generate full conformational ensembles. An encoder represents IDP conformations as vectors in a reduced-dimensional latent space. The mean vector and covariance matrix of the training dataset are calculated to define a multivariate Gaussian distribution, from which vectors are sampled and fed to a decoder to generate new conformations. The ensembles of generated conformations cover those sampled by long MD simulations and are validated by small-angle X-ray scattering profile and NMR chemical shifts. This work illustrates the vast potential of artificial intelligence in conformational mining of IDPs.

Practical Protocols for Efficient Sampling of Kinase-Inhibitor Binding Pathways Using Two-Dimensional Replica-Exchange Molecular Dynamics.

Molecular dynamics (MD) simulations are increasingly used to study various biological processes such as protein folding, conformational changes, and ligand binding. These processes generally involve slow dynamics that occur on the millisecond or longer timescale, which are difficult to simulate by conventional atomistic MD. Recently, we applied a two-dimensional (2D) replica-exchange MD (REMD) method, which combines the generalized replica exchange with solute tempering (gREST) with the replica-exchange umbrella sampling (REUS) in kinase-inhibitor binding simulations, and successfully observed multiple ligand binding/unbinding events. To efficiently apply the gREST/REUS method to other kinase-inhibitor systems, we establish modified, practical protocols with non-trivial simulation parameter tuning. The current gREST/REUS simulation protocols are tested for three kinase-inhibitor systems: c-Src kinase with PP1, c-Src kinase with Dasatinib, and c-Abl kinase with Imatinib. We optimized the definition of kinase-ligand distance as a collective variable (CV), the solute temperatures in gREST, and replica distributions and umbrella forces in the REUS simulations. Also, the initial structures of each replica in the 2D replica space were prepared carefully by pulling each ligand from and toward the protein binding sites for keeping stable kinase conformations. These optimizations were carried out individually in multiple short MD simulations. The current gREST/REUS simulation protocol ensures good random walks in 2D replica spaces, which are required for enhanced sampling of inhibitor dynamics around a target kinase.