SEARCH · Search PubMed

Results for “Molecular Dynamic 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.

SourceResearch connectionFrom yearThrough year

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

At least 1,405 records · Page 78Linked to original sources

PubMed2006

Structural disorder of the CD3zeta transmembrane domain studied with 2D IR spectroscopy and molecular dynamics simulations.

In a recently reported study [Mukherjee, et al. Proc. Natl. Acad. Sci. U.S.A. 2006, 103, 3528] we used 2D IR spectroscopy and 1-(13)C=(18)O isotope labeling to measure the vibrational dynamics of 11 amide I modes in the CD3zeta transmembrane domain. We found that the homogeneous line widths and population relaxation times were all nearly identical, but that the amount of inhomogeneous broadening correlated with the position of the amide group inside the membrane. In this study, we use molecular dynamics simulations to investigate the structural and dynamical origins of these experimental observations. We use two models to convert the simulations to frequency trajectories from which the mean frequencies, standard deviations, frequency correlation functions, and 2D IR spectra are calculated. Model 1 correlates the hydrogen-bond length to the amide I frequency, whereas model 2 uses an ab initio-based electrostatic model. We find that the structural distributions of the peptidic groups and their environment are reflected in the vibrational dynamics of the amide I modes. Environmental forces from the water and lipid headgroups partially denature the helices, shifting the infrared frequencies and creating larger inhomogeneous distributions for residues near the ends. The least inhomogeneously broadened residues are those located in the middle of the membrane where environmental electrostatic forces are weakest and the helices are most ordered. Comparison of the simulations to experiment confirms that the amide I modes near the C-terminal are larger than at the N-terminal because of the asymmetric structure of the peptide bundle in the membrane. The comparison also reveals that residues at a kink in the alpha-helices have broader line widths than more helical parts of the peptide because the peptide backbone at the kink exhibits a larger amount of structural disorder. Taken together, the simulations and experiments reveal that infrared line shapes are sensitive probes of membrane protein structural and environmental heterogeneity.

CD3 Complex↗

PubMed1996

An efficient mean solvation force model for use in molecular dynamics simulations of proteins in aqueous solution.

An empirical solvation model that allows for the elimination of solvent degrees of freedom in molecular dynamics (MD) simulations of biomolecules is proposed. The potential of mean force due to the first solvation shell is approximated by means of a simple, easily derivable analytic function of the solvent-accessible surface area of the molecule. The solvent contribution to the free energy is evaluated by means of only two atomic solvation parameters. This approach requires about 30% more computational effort than an in vacuo simulation, but a factor of 10 to 50 less than a MD simulation involving solvation by explicit water molecules. The implicit solvation model is assessed by application to proteins of different size. Average structural properties are calculated and compared to values obtained from X-ray structures and from MD simulations using explicit water molecules. The complementarity of the implicit solvation force and the intra-solute force field has been checked. The artefacts induced by the use of a vacuum boundary condition without solvation force in a MD simulation are considerably reduced.

Models, Chemical↗

PubMed1999

Molecular dynamics simulation of alpha-lactalbumin and calcium binding c-type lysozyme.

Alpha-lactalbumins (LAs) and c-type lysozymes (LYZs) are two classes of proteins which have a 35-40% sequence homology and share a common three dimensional fold but perform different functions. Lysozymes bind and cleave the glycosidic bond linkage in sugars, where as, alpha-lactalbumin does not bind sugar but participates in the synthesis of lactose. Alpha-lactalbumin is a metallo-protein and binds calcium, where as, only a few of the LYZs bind calcium. These proteins consist of two domains, an alpha-helical and a beta-strand domain, separated by a cleft. Calcium is bound at a loop situated at the bottom of the cleft and is important for the structural integrity of the protein. Calcium is an ubiquitous intracellular signal in higher eukaryotes and structural changes induced on calcium binding have been observed in a number of proteins. In the present study, molecular dynamics simulations of equine LYZ and human LA, with and without calcium, were carried out. We detail the differences in the dynamics of equine LYZ and human LA, and discuss it in the light of experimental data already available and relate it to the behavior of the functionally important regions of both the proteins. These simulations bring out the role of calcium in the conformation and dynamics of these metallo-proteins. In the calcium bound LA, the region of the protein around the calcium binding site is not only frozen but the atomic fluctuations are found to increase away from the binding site and peak at the exposed sites of the protein. This channeling of fluctuations away from the metal binding site could serve as a general mechanism by which the effect of metal binding at a site is transduced to other parts of the protein and could play a key role in protein-ligand and/or protein-protein interaction.

Animals↗

PubMed2001

Extracellular blockade of K(+) channels by TEA: results from molecular dynamics simulations of the KcsA channel.

TEA is a classical blocker of K(+) channels. From mutagenesis studies, it has been shown that external blockade by TEA is strongly dependent upon the presence of aromatic residue at Shaker position 449 which is located near the extracellular entrance to the pore (Heginbotham, L., and R. MacKinnon. 1992. Neuron. 8:483-491). The data suggest that TEA interacts simultaneously with the aromatic residues of the four monomers. The determination of the 3-D structure of the KcsA channel using X-ray crystallography (Doyle, D.A., J.M. Cabral, R.A. Pfuetzner, A. Kuo, J.M. Gulbis, S.L. Cohen, B.T. Chait, and R. MacKinnon. 1998. Science. 280:69-77) has raised some issues that remain currently unresolved concerning the interpretation of these observations. In particular, the center of the Tyr82 side chains in KcsA (corresponding to position 449 in Shaker) forms a square of 11.8-A side, a distance which is too large to allow simultaneous interactions of a TEA molecule with the four aromatic side chains. In this paper, the external blockade by TEA is explored by molecular dynamics simulations of an atomic model of KcsA in an explicit phospholipid bilayer with aqueous salt solution. It is observed, in qualitative accord with the experimental results, that TEA is stable when bound to the external side of the wild-type KcsA channel (with Tyr82), but is unstable when bound to a mutant channel in which the tyrosine residue has been substituted by a threonine. The free energy profile of TEA relative to the pore is calculated using umbrella sampling simulations to characterize quantitatively the extracellular blockade. It is found, in remarkable agreement with the experiment, that the TEA is more stably bound by 2.3 kcal/mol to the channel with four tyrosine residues. In the case of the wild-type KcsA channel, TEA (which has the shape of a flattened oblate spheroid) acts as an ideal plug blocking the pore. In contrast, it is considerably more off-centered and tilted in the case of the mutant channel. The enhanced stability conferred by the tyrosine residues does not arise from Pi-cation interactions, but appears to be due to differences in the hydration structure of the TEA. Finally, it is shown that the experimentally observed voltage dependence of TEA block, which is traditionally interpreted in terms of the physical position of the TEA along the axis of the pore, must arise indirectly via coupling with the ions in the pore.

Bacterial Proteins↗

PubMed1999

Design of an adenosine phosphorylase by active-site modification of murine purine nucleoside phosphorylase. Enzyme kinetics and molecular dynamics simulation of Asn-243 and Lys-244 substitutions of purine nucleoside phosphorylase.

Our objective was to alter the substrate specificity of purine nucleoside phosphorylase such that it would catalyse the phosphorolysis of 6-aminopurine nucleosides. We modified both Asn-243 and Lys-244 in order to promote the acceptance of the C6-amino group of adenosine. The Asn-243-Asp substitution resulted in an 8-fold increase in K(m) for inosine from 58 to 484 microM and a 1000-fold decrease in k(cat)/K(m). The Asn-243-Asp construct catalysed the phosphorolysis of adenosine with a K(m) of 45 microM and a k(cat)/K(m) 8-fold that with inosine. The Lys-244-Gln construct showed only marginal reduction in k(cat)/K(m), 83% of wild type, but had no activity with adenosine. The Asn-243-Asp;Lys-244-Gln construct had a 14-fold increase in K(m) with inosine and 7-fold decrease in k(cat)/K(m) as compared to wild type. This double substitution catalysed the phosphorolysis of adenosine with a K(m) of 42 microM and a k(cat)/K(m) twice that of the single Asn-243-Asp substitution. Molecular dynamics simulation of the engineered proteins with adenine as substrate revealed favourable hydrogen bond distances between N7 of the purine ring and the Asp-243 carboxylate at 2.93 and 2.88 A, for Asn-243-Asp and the Asn-243-Asp;Lys-244-Gln constructs respectively. Simulation also supported a favourable hydrogen bond distance between the purine C6-amino group and Asp-243 at 2.83 and 2.88 A for each construct respectively. The Asn-243-Thr substitution did not yield activity with adenosine and simulation gave unfavourable hydrogen bond distances between Thr-243 and both the C6-amino group and N7 of the purine ring. The substitutions were not in the region of phosphate binding and the apparent S(0.5) for phosphate with wild type and the Asn-243-Asp enzymes were 1.35+/-0.01 and 1.84+/-0.06 mM, respectively. Both proteins exhibited positive co-operativity with phosphate giving Hill coefficients of 7.9 and 3.8 respectively.

Allosteric Regulation↗

PubMed2005

A cylinder-shaped double ribbon structure formed by an amyloid hairpin peptide derived from the beta-sheet of murine PrP: an X-ray and molecular dynamics simulation study.

A structural model of the murine PrP small beta-sheet was obtained by synthesizing the RGYMLGSADPNGNQVYYRG peptide comprising the two beta-strands 127-133 and 159-164 linked by a four-residue sequence of high turn propensity. The DPNG turn sequence is a "short circuit" replacing the original protein sequence between the two strands. This 19-residue peptide spontaneously forms very long single fibrils as observed by electron microscopy. The X-ray diffraction patterns of a partially oriented sample reveals an average arrangement of the hairpin peptides into a structure which can be geometrically approximated by an empty-core cylinder. The hairpins are oriented perpendicular to the cylinder axis and a 130 A helix period is observed. Based on X-ray diffraction constraints and on more indirect general protein structure considerations, a precise and consistent fibril model was built. The structure consists of two beta-sheet ribbons wound around a cylinder and assembled into a single fibril with a hairpin orientation perpendicular to the fibril axis. Subsequent implicit and explicit solvent molecular dynamics simulations provided the final structure at atomic resolution and further insights into the stabilizing interactions. Particularly important are the zipper-like network of polar interactions between the edges of the two ribbons, including the partially buried water molecules. The hydrophobic core is not optimally compact explaining the low density of this region seen by X-ray diffraction. The present findings provide also a simple model for further investigating the sequence-stability relationship using a mutational approach with a quasi-independent consideration of the polar and apolar interactions.

Amino Acid Sequence↗

PubMed2000

Molecular dynamics simulation of human prion protein including both N-linked oligosaccharides and the GPI anchor.

Although glycosylation appears to protect prion protein (PrP(C)) from the conformational transition to the disease-associated scrapie form (PrP(Sc)), available NMR structures are for non-glycosylated PrP(C), only. To investigate the influence of both the two N-linked glycans, Asn181 and Asn197, and of the GPI anchor attached to Ser230, on the structural, dynamical and electrostatic behavior of PrP, we have undertaken molecular dynamics simulations on the C-terminal region of human prion protein HU:PrP(90-230), with and without the three glycans. The simulations used the AMBER94 force field in a periodic box model with explicit water molecules, considering all long-range electrostatic interactions. The results suggest the structured part of the protein, HU:PrP(127-227) is stabilized overall from addition of the glycans, specifically by extensions of Helix-B and Helix-C and reduced flexibility of the linking turn containing Asn197, although some regions such as residues in the turn (165-170) between Strand-B and Helix-B have increased flexibility. The stabilization appears indirect, by reducing the mobility of the surrounding water molecules, and not from specific interactions such as H bonds or ion pairs. The results are consistent with glycosylation at Asn197 having a stabilizing role, while that at Asn181, in a region with already stable secondary structure, having a more functional role, in agreement with literature suggestions. Due to three negatively charged SiaLe(x) groups per N-glycan, the surface electrostatic properties change to a negative electrostatic field covering most of the C-terminal part, including the surface of Helix-B and Helix-C, while the positively charged N-terminal part PrP(90-126) of undefined structure creates a positive potential. The unusual hydrophilic Helix-A (144-152) is not covered by either of these dominant electrostatic fields, and modeling shows it could readily dimerize in anti parallel fashion. In combination with separate simulations of the GPI anchor in a membrane model, the results show the GPI anchor is highly flexible and would maintain the protein at a distance between 9 and 13 A from the membrane surface, with little influence on its structure or orientational freedom.

Amino Acid Sequence↗

PubMed2006

Ab initio molecular dynamics simulation of NO reactivity on the CaO(001) surface.

Ab initio molecular dynamics (MD) is used to investigate NO reaction processes on the (001) surface of CaO. A novel path is proposed for the first steps of nitrogen oxides reactivity catalyzed by the CaO surface. The mechanism consists of the formation of anionic dimers, adsorbing on the surface cations, at the expense of oxidized NO species adsorbed on surface anions. The complete charge-transfer process takes place in two steps, producing first monovalent anionic dimers (NO)2- and, later on, divalent anionic dimers (NO)2(2-). These redox processes cause spin quenching and are observed in the short time scale of the ab initio MD simulation at 300 K. The results presented provide a rationalization of a recent electron spin resonance (ESR) investigation indicating that the spectroscopy is silent to most of the nitrogen oxide species adsorbed on CaO powders, despite deposition of paramagnetic NO molecules at room temperature.

Journal Article↗

PubMed1999

Solution structures of a duplex containing an adenine opposite a gap (absence of one nucleotide). An NMR study and molecular dynamic simulations with explicit water molecules.

We investigated the behaviour of a 15mer DNA duplex, [5'd(CAGAGTCACTGGCTC)3']. [5'd(GAGCCAG)3' + 5'd(GACTCTG)3'] which contained an adenine opposite the gap. Analysis of the NMR data showed the existence of one major species, which was in equilibrium with two minor species. Their relative concentrations varied as a function of pH with a pKa of approximately 4.5. For the major species, the duplex was globally in B conformation with the central adenine stacked in the helix. The two G.C base pairs adjacent to the central adenine were well formed and a gap was present in front of this adenine. For the minor species, major structural perturbations occurred in the centre of the duplex. At neutral pH, the central adenine was involved in a G.A mismatch with G23 adjacent to the gap. Cytosine C7 was then extrahelical and no gap was observed. Under these conditions, the major neutral species corresponded to 70% of the total and the minor species to 30%. At acidic pH, the central adenine of the minor species was protonated and was involved in a G(syn).A+(anti) mismatch. The difference is that C9 is now extrahelical and G22 is implicated in the mispair. Three-dimensional models were built to initiate molecular dynamic simulations, which were in good agreement with the NMR data. Their structural stability in terms of hydrogen bonding and their flexibility are discussed and the biological significance for the interaction with DNA polymerase is evoked.

Adenine↗

PubMed1997

Kinetic folding and cis/trans prolyl isomerization of staphylococcal nuclease. A study by stopped-flow absorption, stopped-flow circular dichroism, and molecular dynamics simulations.

We studied the urea-induced unfolding transition of staphylococcal nuclease (SNase) and its five proline mutants (P47A, P47T, P117G, P47T/P117G, and P47A/P117G) [corrected] by peptide and aromatic circular dichroism and aromatic absorption spectroscopy at equilibrium and the refolding-unfolding kinetics of the proteins by stopped-flow circular dichroism and stopped-flow absorption techniques. Recent studies have revealed that the cis/trans isomerizations about the Pro47 and Pro117 peptide bonds of SNase occur not only in the unfolded state but also in the native state. The mutational effects on the stability and the refolding-unfolding kinetics of SNase were, however, remarkably different between the two sites. The substitution of Ala or Thr for Pro47 neither changed the stability nor affected the refolding-unfolding kinetics of SNase, whereas the substitution of Gly for Pro117 increased the protein stability by 1.2 kcal/mol (pH 7.0 and 20 degrees C) and affected the kinetics. These results have been attributed to the high flexibility of the loop around Pro47, which has been revealed by molecular dynamics simulations of native SNase. Under every condition studied, cooperative refolding-unfolding kinetics of SNase were observed. Refolding of wild-type SNase was represented by two urea concentration-dependent fast phases and a urea concentration-independent slow phase. The double mutant (P47T/P117G) [corrected] of SNase still showed multiphasic refolding kinetics that involved two urea concentration-independent slow phases, suggesting that isomerization of proline residues other than Pro47 and Pro117 may occur in the unfolded state of the mutant. Two phases were observed in the unfolding of the wild-type and mutant proteins that contained Pro117, a fast phase corresponding to the unfolding of the trans isomer and a slow phase corresponding to that of the cis isomer. On the basis of these results, the folding scheme of SNase is discussed.

Alanine↗

PubMed1992

Quenching of fluorescein-conjugated lipids by antibodies. Quantitative recognition and binding of lipid-bound haptens in biomembrane models, formation of two-dimensional protein domains and molecular dynamics simulations.

Three model biomembrane systems, monolayers, micelles, and vesicles, have been used to study the influence of chemical and physical variables of hapten presentation at membrane interfaces on antibody binding. Hapten recognition and binding were monitored for the anti-fluorescein monoclonal antibody 4-4-20 generated against the hapten, fluorescein, in these membrane models as a function of fluorescein-conjugated lipid architecture. Specific recognition and binding in this system are conveniently monitored by quenching of fluorescein emission upon penetration of fluorescein into the antibody's active site. Lipid structure was shown to play a large role in affecting antibody quenching. Interestingly, the observed degrees of quenching were nearly independent of the lipid membrane model studied, but directly correlated with the chemical structure of the lipids. In all cases, the antibody recognized and quenched most efficiently a lipid based on dioctadecylamine where fluorescein is attached to the headgroup via a long, flexible hydrophilic spacer. Dipalmitoyl phosphatidylethanolamine containing a fluorescein headgroup demonstrated only partial binding/quenching. Egg phosphatidylethanolamine with a fluorescein headgroup showed no susceptibility to antibody recognition, binding, or quenching. Formation of two-dimensional protein domains upon antibody binding to the fluorescein-lipids in monolayers is also presented. Chemical and physical requirements for these antibody-hapten complexes at membrane surfaces have been discussed in terms of molecular dynamics simulations based on recent crystallographic models for this antibody-hapten complex (Herron et al., 1989. Proteins Struct. Funct. Genet. 5:271-280).

Animals↗

PubMed1991

A molecular dynamics simulation of polyalanine: an analysis of equilibrium motions and helix-coil transitions.

An understanding of helix dynamics can aid in interpreting the motions of proteins. The conformational transitions that occur also appear to play a role in protein folding. Structural studies of isolated peptides in solution are just becoming available. However, detailed analysis of the helix-coil transition is still not available and will be difficult to obtain experimentally. For these reasons, we performed a long molecular dynamics simulation of polyalanine at high temperature. Using this approach, we obtain a description of the overall structure and inherent flexibility of the chain as well as a structural picture of the conformational changes that occur. In this way, we can address both equilibrium properties of the peptide and the dynamics and mechanisms of the structural transitions. Our results correlate fairly well with the available experimental data and previous simulations aimed at addressing alpha-helix dynamics. The peptide spends the bulk of its time fluctuating between different conformations with intermediate helix contents. Transitions between highly ordered and highly disordered structures were rare, but they occurred rapidly. Our distribution of conformations favored collapsed states. Hence, our transitions to structures with high helical content were from fluctuating compact structures. The conversion between helix and coil occurred sequentially on a residue-by-residue basis. However, there was local cooperativity; the transition of a residue to the coil state was facilitated after a neighboring group became nonhelical. The relevance of our results to protein folding is also discussed.

Models, Chemical↗

PubMed1996

Using a hydrophobic contact potential to evaluate native and near-native folds generated by molecular dynamics simulations.

There are several knowledge-based energy functions that can distinguish the native fold from a pool of grossly misfolded decoys for a given sequence of amino acids. These decoys, which are typically generated by mounting, or "threading", the sequence onto the backbones of unrelated protein structures, tend to be non-compact and quite different from the native structure: the root-mean-squared (RMS) deviations from the native are commonly in the range of 15 to 20 angstroms. Effective energy functions should also demonstrate a similar recognition capability when presented with compact decoys that depart only slightly in conformation from the correct structure (i.e. those with RMS deviations of approximately 5 angstroms or less). Recently, we developed a simple yet powerful method for native fold recognition based on the tendency for native folds to form hydrophobic cores. Our energy measure, which we call the hydrophobic fitness score, is challenged to recognize the native fold from 2000 near-native structures generated for each of five small monomeric proteins. First, 1000 conformations for each protein were generated by molecular dynamics simulation at room temperature. The average RMS deviation of this set of 5000 was 1.5 angstroms. A total of 323 decoys had energies lower than native; however, none of these had RMS deviations greater than 2 angstroms. Another 1000 structures were generated for each at high temperature, in which a greater range of conformational space was explored (4.3 angstroms RMS deviation). Out of this set, only seven decoys were misrecognized. The hydrophobic fitness energy of a conformation is strongly dependent upon the RMS deviation. On average our potential yields energy values which are lowest for the population of structures generated at room temperature, intermediate for those produced at high temperature and highest for those constructed by threading methods. In general, the lowest energy decoy conformations have backbones very close to native structure. The possible utility of our method for screening backbone candidates for the purpose of modelling by side-chain packing optimization is discussed.

Protein Folding↗

PubMed2005

Molecular dynamics simulations of the endocannabinoid N-arachidonoylethanolamine (anandamide) in a phospholipid bilayer: probing structure and dynamics.

The phospholipid bilayer plays a central role in the lifecycle of the endogenous cannabinoid N-arachidonoylethanolamine (anandamide, 1). Compound 1 has been shown to be synthesized from lipids, to interact with the membrane-embedded cannabinoid CB1 receptor, to be transported to intracellular compartments, possibly via caveolae-related endocytosis, and finally, to be degraded by fatty acid amide hydrolase (FAAH), an integral membrane protein which has an active site that is accessed by 1 possibly via the bilayer. Because the anandamide system is intimately associated with the lipid milieu, information concerning the location of 1 in the phospholipid bilayer and the conformations it can adopt is important to our understanding of the mechanism of cannabinoid action at the molecular level. We report here an exploration of the properties of 1 in a 1,2-dioleoyl-sn-glycero-3-phosphocholine (DOPC) phospholipid bilayer via multi-nanosecond molecular dynamics simulations. Our results suggest that the polar headgroup of 1 resides at the lipid-water interface, specifically in the polar phospholipid headgroup region, whereas the nonpolar acyl chain of 1 extends into the hydrocarbon core of the membrane. Our analysis also indicates that (i) an elongated conformation of 1 is preferred in the DOPC bilayer environment; however, many other conformations of 1 are observed; (ii) hydrogen-bonding between the lipid (DOPC) and the headgroup of 1, although extensive, is quite short-lived; and (iii) the C-H bond order parameters for the acyl chain of 1 are low compared to order parameters typically seen for saturated acyl chains of fatty acids, and these order parameters decrease toward the bilayer center. The bilayer location for 1 revealed by these studies may be important for the interaction of 1 with membrane-embedded proteins such as the cannabinoid CB1 receptor and membrane-associated proteins such as FAAH.

Arachidonic Acids↗

PubMed2001

Molecular dynamics simulation of adrenocorticotropin (1-10) peptide in a solvated dodecylphosphocholine micelle.

Adrenocorticotropin (ACTH) (1-10), an adrenocorticotropin hormone fragment, has been studied by molecular dynamics (MD) simulation in an NPT ensemble in an explicit dodecylphosphocholine (DPC) micelle. Two starting configurations of the peptide/micelle system, corresponding to the insertion and surface-binding modes, were used. A common equilibrated configuration, in which the peptide lies parallel to the micellar surface, was reached from both simulations. In the initial part of the simulations, distance restraints derived from NMR nuclear Overhauser enhancements were incorporated before the peptide reached an equilibrium configuration with respect to the micelle. Analyses of the trajectories from the subsequent free (unrestrained) MD simulation showed that ACTH (1-10) does not conform strictly to a helical structure. The loss of the helical structure is due to decreased intramolecular hydrogen bonding accompanied by an increase of hydrogen bonding between the amide protons of the peptide and the micellar head groups. However, the extent of the latter interaction is less pronounced than in the negatively charged SDS micelle. The final structure enhances the amphipathic nature of the peptide, facilitating better interactions at the water-hydrophobic interface. The primary hydrophobic interactions with the micelle came from the side chains of Met4, Phe7, and Trp9. All peptide bonds were either hydrated or were involved in intramolecular hydrogen bonding. The interactions with the DPC micelle, the conformation of the bound peptide, and the dynamics of the peptide, as revealed by the time correlation functions of the N-H bonds, were compared with those of the ACTH (1-10)/SDS system studied previously by MD simulations.

Adrenocorticotropic Hormone↗

PubMed1999

Molecular dynamics simulations of P450 BM3--examination of substrate-induced conformational change.

Cytochrome P450 BM3, of bacterial origin, is one of only five isozymes of the ubiquitous family of over 400 metabolizing heme proteins with a known crystal structure and only one of two with both substrate-free and substrate-bound forms determined. P450 BM3 is of particular interest since it has a similar function and similar substrates as mammalian P450s particularly of the 4A subfamily. Thus, the extent to which the substrate-free form of P450 BM3 undergoes a conformational change upon binding of a typical fatty acid substrate, palmitoleic acid, has been the subject of recent active experimental effort. Surprisingly, direct examination of the substrate-free (pdb2hpd.ent and pdb2bmh.ent) and substrate-bound (pdb1fag.ent) forms do not provide a clear answer to this question. The main reason for this ambiguity is that the two substrate-free monomers reported in the crystal structures themselves have significantly different conformations from each other, one with a more open substrate-access channel than the other. Since there is no way to tell to which substrate-free form the substrate binds, the effect of substrate binding cannot be deduced directly from comparisons of the experimental substrate-bound and substrate-free forms. The computational studies reported here have been designed to more robustly establish the effect of substrate binding on this isozyme. Specifically, molecular dynamics simulations were performed for each of the two substrate-free forms found in the asymmetric unit of the X-ray structure and for the two corresponding substrate-bound forms, constructed by docking palmitloeic acid into each of them. Comparisons of the results showed that palmitoleic acid binding had little effect on the conformation of the more closed substrate-free form of P450 BM3. By contrast, in the more open substrate-free form, this same substrate induced a closing of the entrance to the substrate-binding channel. The MD averaged structure of these two complexes obtained from docking of pamitoleic acid into the two asymmetric units of the substrate-free form were also compared to that obtained starting with the X-ray structure of the substrate-bound form. These results taken together led to the conclusion that, if indeed the substrate induces conformational changes in P450 BM3, the mouth of the substrate-access channel first closes down in response to the presence of the substrate, followed by rotation of the F-G domain to further optimize the P450 BM3-substrate interaction that would occur at a later stage.

Bacterial Proteins↗

PubMed2002

HIV-1(Lai) genomic RNA: combined used of NMR and molecular dynamics simulation for studying the structure and internal dynamics of a mutated SL1 hairpin.

The genome of all retroviruses consists of two identical copies of an RNA sequence associated in a non-covalent dimer. A region upstream from the splice donor (SL1) comprising a self-complementary sequence is responsible for the initiation of the dimerization. This region is able to dimerize in two conformations: a loop-loop complex or an extended duplex. Here, we solve by 2D NMR techniques the solution structure of a 23-nucleotide sequence corresponding to HIV-1 SL1(Lai) in which the mutation G12-->A12 is included to prevent dimerization. It is shown that this monomer adopts a stem-loop conformation with a seven base pairs stem and a nine nucleotide loop containing the G10 C11 A12 C13 G14 C15 sequence. The stem is well structured in an A-form duplex, while the loop is more flexible even though elements of structure are evident. We show that the structure adopted by the stem can be appreciably different from its relaxed structure when the adenines A8, A9 and A16 in the loop are mechanically constrained. This point could be important for the efficiency of the dimerization. This experimental study is complemented with a 10 ns molecular dynamics simulation in the presence of counterions and explicit water molecules. This simulation brings about information on the flexibility of the loop, such as a hinge motion between the stem and the loop and a labile lattice of hydrogen bonds in the loop. The bases of the nucleotides G10 to C15 were found outside of the loop during a part of the trajectory, which is certainly necessary to initiate the dimerization process of the genuine SL1(Lai) sequence.

Base Sequence↗

PubMed2006

Formation and structural anomaly of the metastable phases in an immiscible Ag-Mo system studied by ion beam mixing and molecular dynamics simulation.

For the equilibrium immiscible Ag-Mo system characterized by a large positive heat of formation, the nanosized Ag-Mo multilayered samples are designed and prepared to include sufficient interfacial free energy to elevate their initial energetic states to be higher than that of either the amorphous phase or solid solution and then subject to 200 keV xenon ion irradiation. The results show that a uniform amorphous alloy can be obtained within a composition range, at least, from 25 to 88 atom % of Mo. Interestingly, in the intermediate stage of ion irradiation, a bcc phase, an amorphous phase, and an order (bcc)-disorder coexisting state appear simultaneously in the Ag12Mo88 multilayered sample and extend over the entire bright field image with unanimously homogeneous composition. In thermodynamic modeling, a Gibbs free energy diagram of the Ag-Mo system is constructed, based on Miedema's model, and suggests that within a narrow composition regime of 85-90 atom % of Mo, the energy difference between the bcc and the amorphous phases is extremely small, which is probably the very reason for the order-disorder coexisting state to appear. In atomistic modeling, an ab initio derived Ag-Mo potential is applied to perform molecular dynamics simulations. The simulations not only determine an intrinsic glass-forming ability/range (GFA/GFR) of the Ag-Mo system to be from 10 to 88 atom % of Mo but also reveal the possibility of the formation/appearance of a crystalline and amorphous mixture in a narrow composition regime of 88-92 atom % of Mo. Apparently, the theoretical results are in excellent agreement and/or compatible with the experimental observations in ion beam mixing.

Journal Article↗