Energetics of Ortho-7 (oxime drug) translocation through the active-site gorge of tabun conjugated acetylcholinesterase

Oxime drugs translocate through the 20 Å active-site gorge of acetylcholinesterase in order to liberate the enzyme from organophosphorus compounds’ (such as tabun) conjugation. Here we report bidirectional steered molecular dynamics simulations of oxime drug (Ortho-7) translocation through the gorge of tabun intoxicated enzyme, in which time dependent external forces accelerate the translocation event. The simulations reveal the participation of drug-enzyme hydrogen bonding, hydrophobic interactions and water bridges between them. Employing nonequilibrium theorems that recovers the free energy from irreversible work done, we reconstruct potential of mean force along the translocation pathway such that the desired quantity represents an unperturbed system. The potential locates the binding sites and barriers for the drug to translocate inside the gorge. Configurational entropic contribution of the protein-drug binding entity and the role of solvent translational mobility in the binding energetics is further assessed.     

Molecular dynamics study of unbinding of the avidin-biotin complex.

We report molecular dynamics simulations that induce, over periods of 40-500 ps, the unbinding of biotin from avidin by means of external harmonic forces with force constants close to those of AFM cantilevers. The applied forces are sufficiently large to reduce the overall binding energy enough to yield unbinding within the measurement time. Our study complements earlier work on biotin-streptavidin that employed a much larger harmonic force constant. The simulations reveal a variety of unbinding pathways, the role of key residues contributing to adhesion as well as the spatial range over which avidin binds biotin. In contrast to the previous studies, the calculated rupture forces exceed by far those observed. We demonstrate, in the framework of models expressed in terms of one-dimensional Langevin equations with a schematic binding potential, the associated Smoluchowski equations, and the theory of first passage times, that picosecond to nanosecond simulation of ligand unbinding requires such strong forces that the resulting protein-ligand motion proceeds far from the thermally activated regime of millisecond AFM experiments, and that simulated unbinding cannot be readily extrapolated to the experimentally observed rupture.     

Binding pathway of retinal to bacterio-opsin: a prediction by molecular dynamics simulations

Formation of bacteriorhodopsin (bR) from apoprotein and retinal has been studied experimentally, but the actual pathway, including the point of entry, is little understood. Molecular dynamics simulations provide a surprisingly clear prediction. A window between bR helices E and F in the transmembrane part of the protein can be identified as an entry point for retinal. Steered molecular dynamics, performed by applying a series of external forces in the range of 200-1000 pN over a period of 0.2 ns to retinal, allows one to extract this chromophore from bR once the Schiff base bond to Lys216 is cleaved. Extraction proceeds until the retinal tail forms a hydrogen bond network with Ala144, Met145, and Ser183 side groups lining the exit/entry window. The manipulation induces a distortion with a fitted root mean square deviation of coordinates (ignoring retinal, water, and hydrogen atoms) of less than 1.9 A by the time the retinal carbonyl reaches the protein surface. The forces needed to extract retinal are due to friction and do not indicate significant potential barriers. The simulations therefore suggest a pathway for the binding of retinal. Water molecules are found to play a crucial role in the binding process.     

Using enhanced sampling and structural restraints to refine atomic structures into low-resolution electron microscopy maps

For a variety of problems in structural biology, low-resolution maps generated by electron microscopy imaging are often interpreted with the help of various flexible-fitting computational algorithms. In this work, we systematically analyze the quality of final models of various proteins obtained via molecular dynamics flexible fitting (MDFF) by varying the map-resolution, strength of structural restraints, and the steering forces. We find that MDFF can be extended to understand conformational changes in lower-resolution maps if larger structural restraints and lower steering forces are used to prevent overfitting. We further show that the capabilities of MDFF can be extended by combining it with an enhanced conformational sampling method, temperature-accelerated molecular dynamics (TAMD). Specifically, either TAMD can be used to generate better starting configurations for MDFF fitting or TAMD-assisted MDFF (TAMDFF) can be performed to accelerate conformational search in atomistic simulations.     

Brownian dynamics theory for predicting internal and external blockages of tetraethylammonium in the KcsA potassium channel

The theory of Brownian dynamics is used to model permeation and the blocking of KcsA potassium channels by tetraethylammonium (TEA). A novel Brownian dynamics simulation algorithm is implemented that comprises two free energy profiles; one profile is seen by the potassium ions and the other by the TEA molecules whose shape is approximated by a sphere. Our simulations reveal that internally applied TEA blocks the passage of K⁺ ions by physically occluding the pore. A TEA molecule in the external reservoir encounters an attractive energy-well created by four tyrosine residues at position 82, in addition to all other attractive and repulsive forces impinging on it. Using Brownian dynamics, we investigate how deep the energy-well needs to be to reproduce the experimentally determined inhibitory constant k_i for the TEA blockade of KcsA or the mutant Shaker T449Y. The one-dimensional free energy profile obtained from molecular dynamics is first converted into a one-dimensional potential energy profile, and is then transformed into a three-dimensional free energy profile in Brownian dynamics by adding the short-range potential from the channel walls. When converted, the free energy profile calculated from molecular dynamics gives a well-depth of ∼10 kT. We systematically alter the depths of the profiles, and then use Brownian dynamics simulations to numerically determine the current versus TEA-concentration curves. We show that the sequence of binding and unbinding events of the TEA molecule to the binding pocket can be modeled by a first-order Markov process. The Brownian dynamics simulations also reveal that the probability of a TEA molecule binding to the binding pocket in KcsA potassium channels increases exponentially with TEA concentration and depends also on the applied potential and the K⁺ concentration in the simulation assembly.     

The Effect of Different Force Applications on the Protein-Protein Complex Barnase-Barstar

Steered molecular dynamics simulations are a tool to examine the energy landscape of protein-protein complexes by applying external forces. Here, we analyze the influence of the velocity and geometry of the probing forces on a protein complex using this tool. With steered molecular dynamics, we probe the stability of the protein-protein complex Barnase-Barstar. The individual proteins are mechanically labile. The Barnase-Barstar binding site is more stable than the folds of the individual proteins. By using different force protocols, we observe a variety of responses of the system to the applied tension.     

Mechanical Properties of the Icosahedral Shell of Southern Bean Mosaic Virus: A Molecular Dynamics Study

The mechanical properties of viral shells are crucial for viral assembly and infection. To study their distribution and heterogeneity on the viral surface, we performed atomistic force-probe molecular dynamics simulations of the complete shell of southern bean mosaic virus, a prototypical T = 3 virus, in explicit solvent. The simulation system comprised more than 4,500,000 atoms. To facilitate direct comparison with atomic-force microscopy (AFM) measurements, a Lennard-Jones sphere was used as a model of the AFM tip, and was pushed with different velocities toward the capsid protein at 19 different positions on the viral surface. A detailed picture of the spatial distribution of elastic constants and yielding forces was obtained that can explain corresponding heterogeneities observed in previous AFM experiments. Our simulations reveal three different deformation regimes: a prelinear regime of outer surface atom rearrangements, a linear regime of elastic capsid deformation, and a rearrangement regime that describes irreversible structural changes and the transition from elastic to plastic deformation. For both yielding forces and elastic constants, a logarithmic velocity dependency is evident over nearly two decades, the explanation for which requires including nonequilibrium effects within the established theory of enforced barrier crossing.     

Computational study of the activity,dynamics, energetics and conformations of insulin analogues using molecular dynamics simulations: Application to hyperinsulinemia and the critical residue B26

Due to the increasing prevalence of diabetes, finding therapeutic analogues for insulin has become an urgent issue. While many experimental studies have been performed towards this end, they have limited scope to examine all aspects of the effect of a mutation. Computational studies can help to overcome these limitations, however, relatively few studies that focus on insulin analogues have been performed to date. Here, we present a comprehensive computational study of insulin analogues—three mutant insulins that have been identified with hyperinsulinemia and three mutations on the critical B26 residue that exhibit similar binding affinity to the insulin receptor—using molecular dynamics simulations with the aim of predicting how mutations of insulin affect its activity, dynamics, energetics and conformations. The time evolution of the conformers is studied in long simulations. The probability density function and potential of mean force calculations are performed on each insulin analogue to unravel the effect of mutations on the dynamics and energetics of insulin activation. Our conformational study can decrypt the key features and molecular mechanisms that are responsible for an enhanced or reduced activity of an insulin analogue. We find two key results: 1) hyperinsulinemia may be due to the drastically reduced activity (and binding affinity) of the mutant insulins. 2) Y26_BS and Y26_BE are promising therapeutic candidates for insulin as they are more active than WT-insulin. The analysis in this work can be readily applied to any set of mutations on insulin to guide development of more effective therapeutic analogues.     

Soft ellipsoid potential for biaxial molecules

A soft ellipsoid contact pair potential model for biaxial molecules is proposed considering the configuration dependent energy anisotropy explicitly alongwith their geometrical aspects. Analytic expressions of forces and torques for this coarse-grained potential are provided to make it convenient for molecular dynamics simulation study. Additionally, some molecular dynamics simulations are performed to generate biaxial smectic and nematic liquid crystal phases.     

Characterization of the binding of angiotensin II receptor blockers to human serum albumin using docking and molecular dynamics simulation

Human serum albumin (HSA), the most abundant protein found in blood plasma, transports many drugs and ligands in the circulatory system. The drug binding ability of HSA strongly influences free drug concentrations in plasma, and is directly related to the effectiveness of clinical therapy. In current work, binding of HSA to angiotensin II receptor blockers (ARBs) are investigated using docking and molecular dynamics (MD) simulations. Docking results demonstrate that the main HSA–ARB binding site is subdomain IIIA of HSA. Simulation results reveal clearly how HSA binds with valsartan and telmisartan. Interestingly, electrostatic interactions appear to be more important than hydrophobic interactions in stabilizing binding of valsartan to HSA, and vice versa for HSA–telmisartan. The molecular distance between HSA Trp214 (donor) and the drug (acceptor) can be measured by fluorescence resonance energy transfer (FRET) in experimental studies. The average distances between Trp-214 and ARBs are estimated here based on our MD simulations, which could be valuable to future FRET studies. This work will be useful in the design of new ARB drugs with desired HSA binding affinity.     

Differences between apo and three holo forms of the intestinal fatty acid binding protein seen by molecular dynamics computer calculations

It is commonly believed that binding affinity can be estimated by consideration of local changes of ligand and protein. This paper discusses a set of molecular dynamics simulations of intestinal fatty acid binding protein addressing the protein's response to presence or absence of different ligands. A 5-ns simulation was performed of the protein without a ligand, and three simulations (one 5-ns and two 2-ns) were performed with different fatty acids bound. The results indicate that, although the basic protein structure is unchanged by the presence of the ligand, other properties are significantly affected by ligand binding. For example, zero-time covariance patterns between protein, bound waters, and ligand vary between the different simulations. Moreover, the interaction energies between ligand and specific residues indicate that different ligands are stabilized in different ways. In sum, the results suggest that binding thermodynamics within this system will need to be calculated not from a subset of nearby protein:ligand interactions, but will depend on a knowledge of the motions coupling together water, protein, and ligand.     

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.     

Protocol for MM/PBSA molecular dynamics simulations of proteins

Continuum solvent models have been employed in past years for understanding processes such as protein folding or biomolecular association. In the last decade, several attempts have been made to merge atomic detail molecular dynamics simulations with solvent continuum models. Among continuum models, the Poisson-Boltzmann solvent accessible surface area model is one of the oldest and most fundamental. Notwithstanding its wide usage for simulation of biomolecular electrostatic potential, the Poisson-Boltzmann equation has been very seldom used to obtain solvation forces for molecular dynamics simulation. We propose here a fast and reliable methodology to implement continuum forces in standard molecular mechanics and dynamics algorithms. Results for a totally unrestrained 1 ns molecular dynamics simulation of a small protein are quantitatively similar to results obtained by explicit solvent molecular dynamics simulations.     

Role of flexibility and polarity as determinants of the hydration of internal cavities and pockets in proteins

Molecular dynamics simulations of Staphylococcal nuclease and of 10 variants with internal polar or ionizable groups were performed to investigate systematically the molecular determinants of hydration of internal cavities and pockets in proteins. In contrast to apolar cavities in rigid carbon structures, such as nanotubes or buckeyballs, internal cavities in proteins that are large enough to house a few water molecules will most likely be dehydrated unless they contain a source of polarity. The water content in the protein interior can be modulated by the flexibility of protein elements that interact with water, which can impart positional disorder to water molecules, or bias the pattern of internal hydration that is stabilized. This might explain differences in the patterns of hydration observed in crystal structures obtained at cryogenic and room temperature conditions. The ability of molecular dynamics simulations to determine the most likely sites of water binding in internal pockets and cavities depends on its efficiency in sampling the hydration of internal sites and alternative protein and water conformations. This can be enhanced significantly by performing multiple molecular dynamics simulations as well as simulations started from different initial hydration states.     

Solvation Effects on the Conformational Behaviour of Gellan and Calcium Ion Binding to Gellan Double Helices

The contribution of the presence of solvent to the conformations adopted by disaccharide fragments within the repeat unit of gellan have been studied by molecular modelling techniques. Initial conformational energy searches, using a dielectric continuum to represent the solvent, provided starting geometries for a series of molecular dynamics simulations. The solution behaviour from these simulations was subsequently compared to fibre diffraction data of the potassium gellan salt. The present calculations indicate considerable flexibility of the glycosidic linkages, and this is discussed in relation to its effect on gel formation. One of the fragments was solvated with explicit water molecules. These calculations showed the same conformational behaviour as those simulations conducted in implicit solvent.Finally, a series of molecular dynamics (MD) simulations were performed to study the calcium binding to gellan. The results from this clearly showed a well defined binding site for this ion.Abbreviations MM molecular mechanics - MD molecular dynamics     

On Calculation of the Electrostatic Potential of a Phosphatidylinositol Phosphate-Containing Phosphatidylcholine Lipid Membrane Accounting for Membrane Dynamics

Many signaling events require the binding of cytoplasmic proteins to cell membranes by recognition of specific charged lipids, such as phosphoinositol-phosphates. As a model for a protein-membrane binding site, we consider one charged phosphoinositol phosphate (PtdIns(3)P) embedded in a phosphatidylcholine bilayer. As the protein-membrane binding is driven by electrostatic interactions, continuum solvent models require an accurate representation of the electrostatic potential of the phosphoinositol phosphate-containing membrane. We computed and analyzed the electrostatic potentials of snapshots taken at regular intervals from molecular dynamics simulations of the bilayer. We observe considerable variation in the electrostatic potential of the bilayer both along a single simulation and between simulations performed with the GAFF or CHARMM c36 force fields. However, we find that the choice of GAFF or CHARMM c36 parameters has little effect on the electrostatic potential of a given configuration of the bilayer with a PtdIns(3)P embedded in it. From our results, we propose a remedian averaging method for calculating the electrostatic potential of a membrane system that is suitable for simulations of protein-membrane binding with a continuum solvent model.     

Predicting large-scale conformational changes in proteins using energy-weighted normal modes

We report the development of a method to improve the sampling of protein conformational space in molecular simulations. It is shown that a principal component analysis of energy-weighted normal modes in Cartesian coordinates can be used to extract vectors suitable for describing the dynamics of protein substructures. The method can operate with either atomistic or user-defined coarse-grained models of protein structure. An implicit reverse coarse-graining allows the dynamics of all-atoms to be recovered when a coarse-grained model is used. For an external test set of four proteins, it is shown that the new method is more successful than normal mode analysis in describing the large-scale conformational changes observed on ligand binding. The method has potential applications in protein-ligand and protein-protein docking and in biasing molecular dynamics simulations.     

Energetic Changes Caused by Antigenic Module Insertion in a Virus-Like Particle Revealed by Experiment and Molecular Dynamics Simulations

The success of recombinant virus-like particles (VLPs) for human papillomavirus and hepatitis B demonstrates the potential of VLPs as safe and efficacious vaccines. With new modular designs emerging, the effects of antigen module insertion on the self-assembly and structural integrity of VLPs should be clarified so as to better enabling improved design. Previous work has revealed insights into the molecular energetics of a VLP subunit, capsomere, comparing energetics within various solution conditions known to drive or inhibit self-assembly. In the present study, molecular dynamics (MD) simulations coupled with the molecular mechanics-Poisson-Boltzmann surface area (MM-PBSA) method were performed to examine the molecular interactions and energetics in a modular capsomere of a murine polyomavirus (MPV) VLP designed to protect against influenza. Insertion of an influenza antigenic module is found to lower the binding energy within the capsomere, and a more active state is observed in Assembly Buffer as compared with that in Stabilization Buffer, which has been experimentally validated through measurements using differential scanning calorimetry. Further in-depth analysis based on free-energy decomposition indicates that destabilized binding can be attributed to electrostatic interaction induced by the chosen antigen module. These results provide molecular insights into the conformational stability of capsomeres and their abilities to be exploited for antigen presentation, and are expected to be beneficial for the biomolecular engineering of VLP vaccines.     

Computational studies of the binding mechanisms of fullerenes to human serum albumin

Fullerene and its derivatives show promising prospects for applications in a vast array of biological systems. A key aspect concerning their biomedical applications is how they interact with proteins from molecular levels, which is still poorly understood. In the current study, we investigated the structural and thermodynamic basis of the interactions between two pharmacologically relevant fullerene derivatives and human serum albumin (HSA) using molecular docking, molecular dynamics simulations, and binding free energy calculations. Our results demonstrate that fullerenes steadily bind with HSA at the interfacial cavity formed by subdomains IIA and IIIA. In agreement with available experimental data, our simulations show that the global structure of HSA becomes more compact in the presence of fullerene, while local structural dynamics of the binding cavity behaves diversely depending on the chemical properties of bound fullerenes. Binding free energy calculations confirmed that the interactions between fullerenes and HSA are dominantly stabilized by van der Waals forces and they further allowed the identification of key residues involved in fullerene binding. The structural and energetic insights obtained from this work may help for the development of fullerene-based drug delivery devices and therapeutic agents with improved biological profile.     

Long range molecular dynamics study of regulation of eukaryotic glucosamine-6-phosphate synthase activity by UDP-GlcNAc

Glucosamine-6-phosphate (GlcN-6-P) synthase catalyses the first and practically irreversible step in hexosamine metabolism. The final product of this pathway, uridine 5’ diphospho N-acetyl-D-glucosamine (UDP-GlcNAc), is an essential substrate for assembly of bacterial and fungal cell walls. Moreover, the enzyme is involved in phenomenon of hexosamine induced insulin resistance in type II diabetes, which makes it a potential target for antifungal, antibacterial and antidiabetic therapy. The crystal structure of the isomerase domain of GlcN-6-P synthase from human pathogenic fungus Candida albicans, in complex with UDP-GlcNAc has been solved recently but it has not revealed the molecular mechanism of inhibition taking place under UDP-GlcNAc influence, the unique feature of the eukaryotic enzyme. UDP-GlcNAc is a physiological inhibitor of GlcN-6-P synthase, binding about 1 nm away from the active site of the enzyme. In the present work, comparative molecular dynamics simulations of the free and UDP-GlcNAc-bounded structures of GlcN-6-P synthase have been performed. The aim was to complete static X-ray structural data and detect possible changes in the dynamics of the two structures. Results of the simulation studies demonstrated higher mobility of the free structure when compared to the liganded one. Several amino acid residues were identified, flexibility of which is strongly affected upon UDP-GlcNAc binding. Importantly, the most fixed residues are those related to the inhibitor binding process and to the catalytic reaction. The obtained results constitute an important step toward understanding of mechanism of GlcN-6-P synthase inhibition by UDP-GlcNAc molecule.