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,261 records · Page 70Linked to original sources
PubMed2003

Role of base sequence context in conformational equilibria and nucleotide excision repair of benzo[a]pyrene diol epoxide-adenine adducts.

We investigate the influence of base sequence context on the conformations of the 10S (+)- and 10R (-)-trans-anti-[BP]-N(6)-dA adducts through molecular dynamics (MD) simulations with free energy calculations, and relate the structural findings to results of nucleotide excision repair (NER) assays in human cell extracts. In previous studies, these adducts were studied in the CA*A sequence context, and here we report results for the CA*C sequence. Our simulations indicate that the base sequence context affects the syn-anti conformational equilibrium in the 10S (+) adduct by modulating the barrier heights between these states on the energy surface, with a higher barrier in the CA*C case. Our nucleotide excision repair assay finds greater NER susceptibilities in the 10S (+) adduct for the CA*C sequence context. A structural rationale ties together these results. A sequence specific hydrogen bond, accompanied by a significantly increased roll and consequent bending in the 10S (+) adduct, has been found in our simulations for the CA*C sequence, which could account for the enhanced nucleotide excision repair as well as the syn-anti equilibrium difference we observe in this isomer and sequence. Such sequence specific differential repair could contribute to the existence of mutational hotspots and thereby contribute to the complexity of cancer initiation.

7,8-Dihydro-7,8-dihydroxybenzo(a)pyrene 9,10-oxide↗
PubMed2004

Conformational and molecular modeling studies of beta-cyclodextrin-heptagastrin and the third extracellular loop of the cholecystokinin 2 receptor.

The conformational features of a conjugate of the C-terminus of human gastrin (HG[11-17]), the shortest gastrin sequence retaining biological function, with beta-cyclodextrin ([Nle(15)]-HG[11-17]-betaCD) were determined by NMR spectroscopy in an aqueous solution of dodecylphosphocholine (DPC) micelles. The peptide-betaCD conjugate displays a binding affinity and activation profile comparable to those of HG[11-17] at the cholecysokinin 2 (CCK(2)) receptor, the G protein-coupled receptor responsible for the gastrointestinal function of gastrin. The structure of the peptide consisted of a well-defined beta-turn between Gly(13) and Asp(16) of gastrin. The structural preferences of [Nle(15)]-HG[11-17]-betaCD in DPC micelles and the 5-doxylstearate-induced relaxation of the (1)H NMR resonances support a membrane-associated receptor recognition mechanism. Addition of [Nle(15)]-HG[11-17]-betaCD to the third extracellular loop domain of the CCK(2) receptor, CCK(2)-R(352-379), generated a number of intermolecular nuclear Overhauser enhancements (NOEs) and chemical shift perturbations. NOE-restrained MD simulations of the [Nle(15)]-HG[11-17]-betaCD-CCK(2)-R complex produced a topological orientation in which the C-terminus was located in a shallow hydrophobic pocket near the confluence of TM2 and -3. Despite the steric bulk and physicochemical properties of betaCD, the [Nle(15)]-HG[11-17]-betaCD-CCK(2)-R complex is similar to the CCK-8-CCK(2)-R complex determined previously, providing insight into the mode of ligand binding and the role of electrostatic interactions.

Amino Acid Sequence↗
PubMed2004

How C-terminal carboxyamidation alters the biological activity of peptides from the venom of the eumenine solitary wasp.

Inflammatory peptides display different types of post-transcriptional modifications, such as C-terminal amidation, that alter their biological activity. Here we describe the structural and molecular dynamics features of the mast cell degranulating peptide, eumenine mastoparan-AF (EMP-AF-NH(2)), found in the venom of the solitary wasp, and of its carboxyl-free C-terminal form (EMP-AF-COO(-)) characterized by a reduced activity. Circular dichroism indicates that both peptides switch from a random coil conformation in water to a helical structure in TFE and SDS micelles. NMR data, in 30% TFE, reveal that the two peptides fold into an alpha-helix spanning most of their length, while they differ in terms of molecular rigidity. To understand the origins of the conformational flexibility observed in the case of EMP-AF-COO(-), a 5 ns MD simulation was carried out for each peptide, in an explicit water/TFE environment. The results show that the two peptides differ in an H-bond between Leu14 NH(2) and the backbone carbonyl of Ile11. The loss of that H-bond in EMP-AF-COO(-) leads to a significant modification of its structural dynamics. In fact, as evidenced by essential dynamics analysis, while EMP-AF-NH(2) exists mainly as a rigid structure, EMP-AF-COO(-) presents two helical stretches that fluctuate in some sort of independent fashion. We conclude that the diverse biological activity of the two peptides is not simply due to the reduction of the net positive charge, as generally suggested, but also to a structural perturbation of the amphipathic alpha-helix that affects their ability to perturb the cell membrane.

Amides↗
PubMed2004

Multiple substrate binding states and chiral recognition in cofactor-independent glutamate racemase: a molecular dynamics study.

Glutamate racemase (MurI) catalyzes the racemization of glutamate; two cysteine residues serve as catalytic acid and base. On the basis of the crystal structure of MurI from the hyperthermophilic bacterium Aquifex pyrophilus, we performed molecular dynamics (MD) simulations of six different systems to investigate stereochemistry, substrate ligation, and active site protonation state. The catalytic competence of individual systems was assessed by the abundance of reactive conformers. Only systems in which Cys70 is poised to deprotonate d-Glu were found to be catalytically competent (idem Cys178/l-Glu), in agreement with the experimentally observed stereochemistry of Lactobacillus fermentii MurI [Tanner, M. E. et al. (1993) Biochemistry 32, 3998-4006]. Only systems in which the alpha-amino group of l/d-Glu and the imidazole moiety of His are deprotonated are catalytically competent. The active site of MurI displays an unusual flexibility in substrate ligation, and several transitions between stable binding patterns were observed. In catalytically competent binding states, the conserved threonine residues 72, 114, and 117 ligate the alpha-carboxylate of Glu and the Asn71 amides ligate the alpha-amino group of Glu, whereas the delta-carboxylate of Glu is steered by electrostatic repulsion from the Asp7 and Glu147 side chain carboxylates. A network of hydrogen bonds controls the positioning of each thiol/thiolate. In what we term substrate flipping, Glu suddenly rotates into a binding pattern that resembles the post-racemization state of the other enantiomer, i.e., each enantiomer can be bound in two distinct states. Substrate flipping and unfavorable substrate binding successively trigger dissociation of the substrate, accompanied by an opening of the active site channel. We explain how the weak binding of Glu contributes to catalysis and suggest a mechanism by which binding mismatches are propagated into an opening of the active site.

Amino Acid Isomerases↗
PubMed2005

Effect of lipid composition on the "membrane response" induced by a fusion peptide.

To understand the initial stages of membrane destabilization induced by viral proteins, the factors important for binding of fusion peptides to cell membranes must be identified. In this study, effects of lipid composition on the mode of peptides' binding to membranes are explored via molecular dynamics (MD) simulations of the peptide E5, a water-soluble analogue of influenza hemagglutinin fusion peptide, in two full-atom hydrated lipid bilayers composed of dimyristoyl- and dipalmitoylphosphatidylcholine (DMPC and DPPC, respectively). The results show that, although the peptide has a common folding motif in both systems, it possesses different modes of binding. The peptide inserts obliquely into the DMPC membrane mainly with its N-terminal alpha helix, while in DPPC, the helix lies on the lipid/water interface, almost parallel to the membrane surface. The peptide seriously affects structural and dynamical parameters of surrounding lipids. Thus, it induces local thinning of both bilayers and disordering of acyl chains of lipids in close proximity to the binding site. The "membrane response" significantly depends upon lipid composition: distortions of DMPC bilayer are more pronounced than those in DPPC. Implications of the observed effects to molecular events on initial stages of membrane destabilization induced by fusion peptides are discussed.

Cell Membrane↗
PubMed2006

Temperature dependence of the structure of the substrate and active site of the Thermus thermophilus chorismate mutase E x S complex.

Molecular dynamics (MD) simulations of Thermus thermophilus chorismate mutase substrate complex (TtCM x S) have been carried out at 298 K, 333 K, and the temperature of optimum activity: 343 K. The enzyme exists as trimeric subunits with active sites shared between two neighboring subunits. Two features distinguish intersubunit linkages of the thermophilic and mesophilic enzyme Bacillus subtilis chorismate mutase substrate complex (BsCM x S): (i) electrostatic interactions by intersubunit ion pairs (Arg3-Glu40*/41, Arg76-Glu51* and Arg69*-Asp101, residues labeled with an asterisk are from the neighboring subunit) in the TtCM x S are not present in the structure of the BsCM x S; and (ii) replacement of polar residues with short and nonpolar residues in the interstices of the TtCM x S tighten the intersubunit hydrophobic interactions compared to BsCM x S. Concerning the active site, electrostatic interactions of the critically placed Arg6 and Arg63* with the two carboxylates of chorismate place the latter in a reactive conformation to spontaneously undergo a Claisen rearrangement. The optimum geometry at the active site has the CZ atoms of the two arginines 11 A apart. With a decrease in temperature, Arg63* moves toward Arg6 and the average conformation structure of chorismate moves further away from the reactive ground state conformation. This movement is due to the decrease in distance separating the electrostatic (in the main) and hydrophobic interacting pairs holding the two subunits together.

Amino Acid Substitution↗
PubMed2006

Inhibitor binding to the plasmepsin IV aspartic protease from Plasmodium falciparum.

Plasmepsin IV (Plm IV) is one of the aspartic proteases present in the food vacuole of the malaria parasite Plasmodium falciparum involved in host hemoglobin degradation by the parasite. Using a series of previously synthesized plasmepsin inhibitors [Ersmark, K., et al. (2005) J. Med. Chem. 48, 6090-106], we report here experimental data and theoretical analysis of their inhibitory activity toward Plm IV. All compounds share a 1,2-dihydroxyethylene unit as the transition state mimic. They possess symmetric P1 and P1' side chains and either a diacylhydrazine, a five-membered oxadiazole ring, or a retroamide at the P2 and P2' positions. Experimental binding affinities are compared to those predicted by the linear interaction energy (LIE) method and an empirical scoring function, using both a crystal structure and a homology model for the enzyme. Molecular dynamics (MD) simulations of the modeled complexes allow a rational interpretation of the structural determinants for inhibitor binding. A ligand bearing a P2 and P2' symmetric oxadiazole which is devoid of amide bonds is identified both experimentally and theoretically as the most potent inhibitor of Plm IV. For the P2 and P2' asymmetric compounds, the results are consistent with earlier predictions regarding the mode of binding of this class of inhibitors to Plm II. Theoretical estimation of selectivity for some compounds is also reported. Significant features of the Plm IV binding pocket are discussed in comparison to related enzymes, and the results obtained here should be helpful for further optimization of inhibitors.

Animals↗
PubMed1999

Engineering out motion: a surface disulfide bond alters the mobility of tryptophan 22 in cytochrome b5 as probed by time-resolved fluorescence and 1H NMR experiments.

In the accompanying paper [Storch et al. (1999) Biochemistry 38, 5054-5064] equilibrium denaturation studies and molecular dynamics (MD) simulations were used to investigate localized dynamics on the surface of cytochrome b5 (cyt b5) that result in the formation of a cleft. In those studies, an S18C:R47C disulfide mutant was engineered to inhibit cleft mobility. Temperature- and urea-induced denaturation studies revealed significant differences in Trp 22 fluorescence between the wild-type and mutant proteins. On the basis of the results, it was proposed that wild type populates a conformational ensemble that is unavailable to the disulfide mutant and is mediated by cleft mobility. As a result, the solvent accessibility of Trp 22 is decreased in S18C:R47C, suggesting that the local environment of this residue is less mobile due to the constraining effects of the disulfide on cleft dynamics. To further probe the structural effects on the local environment of Trp 22 caused by inhibition of cleft formation, we report here the results of steady-state and time-resolved fluorescence quenching, differential phase/modulation fluorescence anisotropy, and 1H NMR studies. In Trp fluorescence experiments, the Stern-Volmer quenching constant increases in wild type versus the oxidized disulfide mutant with increasing temperature. At 50 degrees C, KSV is nearly 1.5-fold greater in wild type compared to the oxidized disulfide mutant. In the reduced disulfide mutant, KSV was the same as wild type. The bimolecular collisional quenching constant, kq, for acrylamide quenching of Trp 22 increases 2.7-fold for wild type and only 1.8-fold for S18C:R47C, upon increasing the temperature from 25 to 50 degrees C. The time-resolved anisotropy decay at 25 degrees C was fit to a double-exponential decay for both the wild type and S18C:R47C. Both proteins exhibited a minor contribution from a low-amplitude fast decay, consistent with local motion of Trp 22. This component was more prevalent in the wild type, and the fractional contribution increased significantly upon raising the temperature. The fast rotational component of the S18C:R47C mutant was less sensitive to increasing temperature. A comparison of the 1H NMR monitored temperature titration of the delta-methyl protons of Ile 76 for wild type and oxidized disulfide mutant, S18C:R47C, showed a significantly smaller downfield shift for the mutant protein, suggesting that Trp 22 in the mutant protein experiences comparatively decreased cleft dynamics in core 2 at higher temperatures. Furthermore, comparison of the delta'-methyl protons of Leu 25 in the two proteins revealed a difference in the ratio of the equilibrium heme conformers of 1.2:1 for S18C:R47C versus 1.5:1 for wild type at 40 degrees C. The difference in equilibrium heme orientations between wild type and S18C:R47C suggests that the disulfide bond affects heme binding within core 1, possibly through damped cleft fluctuations. Taken together, the NMR and fluorescence studies support the proposal that an engineered disulfide bond inhibits the formation of a dynamic cleft on the surface of cyt b5.

Amino Acid Substitution↗
PubMed1999

Understanding substrate specificity in human and parasite phosphoribosyltransferases through calculation and experiment.

We present molecular dynamics (MD) simulations on two enzymes: a human hypoxanthine-guanine-phosphoribosyltransferase (HGPRTase) and its analogue in the protozoan parasite Tritrichomonas foetus. The parasite enzyme has an additional ability to process xanthine as a substrate, making it a hypoxanthine-guanine-xanthine phosphoribosyltransferase (HGXPRTase) [Chin, M. S., and Wang, C. C. (1994) Mol. Biochem. Parasitol. 63 (2), 221-229 (1)]. X-ray crystal structures of both enzymes complexed to guanine monoribosyl phosphate (GMP) have been solved, and show only subtle differences in the two active sites [Eads et al. (1994) Cell 78 (2), 325-334 (2); Somoza et al. (1996) Biochemistry 35 (22), 7032-7040 (3)]. Most of the direct contacts with the base region of the substrate are made by the protein backbone, complicating the identification of residues significantly associated with xanthine recognition. Our calculations suggest that the broader specificity of the parasite enzyme is due to a significantly more flexible base-binding region, and rationalize the effect of two mutations, R155E and D163N, that alter substrate specificity [Munagala, N. R., and Wang, C. C. (1998) Biochemistry 37 (47), 16612-16619 (4)]. In addition, our simulations suggested a double mutant (D106E/D163N) that might rescue the D163N mutant. This double mutant was expressed and assayed, and its catalytic activity was confirmed. Our molecular dynamics trajectories were also used with a structure-based design program, Pictorial Representation Of Free Energy Changes (PROFEC), to suggest parasite-selective derivatives of GMP. Our calculations here successfully rationalize the parasite-selectivity of two novel inhibitors derived from the computer-aided design of Somoza et al. (5) and demonstrate the utility of PROFEC in the design of species-selective inhibitors.

Animals↗
PubMedDate not supplied

Mapping the binding site of a large set of quinazoline type EGF-R inhibitors using molecular field analyses and molecular docking studies.

In the current work, three-dimensional QSAR studies for one large set of quinazoline type epidermal growth factor receptor (EGF-R) inhibitors were conducted using two types of molecular field analysis techniques: comparative molecular field analysis (CoMFA) and comparative molecular similarity indices analysis (CoMSIA). These compounds belonging to six different structural classes were randomly divided into a training set of 122 compounds and a test set of 13 compounds. The statistical results showed that the 3D-QSAR models derived from CoMFA were superior to those generated from CoMSIA. The most optimal CoMFA model after region focusing bears significant cross-validated r(2)(cv) of 0.60 and conventional r(2) of 0.92. The predictive power of the best CoMFA model was further validated by the accurate estimation to these compounds in the external test set, and the mean agreement of experimental and predicted log(IC(50)) values of the inhibitors is 0.6 log unit. Separate CoMFA models were conducted to evaluate the influence of different partial charges (Gasteiger-Marsili, Gasteiger-Hückel, MMFF94, ESP-AM1, and MPA-AM1) on the statistical quality of the models. The resulting CoMFA field map provides information on the geometry of the binding site cavity and the relative weights of various properties in different site pockets for each of the substrates considered. Moreover, in the current work, we applied MD simulations combined with MM/PBSA (Molecular mechanics/Possion-Boltzmann Surface Area) to determine the correct binding mode of the best inhibitor for which no ligand-protein crystal structure was present. To proceed, we define the following procedure: three hundred picosecond molecular dynamics simulations were first performed for the four binding modes suggested by DOCK 4.0 and manual docking, and then MM/PBSA was carried out for the collected snapshots. The most favorable binding mode identified by MM/PBSA has a binding free energy about 10 kcal/mol more favorable than the second best one. The most favorable binding mode identified by MM/PBSA can give satisfactory explanation of the SAR data of the studied molecules and is in good agreement with the contour maps of CoMFA. The most favorable binding mode suggests that with the quinazoline-based inhibitor, the N3 atom is hydrogen-bonded to a water molecule which, in turn, interacts with Thr 766, not Thr 830 as proposed by Wissner et al. (J. Med. Chem. 2000, 43, 3244). The predicted complex structure of quinazoline type inhibitor with EGF-R as well as the pharmacophore mapping from CoMFA can interpret the structure activities of the inhibitors well and afford us important information for structure-based drug design.

Binding Sites↗
PubMedDate not supplied

Molecular dynamics integration meets standard theory of molecular vibrations.

An iterative SISM (split integration symplectic method) for molecular dynamics (MD) integration is described. This work explores an alternative for the internal coordinate system prediction in the SISM introduced by JaneZic et al. (J. Chem. Phys. 2005, 122, 174101). The SISM, which employs a standard theory of molecular vibrations, analytically resolves the internal high-frequency molecular vibrations. This is accomplished by introducing a translating and rotating internal coordinate system of a molecule and calculating normal modes of an isolated molecule only. The Eckart frame, which is usually used in the standard theory of molecular vibrations as an internal coordinate system of a molecule, is adopted to be used within the framework of the second order generalized leapfrog scheme. In the presented MD integrator the internal coordinate frame at the end of the integration step is predicted halfway through the integration step using a predictor-corrector type iterative approach thus ensuring the method's time reversibility. The iterative SISM, which is applicable to any system of molecules with one equilibrium configuration, was applied here to perform all-atom MD simulations of liquid CO2 and SO2. The simulation results indicate that for the same level of accuracy, this algorithm allows significantly longer integration time steps than the standard second-order leapfrog Verlet (LFV) method.

Journal Article↗
PubMedDate not supplied

WIN55212-2 docking to the CB1 cannabinoid receptor and multiple pathways for conformational induction.

Key pharmacophoric elements for the (aminoalkyl)indole (AAI) CB1 cannabinoid receptor agonists are the aminoalkyl moiety, the lipophilic aroyl group, and the heterocyclic indole ring. In the present study, the docking space allowed for (R)-[2,3-dihydro-5-methyl-3-[(4-morpholinyl)methyl]pyrrolo[1,2,3-de]-1,4-benzoxazin-6-yl](1-naphthalenyl)methanone (WIN55212-2; 1) within the CB1 receptor was extensively explored by a docking approach that combines Monte Carlo (MC) and molecular dynamics (MD) simulations. The goals were to understand the key binding interactions of AAIs within the CB1 receptor and to examine the role of the ligand in inducing a receptor conformational change. From the findings of extensive SAR studies on the cannabinoid compounds and correlation between AAI binding affinity data and calculated binding energies, we proposed two alternative binding conformations, aroyl-up1 and aroyl-up2. These denote the directionality of the ligand naphthyl ring within the receptor upward with respect to the extracellular side. A comprehensive structural analysis of 1 demonstrated that the aroyl ring moiety could be important as the steric trigger for inducing CB1 receptor conformational change. Thus, it appears that aromatic-aromatic interactions are important not only for the binding of 1 but also for inducing receptor conformational change. It is possible that differences in the nature of the ligand binding could contribute to ligand-specific conformational changes in the receptor.

Benzoxazines↗
PubMed2001

Structural, computational, and (59)Co NMR studies of primary and secondary amine complexes of Co(III) porphyrins.

Four novel low-spin bis(amine) Co(III) porphyrins [Co(TPP)(BzNH(2))(2)](SbF(6)), 1, [Co(TPP)(1-BuNH(2))(2)](SbF(6)), 2, [Co(TPP)(PhCH(2)CH(2)NH(2))(2)](SbF(6)), 3, and [Co(TPP)(1-MePipz)(2)](SbF(6)), 4, have been synthesized and characterized by low-temperature X-ray crystallography, IR, electronic, and NMR ((1)H, (13)C, and (59)Co) spectroscopy. The mean Co-N(p) distance for the four structures is 1.986(1) A. The Co-N(ax) distances for the 1 degrees amine derivatives average to 1.980(5) A; the axial bonds of the 2 degrees amine derivative are significantly longer, averaging 2.040(1) A. The porphyrin core conformation of 4 is significantly nonplanar (mixture of S(4)-ruf and D(2d)-sad distortions) due to a staggered arrangement of the axial ligands over the porphyrin core and meso-phenyl group orientations < 90 degrees. The X-ray structures have been used with the coordinates for [Co(TPP)(Pip)(2)](NO(3)) (Scheidt et al. J. Am. Chem. Soc. 1973, 95, 8289-8294.) to parametrize a molecular mechanics (MM) force field for bis(amine) complexes of Co(III) porphyrins. The calculations show that two types of crystal packing interactions (van der Waals and hydrogen bonding) largely control the crystallographically observed conformations. Gas phase conformational energy surfaces have been computed for these complexes by dihedral angle driving methods and augmented with population distributions calculated by MD simulations at 298 K; the calculations demonstrate that the bis(1 degrees amine) complexes are significantly more flexible than the bis(2 degrees amine) analogues. (59)Co NMR spectra have been acquired for a range of [Co(TPP)(amine)(2)]Cl derivatives as a function of temperature. The (59)Co chemical shifts increase linearly with increasing temperature due to population of thermally excited vibrational levels of the (1)A(1) ground state. Activation energies for molecular reorientation (tumbling) have been determined from an analysis of the (59)Co NMR line widths as a function of 1/T; lower barriers exist for the conformationally rigid 2 degrees amine derivatives (2.6-3.8 kJ mol(-1)). The (59)Co chemical shifts vary linearly with the DFT-calculated radial expectation values (3d) for the Co(III) ion. The correlation leads to the following order for the sigma-donor strengths of the axial ligands: BzNH(2) > or = Cl(-) > 1-BuNH(2) > PhCH(2)CH(2)NH(2) > 1-Bu(2)NH > Et(2)NH. The (59)Co NMR line widths are proportional to the square of the DFT-calculated valence electric field gradient at the Co nucleus. Importantly, this is the first computational rationalization of the (59)Co NMR spectra of Co(III) porphyrins.

Amines↗
PubMed2001

Theoretical studies of d(A:T)-based parallel-stranded DNA duplexes.

Poly d(A:T) parallel-stranded DNA duplexes based on the Hoogsteen and reverse Watson-Crick hydrogen bond pairing are studied by means of extensive molecular dynamics (MD) simulations and molecular mechanics coupled to Poisson-Boltzmann (MM-PB/SA) calculations. The structural, flexibility, and reactivity characteristics of Hoogsteen and reverse Watson-Crick parallel duplexes are described from the analysis of the trajectories. Theoretical calculations show that the two parallel duplexes are less stable than the antiparallel Watson-Crick duplex. The difference in stability between antiparallel and parallel duplexes increases steadily as the length of the duplex increases. The reverse Watson-Crick arrangement is slightly more stable than the Hoogsteen duplex, the difference being also increased linearly with the length of the duplex. A subtle balance of intramolecular and solvation terms is responsible for the preference of a given helical structure.

DNA↗
PubMed2002

Ab initio studies on Al(+)(H(2)O)(n), HAlOH(+)(H(2)O)(n-1), and the size-dependent H(2) elimination reaction.

We report computational studies on Al(+)(H(2)O)(n), and HAlOH(+)(H(2)O)(n-1), n = 6-14, by the density functional theory based ab initio molecular dynamics method, employing a planewave basis set with pseudopotentials, and also by conventional methods with Gaussian basis sets. The mechanism for the intracluster H(2) elimination reaction is explored. First, a new size-dependent insertion reaction for the transformation of Al(+)(H(2)O)(n), into HAlOH(+)(H(2)O)(n-1) is discovered for n > or = 8. This is because of the presence of a fairly stable six-water-ring structure in Al(+)(H(2)O)(n) with 12 members, including the Al(+). This structure promotes acidic dissociation and, for n > or = 8, leads to the insertion reaction. Gaussian based BPW91 and MP2 calculations with 6-31G* and 6-31G** basis sets confirmed the existence of such structures and located the transition structures for the insertion reaction. The calculated transition barrier is 10.0 kcal/mol for n = 9 and 7.1 kcal/mol for n = 8 at the MP2/6-31G** level, with zero-point energy corrections. Second, the experimentally observed size-dependent H(2) elimination reaction is related to the conformation of HAlOH(+)(H(2)O)(n-1), instead of Al(+)(H(2)O)(n). As n increases from 6 to 14, the structure of the HAlOH(+)(H(2)O)(n-1) cluster changes into a caged structure, with the Al-H bond buried inside, and protons produced in acidic dissociation could then travel through the H(2)O network to the vicinity of the Al-H bond and react with the hydride H to produce H(2). The structural transformation is completed at n = 13, coincident approximately with the onset of the H(2) elimination reaction. From constrained ab initio MD simulations, we estimated the free energy barrier for the H(2) elimination reaction to be 0.7 eV (16 kcal/mol) at n = 13, 1.5 eV (35 kcal/mol) at n = 12, and 4.5 eV (100 kcal/mol) at n = 8. The existence of transition structures for the H(2) elimination has also been verified by ab initio calculations at the MP2/6-31G** level. Finally, the switch-off of the H(2) elimination for n > 24 is explored and attributed to the diffusion of protons through enlarged hydrogen bonded H(2)O networks, which reduces the probability of finding a proton near the Al-H bond.

Journal Article↗
PubMed2002

Polyunsaturated fatty acids in lipid bilayers: intrinsic and environmental contributions to their unique physical properties.

Polyunsaturated lipids are an essential component of biological membranes, influencing order and dynamics of lipids, protein-lipid interaction, and membrane transport properties. To gain an atomic level picture of the impact of polyunsaturation on membrane properties, quantum mechanical (QM) and empirical force field based calculations have been undertaken. The QM calculations of the torsional energy surface for rotation about vinyl-methylene bonds reveal low barriers to rotation, indicating an intrinsic propensity toward flexibility. Based on QM and experimental data, empirical force field parameters were developed for polyunsaturated lipids and applied in a 16 ns molecular dynamics (MD) simulation of a 1-stearoyl-2-docosahexaenoyl-sn-glyerco-3-phosphocholine (SDPC) lipid bilayer. The simulation results are in good agreement with experimental data, suggesting an unusually high degree of conformational flexibility of polyunsaturated hydrocarbon chains in membranes. The detailed analysis of chain conformation and dynamics by simulations is aiding the interpretation of experimental data and is useful for understanding the unique role of polyunsaturated lipids in biological membranes. The complete force field is included as Supporting Information and is available from http://www.pharmacy.umaryland.edu/faculty/amackere/research.html.

Computer Simulation↗
PubMed2002

DNA dynamics in a water drop without counterions.

Because of its polyionic character, the DNA double helix is stable and biologically active only in salty aqueous media where its charge is compensated by solvent counterions. Monovalent metal ions are ubiquitous in DNA environment, and they are usually considered as the possible driving force of sequence-dependent modulations of DNA structure that make it recognizable by proteins. In an effort to directly examine this hypothesis, MD simulations of DNA in a water drop surrounded by vacuum were carried out, which relieves the requirement of charge neutrality. Surprisingly, with zero concentration of counterions, a dodecamer DNA duplex appears metastable, and its structure remains similar to that observed in experiment, including the minor groove narrowing in the dodecamer d(CGCGAATTCGCG)(2) often considered as the most evident cation effect. It is suggested that the same computational approach will allow one to simulate dynamics of long DNA chains more efficiently than with periodical boundary conditions.

Computer Simulation↗
PubMed2003

Fundamental reaction mechanism for cocaine hydrolysis in human butyrylcholinesterase.

Butyrylcholinesterase (BChE)-cocaine binding and the fundamental pathway for BChE-catalyzed hydrolysis of cocaine have been studied by molecular modeling, molecular dynamics (MD) simulations, and ab initio calculations. Modeling and simulations indicate that the structures of the prereactive BChE/substrate complexes for (-)-cocaine and (+)-cocaine are all similar to that of the corresponding prereactive BChE/butyrylcholine (BCh) complex. The overall binding of BChE with (-)-cocaine and (+)-cocaine is also similar to that proposed with butyrylthiocholine and succinyldithiocholine, i.e., (-)- or (+)-cocaine first slides down the substrate-binding gorge to bind to Trp-82 and stands vertically in the gorge between Asp-70 and Trp-82 (nonprereactive complex) and then rotates to a position in the catalytic site within a favorable distance for nucleophilic attack and hydrolysis by Ser-198 (prereactive complex). In the prereactive complex, cocaine lies horizontally at the bottom of the gorge. The fundamental catalytic hydrolysis pathway, consisting of acylation and deacylation stages similar to those for ester hydrolysis by other serine hydrolases, was proposed on the basis of the simulated prereactive complex and confirmed theoretically by ab initio reaction coordinate calculations. Both the acylation and deacylation follow a double-proton-transfer mechanism. The calculated energetic results show that within the chemical reaction process the highest energy barrier and Gibbs free energy barrier are all associated with the first step of deacylation. The calculated ratio of the rate constant (k(cat)) for the catalytic hydrolysis to that (k(0)) for the spontaneous hydrolysis is approximately 9.0 x 10(7). The estimated k(cat)/k(0) value of approximately 9.0 x 10(7) is in excellent agreement with the experimentally derived k(cat)/k(0) value of approximately 7.2 x 10(7) for (+)-cocaine, whereas it is approximately 2000 times larger than the experimentally derived k(cat)/k(0) value of approximately 4.4 x 10(4) for (-)-cocaine. All of the results suggest that the rate-determining step of the BChE-catalyzed hydrolysis of (+)-cocaine is the first step of deacylation, whereas for (-)-cocaine the change from the nonprereactive complex to the prereactive complex is rate-determining and has a Gibbs free energy barrier higher than that for the first step of deacylation by approximately 4 kcal/mol. A further analysis of the structural changes from the nonprereactive complex to the prereactive complex reveals specific amino acid residues hindering the structural changes, providing initial clues for the rational design of BChE mutants with improved catalytic activity for (-)-cocaine.

Butyrylcholinesterase↗