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,441 records · Page 80Linked to original sources
PubMed2001

Ligand diffusion in the catalase from Proteus mirabilis: a molecular dynamics study.

The role of the channels and cavities present in the catalase from Proteus mirabilis (PMC) was investigated using molecular dynamics (MD) simulations. The reactant and products of the reaction, H(2)O(2) -->1/2 O(2) + H(2)O, catalyzed by the enzyme were allowed to diffuse to and from the active site. Dynamic fluctuations in the structure are found necessary for the opening of the major channel, identified in the X-ray model, which allows access to the active site. This channel is the only pathway to the active site observed during the dynamics, and both the products and reactant use it. H(2)O and O(2) are also detected in a cavity defined by the heme and Ser196, which could play an important role during the reaction. Free energy profiles of the ligands diffusing through the major channel indicate that the barriers to ligand diffusion are less than 20 kJ mol(-1) for each of the species. It is not clear from our study that minor channels play a role for access to the protein active site or to the protein surface.

Binding Sites↗
PubMed2026

A Novel De Novo STAG1 Variant at the RAD21 Binding Interface Is Associated With Hypoglycemia, Recurrent Fever, Immunodeficiency and Features of Classical Cohesinopathies.

The cohesin complex, composed of SMC1, SMC3, RAD21, and STAG1/STAG2, is essential for chromosome cohesion, DNA repair, and transcriptional regulation. Pathogenic variants in cohesin components cause cohesinopathies. The classical characteristics of cohesinopathies include developmental delay (DD), intellectual disability (ID), feeding difficulties, hypotonia, short stature, hearing loss, and dysmorphic features. Here, we present a 5-year-old boy with classical cohesinopathy features, including DD/ID and feeding difficulties, along with non-classical features such as hypoglycemia, recurrent fever, and immunodeficiency. Trio exome sequencing identified a novel de novo missense variant of uncertain significance (NM_005862.3:c.643G>A(p.Val215Ile)) in the STAG1 gene. The variant localizes to the RAD21 interaction interface, and molecular dynamics (MD) simulations revealed conformational changes comparable to other STAG1 variants reported as likely pathogenic in patients, supporting a deleterious effect which may disrupt the STAG1-RAD21 interaction interface. This case expands the phenotypic and molecular spectrum of STAG1-related cohesinopathy and advances our understanding of the disease mechanism.

STAG1↗
PubMed1994

Synthesis and solution structure of an S-glycosylated cyclic hexapeptide. Evidence for conformational change induced by glycosylation.

Synthesis and conformational analysis of the S-glycosylated cyclic hexapeptide cyclo(-D-Pro1-Phe2-Cys3(tetra-O-acetyl-beta-D-galactopyranosyl)-++ +Trp4-Lys(Z)5- Phe6-) I was carried out to examine the influence of a saccharide residue in position i of a standard beta-turn on the formation of reverse turns and on the biological activity. Synthesis was carried out in the liquid phase employing a galactosylated cysteine building block. The cyclization reagents DPPA/NaHCO3 avoided high dilution conditions. Spectroscopic data were extracted from homo- and heteronuclear 2D-NMR techniques (TOCSY, NOESY, HMQC, HMQC-TOCSY, HMBCS-270). For structural refinement restrained molecular dynamics (MD) simulations in vacuo and with explicit DMSO as solvent were performed. Finally, simulations in DMSO without experimental restraints provided insight in stability and dynamics of the structural model. A comparison of the S-glycosylated Cys3 peptide with the analogous Thr3 peptide exhibits a similar overall conformation of the hexapeptide [beta II' D-Pro-Phe and another beta-turn about Trp4-Lys5(Z)]. However, the latter shows a distinct dynamic flip beta I, beta II in the glycopeptide, whereas the Thr-analogue only populates beta I. This influence is attributed to a beta I stabilizing effect of a hydrogen bridge of Thr-O gamma in position i to the NH of the amino acid in position i + 2, which is lacking in the glycosylated compound.

Amino Acid Sequence↗
PubMed1995

Influence of serine in position i on conformation and dynamics of reverse turns.

NMR spectroscopy has been employed for the conformational analysis of the cyclic hexapeptide cyclo(-D-Pro1-Ala2-Ser3(Bzl)-Trp4-Orn5(Z)-Tyr 6-) with and without protecting groups on Ser3 and Orn5. This peptide sequence was derived from the active loop sequence of the alpha-amylase inhibitor Tendamistat (HOE 467). The aim was to investigate the role of serine in position i of a standard beta-turn on the conformation and stabilization of this turn. Based on distance and torsion constraints from 2D NMR spectroscopic measurements in DMSO-d6 solution, structure refinement was accomplished by restrained molecular dynamics (MD) simulations in vacuo and in DMSO. The analysis of both structures in solution reveals a considerable effect of the unprotected serine sidechain on the adjacent beta-turn conformation. While in the protected peptide with Ser3(Bzl) a beta II-turn is observed between Trp4 and Orn5, the deprotected compound reveals a beta I-turn in this region. The beta I-turn is stabilized by a backbone-sidechain hydrogen bond from Orn5N alpha H to Ser3O gamma. Comparisons with other NMR-derived solution structures of cyclic model peptides and in some protein structures from literature reveal a general structural motif in the stabilization of beta I-turns by serine in the i position through backbone-sidechain interactions.

Amino Acid Sequence↗
PubMed2006

The natural mutation by deletion of Lys9 in the thrombin A-chain affects the pKa value of catalytic residues, the overall enzyme's stability and conformational transitions linked to Na+ binding.

The catalytic competence of the natural thrombin mutant with deletion of the Lys9 residue in the A-chain (deltaK9) was found to be severely impaired, most likely due to modification of the 60-loop conformation and catalytic triad geometry, as supported by long molecular dynamics (MD) simulations in explicit water solvent. In this study, the pH dependence of the catalytic activity and binding of the low-molecular mass inhibitor N-alpha-(2-naphthylsulfonyl-glycyl)-4-amidinophenylalanine-piperidine (alpha-NAPAP) to the wild-type (WT) and deltaK9 thrombin forms were investigated, along with their overall structural stabilities and conformational properties. Two ionizable groups were found to similarly affect the activity of both thrombins. The pKa value of the first ionizable group, assigned to the catalytic His57 residue, was found to be 7.5 and 6.9 in ligand-free deltaK9 and WT thrombin, respectively. Urea-induced denaturation studies showed higher instability of the deltaK9 mutant compared with WT thrombin, and disulfide scrambling experiments proved weakening of the interchain interactions, causing faster release of the reduced A-chain in the mutant enzyme. The sodium ion binding affinity was not significantly perturbed by Lys9 deletion, although the linked increase in intrinsic fluorescence was lower in the mutant. Essential dynamics (ED) analysis highlighted different conformational properties of the two thrombins in agreement with the experimental conformational stability data. Globally, these findings enhanced our understanding of the perturbations triggered by Lys9 deletion, which reduces the overall stability of the molecule, weakens the A-B interchain interactions, and allosterically perturbs the geometry and protonation state of catalytic residues of the enzyme.

Catalytic Domain↗
PubMed2006

Construction of a small peptide library related to inhibitor OM99-2 and its structure-activity relationship to beta-secretase.

AIM: To develop probes for detecting the binding specificity between beta-secretase and substrate, and provide reliable biological activity data for further researching encircling substrate-based inhibitors. METHODS: To prepare the inhibitors, the hydroxyethylene (HE) segment including P1 and P1'was synthesized after multi-step reactions; the combination of all segments was then completed through solid phase synthesis. Recombinant human beta-secretase ectodomain (amino acid residues 1-460) was expressed as a secreted protein with a C-terminal His tag in insect cells using baculovirus infection, and all compounds were evaluated in this beta-secretase enzyme assay. In order to understand the interaction in detail, the theoretical methods, namely molecular dynamics (MD) simulation and molecular mechanics-generalized-born surface area (MM-GBSA) analysis, were performed on the complex of beta-secretase and OM99-2 to obtain the geometrical and energetical information. RESULTS: We designed and constructed a positional scanning combinatorial library including 16 compounds; all members of the library were synthesized based on HE dipeptide isostere. Structure-activity relationship studies at the P4-P1 and P1' -P4'positions led to the discoveries of P and P'sides binding specificity and potent inhibitors 14, 18, and 22. The binding free energy on the whole system and every residue were compared to the biological assay result. CONCLUSION: The removal of P4' yielded inhibitor 22 (A3 *B2) with high potency; further truncation of P3'gave inhibitor 18 (A3 *B1) with equal activity, implying that the right side of the inhibitors play a less important role and could be easily simplified, while change on the P side may cause substantial results.

Amyloid Precursor Protein Secretases↗
PubMed2002

Static and dynamic properties of supercooled thin polymer films.

The dynamic and static properties of a supercooled (non-entangled) polymer melt are investigated via molecular-dynamics (MD) simulations. The system is confined between two completely smooth and purely repulsive walls. The wall-to-wall separation (film thickness), D, is varied from about 3 to about 14 times the bulk radius of gyration. Despite the geometric confinement, the supercooled films exhibit many qualitative features which were also observed in the bulk and could be analyzed in terms of mode-coupling theory (MCT). Examples are the two-step relaxation of the incoherent intermediate scattering function, the time-temperature superposition property of the late time alpha-process and the space-time factorization of the scattering function on the intermediate time scale of the MCT beta-process. An analysis of the temperature dependence of the alpha-relaxation time suggests that the critical temperature, T(c), of MCT decreases with D. If the confinement is not too strong ( D>or=10monomer diameter), the static structure factor of the film coincides with that of the bulk when compared for the same distance, T - T(c)( D), to the critical temperature. This suggests that T - T(c)( D) is an important temperature scale of our model both in the bulk and in the films.

Journal Article↗
PubMed2002

Water structure around enkephalin near a GeO2 surface: a molecular dynamics study.

A molecular model was created consisting of leucine enkephalin (a pentapeptide protein) near a germanium dioxide (GeO2) surface surrounded by water molecules. A molecular dynamic (MD) simulation of the system was conducted, and forces exerted by the water on both the surface and the protein were calculated. Orientational and spatial distribution functions of the water were calculated and were examined to determine if structured water existed and how it contributed to the intermolecular forces. Results from the study demonstrated that there is a strong spatial and orientational structuring of water near the GeO2 surface and that it is dependent on the proximity of the protein to the surface. Only minimal structuring occurred near the protein. A linear correlation was observed between the force of water on the protein and the angular distribution of water molecules in the region with the highest spatial density. When the protein is oriented with its polar side toward the GeO2 surface, the water structure near the protein disrupts the typical structure of the water near the surface, causing a force that pushes the protein away from the surface. When the protein is oriented with its non-polar side toward the surface, there is only minimal disruption of the typical water structure, and the resulting net forces between surface and protein are attractive.

Adsorption↗
PubMed2001

Insight to structural subsite recognition in plant thiol protease-inhibitor complexes : understanding the basis of differential inhibition and the role of water.

BACKGROUND: This work represents an extensive MD simulation / water-dynamics studies on a series of complexes of inhibitors (leupeptin, E-64, E-64-C, ZPACK) and plant cysteine proteases (actinidin, caricain, chymopapain, calotropin DI) of papain family to understand the various interactions, water binding mode, factors influencing it and the structural basis of differential inhibition. RESULTS: The tertiary structure of the enzyme-inhibitor complexes were built by visual interactive modeling and energy minimization followed by dynamic simulation of 120 ps in water environment. DASA study with and without the inhibitor revealed the potential subsite residues involved in inhibition. Though the interaction involving main chain atoms are similar, critical inspection of the complexes reveal significant differences in the side chain interactions in S2-P2 and S3-P3 pairs due to sequence differences in the equivalent positions of respective subsites leading to differential inhibition. CONCLUSION: The key finding of the study is a conserved site of a water molecule near oxyanion hole of the enzyme active site, which is found in all the modeled complexes and in most crystal structures of papain family either native or complexed. Conserved water molecules at the ligand binding sites of these homologous proteins suggest the structural importance of the water, which changes the conventional definition of chemical geometry of inhibitor binding domain, its shape and complimentarity. The water mediated recognition of inhibitor to enzyme subsites (Pn.H2O.Sn) of leupeptin acetyl oxygen to caricain, chymopapain and calotropinDI is an additional information and offer valuable insight to potent inhibitor design.

Amino Acid Chloromethyl Ketones↗
PubMed2003

Mutational analysis of amino acid residues involved in catalytic activity of a family 18 chitinase from tulip bulbs.

We expressed chitinase-1 (TBC-1) from tulip bulbs (Tulipa bakeri) in E. coli cells and used site-directed mutagenesis to identify amino acid residues essential for catalytic activity. Mutations at Glu-125 and Trp-251 completely abolished enzyme activity, and activity decreased with mutations at Asp-123 and Trp-172 when glycolchitin was the substrate. Activity changed with the mutations of Trp-251 to one of several amino acids with side-chains of little hydrophobicity, suggesting that hydrophobic interaction of Trp-251 is important for the activity. Molecular dynamics (MD) simulation analysis with hevamine as the model compound showed that the distance between Asp-123 and Glu-125 was extended by mutation of Trp-251. Kinetic studies of Trp-251-mutated chitinases confirmed these various phenomena. The results suggested that Glu-125 and Trp-251 are essential for enzyme activity and that Trp-251 had a direct role in ligand binding.

Amino Acid Substitution↗
PubMed2005

Acidic sugar degradation pathways: an ab initio molecular dynamics study.

Ab initio molecular dynamics (MD) simulations were employed to elucidate xylose and glucose degradation pathways. In the case of xylose, a 2,5-anhydride intermediate was observed leading to the formation of furfural through elimination of water. This pathway agrees with one of the mechanisms proposed in the literature in that no open chain intermediates were found. In the case of glucose, a series of intermediates were observed before forming the 2,5-anhydride intermediate that eventually leads to hydroxymethylfurfural (HMF). One of these intermediates was a very short-lived open-chain form. Furthermore, two novel side-reaction pathways were identified, which lead to degradation products other than HMF.

Anhydrides↗
PubMed1992

Conformation and fluctuations of free stefin B: a molecular dynamics study.

Molecular dynamics study was performed on the cysteine proteinase inhibitor stefin B. Structure of inhibitor from the complex with papain was used as a starting point. Amino terminal "trunk" of the inhibitor which lies extended along the cleft of the enzyme in the complex, folded onto the body of inhibitor during MD simulation, thereby reducing the total and particularly hydrophobic surface exposed to the solvent. This effect counterbalances hydrophobic contribution of the "trunk" and explains why its deletion in stefin B and related inhibitors doesn't reduce the dissociation constant. The rest of stefin B conformation is conserved together with main chain hydrogen bonds. Fluctuations of C alpha atoms resembles crystallographic B factors with exception of residues in contact with enzyme.

Chemical Phenomena↗
PubMedDate not supplied

Ion transport through gramicidin A. Water structure and functionality.

Molecular Dynamics (MD) simulations were performed on a gramicidin A dimer model representing a transmembrane channel. Different from previous simulations the peptide was in contact with bulk water at both ends of the dimer to guarantee a realistic description of the hydration of the biomolecule. The flexible BJH model for water was employed in the simulations and the gramicidin-water, gramicidin-ion and ion-water potentials used are based on molecular orbital calculations. The water structure near the gramicidin was investigated first by a simulation without ions, while for the energy profiles of the ion transport through the channel a potassium or a sodium ion was added. These investigations provide a detailed and conclusive picture on a molecular level of the role of water in the ion transport through a gramicidin A channel and can explain the experimental results on the selectivity between alkali ions, their double or even triple occupancy, the exclusion or permeability of anions depending upon cation concentration and the consequences of differences in the ionic charge. The investigation demonstrate that the water molecules around the gramicidin behave as an integral part of the peptide and the functionality is the result of the whole complex biomolecule-water.

Anions↗
PubMed2005

Interfacial tryptophan residues: a role for the cation-pi effect?

Integral membrane proteins are characterized by having a preference for aromatic residues, e.g., tryptophan (W), at the interface between the lipid bilayer core and the aqueous phase. The reason for this is not clear, but it seems that the preference is related to a complex interplay between steric and electrostatic forces. The flat rigid paddle-like structure of tryptophan, associated with a quadrupolar moment (aromaticity) arising from the pi-electron cloud of the indole, interacts primarily with moieties in the lipid headgroup region hardly penetrating into the bilayer core. We have studied the interaction between the nitrogen moiety of lipid molecule headgroups and the pi-electron distribution of gramicidin (gA) tryptophan residues (W9, W11, W13, and W15) using molecular dynamics (MD) simulations of gA embedded in two hydrated lipid bilayers composed of 1-palmitoyl-2-oleoylphosphatidylethanolamine (POPE) and 1-palmitoyl-2-oleoylphosphatidyl-choline (POPC), respectively. We use a force field model for tryptophan in which polarizability is only implicit, but we believe that classical molecular dynamics force fields are sufficient to capture the most prominent features of the cation-pi interaction. Our criteria for cation-pi interactions are based on distance and angular requirements, and the results from our model suggest that cation-pi interactions are relevant for W(PE)1), W(PE)13, W(PE)15, and, to some extent, W(PC)11 and W(PC)13. In our model, W9 does not seem to engage in cation-pi interactions with lipids, neither in POPE nor POPC. The criteria for the cation-pi effect are satisfied more often in POPE than in POPC, whereas the H-bonding ability between the indole donor and the carbonyl acceptor is similar in POPE and POPC. This suggests an increased affinity for lipids with ethanolamine headgroups to transmembrane proteins enriched in interfacial tryptophans.

Amino Acids↗
PubMed2006

Role of protein flexibility in ion permeation: a case study in gramicidin A.

Proteins have a flexible structure, and their atoms exhibit considerable fluctuations under normal operating conditions. However, apart from some enzyme reactions involving ligand binding, our understanding of the role of flexibility in protein function remains mostly incomplete. Here we investigate this question in the realm of membrane proteins that form ion channels. Specifically, we consider ion permeation in the gramicidin A channel, and study how the energetics of ion conduction changes as the channel structure is progressively changed from completely flexible to a fixed one. For each channel structure, the potential of mean force for a permeating potassium ion is determined from molecular dynamics (MD) simulations. Using the same molecular dynamics data for completely flexible gramicidin A, we also calculate the average densities and fluctuations of the peptide atoms and investigate the correlations between these fluctuations and the motion of a permeating ion. Our results show conclusively that peptide flexibility plays an important role in ion permeation in the gramicidin A channel, thus providing another reason--besides the well-known problem with the description of single file pore water--why this channel cannot be modeled using continuum electrostatics with a fixed structure. The new method developed here for studying the role of protein flexibility on its function clarifies the contributions of the fluctuations to energy and entropy, and places limits on the level of detail required in a coarse-grained model.

Anti-Bacterial Agents↗
PubMed2026

EWS-RNA Binding Protein 1: Structural Insights into Ewing Sarcoma by Conformational Dynamics Investigations.

BACKGROUND: Prior research has demonstrated that proteins play a significant role in the prognosis and treatments of various sarcomas, including Ewing sarcoma through the interplay of downstream signaling cascades. However, there is limited understanding about the strcucture conformation of EWSR1 and its structural implication in the prognosis of Ewsing Sarcoma by interaction with RNA molecules. AIMS: The primary goal of ongoing research is to determine how EWSR1 contributes to Ewing sarcoma. OBJECTIVE: The current study explores the complexity of EWSR1 structure and its conformational interactions with RNA in relation to Ewing sarcoma. METHODS: Here, we employed a comparative modeling approach to predict EWSR1 domains separately and assembled them into one structural unit using a DEMO server. Additionally, the RNA motifs interacting with EWSR1 were predicted, and the 3D model was built using RNAComposer. Protein-RNA docking and MD simulation studies were carried out to check the intermolecular interactions and stability behavior of docked EWSR1-RNA complexes. RESULTS: The overall results explore the structural insights into EWSR1 and their interactions with RNA, which may play a momentous role in co- and post-transcriptional regulation to control gene expression. CONCLUSION: Taken togather, our findings suggest that EWSR1 may be a useful therapeutic target for the diagnosis and management of Ewing sarcoma.

Sarcoma, Ewing↗
PubMed2005

Definition of an uptake pharmacophore of the serotonin transporter through 3D-QSAR analysis.

The serotonergic system plays a critical role in a wide variety of physiological and behavioral processes. Dysregulation of the tightly controlled extracellular concentration of serotonin (5-hydroxytryptamine, 5-HT) appears to be at the origin of a host of metabolic and psychiatric disorders. Since the plasma membrane 5-HT transporter (SERT) is the major protagonist in regulating extracellular 5-HT concentration, SERT is the target of most drugs interacting with the serotonergic system. Unfortunately, some of the drugs towards SERT (e.g. amphetamine derivatives) interfere with cell homeostasis leading to cell toxicity. Developing new SERT ligands devoid of any side-effect represents a major priority in the treatment of 5-HT-associated pathologies. Here, we report structure-activity relationships (SAR) and three-dimensional QSAR (3D-QSAR) studies of a library of 121 compounds including 5-HT analogs, harmanes, benzothiazoles, indanones, amphetamine derivatives and substrate-type 5-HT releasers, with the goal of identifying the structural determinants crucial for SERT uptake. In the absence of data about the bioactive form of 5-HT, conformational analysis of 5-HT was performed using quantum chemistry calculations. This led to three 5-HT stable conformers with anti, -gauche and +gauche side-chain conformation. These conformers, used as templates for superimposition with all the library compounds, enabled the design of a reliable 6-points pharmacophore representative of SERT uptake activity. Molecular dynamics (MD) simulations performed with compounds that are efficiently, moderately, poorly or not transported by SERT allowed to assess the validity of our pharmacophore. Altogether, our data provide for the first time a reliable pharmacophore of SERT uptake activity, which may help to the design of new drugs targeting SERT.

Adult↗
PubMed2006

Recent advances in molecular modeling and medicinal chemistry aspects of phospho-glycoprotein.

Phospho-glycoprotein (P-gp) is an efflux transporter expressed in many organs (ex: kidney, lung, liver and spleen) and in hormone producing or responsive tissues (ex: adrenal cortex, testis and placenta). It is involved in many important physiological functions. Among them the major one is extrusion of xenobiotics in order to detoxify the cells. This property of P-gp is associated with multidrug resistance (MDR) for many pathological conditions. While the experimental determination of three-dimensional structure is not yet successful, the transmembrane (TM) 5, 6, 11 and 12 are sensitive to mutations and contain substrate binding sites. Designing of potential and selective inhibitors of P-gp is still hampered by a lack of information upon the three dimensional structure of P-gp. The design of P-gp inhibitors was traditionally driven by quantitative structure activity relationship studies, which is complicated by factors such as different types of assays, multiple drug binding sites and diverse chemical structures. Clearly a conclusive and predictive SAR does not seem to be practical, despite progress in the last few years towards more specific SAR suggesting well defined structural features responsible for activity. Advances made recently in solving the crystal structure of prokaryotic ATP binding cassette proteins (ABC) transporters, Ec-MsbA, Vc-MsbA and BtuCD yielded suitable templates for construction of homology models of P-gp. Few molecular dynamics (MD) simulations aimed at elucidating the functional dynamics of ABC transporters have provided useful insights to their mechanism and structure. The present review aims at the general overview of importance, expression, structure, organization and drug binding sites of P-gp. This review also highlights recent developments in the homology modeling, molecular dynamics simulations of P-gp and progress in QSAR, pharmacophore modeling of P-gp modulators.

ATP Binding Cassette Transporter, Subfamily B, Mem↗