A steered molecular dynamics method with adaptive direction adjustments

In this paper, a new steered molecular dynamics (SMD) method with adjusting pulling direction is proposed to search an optimum trajectory of ligand dissociation. A multiobjective model and a searching technique based on information entropy with multi-population are developed to optimize the pulling direction. The improved method has been used to dissociate the substrate-bound complex structure of cytochrome P450 3A4-metyrapone. A more favorable dissociation pathway can be gained. The results show that the new pathway obtained by the proposed method has less dissociation time, smaller rupture force and lower energy barrier than that by the conventional SMD.     

Ca(2+) dissociation from the C-terminal EF-hand pair in calmodulin: a steered molecular dynamics study

We used steered molecular dynamics (SMD) to simulate the process of Ca²⁺ dissociation from the EF-hand motifs of the C-terminal lobe of calmodulin. Based on an analysis of the pulling forces, the dissociation sequences and the structural changes, we show that the Ca²⁺-coordinating residues lose their binding to Ca²⁺ in a stepwise fashion. The two Ca²⁺ ions dissociate from the two EF-hands simultaneously, with two distinct groups among the five Ca²⁺-coordinating residues affecting the EF-hand conformational changes differently. These results provide new insights into the effects of Ca²⁺ on calmodulin conformation, from which a novel sequential mechanism of Ca²⁺-calmodulin dissociation is proposed.     

Carbon Nanotube-Encapsulated Drug Penetration Through the Cell Membrane: An Investigation Based on Steered Molecular Dynamics Simulation

Understanding the penetration mechanisms of carbon nanotube (CNTs)-encapsulated drugs through the phospholipid bilayer cell membrane is an important issue for the development of intracellular drug delivery systems. In the present work, steered molecular dynamics (SMD) simulation was used to explore the possibility of penetration of a polar drug, paclitaxel (PTX), encapsulated inside the CNT, through a dipalmitoylphosphatidylcholine bilayer membrane. The interactions between PTX and CNT and between PTX and the confined water molecules inside the CNT had a significant effect on the penetration process of PTX. The results reveal that the presence of a PTX molecule increases the magnitude of the pulling force. The effect of pulling velocity on the penetration mechanism was also investigated by a series of SMD simulations, and it is shown that the pulling velocity had a significant effect on pulling force and the interaction between lipid bilayer and drug molecule.     

Force-induced unfolding of human telomeric G-quadruplex: A steered molecular dynamics simulation study

We study the unfolding of a parallel G-quadruplex from human telomeric DNA by mechanical stretching using steered molecular dynamics (MD) simulation. We find that the force curves and unfolding processes strongly depend on the pulling sites. With pulling sites located on the sugar-phosphate backbone, the force-extension curve shows a single peak and the unfolding proceeds sequentially. Pulling sites located on the terminal nucleobases lead to a force-extension curve with two peaks and the unfolding is more cooperative. Simulations of the refolding of partially unfolded quadruplexes show very different behavior for the two different pulling modalities. In particular, starting from an unfolded state prepared by nucleobase pulling leads to a long-lived intermediate state whose existence is also corroborated by the free energy profile computed with the Jarzynski equation. Based on this observation, we propose a novel folding pathway for parallel G-quadruplexes with the human telomere sequence.     

Evaluation of protein–ligand affinity prediction using steered molecular dynamics simulations

In computational drug design, ranking a series of compound analogs in a manner that is consistent with experimental affinities remains a challenge. In this study, we evaluated the prediction of protein–ligand binding affinities using steered molecular dynamics simulations. First, we investigated the appropriate conditions for accurate predictions in these simulations. A conic harmonic restraint was applied to the system for efficient sampling of work values on the ligand unbinding pathway. We found that pulling velocity significantly influenced affinity predictions, but that the number of collectable trajectories was less influential. We identified the appropriate pulling velocity and collectable trajectories for binding affinity predictions as 1.25 Å/ns and 100, respectively, and these parameters were used to evaluate three target proteins (FK506 binding protein, trypsin, and cyclin-dependent kinase 2). For these proteins using our parameters, the accuracy of affinity prediction was higher and more stable when Jarzynski’s equality was employed compared with the second-order cumulant expansion equation of Jarzynski’s equality. Our results showed that steered molecular dynamics simulations are effective for predicting the rank order of ligands; thus, they are a potential tool for compound selection in hit-to-lead and lead optimization processes.     

Free-energy landscape of glycerol permeation through aquaglyceroporin GlpF determined from steered molecular dynamics simulations

The free-energy landscape of glycerol permeation through the aquaglyceroporin GlpF has been estimated in the literature by the nonequilibrium method of steered molecular dynamics (SMD) simulations and by the equilibrium method of adaptive biasing force (ABF) simulations. However, the ABF results qualitatively disagree with the SMD results that were based on the Jarzynski equality (JE) relating the equilibrium free-energy difference to the nonequilibrium work of the irreversible pulling experiments. In this paper, I present a new SMD study of the glycerol permeation through GlpF to explore the free-energy profile of glycerol along the permeation channel. Instead of the JE in terms of thermodynamic work, I use the fluctuation-dissipation theorem (FDT) of Brownian dynamics (BD), in terms of mechanical work, for extracting the free-energy difference from the nonequilibrium work of irreversible pulling experiments. The results of this new SMD-BD-FDT study are in agreement with the experimental data and with the ABF results.     

SARS-CoV-2 spike binding to ACE2 is stronger and longer ranged due to glycan interaction

Highly detailed steered molecular dynamics simulations are performed on differently glycosylated receptor binding domains of the severe acute respiratory syndrome coronavirus-2 spike protein. The binding strength and the binding range increase with glycosylation. The interaction energy rises very quickly when pulling the proteins apart and only slowly drops at larger distances. We see a catch-slip-type behavior whereby interactions during pulling break and are taken over by new interactions forming. The dominant interaction mode is hydrogen bonds, but Lennard-Jones and electrostatic interactions are relevant as well.     

Force-based analysis of multidimensional energy landscapes: application of dynamic force spectroscopy and steered molecular dynamics simulations to an antibody fragment-peptide complex

Multidimensional energy landscapes are an intrinsic property of proteins and define their dynamic behavior as well as their response to external stimuli. In order to explore the energy landscape and its implications on the dynamic function of proteins dynamic force spectroscopy and steered molecular dynamics (SMD) simulations have proved to be important tools. In this study, these techniques have been employed to analyze the influence of the direction of the probing forces on the complex of an antibody fragment with its peptide antigen. Using an atomic force microscope, experiments were performed where the attachment points of the 12 amino acid long peptide antigen were varied. These measurements yielded clearly distinguishable basal dissociation rates and potential widths, proving that the direction of the applied force determines the unbinding pathway. Complementary atomistic SMD simulations were performed, which also show that the unbinding pathways of the system are dependent on the pulling direction. However, the main barrier to be crossed was independent of the pulling direction and is represented by a backbone hydrogen bond between Gly^H-H40 of the antibody fragment and Glu^Oε-6_peptide of the peptide. For each pulling direction, the observed barriers can be correlated with the rupture of specific interactions, which stabilize the bound complex. Furthermore, although the SMD simulations were performed at loading rates exceeding the experimental rates by orders of magnitude due to computational limitations, a detailed comparison of the barriers that were overcome in the SMD simulations with the data obtained from the atomic force microscope unbinding experiments show excellent agreement.     

Mechanism of CDK5 activation revealed by steered molecular dynamics simulations and energy calculations

In the current work, CDK5/p25 complexes were pulled apart by applying external forces with steered molecular dynamics (SMD) simulations. The crucial interactions between the kinase and the activation protein were investigated and the SMD simulations showed that several activation-relevant motifs of CDK5 leave p25 in sequence during the pulling and lead to an apo-CDK2 like CDK5 structure after separation. Based on systematic examination of hydrogen bond breaking and classical MD/molecular mechanics-generalized Born/surface area) (MM-GBSA) calculations, a CDK5 activation mechanism by p25 is suggested. This is the first step towards the systemic development of CDK inhibitors and the mechanism proposed could lead to a better understanding of the protein–protein recognition characteristics between the kinase and its activator.     

The Molecular Mechanism Underlying Mechanical Anisotropy of the Protein GB1

Mechanical responses of elastic proteins are crucial for their biological function and nanotechnological use. Loading direction has been identified as one key determinant for the mechanical responses of proteins. However, it is not clear how a change in pulling direction changes the mechanical unfolding mechanism of the protein. Here, we combine protein engineering, single-molecule force spectroscopy, and steered molecular dynamics simulations to systematically investigate the mechanical response of a small globular protein GB1. Force versus extension profiles from both experiments and simulations reveal marked mechanical anisotropy of GB1. Using native contact analysis, we relate the mechanically robust shearing geometry with concurrent rupture of native contacts. This clearly contrasts the sequential rupture observed in simulations for the mechanically labile peeling geometry. Moreover, we identify multiple distinct mechanical unfolding pathways in two loading directions. Implications of such diverse unfolding mechanisms are discussed. Our results may also provide some insights for designing elastomeric proteins with tailored mechanical properties.     

Molecular dynamics simulations to the bidirectional adhesion signaling pathway of integrin αVβ3

The bidirectional force transmission process of integrin through the cell membrane is still not well understood. Several possible mechanisms have been discussed in literature on the basis of experimental data, and in this study, we investigate these mechanisms by free and steered molecular dynamics simulations. For the first time, constant velocity pulling on the complete integrin molecule inside a dipalmitoyl-phosphatidylcholine membrane is conducted. From the results, the most likely mechanism for inside-out and outside-in signaling is the switchblade model with further separation of the transmembrane helices.     

Investigations of Takeout proteins’ ligand binding and release mechanism using molecular dynamics simulation

Takeout (To) proteins exist in a diverse range of insect species. They are involved in many important processes of insect physiology and behaviors. As the ligand carriers, To proteins can transport the small molecule to the target tissues. However, ligand release mechanism of To proteins is unclear so far. In this contribution, the process and pathway of the ligand binding and release are revealed by conventional molecular dynamics simulation, steered molecular dynamics simulation and umbrella sampling methods. Our results show that the α4-side of the protein is the unique gate for the ligand binding and release. The structural analysis confirms that the internal cavity of the protein has high rigidity, which is in accordance with the recent experimental results. By using the potential of mean force calculations in combination with residue cross correlation calculation, we concluded that the binding between the ligand and To proteins is a process of conformational selection. Furthermore, the conformational changes of To proteins and the hydrophobic interactions both are the key factors for ligand binding and release.     

Predicting doxorubicin drug delivery by single-walled carbon nanotube through cell membrane in the absence and presence of nicotine molecules: a molecular dynamics simulation study

Abstract

In order to study the interaction of the anticancer agent Doxorubicin with the single-walled carbon nanotubes with different diameters as drug delivery systems, the molecular dynamics (MD) simulations have been used. Also, for design and development of intracellular Doxorubicin drug delivery systems, a series of steered MD simulations are applied to explore the possibility of encapsulated Doxorubicin–carbon nanotube penetration through a lipid bilayer in presence and absence of Nicotine molecules at different pulling rates. Our simulation results showed that in spite of the adsorption of drug molecules on the outer sidewall of the nanotubes, the spontaneous localization of one Doxorubicin molecule into the cavity of the nanovectors with larger diameters is observed. It is found that the presence of Nicotine molecules in extracellular medium increases the required force for pulling nanotube-encapsulated drug as well as the required time for penetration process, especially at higher velocity. Also, the entering process of the Nicotine molecules into the carbon nanotube causes that the encapsulated drug molecule is fully released in the hydrophobic phase of the lipid bilayer.

Communicated by Ramaswamy H. Sarma     

On the impact of nanotube diameter on biomembrane indentation – Computer simulations study

The influence of the single-walled carbon nanotubes on the phospholipid bilayer has been studied using steered molecular dynamics (SMD) simulations. The impact of different nanotubes on the phospholipid bilayer structure is discussed as well as the speed of indentation. Additionally, a series of simulations with pulling out of the nanotubes from the membrane were performed. The deflection of the membrane in both nanoindenation and extraction processes is also discussed. The self-sealing ability of membrane during this process is examined. Complete degradation of the bilayer was not observed even for the most invasive nanoindentation process studied. The obtained results show that carbon nanotubes can be regarded as potential drug carriers for targeted therapy.     

Can an Atomic Force Microscope Sequence DNA Using a Nanopore?

R. Bension has proposed that single molecules of DNA could be sequenced rapidly, in long sequential reads, by reading off the force required to pull a tightly fitting molecular ring over each base in turn using an atomic force microscope (AFM). We present molecular dynamics simulations that indicate that pulling DNA very rapidly (m/s) could generate large force peaks as each base is passed (∼1 nN) with significant differences (∼0.5 nN) between purine and pyrimidine. These speeds are six orders of magnitude faster than could be read out by a conventional AFM, and extending the calculations to accessible speeds using Kramers’ theory shows that thermal fluctuations dominate the process with the result that purine and pyrimidine cannot be distinguished with the pulling speeds attained by current AFM technology.     

The Ca2+ Influence on Calmodulin Unfolding Pathway: A Steered Molecular Dynamics Simulation Study

The force-induced unfolding of calmodulin (CaM) was investigated at atomistic details with steered molecular dynamics. The two isolated CaM domains as well as the full-length CaM were simulated in N-C-terminal pulling scheme, and the isolated N-lobe of CaM was studied specially in two other pulling schemes to test the effect of pulling direction and compare with relevant experiments. Both Ca²⁺-loaded CaM and Ca²⁺-free CaM were considered in order to define the Ca²⁺ influence to the CaM unfolding. The results reveal that the Ca²⁺ significantly affects the stability and unfolding behaviors of both the isolated CaM domains and the full-length CaM. In Ca²⁺-loaded CaM, N-terminal domain unfolds in priori to the C-terminal domain. But in Ca²⁺-free CaM, the unfolding order changes, and C-terminal domain unfolds first. The force-extension curves of CaM unfolding indicate that the major unfolding barrier comes from conquering the interaction of two EF-hand motifs in both N- and C- terminal domains. Our results provide the atomistic-level insights in the force-induced CaM unfolding and explain the observation in recent AFM experiments.     

Calculation of the Generalized Susceptibility for a Highly Supercooled Fluid Through Molecular-Dynamics Simulation

Abstract

A new method of computation of generalized susceptibility and dynamical structure factor through molecular dynamics (MD) simulation is proposed. This gives rise to a reliable and accurate result more than that calculated from a conventional method with a direct Fourier transformation. Computational results are presented for the imaginary part of the generalized susceptibility, X″ (ω), for a binary soft-sphere fluid with a super-long-time molecular dynamics (MD) simulation. Both α- and β-peaks in X″ (ω) in a supercooled fluid is shown for the first time through the present MD computation. The MD result obtained is in a good agreement with that obtained by the trapping diffusion model, which we have previously proposed for the glass transition.