Molecular dynamics simulation of a synthetic four-alpha-helix bundle that binds the anesthetic halothane.

The structural features of binding sites for volatile anesthetics are examined by performing a molecular dynamics simulation study of the synthetic four-alpha-helix bundles (Aalpha2)2, which are formed by association of two 62-residue di-alpha-helical peptides. The peptide bundle (Aalpha2)2 was designed by Johansson et al. [Biochemistry 37 (1998) 1421-1429] and was shown experimentally to have a high affinity for the binding of the anesthetic halothane (CF3CBrCIH) in a hydrophobic cavity. Since (Aalpha2)2 can exhibit either the anti or syn topologies, the two distinct bundles are simulated both in the presence and in the absence of halothane. Nanosecond length molecular dynamics trajectories were generated for each system at room temperature (T = 298 K). The structural and dynamic effects of the inclusion of halothane are compared, illustrating that the structures are stable over the course of the simulation; that the (Aalpha2)2 bundles have suitable pockets that can accommodate halothane; that the halothane remains in the designed hydrophobic cavity in close proximity to the Trp residues with a preferred orientation; and that the dimensions of the peptide are perturbed by the inclusion of an anesthetic molecule.     

A designed four-alpha-helix bundle that binds the volatile general anesthetic halothane with high affinity

The structural features of volatile anesthetic binding sites on proteins are being examined with the use of a defined model system consisting of a four-alpha-helix bundle scaffold with a hydrophobic core. Previous work has suggested that introducing a cavity into the hydrophobic core improves anesthetic binding affinity. The more polarizable methionine side chain was substituted for a leucine, in an attempt to enhance the dispersion forces between the ligand and the protein. The resulting bundle variant has an improved affinity (K(d) = 0.20 +/- 0.01 mM) for halothane binding, compared with the leucine-containing bundle (K(d) = 0.69 +/- 0.06 mM). Photoaffinity labeling with (14)C-halothane reveals preferential labeling of the W15 residue in both peptides, supporting the view that fluorescence quenching by bound anesthetic reports both the binding energetics and the location of the ligand in the hydrophobic core. The rates of amide hydrogen exchange were similar for the two bundles, suggesting that differences in binding affinity were not due to changes in protein stability. Binding of halothane to both four-alpha-helix bundle proteins stabilized the native folded conformations. Molecular dynamics simulations of the bundles illustrate the existence of the hydrophobic core, containing both W15 residues. These results suggest that in addition to packing defects, enhanced dispersion forces may be important in providing higher affinity anesthetic binding sites. Alternatively, the effect of the methionine substitution on halothane binding energetics may reflect either improved access to the binding site or allosteric optimization of the dimensions of the binding pocket. Finally, preferential stabilization of folded protein conformations may represent a fundamental mechanism of inhaled anesthetic action.     

Expression and characterization of a four-alpha-helix bundle protein that binds the volatile general anesthetic halothane

The structural features of volatile anesthetic binding sites on proteins are being investigated with the use of a defined model system consisting of a four-alpha-helix bundle scaffold with a hydrophobic core. The current study describes the bacterial expression, purification, and initial characterization of the four-alpha-helix bundle (Aalpha(2)-L1M/L38M)(2). The alpha-helical content and stability of the expressed protein are comparable to that of the chemically synthesized four-alpha-helix bundle (Aalpha(2)-L38M)(2) reported earlier. The affinity for binding halothane is somewhat improved with a K(d) = 120 +/- 20 microM as determined by W15 fluorescence quenching, attributed to the L1M substitution. Near-UV circular dichroism spectroscopy demonstrated that halothane binding changes the orientation of the aromatic residues in the four-alpha-helix bundle. Nuclear magnetic resonance experiments reveal that halothane binding results in narrowing of the peaks in the amide region of the one-dimensional proton spectrum, indicating that bound anesthetic limits protein dynamics. This expressed protein should prove to be amenable to nuclear magnetic resonance structural studies on the anesthetic complexes, because of its relatively small size (124 residues) and the high affinities for binding volatile anesthetics. Such studies will provide much needed insight into how volatile anesthetics interact with biological macromolecules and will provide guidelines regarding the general architecture of binding sites on central nervous system proteins.     

Effect of four-alpha-helix bundle cavity size on volatile anesthetic binding energetics

Currently, it is thought that inhalational anesthetics cause anesthesia by binding to ligand-gated ion channels. This is being investigated using four-alpha-helix bundles, small water-soluble analogues of the transmembrane domains of the "natural" receptor proteins. The study presented here specifically investigates how multiple alanine-to-valine substitutions (which each decrease the volume of the internal binding cavity by 38 A(3)) affect structure, stability, and anesthetic binding affinity of the four-alpha-helix bundles. Structure remains essentially unchanged when up to four alanine residues are changed to valine. However, stability increases as the number of these substitutions is increased. Anesthetic binding affinities are also affected. Halothane binds to the four-alpha-helix bundle variants with 0, 1, and 2 substitutions with equivalent affinities but binds to the variants with 3 and 4 more tightly. The same order of binding affinities was observed for chloroform, although for a particular variant, chloroform was bound less tightly. The observed differences in binding affinities may be explained in terms of a modulation of van der Waals and hydrophobic interactions between ligand and receptor. These, in turn, could result from increased four-alpha-helix bundle binding cavity hydrophobicity, a decrease in cavity size, or improved ligand/receptor shape complementarity.     

Molecular dynamics of the local anesthetic tetracaine in phospholipid vesicles

Upon introduction into phosphatidylcholine vesicles, the 13C magnetic resonance peaks of the aromatic resonances of tetracaine are broadened while the T1 relaxation times show little change. Addition of tetracaine to vesicles containing 30% cholesterol produces a similar broadening in the 13C NMR spectrum of tetracaine. Nuclear magnetic resonance parameters of phosphatidylcholine in vesicles which are unchanged by the addition of equimolar tetracaine include 13C T1 relaxation time and 31P linewidth, T1 relaxation time, and nuclear Overhauser effect enhancement. These results are interpreted as indicating a hydrophobic interaction between hydrocarbon portions of the anesthetic and phospholipid bilayer. The rotational correlation time of tetracaine about its long axis in the vesicles has been calculated from the 13C NMR spin lattice relaxation times to be about 10(-10.3) s and is unchanged by incorporation into the phospholipid bilayer. The positively charged ammonium group of tetracaine interacts with the negatively charged phosphate group of the vesicle lipids. Using shift reagents and 31P NMR, tetracaine has been shown to displace cations from the bilayer surface, and does not undergo fast flip-flop across the vesicle bilayer.     

Inhibition of sperm motility by the volatile anesthetic halothane

The inhalational anesthetic halothane reversibly inhibits the motility of sea urchin sperm dose-dependently at concentrations up to 5 mM. Experiments with Triton X-100 extracted, trypsinized axonemes showed that halothane has no effect on the rate of axonemal disintegration in the presence of ATP. These results suggest that halothane inhibits flagellar activity by acting at a site other than the dynein ATPase component of the flagellum.     

Molecular dynamics simulation of silicate glasses

Summary Although the structure of glasses is not really accessible by experimental methods, molecular dynamics is a very useful alternative, as we have tried to demonstrate in this chapter. The simulations reproduce the broad macroscopic features found in these glasses, both structural and transport-related, providing a basis for the more detailed atomic scale features found in the simulated structures. An understanding of important aspects of alkali ion transport, such as the mixed alkali effect and anomalous behaviour in some alumino-silicates, can thus be approached from the atomistic pictures of the glasses produced by the simulations. Although there is room for improvements to the potential models available, it should be clear that the further application of computer simulation methods, such as molecular dynamics, promises to provide much needed advances in glass science and engineering.     

Molecular dynamics simulation and steered molecular dynamics simulation on irisin dimers

Irisin is found closely associated with promoting the browning of beige fat cells in white adipose tissue. The crystal structure reveals that irisin forms a continuous inter-subunit β-sheet dimer. Here, molecular dynamics (MD) simulation and steered molecular dynamics (SMD) simulation were performed to investigate the dissociation process and the intricate interactions between the two irisin monomers. In the process of MD, the interactions between the monomers were roughly analyzed through the average numbers of both hydrophobic contacts and H-bonds. Then, SMD was performed to investigate the accurate interaction energy between the monomers. By the analysis of dissociation energy, the van der Waals (vdW) force was identified as the major energy to maintain the dimer structure, which also verified the results of MD simulation. Meanwhile, 11 essential residues were discovered by the magnitude of rupture force during dissociation. Among them, residues Arg75, Glu79, Ile77, Ala88, and Trp90 were reported in a previous study using the method of mutagenesis and size exclusion chromatography, and several new important residues (Arg72, Leu74, Phe76, Gln78, Val80, and Asp91) were also identified. Interestingly, the new important residues that we discovered and the important residues that were reported are located in the opposite side of the β-sheet of the dimer.     

Molecular dynamics simulations of local anesthetic articaine in a lipid bilayer

In order to investigate structural and dynamical properties of local anesthetic articaine in a model lipid bilayer, a series of molecular dynamics simulations have been performed. Simulations were carried out for neutral and charged (protonated) forms of articaine inserted in fully hydrated dimyristoylphosphatidylcholine (DMPC) lipid bilayer. For comparison purpose, a fully hydrated DMPC bilayer without articaine was also simulated. The length of each simulation was 200 ns. Various properties of the lipid bilayer systems in the presence of both charged and uncharged forms of articaine taken at two different concentrations have been examined: membrane area per lipid, mass density distributions, order parameters, radial distribution functions, head group tilt, diffusion coefficients, electrostatic potential, etc, and compared with results of previous simulations of DMPC bilayer in the presence of lidocaine. It was shown that addition of both charged and neutral forms of articaine causes increase of the dipole electrostatic potential in the membrane interior.     

Molecular dynamics simulation of a phosphatidylglycerol membrane

Although molecular dynamics simulations are an important tool for studying membrane systems, relatively few simulations have used anionic lipids. This paper reports the first simulation of a pure phosphatidylglycerol (PG) bilayer. The properties of this equilibrated palmitoyloleoylphosphatidylglycerol membrane agree with experimental observations of PG membranes and with previous simulations of monolayers and mixed bilayers containing PG lipids. These simulations also provide interesting insights into hydrogen bonding interactions in PG membranes. This equilibrated membrane will be a useful starting point for simulations of membrane proteins interacting with PG lipids.     

Molecular dynamics simulation of a phospholipid membrane

We present the results of molecular dynamics (MD) simulations of a phospholipid membrane in water, including full atomic detail. The goal of the simulations was twofold: first we wanted to set up a simulation system which is able to reproduce experimental results and can serve as a model membrane in future simulations. This goal being reached it is then further possible to gain insight in to those properties that are experimentally more difficult to access. The system studied is dipalmitoylphosphatidylcholine/water, consisting of 5408 atoms. Using original force field parameters the membrane turned out to approach a gel-like state. With slight changes of the parameters, the system adopted a liquid-crystalline state. Separate 80 ps runs were performed on both the gel and liquid-crystalline systems. Comparison of MD results with reliable experimental data (bilayer repeat distance, surface area per lipid, tail order parameters, atom distributions) showed that our simulations, especially the one in the liquid-crystalline phase, can serve as a realistic model for a phospholipid membrane. Further analysis of the trajectories revealed valuable information on various properties. In the liquid-crystalline phase, the interface turns out to be quite diffuse, with water molecules penetrating into the bilayer to the position of the carbonyl groups. The 10–90% width of the interface turns out to be 1.3 nm and the width of the hydrocarbon interior 3.0 nm. The headgroup dipoles are oriented at a small angle with respect to the bilayer plane. The resulting charge distribution is almost completely cancelled by the water molecules. The electron density distribution shows a large dip in the middle of the membrane. In this part the tails are more flexible. The mean life time between dihedral transitions is 20 ps. The average number of gauche angles per tail is 3.5. The occurrence of kinks is not a significant feature.Abbreviations MD molecular dynamics - DPPC dipalmitoylphosphatidylcholine - SPC simple point charges - DPPE dipalmitoylphosphatidylethanolamine Correspondence to: H. J. C. Berendsen     

Molecular dynamics simulation of shell-symmetric Pd nanoclusters

The melting behaviour of Palladium (Pd) isolated shell-symmetric cubooctahedron and icosahedron nanoclusters, both consisting of 309 atoms, were simulated by Molecular Dynamics (MD) simulation, using the Sutton-Chen many-body potential (SC) for the interaction between the Pd atoms. The thermal, structural and dynamic properties were calculated for the Pd nanoclusters along the heating process. The cubooctahedron nanocluster melts around 1040 K, much lower than the melting point of bulk Pd system (1828.05 K). The icosahedron nanocluster melts around 1070 K. Furthermore, structural and dynamic evidence is found for the pre-melting of the nanoclusters. The outer two shells of the shell-symmetric nanocluster melt prior to their homogeneous melting of the whole nanoclusters.     

Molecular dynamics simulation of truncated bovine adrenodoxin

The 128 amino acid long soluble protein adrenodoxin (Adx) is a typical member of the ferredoxin protein family that are electron carrier proteins with an iron-sulfur cofactor. Adx carries electrons from adrenodoxin reductase (AdR) to cytochrome P450s. Its binding modes to these proteins were previously characterized by site-directed mutagenesis, by X-ray crystallography for the complex Adx:AdR, and by NMR. However, no clear evidence has been provided for the driving force that promotes Adx detachment from AdR upon reduction. Here, we characterized the conformational dynamics of unbound Adx in the oxidized and reduced forms using 2-20 ns long molecular dynamics simulations. The most noticeable difference between both forms is the enhanced flexibility of the loop (47-51) surrounding the iron-sulfur cluster in the reduced form. Together with several structural displacements at the binding interface, this increased flexibility may be the key factor promoting unbinding of reduced Adx from AdR. This points to an intrinsic property of reduced Adx that drives dissociation.     

Molecular dynamics simulation of hydration in myoglobin

This study was carried out to evaluate the stability of the 89 bound water molecules that were observed in the neutron diffraction study of CO myoglobin. The myoglobin structure derived from the neutron analysis was used as the starting point in the molecular dynamics simulation using the software package CHARMM. After solvation of the protein, energy minimization and equilibration of the system, 50 ps of Newtonian dynamics was performed. This data showed that only 4 water molecules are continously bound during the length of this simulation while the other solvent molecules exhibit considerable mobility and are breaking and reforming hydrogen bonds with the protein. At any instant during the simulation, 73 of the hydration sites observed in the neutron structure are occupied by water. © 1995 Wiley-Liss, Inc.     

Molecular dynamics simulation of intrinsically disordered proteins

Intrinsically disordered proteins are biomolecules that do not have a definite 3D structure; therefore, their dynamical simulation cannot start from a known list of atomistic positions, such as a Protein Data Bank file. We describe a method to start a computer simulation of these proteins. The first step of the procedure is the creation of a multi-rod configuration of the molecule, derived from its primary sequence. This structure is dynamically evolved in vacuo until its gyration radius reaches the experimental average value; at this point solvent molecules, in explicit or implicit implementation, are added to the protein and a regular molecular dynamics simulation follows. We have applied this procedure to the simulation of tau, one of the largest totally disordered proteins.     

Molecular dynamics simulation of the tubulin dimer with cytostatics

Colchicin, podophylotoxin and indibulin are natural cytostatics that are used in the treatment of neoplasms. But applications of those compounds are rather restricted due to the high toxicity and low specificity. It seems very promising to investigate possibility to design new analogs of the above mentioned drugs that will possess higher cytostatic activity and less toxicity. For this purpose we see computer modeling experiments of tubulin and above mentioned compounds interaction as a powerful approach to design new artificial cytostatics with desired properties. In the current study the CHARMM software of protein-ligand interaction molecular dynamics method has been used. Particularly the following strategy has been applied. Molecules of the cytostatits have been positioned at several random positions around binding sites of tubulin and after energy minimization several binding sites have been identified on the tubulin macromolecule. In these binding sites structural changes that may be responsible for tubulin polymerization have been detected.     

Molecular dynamics simulation by atomic mass weighting

A molecular dynamics-based simulation method in which atomic masses are weighted is described. Results from this method showed that the capability for conformation search in molecular dynamics simulation of a short peptide (FMRF-amide) is significantly increased by mass weighting.     

Molecular dynamics simulation of the docking of substrates to proteins

A simple method is described to perform docking of subtrates to proteins or probes to receptor molecules by a modification of molecular dynamics simulations. The method consists of a separation of the center-of-mass motion of the substrate from its internal and rotational motions, and a separate coupling to different thermal baths for both types of motion of the substrate and for the motion of the receptor. Thus the temperatures and the time constants of coupling to the baths can be arbitrarily varied for these three types of motion, allowing either a frozen or a flexible receptor and allowing control of search rate without disturbance of internal structure. In addition, an extra repulsive term between substrate and protein was applied to smooth the interaction. The method was applied to a model substrate docking onto a model surface, and to the docking of phosphocholine onto immunoglobulin McPC603, in both cases with a frozen receptor. Using transrational temperatures of the substrate in the range of 1300–1700 K and room temperature for the internal degrees of freedom of the substrate, an efficient nontrapping exploratory search (“helicopter view”) is obtained, which visits the correct binding sites. Low energy conformations can then be further investigated by separate search or by dynamic simulated annealing. In both cases the correct minima were identified. The possibility to work with flexible receptors is discussed. © 1994 Wiley-Liss, Inc.     

Molecular dynamics simulation of amyloid beta dimer formation

Recent experiments with amyloid beta (Abeta) peptide indicate that formation of toxic oligomers may be an important contribution to the onset of Alzheimer's disease. The toxicity of Abeta oligomers depends on their structure, which is governed by assembly dynamics. Due to limitations of current experimental techniques, a detailed knowledge of oligomer structure at the atomic level is missing. We introduce a molecular dynamics approach to study Abeta dimer formation. 1), We use discrete molecular dynamics simulations of a coarse-grained model to identify a variety of dimer conformations; and 2), we employ all-atom molecular mechanics simulations to estimate thermodynamic stability of all dimer conformations. Our simulations of a coarse-grained Abeta peptide model predicts 10 different planar beta-strand dimer conformations. We then estimate the free energies of all dimer conformations in all-atom molecular mechanics simulations with explicit water. We compare the free energies of Abeta(1-42) and Abeta(1-40) dimers. We find that 1), dimer conformations have higher free energies compared to their corresponding monomeric states; and 2), the free-energy difference between the Abeta(1-42) and the corresponding Abeta(1-40) dimer conformation is not significant. Our results suggest that Abeta oligomerization is not accompanied by the formation of thermodynamically stable planar beta-strand dimers.