Adaptive steered molecular dynamics of biomolecules

Mechanistic Study of Protein Interaction with Natto Inhibitory Peptides Targeting Xanthine Oxidase: Insights from Machine Learning and Molecular Dynamics Simulations

Bioactive peptides from food sources offer a safe and biocompatible approach to enzyme inhibition, with potential applications in managing metabolic disorders such as hyperuricemia and gout, conditions linked to excessive xanthine oxidase activity. Using a machine learning-based screening approach inspired by the bioactivity of natto, two peptides, ECFK and FECK, were identified from the Bacillus subtilis proteome and validated as xanthine oxidase inhibitors with IC50 values of 37.36 and 71.57 mM, respectively. Further experiments confirmed their safety through cytotoxicity assays, and electronic tongue analysis demonstrated their mild sensory properties, supporting their edibility. Molecular dynamics simulations revealed that these peptides stabilize critical enzyme regions, with ECFK showing a higher dissociation energy barrier (52.08 kcal/mol) than FECK (46.39 kcal/mol), indicating strong, stable interactions. This study highlights food-derived peptides as safe and natural inhibitors of xanthine oxidase, offering promising therapeutic potential for metabolic disorder management.

A discard-and-restart MD algorithm for the sampling of protein intermediate states

We introduce a discard-and-restart molecular dynamics (MD) algorithm tailored for the sampling of realistic protein intermediate states. It aids computational structure-based drug discovery by reducing the simulation times to compute a "quick sketch" of folding pathways by up to 2000×. The algorithm iteratively performs short MD simulations and measures their proximity to a target state via a collective variable loss, which can be defined in a flexible fashion, locally or globally. Using the loss, if the trajectory proceeds toward the target, the MD simulation continues. Otherwise, it is discarded, and a new MD simulation is restarted, with new initial velocities randomly drawn from a Maxwell-Boltzmann distribution. The discard-and-restart algorithm demonstrates efficacy and atomistic accuracy in capturing the folding pathways in several contexts: 1) fast-folding small protein domains, 2) the folding intermediate of the prion protein PrP, and 3) the spontaneous partial unfolding of α-tubulin, a crucial event for microtubule severing. During each iteration of the algorithm, we can perform AI-based analysis of the transitory conformations to find potential binding pockets, which could represent druggable sites. Overall, our algorithm enables systematic and computationally efficient exploration of conformational landscapes, enhancing the design of ligands targeting dynamic protein states.

A study of alpha-synuclein and poly(N-isopropylacrylamide) complex formation through detailed atomistic simulations

This work presents an investigation of the influence of poly(N-isopropylacrylamide) (PNIPAM) polymer on the structural dynamics of intrinsically disordered alpha-synuclein (α-syn) protein, exploring the formation and intricate features of the resulting α-syn/PNIPAM complexes. Using atomistic molecular dynamics (MD) simulations, our study analyzes the impact of initial configuration, polymer molecular weight, and protein mutations on the α-syn and the α-syn/PNIPAM complex. Atomistic simulations, of a few μs, of the protein/polymer complex reveal crucial insights into molecular interactions within the complex, emphasizing a delicate balance of forces governing its stability and structural evolution. Our findings indicate that PNIPAM polymer engages in significant non-polar interactions with the non-amyloid component (NAC) region of α-syn, which plays a crucial role in fibril formation, under various conditions such as the mutations in the protein structure and polymer chain length. Especially the PNIPAM polymer with a 40mer monomer exhibits a stabilizing effect on the structural properties of the protein, reducing intramolecular interactions that contribute to misfolding. These findings, which delve into protein/polymer interactions, hold promise as potential guidance for therapeutic strategies in various neurodegenerative disorders.

Computational Insights into Membrane Disruption by Cell-Penetrating Peptides

Cell-penetrating peptides (CPPs) can translocate into cells without inducing cytotoxicity. The internalization process implies several steps at different time scales ranging from microseconds to minutes. We combine adaptive Steered Molecular Dynamics (aSMD) with conventional Molecular Dynamics (cMD) to observe nonequilibrium and equilibrium states to study the early mechanisms of peptide-bilayer interaction leading to CPPs internalization. We define three membrane compositions representing bilayer sections, neutral lipids (i.e., upper leaflet), neutral lipids with cholesterol (i.e., hydrophobic core), and neutral/negatively charged lipids with cholesterol (i.e., lower leaflet) to study the energy barriers and disruption mechanisms of Arg9, MAP, and TP2, representing cationic, amphiphilic, and hydrophobic CPPs, respectively. Cholesterol and negatively charged lipids increase the energetic barriers for the peptide-bilayer crossing. TP2 interacts with the bilayer by hydrophobic insertion, while Arg9 disrupts the bilayer by forming transient or stable pores. MAP has shown both behaviors. Collectively, these findings underscore the significance of innovative computational approaches in studying membrane-disruptive peptides and, more specifically, in harnessing their potential for cell penetration.

Elucidating Protein Structures in the Gas Phase: Traversing Configuration Space with Biasing Methods

Achieving accurate characterization of protein structures in the gas phase continues to be a formidable challenge. To tackle this issue, the present study employs Molecular Dynamics (MD) simulations in tandem with enhanced sampling techniques (methods designed to efficiently explore protein conformations). The objective is to identify suitable structures of proteins by contrasting their calculated Collision Cross-Section (CCS) with those observed experimentally. Significant discrepancies were observed between the initial MD-simulated and experimentally measured CCS values through Ion Mobility-Mass Spectrometry (IMS-MS). To bridge this gap, we employed two distinct enhanced sampling methods, Harmonic Biasing Potential and Adaptive Biasing Force, which help the proteins overcome energy barriers to adopt more compact configurations. These techniques leverage the radius of gyration as a reaction coordinate (guiding parameter), guiding the system toward compressed states that potentially match experimental configurations more closely. The guiding forces are only employed to overcome existing barriers and are removed to allow the protein to naturally arrive at a potential gas phase configuration. The results demonstrated close alignment (within ∼4%) between simulated and experimental CCS values despite using different strengths and/or methods, validating their efficacy. This work lays the groundwork for future studies aimed at optimizing biasing methods and expanding the collective variables used for more accurate gas-phase structural predictions.

Molecular insights and inhibitory dynamics of flavonoids in targeting Pim-1 kinase for cancer therapy

Pim-1 kinase, a serine/threonine kinase, is often overexpressed in various cancers, contributing to disease progression and poor prognosis. In this study, we explored the potential of flavonoids as inhibitors of Pim-1 kinase using a combination of molecular docking and steered molecular dynamics (SMD) simulations. Our docking studies revealed two main binding orientations for the flavonoid molecules. The SMD simulations showed that the binding mode with higher pulling forces was linked to stronger inhibitory activity, with a strong positive correlation (R 2 ≈ 0.92) between pulling forces and IC50 values. Quercetin stood out as the most potent inhibitor, showing a pulling force of about 820 pN and an IC_(5) 0 of less than 6 µM. Further dynamic simulations indicated that quercetin's hydroxyl groups at the C3, C-5 and C-7 positions formed stable hydrogen bonds with key residues GLU-121, Leu-44 and Val-126, respectively enhancing its binding stability and effectiveness. Our results emphasized the critical role of the hydroxyl group at the C-3 position, which plays a pivotal function in effectively anchoring these molecules in the active site of Pim-1 kinase. Principal component analysis (PCA) of Pim-1 kinase's conformational changes revealed that potent inhibitors like quercetin, galangin, and kaempferol significantly restricted the enzyme's flexibility, suggesting potential inhibitory effect. These findings provide insights into the structural interactions between flavonoids and Pim-1 kinase, offering a foundation for future experimental investigations. However, further studies, including in vitro and in vivo validation, are necessary to assess the pharmacological relevance and specificity of flavonoids in cancer therapy.

Tertiary Plasticity Drives the Efficiency of Enterocin 7B Interactions with Lipid Membranes

The ability of antimicrobial peptides to efficiently kill their bacterial targets depends on the efficiency of their binding to the microbial membrane. In the case of enterocins, there is a three-part interaction: initial binding, unpacking of helices on the membrane surface, and permeation of the lipid bilayer. Helical unpacking is driven by disruption of the peptide hydrophobic core when in contact with membranes. Enterocin 7B is a leaderless enterocin antimicrobial peptide produced from Enterococcus faecalis that functions alone, or with its cognate partner enterocin 7A, to efficiently kill a wide variety of Gram-stain positive bacteria. To better characterize the role that tertiary structural plasticity plays in the ability of enterocin 7B to interact with the membranes, a series of arginine single-site mutants were constructed that destabilize the hydrophobic core to varying degrees. A series of experimental measures of structure, stability, and function, including CD spectra, far UV CD melting profiles, minimal inhibitory concentrations analysis, and release kinetics of calcein, show that decreased stabilization of the hydrophobic core is correlated with increased efficiency of a peptide to permeate membranes and in killing bacteria. Finally, using the computational technique of adaptive steered molecular dynamics, we found that the atomistic/energetic landscape of peptide mechanical unfolding leads to free energy differences between the wild type and its mutants, whose trends correlate well with our experiment.

Telescoping-Solvation-Box Protocol-Based Umbrella Sampling Method for Potential Mean Force Estimation

Previously, it was shown that the telescoping box scheme, in combination with adaptive steered molecule dynamics (ASMD), can be used to estimate the potential of mean force (PMF) with a decrease in computational cost associated with large solvation boxes. Since ASMD reduces to umbrella sampling (US) in the limit of infinitely slow pulling velocity, a hypothesis was made that the telescoping box scheme can be extended to include the US method. The hypothesis was tested using the unfolding pathway of a polyalanine peptide in a water box and translocation of α-tocopherol through the human membrane. Two different approaches were tried: telescoping US (TELUS), in which the number of solvent molecules was linearly coupled to the reaction coordinate, and Block-TELUS, which was a compromise between the fixed solvation box of the US and the window-wise variable solvation box of TELUS. In the case of polyalanine peptide in a water box, both approaches gave overlapping potential of mean force (PMF) concerning the benchmark US-PMF. Window-wise comparison of properties like root-mean-square inner product, Ramachandran plot, α-helix content, and hydrogen bond formation was used to verify that both approaches captured the same dynamics as the US method. Applying the Block-TELUS protocol in the system with diffusing α-tocopherol through the bilayer resulted in overlapping PMF to its US benchmark. A comparison between the window-wise orientation of the chromanol headgroup also produced almost identical results. This study concluded that with the careful application of telescoping solvation boxes, a less computationally expensive US could be performed for systems where the effect of distant solvent molecules on the configurational space sampled in the window depends on the value of the reaction coordinate.

Correlation between chemical denaturation and the unfolding energetics of Acanthamoeba actophorin

The actin filament network is in part remodeled by the action of a family of filament severing proteins that are responsible for modulating the ratio between monomeric and filamentous actin. Recent work on the protein actophorin from the amoeba Acanthamoeba castellani identified a series of site-directed mutations that increase the thermal stability of the protein by 22°C. Here, we expand this observation by showing that the mutant protein is also significantly stable to both equilibrium and kinetic chemical denaturation, and employ computer simulations to account for the increase in thermal or chemical stability through an accounting of atomic-level interactions. Specifically, the potential of mean force (PMF) can be obtained from steered molecular dynamics (SMD) simulations in which a protein is unfolded. However, SMD can be inefficient for large proteins as they require large solvent boxes, and computationally expensive as they require increasingly many SMD trajectories to converge the PMF. Adaptive steered molecular dynamics (ASMD) overcomes the second of these limitations by steering the particle in stages, which allows for convergence of the PMF using fewer trajectories compared with SMD. Use of the telescoping water scheme within ASMD partially overcomes the first of these limitations by reducing the number of waters at each stage to only those needed to solvate the structure within a given stage. In the PMFs obtained from ASMD, the work of unfolding Acto-2 was found to be higher than the Acto-WT by approximately 120 kCal/mol and reflects the increased stability seen in the chemical denaturation experiments. The evolution of the average number of hydrogen bonds and number of salt bridges during the pulling process provides a mechanistic view of the structural changes of the actophorin protein as it is unfolded, and how it is affected by the mutation in concert with the energetics reported through the PMF.

Conservation and divergence of the G-interfaces of Drosophila melanogaster septins

Septins possess a conserved guanine nucleotide-binding (G) domain that participates in the stabilization of organized hetero-oligomeric complexes which assemble into filaments, rings and network-like structures. The fruit fly, Drosophila melanogaster, has five such septin genes encoding Sep1, Sep2, Sep4, Sep5 and Pnut. Here, we report the crystal structure of the heterodimer formed between the G-domains of Sep1 and Sep2, the first from an insect to be described to date. A G-interface stabilizes the dimer (in agreement with the expected arrangement for the Drosophila hexameric particle) and this bears significant resemblance to its human counterparts, even down to the level of individual amino acid interactions. On the other hand, a model for the G-interface formed between the two copies of Pnut which occupy the centre of the hexamer, shows important structural differences, including the loss of a highly favourable bifurcated salt-bridge network. Whereas wild-type Pnut purifies as a monomer, the reintroduction of the salt-bridge network results in stabilizing the dimeric interface in solution as shown by size exclusion chromatography and thermal stability measurements. Adaptive steered molecular dynamics reveals an unzipping mechanism for dimer dissociation which initiates at a point of electrostatic repulsion within the switch II region. Overall, the data contribute to a better understanding of the molecular interactions involved in septin assembly/disassembly.

Dissociation Rate Calculation via Constant-Force Steered Molecular Dynamics Simulation

Steered molecular dynamics (SMD) simulations are used to study molecular dissociation events by applying a harmonic force to the molecules and pulling them at a constant velocity. Instead of constant-velocity pulling, we use a constant force: the constant-force SMD (CF-SMD) simulation. The CF-SMD simulation employs a constant force to reduce the activation barrier of molecular dissociation, thereby enhancing the dissociation event. Here, we present the capability of the CF-SMD simulation to estimate the dissociation time at equilibrium. We performed all-atom CF-SMD simulations for NaCl and protein-ligand systems, producing dissociation time at various forces. We extrapolated these values to the dissociation rate without a constant force using Bell's model or the Dudko-Hummer-Szabo model. We demonstrate that the CF-SMD simulations with the models predicted the dissociation time in equilibrium. A CF-SMD simulation is a powerful tool for estimating the dissociation rate in a direct and computationally efficient manner.

Free energy and kinetic rate calculation via non-equilibrium molecular simulation: application to biomolecules

Non-equilibrium molecular dynamics (NEMD) simulation has been recognized as a powerful tool for examining biomolecules and provides fruitful insights into not only non-equilibrium but also equilibrium processes. We review recent advances in NEMD simulation and relevant, fundamental results of non-equilibrium statistical mechanics. We first introduce Crooks fluctuation theorem and Jarzynski equality that relate free energy difference to work done on a physical system during a non-equilibrium process. The theorems are beneficial for the analysis of NEMD trajectories. We then describe rate theory, a framework to calculate molecular kinetics from a non-equilibrium process; this theoretical framework enables us to calculate a reaction time-mean-first passage time-from NEMD trajectories. We, in turn, present recent NEMD techniques that apply an external force to a system to enhance molecular dissociation and introduce their application to biomolecules. Lastly, we show the current status of an appropriate selection of reaction coordinates for NEMD simulation.

Benchmarking Adaptive Steered Molecular Dynamics (ASMD) on CHARMM Force Fields

The potentials of mean force (PMFs) along the end-to-end distance of two different helical peptides have been obtained and benchmarked using the adaptive steered molecular dynamics (ASMD) method. The results depend strongly on the choice of force field driving the underlying all-atom molecular dynamics, and are reported with respect to the three most popular CHARMM force field versions: c22, c27 and c36. Two small peptides, ALA 10 and 1PEF, serve as the particular case studies. The comparisons between the versions of the CHARMM force fields provides both a qualitative and quantitative look at their performance in forced unfolding simulations in which peptides undergo large changes in structural conformations. We find that ASMD with the underlying c36 force field provides the most robust results for the selected benchmark peptides.

Implementation of Telescoping Boxes in Adaptive Steered Molecular Dynamics

Long-time dynamical processes, such as those involving protein unfolding and ligand interactions, can be accelerated and realized through steered molecular dynamics (SMD). The challenge has been the extraction of information from such simulations that generalize for complex nonequilibrium processes. The use of Jarzynski's equality opened the possibility of determining the free energy along the steered coordinate, but sampling over the nonequilibrium trajectories is slow to converge. Adaptive steered molecular dynamics (ASMD) and other related techniques have been introduced to overcome this challenge through the use of stages. Here, we take advantage of these stages to address the numerical cost that arises from the required use of very large solvent boxes. We introduce telescoping box schemes within adaptive steered molecular dynamics (ASMD) in which we adjust the solvent box between stages and thereby vary (and optimize) the required number of solvent molecules. We have benchmarked the method on a relatively long α-helical peptide, Ala₃₀, with respect to the potential of mean force and hydrogen bonds. We show that the use of telescoping boxes introduces little numerical error while significantly reducing the computational cost.

Picking on Carbonate: Kinetic Selectivity in the Encapsulation of Anions

Supramolecular hosts bind to inorganic anions at a fast rate and select them in proportion with thermodynamic stability of the corresponding [anion⊂host] complexes, forming in a reversible manner. In this study, we describe the action of hexapodal capsule 1 and its remarkable ability to select anions based on a large span of rates by which they enter this host. The thermodynamic affinity of 1 toward eighteen anions extends over eight orders of magnitude (0<K_a <10⁸  M^-1 ; ¹ H NMR spectroscopy). The capsule would retain CO₃ ^2- (K_a =10⁷  M^-1 ) for hours in the presence of eleven competing anions, including stronger binding SO₄ ^2- , HAsO₄ ^2- and HPO₄ ^2- (K_a =10⁷ -10⁸  M^-1 ). The observed selection resulted from 1 possessing narrow apertures (ca. 3×6 Å) comparable in size to anions (d=3.5-7.1 Å) slowing down the encapsulation to last from seconds to days. The unorthodox mode of action of 1 sets the stage for creating hosts that pick anions by their ability to access the host.

Adaptative Steered Molecular Dynamics Study of Mutagenesis Effects on Calcium Affinity in the Regulatory Domain of Cardiac Troponin C

Calcium-dependent cardiac muscle contraction is regulated by the protein complex troponin (cTn) and specifically by the regulatory N-terminal domain (N-cTnC) which contains one active Ca²⁺ binding site (site II). It has been previously shown that cardiac muscle contractility and functionality is affected by mutations in N-cTnC which alter calcium binding affinity. Here, we describe the application of adaptive steered molecular dynamics to characterize the influence of N-cTnC mutations on site II calcium binding affinity. We observed the correct trends for all of the studied calcium sensitizing and desensitizing mutants, in conjunction with loop II perturbations. Additionally, the potential of mean force accuracy was shown to increase substantially with increasingly slower speeds and using fewer trajectories. This study presents a novel approach to computationally estimate the Ca²⁺ binding affinity of N-cTnC structures and is a valuable potential tool to support the design and characterization of novel mutations with potential therapeutic benefits.