Search PubMed⌕ Search

SEARCH · Search PubMed

Results for “MD simulation”

Search indexed PubMed citations on genomics, clinical trials, systematic reviews and public health. Explore titles, authors and supplied subject terms, then open the PubMed record.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 1,063 records · Page 59Linked to original sources
PubMed1999

Comparison of the solution conformations of a cell-adhesive peptide LBE and its reverse sequence EBL.

T-cell adhesion is mediated by an ICAM-1/LFA-1 interaction; this interaction plays a crucial role in T-cell activation during immune response. LBE peptide, which is derived from the beta-subunit of LFA-1, has been shown to inhibit ICAM-1/LFA-1-mediated T-cell adhesion. In this work, we studied the solution conformations of LBE peptide and its reverse sequence (EBL) by NMR, CD and molecular dynamics simulations. Reverse peptides have been used as controls in biological studies. The effect of reversing the sequence of LBE to EBL peptides on their respective conformations is important in understanding their biological properties in vitro or in vivo. The NMR studies for these peptides were carried out in water and in TFE/water solvent systems. In 40% TFE/water, both peptides exhibited helical conformation. CD studies suggested that the LBE exhibits 30% helical conformation, while the EBL exhibits 20% helical conformation. From the NMR and MD simulation studies, it was evident that the peptides exhibited a stable helical conformation; a stable helical structure was found at Leu6 to Leu15 for LBE and at Gly9 to Leu17 for EBL. The helical conformations of LBE and EBL may be in equilibrium with other possible conformers; the other conformers contain loop and turn structures. Both peptides bind to divalent cations because the LBE is derived from the cation-binding region of the LFA-1. This study shows that reversing the peptide sequence did not alter the secondary structure of the corresponding sequence. Hence, caution must be exercised when using reverse peptides as controls in biological studies. This report will improve our ability to design a better inhibitor of ICAM-1/LFA-1 interaction.

Amino Acid Sequence↗
PubMed2006

Molecular dynamics simulation studies of a protein-RNA complex with a selectively modified binding interface.

The RNA recognition motif (RRM) is one of the most common RNA binding domains. We have investigated the contribution of three highly conserved aromatic amino acids to RNA binding by the N-terminal RRM of the U1A protein. Recently, we synthesized a modified base (A-4CPh) in which a phenyl group is tethered to adenine using a linker of 4 methylene groups. The substitution of this base for adenine in the target RNA selectively stabilizes the complex formed with a U1A protein in which one of the conserved aromatic amino acids is replaced with Ala (Phe56Ala). In this article, we report molecular dynamics (MD) simulations that probe the structural consequences of the substitution of A-4CPh for adenine in the wild type and Phe56Ala U1A-RNA complexes and in the free RNA. The simulations suggest that A-4CPh stabilizes the complex formed with Phe56Ala by adopting a folded conformation in which the tethered phenyl group fills the site occupied by the phenyl group of Phe56 in the wild-type complex. In contrast, an extended conformation of A-4CPh is predicted to be most stable in the complex formed with the wild-type protein. The calculations indicate A-4CPh is in an extended conformation in the free RNA. Therefore, preorganizing the structure of the phenyl-tethered base for binding may improve both the affinity and specificity of the RNA containing A-4CPh for the Phe56Ala U1A protein. Taken together, the previous experimental work and the calculations reported here suggest a general design strategy for altering RRM-RNA complex stability.

Computer Simulation↗
PubMed2001

Molecular dynamics and binding specificity analysis of the bovine immunodeficiency virus BIV Tat-TAR complex.

We have performed molecular dynamics (MD) simulations, with particle-mesh Ewald, explicit waters, and counterions, and binding specificity analyses using combined molecular mechanics and continuum solvent (MM-PBSA) on the bovine immunodeficiency virus (BIV) Tat peptide-TAR RNA complex. The solution structure for the complex was solved independently by Patel and co-workers and Puglisi and co-workers. We investigated the differences in both structures and trajectories, particularly in the formation of the U-A-U base triple, the dynamic flexibility of the Tat peptide, and the interactions at the binding interface. We observed a decrease in RMSD in comparing the final average RNA structures and initial RNA structures of both trajectories, which suggests the convergence of the RNA structures to a MD equilibrated RNA structure. We also calculated the relative binding of different Tat peptide mutants to TAR RNA and found qualitative agreement with experimental studies.

Amino Acid Sequence↗
PubMed2002

A novel class of peptides with facilitating action on neuronal nicotinic receptors of rat chromaffin cells in vitro: functional and molecular dynamics studies.

Peptides related to the N-terminal region of calcitonin gene-related peptide (CGRP) were tested for their ability to modulate neuronal nicotinic acetylcholine receptors (nAChRs) of rat cultured chromaffin cells under whole cell patch-clamp conditions. Although CGRP(1-7) and CGRP(2-7) depressed responses mediated by nAChRs, CGRP(1-6), CGRP(1-5), or CGRP(1-4) rapidly and reversibly potentiated submaximal nicotine currents while sparing maximal currents. CGRP(1-3) was inactive. The threshold concentration for the enhancing effect of CGRP(1-6) was 0.1 microM. CGRP(1-5) or CGRP(1-4) were less effective than CGRP(1-6). Coapplication of CGRP(1-6) and of the allosteric potentiator physostigmine (0.5 microM) gave additive effects on nicotine currents. CGRP(1-6) did not enhance responses generated by muscle-type nicotinic receptors of cultured myoblasts or by gamma-aminobutyric acid(A) receptors expressed by human embryonic kidney cells. Molecular dynamics (MD) simulations suggested that CGRP(1-7) exhibited a relatively rigid ring structure imparted by the disulfide bridge between Cys(2) and Cys(7). The circular dichroism (CD) spectrum recorded from the same peptide was in agreement with this result. Shorter peptides, missing such a bridge, exhibited propensity for alpha-helix configuration. Replacing Cys(7) with Ala yielded CGRP(1-7A), a fragment with partial alpha-helix structure and ability to enhance nicotine currents. CD measurements on CGRP(1-6) were compatible with these MD structural findings. Short terminal fragments of CGRP represent a novel class of substances with selective, rapid, and reversible potentiation of nAChRs.

Allosteric Regulation↗
PubMed2000

Transmembrane peptide NB of influenza B: a simulation, structure, and conductance study.

The putative transmembrane segment of the ion channel forming peptide NB from influenza B was synthesized by standard solid-phase peptide synthesis. Insertion into the planar lipid bilayer revealed ion channel activity with conductance levels of 20, 61, 107, and 142 pS in a 0.5 M KCl buffer solution. In addition, levels at -100 mV show conductances of 251 and 413 pS. A linear current-voltage relation reveals a voltage-independent channel formation. In methanol and in vesicles the peptide appears to adopt an alpha-helical-like structure. Computational models of alpha-helix bundles using N = 4, 5, and 6 NB peptides per bundle revealed water-filled pores after 1 ns of MD simulation in a solvated lipid bilayer. Calculated conductance values [using HOLE (Smart et al. (1997) Biophys. J. 72, 1109-1126)] of ca. 20, 60, and 90 pS, respectively, suggested that the multiple conductance levels seen experimentally must correspond to different degrees of oligomerization of the peptide to form channels.

Amino Acid Sequence↗
PubMed2005

Induced fit and the entropy of structural adaptation in the complexation of CAP and lambda-repressor with cognate DNA sequences.

Molecular dynamics (MD) simulations of 5 ns on protein-DNA complexes of catabolite-activator protein (CAP), lambda-repressor, and their corresponding uncomplexed protein and DNA, are reported. These cases represent two extremes of DNA bending, with CAP DNA bent severely and the lambda-operator nearly straight when complexed with protein. The calculations were performed using the AMBER suite of programs and the parm94 force field, validated for these studies by good agreement with experimental nuclear magnetic resonance data on DNA. An explicit computational model of structural adaptation and computation of the quasiharmonic entropy of association were obtained from the MD. The results indicate that, with respect to canonical B-form DNA, the extreme bending of the DNA in the complex with CAP is approximately 60% protein-induced and 40% intrinsic to the sequence-dependent structure of the free oligomer. The DNA in the complex is an energetically strained form, and the MD results are consistent with a conformational-capture mechanism. The calculated quasiharmonic entropy change accounts for the entropy difference between the two cases. The calculated entropy was decomposed into contributions from protein adaptation, DNA adaptation, and protein-DNA structural correlations. The origin of the entropy difference between CAP and lambda-repressor complexation arises more from the additional protein adaptation in the case of lambda, than to DNA bending and entropy contribution from DNA bending. The entropy arising from protein DNA cross-correlations, a contribution not previously discussed, is surprisingly large.

Allosteric Site↗
PubMed2002

Electrophoresis of an end-labeled polymer chain: a molecular dynamics study.

We study the conformational and the dynamic properties of an end-labeled (telechelic) polymer chain embedded in a porous medium made of randomly distributed immobile spherical obstacles using a stochastic molecular dynamics (MD) simulation method for several obstacle densities rho(imp) and for various field intensities Fx applied only to one end of the chain. For F(x)=0, the chain initially shrinks with increasing density of the obstacles rho(imp). In general, for small Fx and low rho(imp), the chain elongates along the direction of the force and shrinks in the transverse direction whereby this effect becomes more pronounced at larger chain lengths. However, we notice that for moderate values of rho(imp) and Fx, the conformational properties exhibit extrema before reaching a saturation at larger values of rho(imp) and Fx. Likewise, we also find that the drift velocity V(d) of the center of mass of the chain is a nonmonotonic function of the field intensity in the sense that V(d) also exhibits a maximum at a critical value of the field intensity F(crit)(x) beyond which it decreases. Our MD results indicate that for large rho(imp) the chain still can be described by a self-avoiding random walk, which contradicts the prediction of variational calculation using the replica trick, but supports a more recent analytical result using the optimal fluctuation method, as well as a Monte Carlo simulation result for a slightly different disordered medium.

Journal Article↗
PubMed2005

Quantum chemical investigation of hyperfine coupling constants on first coordination sphere water molecule of gadolinium(III) aqua complexes.

Hyperfine interactions (HFI) on the nuclei of the first coordination sphere water molecules in a model [Gd(H(2)O)(8)](3+) aqua complex and in the magnetic resonance imaging contrast agent [Gd(DOTA)(H(2)O)](-) were studied theoretically. Density functional theory (DFT) calculations combined with classical molecular dynamics (MD) simulations have been used in order to take into account dynamic effects in aqueous solution. DFT relativistic calculations show a strong spin-polarization of the first coordination sphere water molecules. This spin-polarization leads to a positive (17)O isotropic hyperfine coupling constant (A(iso)((17)O) = 0.58 +/- 0.11 MHz) and to a significant increase of the effective distance (r(eff)(Gd-O) = 2.72 +/- 0.06 A) of dipolar interaction compared to the mean internuclear distance (r(Gd-O) = 2.56 +/- 0.06 A) obtained from the MD trajectory of [Gd(DOTA)(H(2)O)](-) in aqueous solution. The point-dipole model for anisotropic hyperfine interaction overestimates therefore the longitudinal relaxation rate of the (17)O nucleus by approximately 45%. The (1)H isotropic hyperfine coupling constant of the bound water molecule is predicted to be very small (A(iso)((1)H) = 0.03 +/- 0.02 MHz), and the point-dipole approximation for first coordination sphere water protons holds. The calculated hyperfine parameters are in good agreement with available experimental data.

Journal Article↗
PubMed2005

Ion solvation in water from molecular dynamics simulation with the ABEEM/MM force field.

A systematic study on monovalent ions in water clusters and in aqueous solution is presented for providing insight into their solvation structures, charge distributions, binding energies, as well as dynamic and thermodynamic properties in terms of the atom-bond electronegativity equalization method fused into molecular mechanics (ABEEM/MM) that is to take ABEEM charges into the Coulomb term in MM. For hydrated systems of M+(H2O)n, M+ being Li+, Na+, and K+, as well as X-(H2O)n, X- being F-, Cl-, and Br-, with n = 1-6, parameters for the effective interaction between the ion and the water molecules were determined, so as to reproduce the experimental or ab initio results. The corresponding parameters were tested with molecular dynamics (MD) simulations of these ions in liquid water and with solvation free energy calculations using the perturbation technique. The results of aqueous ionic solution simulations with the ABEEM/MM force field provide a reasonable description of many important properties, which are in good agreement with the experimental measurements. This work demonstrates that the combination of ABEEM/MM-MD provides a powerful tool in analyzing solvation processes of monovalent ions in water.

Chemistry, Physical↗
PubMed2002

Comparative molecular dynamics simulations of amphotericin B-cholesterol/ergosterol membrane channels.

Amphotericin B (AmB) is a very effective anti-fungal polyene macrolide antibiotic whose usage is limited by its toxicity. Lack of a complete understanding of AmB's molecular mechanism has impeded attempts to design less toxic AmB derivatives. The antibiotic is known to interact with sterols present in the cell membrane to form ion channels that disrupt membrane function. The slightly higher affinity of AmB toward ergosterol (dominant sterol in fungal cells) than cholesterol (mammalian sterol) is regarded as the most essential factor on which antifungal chemotherapy is based. To study these differences at the molecular level, two realistic model membrane channels containing molecules of AmB, sterol (cholesterol or ergosterol), phospholipid, and water were studied by molecular dynamics (MD) simulations. Comparative analysis of the simulation data revealed that the sterol type has noticeable effect on the properties of AmB membrane channels. In addition to having a larger size, the AmB channel in the ergosterol-containing membrane has a more pronounced pattern of intermolecular hydrogen bonds. The interaction between the antibiotic and ergosterol is more specific than between the antibiotic and cholesterol. These observed differences suggest that the channel in the ergosterol-containing membrane is more stable and, due to its larger size, would have a higher ion conductance. These observations are in agreement with experiments.

Amphotericin B↗
PubMed2006

Roles of conserved basic amino acid residues and activation mechanism of the hyperthermophilic aspartate racemase at high temperature.

X-ray crystallography has revealed two similar alpha/beta domains of the aspartate racemase from the hyperthermophilic archaeon, Pyrococcus horikoshii OT3. The active site is located in the cleft between the two domains where two cysteine residues face each other. This arrangement allows the substrate to enter the cleft and enables the two cysteine residues to act synergistically. However, the distance between their thiolates was estimated to be 9.6 angstroms, which is beyond the distance for cooperative action of them. We examined the molecular mechanism for the racemization reaction of this hyperthermophilic aspartate racemase by mutational analyses and molecular dynamics simulations. The mutational analyses revealed that Arg48 and Lys164 were essential for catalysis in addition to the putative catalytic cysteine residues. The molecular dynamics simulations revealed that the distance between the two active gamma-sulfur atoms of cysteine residues oscillate to periodically become shorter than the predicted cooperative distance at high temperature. In addition, the conformation of Tyr160, which is located at the entrance of the cleft and inhibits the entry of a substrate, changes periodically to open the entrance at 375 K. The opening of the gate is likely to be induced by the motion of the adjacent amino acid, Lys164. The entrance of an aspartate molecule was observed by molecular dynamics (MD) simulations driven by the force of the electrostatic interaction with Arg48, Lys164, and also Asp47. These results provide insights into the roles of amino acid residues at the catalytic site and also the activation mechanism of a hyperthermophilic aspartate racemase at high temperature.

Amino Acid Isomerases↗
PubMed2005

Computational studies of the farnesyltransferase ternary complex part I: substrate binding.

Farnesyltransferase (FTase) catalyzes the transfer of farnesyl from farnesyl diphosphate (FPP) to cysteine residues at or near the C-terminus of protein acceptors with a CaaX motif (a, aliphatic; X, Met). Farnesylation is a critical modification to many switch proteins involved in the extracellular signal transduction pathway, which facilitates their fixation on the cell membrane where the extracellular signal is processed. The malfunction caused by mutations in these proteins often causes uncontrolled cell reproduction and leads to tumor formation. FTase inhibitors have been extensively studied as potential anticancer agents in recent years, several of which have advanced to different phases of clinical trials. However, the mechanism of this biologically important enzyme has not been firmly established. Understanding how FTase recruits the FPP substrate is the first and foremost step toward further mechanistic investigations and the design of effective FTase inhibitors. Molecular dynamic simulations were carried out on the ternary structure of FTase complexed with the FPP substrate and an acetyl-capped tetrapeptide (acetyl-CVIM), which revealed that the FPP substrate maintains an inactive conformation and the binding of the diphosphate group can be largely attributed to residues R291beta, K164alpha, K294beta, and H248beta. The FPP substrate assumes an extended conformation in the binding site with restricted rotation of the backbone dihedral angles; however, it does not have a well-defined conformation when unbound in solution. This is evident from multinanosecond MD simulations of the FPP substrate in a vacuum and solution. Our conclusion is further supported by theoretical J coupling calculations. Our results on the FPP binding are in good agreement with previous experimental kinetic studies on FTase mutants. The hypothetical conformational activation of the FPP substrate is currently under investigation.

Farnesyltranstransferase↗
PubMed1999

Peptide-bilayer interactions: simulations of dermaseptin B, an antimicrobial peptide.

Dermaseptins, a family of antimicrobial peptides, are believed to act by forming amphipathic alpha-helices which associate with the cell membrane, leading to its permeabilisation and disruption. A simple mean field method is described for simulation of the interactions of peptides with lipid bilayers which includes an approximate representation of the electrostatic effects of the head-group region of the bilayer. Starting from an atomistic model of a PC phospholipid bilayer we calculate an average electrostatic potential along the bilayer normal. By combining the interaction of the peptide with this electrostatic potential and with the hydrophobic core of the membrane we arrive at a more complete description of peptide-bilayer energetics than would be obtained using sidechain hydrophobicities alone. Using this interaction potential in MD simulations of the frog skin peptide dermaseptin B reveals that the lipid bilayer stabilises the alpha-helical conformation of the peptide. This is in agreement with FTIR data. A surface associated orientation thus appears to be the most stable arrangement of the peptide, at least at zero ionic strength and without taking account of possible peptide-peptide interactions.

Algorithms↗
PubMed2006

Shear dynamics of hydration layers.

Molecular dynamics (MD) simulations have been performed to investigate the shear dynamics of hydration layers of the thickness of D=0.61-2.44 nm confined between two mica surfaces. Emphases are placed on the external shear response and internal relaxation properties of aqueous films. For D=0.92-2.44 nm liquid phase, the shear responses are fluidic and similar to those observed in surface force balance experiments [U. Raviv and J. Klein, Science 297, 1540 (2002)]. However, for the bilayer ice (D=0.61 nm) [Y. S. Leng and P. T. Cummings, J. Chem. Phys. 124, 74711 (2006)] significant shear enhancement and shear thinning over a wide range of shear rates in MD regime are observed. The rotational relaxation time of water molecules in this bilayer ice is found to be as high as 0.017 ms (10(-5) s). Extrapolating the shear rate to the inverse of this longest relaxation time, we obtain a very high shear viscosity for the bilayer ice, which is also observed quite recently for D< or =0.6+/-0.3 nm hydration layers [H. Sakuma et al., Phys. Rev. Lett. 96, 46104 (2006)]. We further investigate the boundary slip of water molecules and hydrated K(+) ions and concluded that no-slip boundary condition should hold for aqueous salt solution under extreme confinement between hydrophilic mica surfaces, provided that the confined film is of Newtonian fluid.

Aluminum Silicates↗
PubMed2007

Conformational changes upon calcium binding and phosphorylation in a synthetic fragment of calmodulin.

We have recently investigated by far-UV circular dichroism (CD) the effects of Ca(2+) binding and the phosphorylation of Ser 81 for the synthetic peptide CaM [54-106] encompassing the Ca(2+)-binding loops II and III and the central alpha helix of calmodulin (CaM) (Arrigoni et al., Biochemistry 2004, 43, 12788-12798). Using computational methods, we studied the changes in the secondary structure implied by these spectra with the aim to investigate the effect of Ca(2+) binding and the functional role of the phosphorylation of Ser 81 in the action of the full-length CaM. Ca(2+) binding induces the nucleation of helical structure by inducing side chain stacking of hydrophobic residues. We further investigated the effect of Ca(2+) binding by using near-UV CD spectroscopy. Molecular dynamics simulations of different fragments containing the central alpha-helix of CaM using various experimentally determined structures of CaM with bound Ca(2+) disclose the structural effects provided by the phosphorylation of Ser 81. This post-translational modification is predicted to alter the secondary structure in its surrounding and also to hinder the physiological bending of the central helix of CaM through an alteration of the hydrogen bond network established by the side chain of residue 81. Using quantum mechanical methods to predict the CD spectra for the frames obtained during the MD simulations, we are able to reproduce the relative experimental intensities in the far-UV CD spectra for our peptides. Similar conformational changes that take place in CaM [54-106] upon Ca(2+) binding and phosphorylation may occur in the full-length CaM.

Calcium↗
PubMed2001

Molecular dynamics simulation of n-dodecyl phosphate aggregate structures.

Aggregates of n-dodecyl phosphate present an attractive model system of simple phospholipid amphiphile supramolecular structures for study by molecular dynamics simulation, since these systems have previously been studied experimentally under various conditions. A detailed molecular dynamics description of the properties of planar bilayer membranes (as a model for unilamellar vesicular membranes) and spherical micelles under various simulated conditions is presented. It is shown that the united-atom model of GROMOS96 applying the force-field parameter set 43A2 for biomolecular systems yields properties in agreement with experimental ones in most cases. Hydrogen bonding plays a role in stabilizing the bilayer aggregates at low pH, but not for the micelles, which are energetically favoured at high pH. NMR -S(CD) order parameters for a lipid bilayer system, the diffusion of amphiphiles within aggregates and of counterions, and lifetimes of hydrogen bonds between amphiphiles and to water are estimated from the MD simulations.

Hydrogen Bonding↗
PubMed2003

Dynamical structure functions, collective modes, and energy gap in charged-particle bilayers.

The dynamical properties of strongly coupled charged-particle bilayers are investigated by molecular dynamics (MD) simulation and theoretical analysis. The spectra of the current correlation functions show the existence of two (in-phase and out-of-phase) longitudinal and two (in-phase and out-of-phase) transverse collective modes. The out-of-phase modes possess finite frequencies at wave numbers k-->0, confirming the existence of the predicted long-wavelength energy gap in the bilayer system. A theoretical model based on an extended Feynman ansatz for the dynamical structure functions provides predictions on the strength of the collective modes that are verified by the MD experiment.

Journal Article↗
PubMed2005

The influence of aqueous versus glassy solvents on protein dynamics: vibrational echo experiments and molecular dynamics simulations.

Spectrally resolved infrared stimulated vibrational echo measurements are used to measure the vibrational dephasing of the CO stretching mode of carbonmonoxy-hemoglobin (HbCO), a myoglobin mutant (H64V), and a bacterial cytochrome c(552) mutant (Ht-M61A) in aqueous solution and trehalose glasses. The vibrational dephasing of the heme-bound CO is significantly slower for all three proteins embedded in trehalose glasses compared to that of aqueous protein solutions. All three proteins exhibit persistent but notably slower spectral diffusion when the protein surface is fixed by the glassy solvent. Frequency-frequency correlation functions (FFCFs) of the CO are extracted from the vibrational echo data to reveal that the structural dynamics, as sensed by the CO, of the three proteins in trehalose and aqueous solution are dominated by fast (tens of femtoseconds), motionally narrowed fluctuations. MD simulations of H64V in dynamic and "static" water are presented as models of the aqueous and glassy environments. FFCFs are calculated from the H64V simulations and qualitatively reproduce the important features of the experimentally extracted FFCFs. The suppression of long time scale (picoseconds to tens of picoseconds) frequency fluctuations (spectral diffusion) in the glassy solvent is the result of a damping of atomic displacements throughout the protein structure and is not limited to structural dynamics that occur only at the protein surface. The analysis provides evidence that some dynamics are coupled to the hydration shell of water, supporting the idea that the bioprotection offered by trehalose is due to its ability to immobilize the protein surface through a thin, constrained layer of water.

Computer Simulation↗