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,207 records · Page 67Linked to original sources
PubMed2005

Protonation state of Asp30 exerts crucial influence over surface loop rearrangements responsible for NO release in nitrophorin 4.

pK(a) values of ionizable residues were calculated for the crystal structures describing the pH and NO binding dependant conformations of nitrophorin 4, a pH sensitive NO carrier heme protein. Comparison of resultant H-bonding patterns allowed the identification of the amino acids that take part in signaling pH change. We carried out MD simulations to show that the protonation state of Asp30, buried in the closed conformation, is crucial for maintaining the tight packed conformation of the closed form of the complex - presenting a model for the functional decrease of NO binding affinity of nitrophorins at physiological pH.

Asparagine↗
PubMed2026

Inhibitory mechanism of phloretin on the AgrA LytTR domain-agr operon complex formation and its application in beef.

Staphylococcus aureus (S. aureus) represents a major foodborne pathogen whose enterotoxin production poses significant challenges to food safety due to its high environmental resistance and limited efficacy of conventional sterilization. Since the expression of enterotoxins is predominantly governed by the agr quorum sensing system, targeting this regulatory pathway has become a strategic choice for virulence control. This study elucidated the mechanism by which phloretin, a potential quorum sensing inhibitor, interferes with the agr system to attenuate virulence. To achieve this, the recombinant AgrA LytTR domain was expressed and purified, and its interaction with phloretin was characterized using thermal shift assays (TSA), electrophoretic mobility shift assays (EMSA), and molecular dynamics (MD) simulations. The results showed that phloretin specifically binds to the AgrA LytTR domain, enhancing its thermal stability and disrupting AgrA LytTR-agr operon binding by reducing the free energy of interaction between them, without causing significant structural rearrangement. Mechanistic analysis indicated that phloretin sterically hinders key β-sheet turn residues (HIS169, ASN201, ARG233), thereby impairing DNA recognition, downregulating RNAIII transcription, and inhibiting agr signaling. In cooked beef, phloretin significantly inhibited the secretion of enterotoxins and α-hemolysin, while delaying lipid oxidation and protein degradation, and maintaining the meat texture. These findings suggested that phloretin is a multifunctional substance with anti-virulence, antioxidant, and preservative properties, demonstrating its potential as a natural food preservative.

Phloretin↗
PubMed2026

Dual solvent cavities and hydrogen-bond networks define the chromophore environment in a far-red/orange-sensing cyanobacteriochrome.

Cyanobacteriochromes (CBCRs) are bilin-binding photoreceptors that exhibit remarkable spectral diversity and mediate light-dependent signaling in cyanobacteria. Far-red/orange-sensing CBCRs (froCBCRs) have attracted interest because of their unusually red-shifted absorption properties, yet structural information for their illuminated states has been lacking. Here, we report the first high-resolution (1.8 Å) crystal structure of the orange-absorbing (Po) state of the froCBCR ToFrO from Tolypothrix sp. PCC 7910. The structure reveals a compact, cyclic bilin configuration and water-mediated hydrogen-bonding networks within two solvent-accessible cavities. Within the GAF domain, the D-ring remains nearly perpendicular to the planar A-to-C ring system through interactions involving a flexible loop region. Comparative analyses of cryogenic synchrotron and room-temperature X-ray free-electron laser (XFEL) structures, together with molecular dynamics (MD) simulations, revealed alternative Met636 conformations associated with dynamic water exchange through the solvent-accessible cavity. Site-directed mutagenesis of cavity-lining and water-interacting residues resulted in modest spectral shifts. By contrast, mutations of two Trp residues, participating in π-π stacking with the D-ring and likely imposing a steric constraint near the A-ring, resulted in substantial blue and red shifts in the dark and illuminated states, respectively. Together with the observed chromophore geometry, these findings indicate that the spectral properties of ToFrO are governed by chromophore conformation and its direct interaction with surrounding residues through hydrogen-bonding, electrostatic, and π-π interactions. These results further suggest that cavity-mediated solvent organization contributes to stabilizing the local structural environment surrounding the chromophore and adjacent protein backbone. Collectively, these findings elucidate the structural basis for photoconversion and spectral tuning in froCBCRs.

Cyanobacteria↗
PubMed2005

Thermophoresis in liquids: a molecular dynamics simulation study.

Thermophoresis in liquids is studied by molecular dynamics simulation (MD). A theory is developed that divides the problem in the way consistent with the characteristic scales. MD is then conducted to obtain the solution of each problem, which is to be all combined for macroscopic predictions. It is shown that when the temperature gradient is applied to the nonconducting liquid bath that contains neutral particles, there occurs a pressure gradient tangential to the particle surface at the particle-liquid interface. This may induce the flow in the interfacial region and eventually the particle to move. This applies to the material system that interacts through van der Waals forces and may be a general source of the thermophoresis phenomenon in liquids. The particle velocity is linearly proportional to the temperature gradient. And, in a large part of the given temperature range, the particle motion is in the direction toward the cold end and decreases with respect to the temperature. It is also shown that the particle velocity decreases or even reverses its sign in the lowest limit of the temperature range or with a particle of relatively weak molecular interactions with the liquid. The characteristics of the phenomenon are analyzed in molecular details.

Journal Article↗
PubMed2005

Topological frustration and the folding of interleukin-1 beta.

The cytokine, interleukin-1beta (IL-1beta), adopts a beta-trefoil fold. It is known to be much slower folding than similarly sized proteins, despite having a low contact order. Proteins are sufficiently well designed that their folding is not dominated by local energetic traps. Therefore, protein models that encode only the folded structure and are energetically unfrustrated (Gō-type), can capture the essentials of the folding routes. We investigate the folding thermodynamics of IL-1beta using such a model and molecular dynamics (MD) simulations. We develop an enhanced sampling technique (a modified multicanonical method) to overcome the sampling problem caused by the slow folding. We find that IL-1beta has a broad and high free energy barrier. In addition, the protein fold causes intermediate unfolding and refolding of some native contacts within the protein along the folding trajectory. This "backtracking" occurs around the barrier region. Complex folds like the beta-trefoil fold and functional loops like the beta-bulge of IL-1beta can make some of the configuration space unavailable to the protein and cause topological frustration.

Computer Simulation↗
PubMed2006

Structural studies of a two-domain chitinase from Streptomyces griseus HUT6037.

Chitinase C (ChiC) from Streptomyces griseus HUT6037 was the first glycoside hydrolase family 19 chitinase that was found in an organism other than higher plants. An N-terminal chitin-binding domain and a C-terminal catalytic domain connected by a linker peptide constitute ChiC. We determined the crystal structure of full-length ChiC, which is the only representative of the two-domain chitinases in the family. The catalytic domain has an alpha-helix-rich fold with a deep cleft containing a catalytic site, and lacks three loops on the domain surface compared with the catalytic domain of plant chitinases. The chitin-binding domain is an all-beta protein with two tryptophan residues (Trp59 and Trp60) aligned on the surface. We suggest the binding mechanism of tri-N-acetylchitotriose onto the chitin-binding domain on the basis of molecular dynamics (MD) simulations. In this mechanism, the ligand molecule binds well on the surface-exposed binding site through two stacking interactions and two hydrogen bonds and only Trp59 and Trp60 are involved in the binding. Furthermore, the flexibility of the Trp60 side-chain, which may be involved in adjusting the binding surface to fit the surface of crystalline chitin by the rotation of chi2 angle, is shown.

Amino Acid Sequence↗
PubMed2004

Application of the linear interaction energy method (LIE) to estimate the binding free energy values of Escherichia coli wild-type and mutant arginine repressor C-terminal domain (ArgRc)-l-arginine and ArgRc-l-citrulline protein-ligand complexes.

Protein-ligand binding free energy values of wild-type and mutant C-terminal domain of Escherichia coli arginine repressor (ArgRc) protein systems bound to L-arginine or L-citrulline molecules were calculated using the linear interaction energy (LIE) method by molecular dynamics (MD) simulation. The binding behaviour predicted by the dissociation constant (K(d)) calculations from the binding free energy values showed preferences for binding of L-arginine to the wild-type ArgRc but not to the mutant ArgRc(D128N). On the other hand, L-citrulline do not favour binding to wild-type ArgRc but prefer binding to mutant ArgRc(D128N). The dissociation constant for the wild-type ArgRc-L-arginine complex obtained in this study is in agreement with reported experimental results. Our results also support the experimental data for the binding of L-citrulline to the mutant ArgRc(D128N). These showed that LIE method for protein-ligand binding free energy calculation could be applied to the wild-type and the mutant E. coli ArgRc-L-arginine and ArgRc-L-citrulline protein-ligand complexes and possibly to other transcriptional repressor-co-repressor systems as well.

Arginine↗
PubMed2004

New perspectives on multiple-copy, mean-field molecular dynamics methods.

Mean-field molecular dynamics (MD) techniques are designed to improve phase-space sampling in MD simulations. Reviewed here are theoretical and practical contributions from our group, and the ideas our contributions are based upon, beginning with the original time-dependent Hartree technique, locally enhanced sampling (LES), and a new, purely classical derivation of multiple-copy, mean-field equations of motion. This view is used to provide new insights into approximations inherent in LES, as well as the basis for a recently proposed method, ensembles eXtracted from atomic coordinate transformations (the EXACT approximation).

Algorithms↗
PubMed2006

Predicting proteinase specificities from free energy calculations.

The role of the primary binding residue (P1) in complexes between three different subtilases (subtilisin Carlsberg, thermitase and proteinase K) and their canonical protein inhibitor eglin c have been studied by free energy calculations. Based on the crystal structures of eglin c in complex with subtilisin Carlsberg and thermitase, and a homology model of the eglin c-proteinase K complex, a total of 57 mutants have been constructed and docked into their host proteins. The binding free energy was then calculated using molecular dynamics (MD) simulations combined with the linear interaction energy (LIE) method for all complexes differing only in the nature of the amino acid at the P1 position. LIE calculations for 19 different complexes for each subtilase were thus carried out excluding proline. The effects of substitutions at the P1 position on the binding free energies are found to be very large, and positively charged residues (Arg, Lys and His) are particularly deleterious for all three enzymes. The charged variants of the acidic side chains are found to bind more favorably as compared to their protonated states in all three subtilases. Furthermore, hydrophobic amino acids are accommodated most favorably at the S1-site in all three enzymes. Comparison of the three series of binding free energies shows only minor differences in the 19 computed relative binding free energies among these subtilases. This is further reflected in the correlation coefficient between the 23 relative binding free energies obtained, including the possible protonation states of ionizable side chains, but excluding the P1 Pro, for subtilisin Carlsberg versus thermitase (0.95), subtilisin versus proteinase K (0.94) and thermitase versus proteinase K (0.96).

Algorithms↗
PubMed2006

Systematic molecular dynamics searching in a lipid bilayer: application to the glycophorin A and oncogenic ErbB-2 transmembrane domains.

Molecular dynamics (MD) simulations of proteins in a lipid bilayer environment are usually undertaken with one or a few starting structures. Here we report a search protocol for systematically exploring the possible interactions in helical bundle transmembrane proteins, a frequently occurring structural motif. The search protocol correctly identifies the experimentally known structure of the dimeric human glycophorin A transmembrane domain as the lowest energy structure among five different models without any prior assumptions, whilst an identical in vacuo search fails to identify the correct structure. The lowest energy structure from the search in a lipid bilayer has a root mean square deviation of 1.1A to the experimental structure. We have applied the same search protocol to the unknown transmembrane structure of the oncogenic mutant ErbB-2 protein, a member of the family of epidermal growth factor receptors. Resulting structures show the role of glutamic acid hydrogen bonding and close helical packing. Water molecules may also play a key role in stabilisation of the transmembrane helix association.

Amino Acid Sequence↗
PubMed2006

Homology modeling and molecular dynamics study of chorismate synthase from Shigella flexneri.

Shigellosis is a major public health problem in many developing countries. Antibiotic therapy can reduce the severity of the dysentery and prevent potentially lethal complication. However, owing to the increased resistance to most of the widely used and inexpensive antibiotics, there is an urgent need for new antibacterial agents, particularly those that act on novel targets. Chorismate synthase (CS) is a key enzyme in the shikimic acid pathway, which is essential for the synthesis of aromatic amino acids in bacteria. As an anti-bacterial drug target, CS has been well validated. A homology model of Shigella-CS with the flavin mononucleotide (FMN) binding was constructed using the crystal structure of CS from other species. The substrate 5-enolpyruvylshikimate 3-phosphate (EPSP) was subsequently docked into the active site based on previous theoretical studies. Molecular dynamics (MD) was used to refine the starting ternary model. The model was well conserved during the 1.8 ns MD simulation with the equilibrium root mean square deviation (RMSD) value of 3.5 angstrom. The substrate binding energy was calculated and the electrostatic energy was found to be the most important term for binding. Decomposition of binding energies revealed that R129, R125, R327, R134 and R48 are important residues involved in substrate binding, which is useful for further site-directed mutagenesis experiments. In the absence of crystal structure, our study provides an early insight into the structure of CS from Shigella flexneri and its binding to the substrate and cofactor, thus facilitating the inhibitor design.

Amino Acid Sequence↗
PubMed2006

Exploring the molecular basis of the enantioselective binding of penicillin G acylase towards a series of 2-aryloxyalkanoic acids: a docking and molecular dynamics study.

In the present paper, molecular modeling studies were undertaken in order to shed light on the molecular basis of the observed enantioselectivity of penicillin G acylase (PGA), a well known enzyme for its industrial applications, towards 16 racemic 2-aryloxyalkanoic acids, which have been reported to affect several biological systems. With this intention docking calculations and MD simulations were performed. Docking results indicated that the (S)-enantiomers establish several electrostatic interactions with SerB1, SerB386 and ArgB263 of PGA. Conversely, the absence of specific polar interactions between the (R)-enantiomers and ArgB263 seems to be the main reason for the different binding affinities observed between the two enantiomers. Results of molecular dynamics simulations demonstrated that polar interactions are responsible for both the ligand affinity and PGA enantiospecificity. Modeling calculations provided possible explanations for the observed enantioselectivity of the enzyme that rationalize available experimental data and could be the basis for future protein engineering efforts.

Alkanes↗
PubMed2006

Optimization of electrostatics as a strategy for cold-adaptation: a case study of cold- and warm-active elastases.

Adaptation to both high and low temperatures requires proteins with special properties. While organisms living at or close to the boiling point of water need to have proteins with increased stability, other properties are required at temperatures close to the freezing point of water. Indeed, it has been shown that enzymes adapted to cold environments are less resistant to heat with a concomitant increased activity as compared to their warm-active counter-parts. Several recent studies have pointed in the direction that electrostatic interactions play a central role in temperature adaptation, and in this study we investigate the role such interactions have in adaptation of elastase from Atlantic salmon and pig. Molecular dynamics (MD) simulations have been used to generate structural ensembles at 283 and 310 K of the psychrophilic and mesophilic elastase, and a total of eight 12 ns simulations have been carried out. Even though the two homologues have a highly similar three-dimensional structure, the location and number of charged amino acids are very different. Based on the simulated structures we find that very few salt-bridges are stable throughout the simulations, and provide little stabilization/destabilization of the proteins as judged by continuum electrostatic calculations. However, the mesophilic elastase is characterized by a greater number of salt-bridges as well as a putative salt-bridge network close to the catalytic site, indicating a higher rigidity of the components involved in the catalytic cycle. In addition, subtle differences are also found in the electrostatic potentials in the vicinity of the catalytic residues, which may explain the increased catalytic efficiency of the cold-adapted elastase.

Acclimatization↗
PubMed2006

Structural, energetic, and dynamical properties of rotaxanes constituted of alpha-cyclodextrins and an azobenzene chain.

The [3]rotaxane synthesised as a single isomer constituted of two cyclodextrins (CDs) and an azobenzene chain [M.R. Craig, T.D.W. Claridge, M.G. Hutchings, H.L. Anderson, Synthesis of a cyclodextrin azo dye [3]rotaxane as a single isomer, Chem. Commun. 16 (1999) 1537-1538] has been investigated using molecular mechanics (MM) and dynamics (MD) with the MM3 force field in order to evaluate the stability of various configurations of the complex in the isolated and solvated states. The influence of the blocking groups and the presence of energy barriers along the azo chain were first investigated through the calculation of an energy profile. It revealed that the CD could translate along the chain at room temperature without any decomplexation. Next, MD simulations of three different types of configurations, i.e., head-to-head, head-to-tail, and tail-to-tail, of the [3]rotaxane were carried out. The non-solvated phase simulations showed structures with the CDs close to each other while the solvated ones showed structures with CDs separated by larger distances. This separation occurs due to the solute-solvent interactions. When the systems are in isolated state, the observed structure of the complexes are less stable due to an unfavourable arrangement of the hydroxyls groups of the adjacent CD faces. When considering solvation, energies of the three configurations are roughly identical due to the large distance between the faces of the CDs.

Azo Compounds↗
PubMed2026

In silico identification of DNMT1 inhibitors from the PlantCyc database through computational approach to assess the anti-cancer potential of nutraceutical compounds in breast cancer.

Breast cancer accounts for a disproportionate share of global cancer-related deaths, with 670,000 fatalities and 2.3 million new diagnoses recorded in women during 2022 alone. Existing treatment modalities carry considerable toxicity burdens, and resistance to available agents remains an unresolved clinical problem. DNA methyltransferase 1 (DNMT1), the enzyme chiefly responsible for maintaining genome-wide methylation patterns during DNA replication, has been mapped out as a high-value target in breast cancer because its dysregulation silences tumour suppressor genes through promoter hypermethylation. The present work involves hierarchical in silico workflow to screen 4549 plant-derived compounds from the PlantCyc database (v16.0.3) against the human DNMT1 catalytic domain (PDB ID: 4WXX). Ten top-scoring compounds were taken forward for molecular docking via AutoDock Vina; Quercetin and Kaempferol both recorded the highest binding affinities at -9.5 kcal/mol, Wogonin (-9.3 kcal/mol) and Xanthohumol (-8.1 kcal/mol) also emerged as strong binders. Pharmacokinetic evaluation using ADMET-AI confirmed that all 10 compounds met Lipinski's rule of five, with human intestinal absorption values at or above 0.98. Wogonin and Xanthohumol were selected for a 100 ns all-atom molecular dynamics (MD) simulation in GROMACS due to their well-rounded ADMET profiles and limited existing data on their specific interactions with DNMT1 in breast cancer. Across all measured trajectory metrics, backbone RMSD, residue fluctuation, radius of gyration, solvent-accessible surface area, and intermolecular hydrogen bond count, Wogonin formed a more stable, compact complex. These findings suggest that Wogonin and Xanthohumol are non-toxic nutraceutical candidates suitable for DNMT1 targeted epigenetic therapy, with computational foundation strong enough to facilitate future in vitro and in vivo validation work.

Humans↗
PubMed2005

Dynamics of prey-flock escaping behavior in response to predator's attack.

The dynamic behavior of prey-flock in response to predator's attack was investigated by using molecular dynamics (MD) simulations in a two-dimensional (2D) continuum model. By locally applying interactive forces between prey individuals (e.g. attraction, repulsion, and alignment), a coherently moving state in the same direction was obtained among individuals in prey-flock. When a single predator was introduced to the prey population, the prey-flock was correspondingly deformed by the predator's continuous attacks towards the prey-flock's center. In response to the predator's attack, three regimes in the flock size (compression (Regime I), expansion (Regime II), compression (Regime III)) were revealed if the predator's attack speed (kappa) was comparatively low to the escape speed of prey-flock. If noise was added to the predator's attacking course, a higher degree of variation was observed in the patterns of compression and expansion in the prey-flock size. However, the scaling behavior in the changes in prey-flock size was present in different levels of noise with the increase in predation risk (R) when kappa takes an appropriately low value. During the procedure of escaping, order breaking in alignment (phi) of prey population was observed, while the degree of alignment was dependent upon the changes in parameters of kappa and R.

Animals↗
PubMed2005

Molecular modeling of SWR-0342SA, a beta3-selective agonist, with beta1- and beta3-adrenoceptor.

The molecular dynamics (MD) simulations study in the formation of the complex between compound SWR-0342SA and beta-ARs suggested that upon binding SWR-0342SA stimulates receptor activation through residues network (Asp104, Leu335 in beta(1)-AR; Asp117, Ser209, Leu303, Ser191 in beta(3)-AR) in an active conformation state. The models suggest that the structural origin of the selectivity of SWR-0342SA to beta(3)-AR vs. beta(1)-AR comes from the following results: (a) the tight interaction between the agonist and the TMs 3, 5, 6 and 2 nd EC loop. Asp117 interacts with the cationic amino group of the agonist molecule. (b) Additional contacts are done with Ser209, Leu303 and Ser191. These results are in good agreement with the binding affinities (pKi values) of SWR-0342SA to beta-AR family expressed in recombinant mammalian cells.

Acetates↗
PubMed2005

Spectroscopy, NMR and DFT studies on molecular recognition of crown ether bridged chiral heterotrinuclear salen Zn(II) complex.

A barium-containing crown ether bridged chiral heterotrinuclear salen Zn(II) complex BaZn2L(ClO4)2, where L is a folded dinuclear chiral (R,R)-salen ligand, has been synthesized and characterized by elemental analysis, 1H NMR, UV-vis, IR, circular dichroism (CD) spectra, and mass spectra. As a folded dinuclear chiral host, its recognition with achiral guests (imidazole derivatives), rigid bidentate guest (1,4-diazobicyclo[2,2,2]octane, DABCO) and chiral guests (amino acid methyl esters) was investigated by means of UV-vis spectrophotometric titration, CD spectra. The association constants of D-amino acid methyl esters are found to be higher than those of their L-enantiomer. The sandwich-type binding of BaZn2L(ClO4)2-DABCO supramolecular assembly was specially studied via 1H NMR titration and 1H ROESY. To understand the recognition on molecular level, density functional theory (DFT) calculations on B3LYP/LanL2DZ were performed on the minimal energy conformations of host, guests, and host-guest complexes. The minimal energy conformations were obtained by molecular mechanics (MM) optimization and molecular dynamics (MD) simulation. The results of single point energy, HOMO energy, and charges transfer were analyzed. The results of theoretical calculations are in good agreement with the experimental data.

Barium↗