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,027 records · Page 57Linked to original sources
PubMed2004

Analysis and simulation of the structure of nanoparticles that undergo a surface-driven structural transformation.

A room temperature solid-state structural transformation was observed in 3 nm ZnS nanoparticles in methanol following the addition of water (Zhang et al., Nature 424, 1025, 2003). Experimental wide angle x-ray scattering (WAXS), x-ray absorption near edge structure (XANES) and extended x-ray absorption fine structure (EXAFS) spectroscopy measurements show a large increase in crystallinity associated with water addition, in agreement with molecular dynamics (MD) predictions. Here we perform first-shell EXAFS and pair distribution function analysis and whole-nanoparticle calculations of WAXS, EXAFS and XANES to compare structural data with the MD predictions. The predicted WAXS patterns give excellent agreement with data, while the predicted EXAFS and XANES spectra give poor agreement. Relative to WAXS, XANES and EXAFS spectra contain additional structural information related to the distribution of disorder. The discrepancy between the x-ray diffraction and x-ray absorption results indicates that structural disorder is partitioned between interior and surface regions more strongly than predicted in the MD simulations.

Journal Article↗
PubMed2026

Structural Basis for Stepwise Substrate Transport and Disease Phenotypic Heterogeneity of the Mitochondrial ADP/ATP Carrier.

The mitochondrial ADP/ATP carrier (AAC) is essential for cellular energy metabolism and responsible for exchanging ADP for ATP across the inner mitochondrial membrane. However, the precise molecular determinants of substrate binding and the mechanisms underlying the phenotypic heterogeneity of AAC-related diseases remain poorly understood. Here, we combined AlphaFold3 predictions, molecular dynamics simulations, and experimental validation to identify and characterize a previously unrecognized ADP-binding site in AAC (site S2), which is distinct from the canonical bottom site (site S1). AlphaFold3 predictions on AAC variants with disrupted site S1 consistently placed ADP at site S2, interacting with residues R188, K92, and K96, a finding that was independently corroborated by our prior MD simulations. Systematic mutagenesis and functional analysis revealed distinct roles for site S2 residues: R188 serves as the primary phosphate-specific anchor; K92 and K96 facilitate initial recruitment and stabilization; and the aromatic ladder (Y187/Y191/F192/Y195) assists the conformational transition of ADP from anti to syn, a critical step enabling downward translocation. Functional characterization demonstrated that mutations in site S2 significantly impaired ADP transport and oxidative phosphorylation. Notably, the spatial distribution of AAC disease mutations correlates with this bipartite architecture: mild PEOA2-associated mutations cluster near site S2 and perturb the local conformation without abolishing binding, whereas lethal mutations cluster near site S1 and disrupt both ADP binding and structural integrity of the m-gate. Our findings provide new mechanistic insights into stepwise substrate transport and potential therapeutic targets for AAC-related diseases.

Mitochondrial ADP, ATP Translocases↗
PubMed1997

Refolding simulations of an isolated fragment of barnase into a native-like beta hairpin: evidence for compactness and hydrogen bonding as concurrent stabilizing factors.

Experimental evidence and theoretical models both suggest that protein folding is initiated within specific fragments intermittently adopting conformations close to that found in the protein native structure. These folding initiation sites encompassing short portions of the protein are ideally suited for study in isolation by computational methods aimed at peering into the very early events of folding. We have used Molecular Dynamics (MD) technique to investigate the behavior of an isolated protein fragment formed by residues 85 to 102 of barnase that folds into a beta hairpin in the protein native structure. Three independent MD simulations of 1.3 to 1.8 ns starting from unfolded conformations of the peptide portrayed with an all-atom model in water were carried out at gradually decreasing temperature. A detailed analysis of the conformational preferences adopted by this peptide in the course of the simulations is presented. Two of the unfolded peptides conformations fold into a hairpin characterized by native and a larger bulk of nonnative interactions. Both refolding simulations substantiate the close relationship between interstrand compactness and hydrogen bonding network involving backbone atoms. Persistent compactness witnessed by side-chain interactions always occurs concomitantly with the formation of backbone hydrogen bonds. No highly populated conformations generated in a third simulation starting from the remotest unfolded conformer relative to the native structure are observed. However, nonnative long-range and medium-range contacts with the aromatic moiety of Trp94 are spotted, which are in fair agreement with a former nuclear magnetic resonance study of a denaturing solution of an isolated barnase fragment encompassing the beta hairpin. All this lends reason to believe that the 85-102 barnase fragment is a strong initiation site for folding.

Bacterial Proteins↗
PubMed2006

A molecular dynamics approach to the structural characterization of amyloid aggregation.

A novel computational approach to the structural analysis of ordered beta-aggregation is presented and validated on three known amyloidogenic polypeptides. The strategy is based on the decomposition of the sequence into overlapping stretches and equilibrium implicit solvent molecular dynamics (MD) simulations of an oligomeric system for each stretch. The structural stability of the in-register parallel aggregates sampled in the implicit solvent runs is further evaluated using explicit water simulations for a subset of the stretches. The beta-aggregation propensity along the sequence of the Alzheimer's amyloid-beta peptide (Abeta(42)) is found to be highly heterogeneous with a maximum in the segment V(12)HHQKLVFFAE(22) and minima at S(8)G(9), G(25)S(26), G(29)A(30), and G(38)V(39), which are turn-like segments. The simulation results suggest that these sites may play a crucial role in determining the aggregation tendency and the fibrillar structure of Abeta(42). Similar findings are obtained for the human amylin, a 37-residue peptide that displays a maximal beta-aggregation propensity at Q(10)RLANFLVHSSNN(22) and two turn-like sites at G(24)A(25) and G(33)S(34). In the third application, the MD approach is used to identify beta-aggregation "hot-spots" within the N-terminal domain of the yeast prion Ure2p (Ure2p(1-94)) and to design a double-point mutant (Ure2p-N4748S(1-94)) with lower beta-aggregation propensity. The change in the aggregation propensity of Ure2p-N4748S(1-94) is verified in vitro using the thioflavin T binding assay.

Amino Acid Sequence↗
PubMed2006

Structural dynamics of the mitochondrial ADP/ATP carrier revealed by molecular dynamics simulation studies.

The mitochondrial adenosine diphosphate/adenosine triphosphate (ADP/ATP) carrier has been recently crystallized in complex with its specific inhibitor carboxyatractyloside (CATR). In the crystal structure, the six-transmembrane helix bundle that defines the nucleotide translocation pathway is closed on the matrix side due to sharp kinks in the odd-numbered helices. The closed conformation is further sealed by the loops protruding into the matrix that interact through an intricate network of charge-pairs. To gain insight into its structural dynamics we performed molecular dynamics (MD) simulation studies of the ADP/ATP carrier with and without its cocrystallized inhibitor. The two trajectories sampled a conformational space around two different configurations characterized by distinct salt-bridge networks with a significant shift from inter- to intrarepeat bonding on the matrix side in the absence of CATR. Analysis of the geometrical parameters defining the transmembrane helices showed that even-numbered helices can undergo a face rotation, whereas odd-numbered helices can undergo a change in the wobble angle with a conserved proline acting as molecular hinge. Our results provide new information on the dynamical properties of the ADP/ATP carrier and for the first time yield a detailed picture of a stable carrier conformation in absence of the inhibitor.

Animals↗
PubMed2005

Water chain formation and possible proton pumping routes in Rhodobacter sphaeroides cytochrome c oxidase: a molecular dynamics comparison of the wild type and R481K mutant.

Cytochrome c oxidase (CcO) converts the energy from redox and oxygen chemistry to support proton translocation and create a transmembrane DeltamuH(+) used for ATP production. Molecular dynamics (MD) simulations were carried out to probe for the formation water chains capable of participating in proton translocation. Attention was focused on the region between and above the a and a(3) hemes where well-defined water chains have not been identified in crystallographic studies. An arginine (R481) (Rhodobacter sphaeroides numbering), positioned between the D-propionates of the hemes, had been mutated in vivo to lysine and showed to have altered activity consistent with an altered proton conductance [Qian, J., Mills, D. A., Geren, L., Wang, K. F., Hoganson, C. W., Schmidt, B., Hiser, C., Babcock, G. T., Durham, B., Millett, F., and Ferguson-Miller, S. (2004) Role of the conserved arginine pair in proton and electron transfer in cytochrome c oxidase, Biochemistry 43, 5748-5756; also see the accompanying paper by Mills et al.]. This mutant was created in silico, and the MD results for the mutant and wild type were compared to explore the effects on the formation of hydrogen-bonded water chains by this mutation. The simulations reveal the presence of hydrogen-bonded water chains that lead from E286 through the region above the hemes to the Mg(2+), and from E286 to the heme a(3) D-propionate and the binuclear center. The R481K mutant does not form as many, or as extensive, water chains as wild-type CcO, due to a new conformation of residues in a large loop between helices III and IV in subunit I, indicating a reduction in the level of water chain formation in the mutant. This loop appears to play a role in controlling the formation of hydrogen-bonded water chains above the hemes. The results suggest a possible gating mechanism for proton movement that includes key residues W172 and Y175 on the loop and F282 on helix VI.

Arginine↗
PubMed1989

Molecular dynamics simulation as a complement to diffraction in the study of disorder in crystals.

It is not uncommon that a crystal structure contains a minority component which does not share the overall space-group of the remainder of the structure, but has a lower (or no) effective space-group symmetry. A conventional diffraction treatment will inevitably be made within the space-group of the majority of the system, and will therefore result in an average structure also with this symmetry. We here indicate the value of molecular dynamics (MD) simulation within a supercell of the structure as a complementary tool in providing information concerning the local structure for the minority component. Such information, though illustrative, is often the most interesting and the physically most significant aspect of the structure. It can also be crucial information for the testing of new theory, and in the design of new materials. The approach is illustrated for the case of the Na+ ion distribution in Na+ beta"-alumina, and for the experimentally more inaccessible mixed-ion system Na+/Ba2+ beta"-alumina.

Aluminum↗
PubMed1992

Unusual thionation of a cyclic hexapeptide. Conformational changes and dynamics.

One carbonyl oxygen of the cyclic hexapeptide cyclo(-Gly1-Pro2-Phe3-Val4-Phe5-Phe6-) (A) can be selectively exchanged with sulphur using Yokoyama's reagent. Surprisingly it was not the C=] of Gly1 but that of Phe5 which was substituted and cyclo(-Gly1-Pro2-Phe3-Val4-Phe5 psi [CS-NH]Phe6-) (B) was obtained. Thionation results in a conformational change of the peptide backbone although the C=O of Phe5 and the corresponding C=S are not involved in internal hydrogen bonds. Two isomers in slow exchange, containing a cis Gly1-Pro2 bond in a beta VIa-turn (minor) and a trans Gly-Pro bond in a beta II'-turn (major), were analyzed by restrained molecular dynamics in vacuo and in DMSO as well as using time dependent distance constraints. It is impossible to fit all experimental data to a static structure of each isomer. Interpreting the conflicting NOEs, local segment flexibility is found. MD simulations lead to a dynamic model for each structure with evidence of an equilibrium between a beta I- and beta II-turn about the Val4-Phe5 amide bond in both the cis and trans isomers. Additionally proton relaxation rates in the rotating frame (R1 rho) were measured to verify the assumption of this fast beta I/beta II equilibrium within each isomer. Significant contributions to R1 rho-rates from intramolecular motions were found for both isomers. Therefore it is possible to distinguish between at least four conformers interconverting on different time scales based on NMR data and MD refinement. This work shows that thionation is a useful modification of peptides for conformation-activity investigations.

Amino Acid Sequence↗
PubMedDate not supplied

Structure based prediction of binding affinity of human immunodeficiency virus-1 protease inhibitors.

A series of computations were performed to derive a strategy for the prediction of binding affinities of non-peptidic human immunodeficiency virus-1 (HIV-1) protease inhibitors. This paper describes the development of a 3D quantitative structure-activity relationship (3D-QSAR) methodology by using receptor information of HIV-1 protease. The docking and molecular dynamics simulations were performed on a model ligand/enzyme complex to optimize the variables involved in the generation of ligand/enzyme models. The protonation scheme of the active site aspartic acid residues of HIV-1 protease was derived from a computational study. The active site aspartate is monoprotonated with a proton placed on the OD1 atom of the ASP B25. This protocol of docking and molecular dynamics (MD) simulation was then used to derive the ligand-enzyme complexes of the molecules used in the present study. The molecular mechanics interaction descriptors were calculated from these ligand/enzyme models. A partial least squares (PLS) method was used to derive a linear correlation between the interaction descriptors and the biological activity. A good correlation was observed when the change in the energy of the ligand upon complex formation and the electrostatic contributions to the solvation energy of the ligand were included in the QSAR analysis. A highest cross-validated q2 value of 0.649 was observed. This model had a conventional r2 of 0.725, and when this model was used to predict the activity of the external test set, it produced a r2pred of 0.761. The total interaction energy was partitioned into interactions in different subsites and interactions with each of the amino acid residues of the enzyme. The PLS analysis using these descriptors helped to identify the important interactions which can be exploited for the design of HIV-1 protease inhibitors.

HIV Protease↗
PubMedDate not supplied

Constructing protein models for ligand-receptor binding thermodynamic simulations: an application to a set of peptidometic renin inhibitors.

Structure-based design is the application of ligand-receptor modeling to predict the activity of a series of molecules that bind to a common receptor for which the molecular geometry is available. Successful structure-based design requires an accurate receptor model which can be economically employed in the design calculations. One goal of the work reported here has been to reduce the size of a model structure of a macromolecular receptor to allow multiple ligand-receptor molecular dynamic (MD) simulations to be computationally economical yet still provide meaningful binding thermodynamic data. A scaled-down 10 A receptor model of the enzyme renin, when subjected to an alternate atomic mass constraint, maintains the structural integrity of the composite parent crystal structure. A second goal of the work has been to develop schemes to explore and characterize the protonation states of receptors and ligand-receptor systems. Application of the charge state characterization schemes to the hydroxyethylene and statine transition state inhibitors of renin in the training set suggests a monoprotonation state of the two active-site aspartate residues, where the lone proton resides on the outer carboxylate oxygen of Asp226 is most likely. For the reduced amide transition state inhibitors an active site consisting of both aspartates in the totally ionized state, and the ligand carrying a net +1.0 charge, is most stable and consistent with experimental data.

Amino Acid Sequence↗
PubMed2002

Hoogsteen-based parallel-stranded duplexes of DNA. Effect of 8-amino-purine derivatives.

The structure of parallel-stranded duplexes of DNA-containing a mixture of guanines (G) and adenines (A) is studied by means of molecular dynamics (MD) simulation, as well as NMR and circular dichroism (CD) spectroscopy. Results demonstrate that the structure is based on the Hoogsteen motif rather than on the reverse Watson-Crick one. Molecular dynamics coupled to thermodynamic integration (MD/TI) calculations and melting experiments allowed us to determine the effect of 8-amino derivatives of A and G and of 8-amino-2'-deoxyinosine on the stability of parallel-stranded duplexes. The large stabilization of the parallel-stranded helix upon 8-amino substitution agrees with a Hoogsteen pairing, confirming MD, NMR, and CD data, and suggests new methods to obtain DNA triplexes for antigene and antisense purposes.

Circular Dichroism↗
PubMed2005

Construction of a 3D model of nattokinase, a novel fibrinolytic enzyme from Bacillus natto. A novel nucleophilic catalytic mechanism for nattokinase.

A three-dimensional structural model of nattokinase (NK) from Bacillus natto was constructed by homology modeling. High-resolution X-ray structures of Subtilisin BPN' (SB), Subtilisin Carlsberg (SC), Subtilisin E (SE) and Subtilisin Savinase (SS), four proteins with sequential, structural and functional homology were used as templates. Initial models of NK were built by MODELLER and analyzed by the PROCHECK programs. The best quality model was chosen for further refinement by constrained molecular dynamics simulations. The overall quality of the refined model was evaluated. The refined model NKC1 was analyzed by different protein analysis programs including PROCHECK for the evaluation of Ramachandran plot quality, PROSA for testing interaction energies and WHATIF for the calculation of packing quality. This structure was found to be satisfactory and also stable at room temperature as demonstrated by a 300ps long unconstrained molecular dynamics (MD) simulation. Further docking analysis promoted the coming of a new nucleophilic catalytic mechanism for NK, which is induced by attacking of hydroxyl rich in catalytic environment and locating of S221.

Amino Acid Sequence↗
PubMed2001

Fourier acceleration of Langevin molecular dynamics.

Fourier acceleration has been successfully applied to the simulation of lattice field theories for more than a decade. In this paper, we extend the method to the dynamics of discrete particles moving in a continuum. Although our method is based on a mapping of the particles' dynamics to a regular grid so that discrete Fourier transforms may be taken, it should be emphasized that the introduction of the grid is a purely algorithmic device and that no smoothing, coarse-graining, or mean-field approximations are made. The method thus can be applied to the equations of motion of molecular dynamics (MD) or its Langevin or Brownian variants. For example, in Langevin MD simulations our acceleration technique permits a straightforward spectral decomposition of forces so that the long-wavelength modes are integrated with a longer time step, thereby reducing the time required to reach equilibrium or to decorrelate the system in equilibrium. Speedup factors of up to 30 are observed relative to pure (unaccelerated) Langevin MD. As with acceleration of critical lattice models, even further gains relative to the unaccelerated method are expected for larger systems. Preliminary results for Fourier-accelerated molecular dynamics are presented in order to illustrate the basic concepts. Possible extensions of the method and further lines of research are discussed.

Journal Article↗
PubMed2001

Mobility of the active site bound paraoxon and sarin in zinc-phosphotriesterase by molecular dynamics simulation and quantum chemical calculation.

The kinetic data published on phosphotriesterase (PTE), with various complexed metals, clearly indicates that the P=O and P=S bonds of phosphotriester and thiophosphotriester substrates, respectively, are strongly polarized by one or both of the active site complexed metal ions. However, this observation is not consistent with the three-dimensional X-ray crystal structure of zinc-substituted PTE with active site bound substrate analogue diethyl 4-methylbenzylphosphonate. In this structure, the distance between the phosphoryl oxygen and the nearest zinc is 3.4 A, a distance too large to afford strong polarization. In the present paper, the geometry and mobility of various PTE active site-substrate complexes are examined by performing both molecular dynamics (MD) simulations and quantum mechanical calculations. Two known substrates are considered, paraoxon and sarin, although their turnover rates vary about 100-fold. The results indicate that PTE forms a complex with either substrate in which the phosphoryl oxygen becomes strongly coordinated with the less buried zinc atom. It is shown that the geometry of the active site is changed when the protein is immersed in a water bath and relaxed by MD. The most substantial conformational change is the opening of the gateway in a pocket where the location of the leaving group is expected. The opening is observed for the pure enzyme as well as for the enzyme/substrate complexes and it ranges from 11 to 18 A. It is also shown that the pockets, in which the substrate substituents are localized, exhibit different flexibility and interact with the substrate with coordinated conformational adjustments.

Aryldialkylphosphatase↗
PubMed2006

Ion permeation through a narrow channel: using gramicidin to ascertain all-atom molecular dynamics potential of mean force methodology and biomolecular force fields.

We investigate methods for extracting the potential of mean force (PMF) governing ion permeation from molecular dynamics simulations (MD) using gramicidin A as a prototypical narrow ion channel. It is possible to obtain well-converged meaningful PMFs using all-atom MD, which predict experimental observables within order-of-magnitude agreement with experimental results. This was possible by careful attention to issues of statistical convergence of the PMF, finite size effects, and lipid hydrocarbon chain polarizability. When comparing the modern all-atom force fields of CHARMM27 and AMBER94, we found that a fairly consistent picture emerges, and that both AMBER94 and CHARMM27 predict observables that are in semiquantitative agreement with both the experimental conductance and dissociation coefficient. Even small changes in the force field, however, result in significant changes in permeation energetics. Furthermore, the full two-dimensional free-energy surface describing permeation reveals the location and magnitude of the central barrier and the location of two binding sites for K(+) ion permeation near the channel entrance--i.e., an inner site on-axis and an outer site off-axis. We conclude that the MD-PMF approach is a powerful tool for understanding and predicting the function of narrow ion channels in a manner that is consistent with the atomic and thermally fluctuating nature of proteins.

Cell Membrane Permeability↗
PubMed2004

Toward the understanding of the structure and dynamics of protein-carbohydrate interactions: molecular dynamics studies of the complexes between hevein and oligosaccharidic ligands.

Herein we study, through all atom molecular dynamics simulations, the complex between hevein and two N-acetylated chitin oligomers, namely N,N(')-diacetylchitobiose and N,N('),N(")-triacetylchitotriose. The results of the simulations for two disaccharide complexes and one trisaccharide complex show that a carbohydrate oligomer is able to move on the surface of the relatively flat binding pocket of hevein, therefore occupying different binding subpockets. Statistical analysis methods were also applied in order to define the principal overall motions in the complexes, showing how the different ligands in the simulations modulate the protein motions. The oligosaccharide binding can be considered as defined by a subtle balance between enthalpic (formation of intermolecular interactions between the ligand and the receptor) and entropic (due mainly to the possibility for the sugar to move on the surface of the protein domain) effects, determining multiple binding conformations. This structural and dynamical view could parallel the results obtained by regularly used restrained MD simulations based on NOE NMR data that provide a well defined structure for both the disaccharide and trisaccharide complexes, and agrees with the observations for longer oligosaccharide chains.

Antimicrobial Cationic Peptides↗
PubMed2006

Solvation pressure in chloroform.

Molecular dynamics (MD) simulations of chloroform vapor and liquid at normal temperature and pressure and liquid under hydrostatic pressure are presented, giving bond lengths and vibrational frequencies as functions of pressure. The change in bond lengths between vapor and liquid at normal temperature and pressure is consistent with a pressure equivalent to the cohesive energy density (CED) of the liquid, supporting the solvation pressure model which predicts that solvated molecules or nanoparticles experience a pressure equal to the CED of the liquid. Experimental data for certain Raman frequencies of chloroform in the vapor phase, in the liquid, and in the liquid under pressure are presented and compared to MD. Results for C-Cl vibrational modes are in general agreement with the solvation pressure model whereas frequencies associated with the C-H bond are not. The results demonstrate that masking interactions exist in the real liquid that can be reduced or eliminated in simplified simulations.

Journal Article↗
PubMed2006

Toward understanding the inactivation mechanism of monooxygenase P450 BM-3 by organic cosolvents: a molecular dynamics simulation study.

Cytochrome P450 BM-3 from Bacillus megaterium is an extensively studied enzyme for industrial applications. A major focus of current protein engineering research is directed to improving the catalytic performance of P450 BM-3 toward nonnatural substrates of industrial importance in the presence of organic solvents or cosolvents. For the latter reason, it is important to study the effect of organic cosolvent molecules on the structure and dynamics of the enzyme, in particular, the effect of cosolvent molecules on the active site's structure and dynamics. In this paper, we have studied, using molecular dynamics (MD) simulations, the F87A mutant of P450 BM-3 in the presence of DMSO as cosolvent, to understand the role of the F87A substitution for its catalytic activity. This mutant exhibits an altered regioselectivity and substrate specificity compared with wild-type; however, it has lower tolerance toward DMSO. The simulation results offer an explanation for the DMSO sensitivity of the F87A mutant. Our simulation results show that the F87 side chain prevents the disturbance of the water molecule bound to the heme iron by DMSO molecules. The absence of the phenyl ring in F87A mutant promotes interactions of the DMSO molecule with the heme iron resulting in water displacement by DMSO at the catalytic heme center.

Bacillus megaterium↗