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.

SourceResearch connectionFrom yearThrough year

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

At least 631 records · Page 35Linked to original sources

PubMed2000

Diffusion of water molecules in crystalline beta-cyclodextrin hydrates.

To understand the rapid diffusion mechanism of water molecules in the crystal lattice of hydrated beta-cyclodextrin (beta-CD), molecular dynamics (MD) simulations of crystalline beta-CD were performed at five different relative humidities corresponding to hydration states ranging from beta-CD-9.4H2O to beta-CD-12.3H2O, and in aqueous solution. The trajectories for the crystalline beta-CD hydrates had lengths of 4 ns each, whereas the simulation in aqueous solution extended to 2 ns. Transport of water molecules in the crystal was characterized in terms of a spatially varying diffusion constant and the main direction of diffusion, which were evaluated using newly developed algorithms. The main diffusion pathway winds through the cavities of adjacent doughnut shaped beta-CD molecules and is slightly slanted with respect to the crystallographic b-axis. Water molecules outside the beta-CD cavities have access to the main diffusion pathway. The diffusion constant for transport of water molecules along the main pathway calculated from the MD simulation data adopts 1/30 of the value in bulk water at room temperature. This is in agreement with estimates that can be made from experimental data on the adjustment of a beta-CD crystal to changes in relative humidity.

Algorithms↗

PubMed2000

Molecular contacts in the transmembrane c-subunit oligomer of F-ATPases identified by tryptophan substitution mutagenesis.

When isolated in its monomeric form, subunit c of the proton transporting ATP synthase of Escherichia coli was shown to fold in a hairpin-like structure consisting of two hydrophobic membrane spanning helices and a short connecting hydrophilic loop. In the plasma membrane of Escherichia coli, however, about 9-12 c-subunit monomers form an oligomeric complex that functions in transmembrane proton conduction and in energy transduction to the catalytic F1 domain. The arrangement of the monomers and the molecular architecture of the complex were studied by tryptophan scanning mutagenesis and restrained MD simulations. Residues 12-24 of the N-terminal transmembrane segment of subunit c were individually substituted by the large and moderately hydrophobic tryptophan side chain. Effects on the activity of the mutant proteins were studied in selective growth experiments and various ATP synthase specific activity assays. The results identify potential intersubunit contacts and structurally non-distorted, accessible residues in the c-oligomer and add constraints to the arrangement of monomers in the oligomeric complex. Results from our mutagenesis experiments were interpreted in structural models of the c-oligomer that have been obtained by restrained MD simulations. Different stoichiometries and monomer orientations were applied in these calculations. A cylindrical complex consisting of 10 monomers that are arranged in two concentric rings with the N-terminal helices of the monomers located at the periphery shows the best match with the experimental data.

Adenosine Triphosphate↗

PubMed2005

Convergence of replica exchange molecular dynamics.

Replica exchange molecular dynamics (REMD) method is one of the generalized-ensemble algorithms which performs random walk in energy space and helps a system to escape from local energy traps. In this work, we studied the accuracy and efficiency of REMD by examining its ability to reproduce the results of multiple extended conventional molecular dynamics (MD) simulations and to enhance conformational sampling. Two sets of REMD simulations with different initial configurations, one from the fully extended and the other from fully helical conformations, were conducted on a fast-folding 21-amino-acid peptide with a continuum solvent model. Remarkably, the two REMD simulation sets started to converge even within 1.0 ns, despite their dramatically different starting conformations. In contrast, the conventional MD within the same time and with identical starting conformations did not show obvious signs of convergence. Excellent convergence between the REMD sets for T>300 K was observed after 14.0 ns REMD simulations as measured by the average helicity and free-energy profiles. We also conducted a set of 45 MD simulations at nine different temperatures with each trajectory simulated to 100.0 and 200.0 ns. An excellent agreement between the REMD and the extended MD simulation results was observed for T>300 K, showing that REMD can accurately reproduce long-time MD results with high efficiency. The autocorrelation times of the calculated helicity demonstrate that REMD can significantly enhance the sampling efficiency by 14.3+/-6.4, 35.1+/-0.2, and 71.5+/-20.4 times at, respectively, approximately 360, approximately 300, and approximately 275 K in comparison to the regular MD. Convergence was less satisfactory at low temperatures (T<300 K) and a slow oscillatory behavior suggests that longer simulation time was needed to reach equilibrium. Other technical issues, including choice of exchange frequency, were also examined.

Algorithms↗

PubMed1998

Conformational analysis and automated receptor docking of selective arylacetamide-based kappa-opioid agonists.

The three-dimensional structure, dynamics, and binding modes of representative kappa-opioid agonists of the arylacetamide class (U50, 488; U69,593; U62,066; CI-977; ICI199,441; ICI197,067; BRL52,537; and BRL52,656) have been investigated using molecular modeling techniques. Systematic exploration of the conformational space of the ligand combined with molecular dynamics (MD) simulations in water revealed consistent conformational preferences for all the kappa-agonists in this series. The results were further compared with available X-ray and 1D- and 2D-NMR data to identify potential "lead" conformers for molecular docking. Ligand binding modes were initially determined using automated docking of two of the ligands (U50,488 and BRL52,537) to the kappa-opioid receptor. Extrapolation of the predicted binding mode to other members in this ligand series revealed similar docking preferences, with each ligand docked along the receptor helical axis. The binding modes were further refined using MD simulations of the receptor-ligand complexes. The results show a that salt bridge is formed between the amino proton of the ligands and the carboxylate group of Asp138 in TM3. This interaction most likely serves as a key anchoring point for the agonist association. Additional ligand contacts were noted with kappa-specific residues Ile294, Leu295, and Ala298, which may, in part, explain the kappa-selectivity in this series. In comparing the arylacetamides with opiate-based ligands, no evidence was found to link these classes through a common binding motif (except for the ion pair). The binding site model was also applied to explain the enantiomeric preference of U50,488 and to provide insight to the mu/kappa-selectivity of representative ligands in this series. Overall, the results provide a structure-based rationale for ligand recognition that is consistent both with site-directed mutagenesis experiments and structure-function relationship data.

Acetamides↗

PubMed2002

Substrate entropy in enzyme enantioselectivity: an experimental and molecular modeling study of a lipase.

The temperature dependence of the enantioselectivity of Candida antarctica lipase B for 3-hexanol, 2-butanol, 3-methyl-2-butanol, 3,3-dimethyl-2-butanol, and 1-bromo-2-butanol revealed that the differential activation entropy, deltaR-SdeltaS, was as significant as the differential activation enthalpy, DeltaR-SdeltaH, to the enantiomeric ratio, E. 1-Bromo-2-butanol, with isosteric substituents, displayed the largest deltaR-SdeltaS. 3-Hexanol displayed, contrary to other sec-alcohols, a positive deltaR-SdeltaS. In other words, for 3-hexanol the preferred R-enantiomer is not only favored by enthalpy but also by entropy. Molecular dynamics (MD) simulations and systematic search calculations of the substrate accessible volume within the active site revealed that the (R)-3-hexanol transition state (TS) accessed a larger volume within the active site than the (S)-3-hexanol TS. This correlates well with the higher TS entropy of (R)-3-hexanol. In addition, this enantiomer did also yield a higher number of allowed conformations, N, from the systematic search routines, than did the S-enantiomer. The substrate accessible volume was greater for the enantiomer preferred by entropy also for 2-butanol. For 3,3-dimethyl-2-butanol, however, neither MD-simulations nor systematic search calculations yielded substrate accessible volumes that correlate to TS entropy. Ambiguous results were achieved for 3-methyl-2-butanol.

Alcohols↗

PubMed2004

Activation and inhibition of cyclin-dependent kinase-2 by phosphorylation; a molecular dynamics study reveals the functional importance of the glycine-rich loop.

Nanoseconds long molecular dynamics (MD) trajectories of differently active complexes of human cyclin-dependent kinase 2 (inactive CDK2/ATP, semiactive CDK2/Cyclin A/ATP, fully active pT160-CDK2/Cyclin A/ATP, inhibited pT14-; pY15-; and pT14,pY15,pT160-CDK2/Cyclin A/ATP) were compared. The MD simulations results of CDK2 inhibition by phosphorylation at T14 and/or Y15 sites provide insight into the structural aspects of CDK2 deactivation. The inhibitory sites are localized in the glycine-rich loop (G-loop) positioned opposite the activation T-loop. Phosphorylation of T14 and both inhibitory sites T14 and Y15 together causes ATP misalignment for phosphorylation and G-loop conformational change. This conformational change leads to the opening of the CDK2 substrate binding box. The phosphorylated Y15 residue negatively affects substrate binding or its correct alignment for ATP terminal phospho-group transfer to the CDK2 substrate. The MD simulations of the CDK2 activation process provide results in agreement with previous X-ray data.

CDC2-CDC28 Kinases↗

PubMed2004

Two-dimensional packing patterns of amino acid surfactant and higher alcohols in an aqueous phase and their associated packing parameters.

An amino acid surfactant, monosodium N-stearoyl-L-glutamate (MSSG), was assumed to associate with higher alcohols (HAs) in 1:1 molar ratio in an aqueous phase, and the packing parameters of the 1/1 MSSG/HA associates, Pa, were calculated by molecular dynamics (MD). Pa is defined as At1(l(t1)+l(t2))/a(h)l(t1), where At1 stands for the cross-sectional area of the MSSG tail, l(t1) for the length of the MSSG tail, l(t2) for the length of the HA tail, and a(h) for the area of the MSSG head group. The Pa value increased from 0.83 to 1.02 as the HA tail length (l(t2)) increased from C14 to C22. Associates of HAs having longer tails with MSSG are likely to fit better into bilayers because Pa is closer to 1 when l(t2) is longer. In the graphical results of the MD simulation for the association, however, a steric hindrance was found between the head groups of MSSG/HA when l(t2) was > or = C19. Based on this result, HA was classified into short-chain HA (SCHA, l(t2) < C19) and long-chain HA (LCHA, l(t2) > or = C19), and several possible packing units, composed of compositional combinations of MSSG, SCHA, and LCHA, were proposed. The packing unit is a building block which could constitute bilayers, and it is composed of variable compositional combinations of MSSG/HA. Assuming that SCHA associates with MSSG, the packing parameter, Punit, was calculated in a water box by MD for each packing unit. Punit is the packing parameter of a packing unit and it is defined as Vu/l(u)a(u), where Vu is the tail volume of the packing unit, l(u) is the chain length, and a(u) is the head area. For the calculation, stearyl alcohol (C18-OH, SA) was chosen as a SCHA and behenyl alcohol (C22-OH, BA) as a LCHA. When the compositional ratio MSSG:SCHA:LCHA was 1:1:1, Punit was around 1. The packing unit having Punit of around 1 formed a colloidally stable suspension for 30 days and its aggregate was a lamellar structure. However, the other packing units, for which Punit deviates significantly from 1, precipitated out in their suspensions and showed no evidence of a lamellar structure. According to the graphical MD simulations for the compositional MSSG/SCHA/LCHA associations in bilayers, vertical steric hindrance was found between LCHAs when Punit deviated significantly from 1. The steric hindrance would prevent the packing units from forming a stable bilayer and induce precipitation in the suspensions. Therefore, a proper combinational ratio of MSSG:LCHA:SCHA would play a major role in forming a lamellar structure.

Alcohols↗

PubMed2025

Oncogenic EME1 promotes tumor progression and immune modulation in human cancers with therapeutic targeting potential.

BACKGROUND: EME1, a critical DNA repair endonuclease, has emerged as a potential oncogene implicated in genome instability and cancer progression. However, its pan-cancer roles, prognostic significance, immune interactions, and therapeutic targeting remain underexplored. METHODS: We conducted a comprehensive pan-cancer analysis integrating multi-omics data from public databases, including TIMER2.0, GEPIA2, TISIDB, and cBioPortal, to evaluate EME1 expression, genetic alterations, and their association with clinical outcomes, immune infiltration, and molecular pathways. Virtual screening of 3180 FDA-approved drugs and molecular dynamics (MD) simulations were employed to identify and validate potential EME1 inhibitors. RESULTS: EME1 was significantly overexpressed in various human cancers and positively associated with advanced tumor grade and stage. High EME1 expression and mutations were linked to poor overall and disease-free survival. Immunogenomic profiling revealed strong positive correlations between EME1 and myeloid-derived suppressor cells (MDSCs), alongside a negative association with endothelial cell function, suggesting immunosuppressive roles. Machine learning models based on EME1-associated genes demonstrated high predictive accuracy for liver hepatocellular carcinoma (AUC > 0.90). Virtual screening identified eight promising drug candidates, including Everolimus and Dioscin, with strong binding affinities. MD simulations confirmed the stability of these interactions, particularly for Dioscin. CONCLUSION: This study reveals the multifaceted oncogenic roles of EME1 in tumor progression, immune evasion, and prognosis. It proposes EME1 as a promising biomarker and therapeutic target across multiple cancer types. The identified drug candidates warrant further in vitro and in vivo validation for potential repurposing in EME1-targeted cancer therapy.

EME1↗

PubMed1991

Conformation and aggregation of M13 coat protein studied by molecular dynamics.

Molecular dynamics (MD) simulations are performed on M13 coat protein, a small membrane protein for which both alpha- and beta-structures have been suggested. The simulations are started from initial conformations that are either monomers or dimers of alpha-helices or U-shaped beta-sheets. The lipid bilayer is represented by a hydrophobic potential. The results are analyzed in terms of stability, energy and secondary structure. The U-shaped beta-structure changes from a planar to a twisted form with larger twist for the monomer than the dimer. The beta-sheet is much more flexible than the alpha-helix as monitored by the root mean square (rms) fluctuations of the C alpha atoms. A comparison of the energies after 100 ps MD simulation shows that of the monomers, the alpha-helix has the lowest energy. The energy difference between alpha- and beta-structures decreases from 266 kJ/mol to 148 kJ/mol, when going from monomers to dimers. It is expected that this difference will decrease with higher aggregation numbers.

Amino Acid Sequence↗

PubMed2006

Modeling the structure of the StART domains of MLN64 and StAR proteins in complex with cholesterol.

Steroidogenic acute regulatory protein-related lipid transfer (StART) domains are ubiquitously involved in intracellular lipid transport and metabolism and other cell-signaling events. In this work, we use a flexible docking algorithm, comparative modeling, and molecular dynamics (MD) simulations to generate plausible three-dimensional atomic models of the StART domains of human metastatic lymph node 64 (MLN64) and steroidogenic acute regulatory protein (StAR) proteins in complex with cholesterol. Our results show that cholesterol can adopt a similar conformation in the binding cavity in both cases and that the main contribution to the protein-ligand interaction energy derives from hydrophobic contacts. However, hydrogen-bonding and water-mediated interactions appear to be important in the fine-tuning of the binding affinity and the position of the ligand. To gain insights into the mechanism of binding, we carried out steered MD simulations in which cholesterol was gradually extracted from within the StAR model. These simulations indicate that a transient opening of loop Omega1 may be sufficient for uptake and release, and they also reveal a pathway of intermediate states involving residues known to be crucial for StAR activity. Based on these observations, we suggest specific mutagenesis targets for binding studies of cholesterol and its derivatives that could improve our understanding of the structural determinants for ligand binding by sterol carrier proteins.

Algorithms↗

PubMed2003

An all-atom solution-equilibrated model for human extrinsic blood coagulation complex (sTF-VIIa-Xa): a protein-protein docking and molecular dynamics refinement study.

Tissue factor (TF)-bound factor (F)VIIa plays a critical role in activating FX, an event that rapidly results in blood coagulation. Despite recent advances in the structural information about soluble TF (sTF)-bound VIIa and Xa individually, the atomic details of the ternary complex are not known. As part of our long-term goal to provide a structural understanding of the extrinsic blood coagulation pathway, we built an all atom solution-equilibrated model of the human sTF-VIIa-Xa ternary complex using protein-protein docking and molecular dynamics (MD) simulations. The starting structural coordinates of sTF-VIIa and Xa were derived from dynamically equilibrated solution structures. Due to the flexible nature of the light-chain of the Xa molecule, a three-stage docking approach was employed in which SP (Arg195-Lys448)/EGF2 (Arg86-Arg139), EGF1 (Asp46-Thr85) and GLA (Ala1-Lys45) domains were docked in a sequential manner. The rigid-body docking approach of the FTDOCK method in conjunction with filtering based on biochemical knowledge from experimental site-specific mutagenesis studies provided the strategy. The best complex obtained from the docking experiments was further refined using MD simulations for 3 ns in explicit water. In addition to explaining most of the known experimental site-specific mutagenesis data pertaining to sTF-VIIa, our model also characterizes likely enzyme-binding exosites on FVIIa and Xa that may be involved in the ternary complex formation. According to the equilibrated model, the 140s loop of VIIa serves as the key recognition motif for complex formation. Stable interactions occur between the FVIIa 140s loop and the FXa -strand B2 region near the sodium-binding domain, the 160 s loop and the N-terminal activation loop regions. The helical-hydrophobic stack region that connects the GLA and EGF1 domains of VIIa and Xa appears to play a potential role in the membrane binding region of the ternary complex. The proposed model may serve as a reasonable structural basis for understanding the exosite-mediated substrate recognition of sTF-VIIa and to advance understanding of the TFPI-mediated regulatory pathway of the extrinsic blood coagulation cascade.

Binding Sites↗

PubMed2003

Structural evaluation of the agonistic action of a vitamin D analog with two side chains binding to the nuclear vitamin D receptor.

The vitamin D receptor (VDR) is one of the endocrine members of the nuclear receptor superfamily and has a characteristic high affinity for its natural ligand 1alpha,25-dihydroxyvitamin D3 [1alpha,25(OH)2D3]. From a mechanistic point of view, the most interesting analog of 1alpha,25(OH)2D3 is the one that carries two side chains, referred to as Gemini. In this study, molecular dynamics (MD) simulations of the Gemini-VDR complex were performed that demonstrated that the binding of a ligand with a 25% increased volume does not disturb the overall structure of the ligand-binding domain (LBD). It was found that one of the two side chains takes exactly the same position as the single side chain of the natural ligand, which suggests that the molecular mechanism of the agonism of Gemini is identical to that of 1alpha,25(OH)2D3. VDR single and double point mutants at L227, A303, I313, and L397 and in vitro and ex vivo assessment of their agonistic action confirmed the predictions of the MD simulations. Moreover, it was found that the second side chain of Gemini can choose between two binding positions within the ligand-binding pocket of the VDR. These two newly identified "corners" were characterized most specifically by the amino acids pairs L227/A303 and I313/L397. Therefore, Gemini is an important model compound that allows further insight into the molecular actions of the VDR but is, in parallel, also a promising precursor for the design of even more potent 1alpha,25(OH)2D3 analogs.

Animals↗

PubMed2001

Structural features and molecular assembly of amorphous phosphazenic materials in the bulk--combined theoretical and experimental techniques: Tris-(2,2'-dioxy-1,1'-binaphthyl)cyclotriphosphazen.

The structure and the assembly of tris-(2,2'-dioxy-binaphthyl)cyclotriphosphazene [(+)-[NP3(O2C20H12)3], DBNP, in the solid amorphous state was studied using molecular dynamics (MD) including ad hoc quantum mechanically derived force field (FF) parameters, in combination with the energy dispersive X-ray diffraction (EDXD) technique. The atom-atom radial distribution function (RDF) curve obtained through the EDXD experiment revealed low intensity peaks not attributable to the intramolecular distances of the single molecule, but clearly featuring a low energy state of long-distance three-dimensional assembly. The radial distribution functions (RDF) were calculated for various models of DBNP submitted to theoretical MD simulations. Based on the comparison of theoretically calculated RDFs and those obtained from the EDXD experiment, the predominant structural motif of the material in the bulk was found to have DBNP molecules laid one upon the other to form tubular nanostructures. These contain eight DBNP units each (length ca. 46 A) with two and three of these units aligned in parallel and held together. The material can be represented as a bulk of tubular snake-like chains undergoing distortions with a step of eight DBNP units. The bending angles, that vary randomly, attain limited values sufficient to induce disorder and thus nonperiodic structure. The present application of MD simulations combined with EDXD data appear to be a general approach to solve for the first time otherwise intractable issues concerning structural features and assembly of molecular materials in the bulk.

Journal Article↗

PubMed2026

Elucidating the Mechanism of Xiaoqinglong Decoction in Chronic Urticaria Treatment: An Integrated Approach of Network Pharmacology, Bioinformatics Analysis, Molecular Docking, and Molecular Dynamics Simulations.

INTRODUCTION: Xiaoqinglong Decoction (XQLD) is a traditional Chinese medicinal formula commonly used to treat chronic urticaria (CU). However, its underlying therapeutic mechanisms remain incompletely characterized. This study employed an integrated approach combining network pharmacology, bioinformatics, molecular docking, and molecular dynamics simulations to identify the active components, potential targets, and related signaling pathways involved in XQLD's therapeutic action against CU, thereby providing a mechanistic foundation for its clinical application. METHODS: The active components of XQLD and their corresponding targets were identified using the Traditional Chinese Medicine Systems Pharmacology (TCMSP) database. CU-related targets were retrieved from the OMIM and GeneCards databases. Subsequently, core components and targets were determined via protein-protein interaction (PPI) network analysis and component-target-pathway network construction. Topological analyses were performed using Cytoscape software to prioritize core nodes within these networks. Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment analyses were conducted via the DAVID database to identify enriched biological processes and signaling pathways. Molecular docking was performed to evaluate binding interactions between key components and core targets, while molecular dynamics (MD) simulations were employed to assess the stability of the component-target complexes with the lowest binding energy. Finally, CU-related targets of XQLD were validated using datasets from the Gene Expression Omnibus (GEO) database. RESULTS: A total of 135 active components and 249 potential targets of XQLD were identified, alongside 1,711 CU-related targets. Core components, such as quercetin, kaempferol, beta-sitosterol, naringenin, stigmasterol, and luteolin, exhibited high degree values in the constructed networks. The core targets identified included AKT1, TNF, IL6, TP53, PTGS2, CASP3, BCL2, ESR1, PPARG, and MAPK3. GO and KEGG pathway enrichment analyses revealed the PI3K-Akt signaling pathway as a central regulatory mechanism. Molecular docking studies demonstrated strong binding affinities between active components and core targets, with the stigmasterol-AKT1 complex exhibiting the lowest binding energy (-11.4 kcal/mol) and high stability in MD simulations. Validation using GEO datasets identified 12 core genes shared between CU-related targets and XQLD-associated targets, including PTGS2 and IL6, which were also prioritized as core targets in the network pharmacology analyses. DISCUSSION: This study comprehensively integrates multidisciplinary approaches to clarify the potential molecular mechanisms of XQLD in treating CU, highlighting its multitarget and multipathway synergistic effects. Molecular docking and dynamics simulations confirm the stable interaction between stigmasterol and the core target AKT1. Additionally, GEO dataset analysis verifies the pathogenic relevance of targets such as PTGS2 and IL6, significantly enhancing the credibility of our findings. These results provide a modern scientific basis for the traditional therapeutic effects of XQLD on CU and have important implications for developing multitarget treatments for this condition. However, this study mainly relies on database mining and computational simulations. Further in vitro and in vivo experimental validations are needed to confirm the predicted component-target-pathway interactions. CONCLUSION: This study identifies the active components, potential targets, and pathways through which XQLD exerts therapeutic effects on CU. These findings provide a theoretical foundation for further mechanistic studies and support their clinical application in the treatment of CU.

Molecular Docking Simulation↗

PubMed2003

Modeling of flap endonuclease interactions with DNA substrate.

Structure-specific 5' nucleases play an important role in DNA replication and repair uniquely recognizing an overlap flap DNA substrate and processing it into a DNA nick. However, in the absence of a high-resolution structure of the enzyme/DNA complex, the mechanism underlying this recognition and substrate specificity, which is key to the enzyme's function, remains unclear. Here, we propose a three-dimensional model of the structure-specific 5' flap endonuclease from Pyrococcus furiosus in its complex with DNA. The model is based on the known X-ray structure of the enzyme and a variety of biochemical and molecular dynamics (MD) data utilized in the form of distance restraints between the enzyme and the DNA. Contacts between the 5' flap endonuclease and the sugar-phosphate backbone of the overlap flap substrate were identified using enzyme activity assays on substrates with methylphosphonate or 2'-O-methyl substitutions. The enzyme footprint extends two to four base-pairs upstream and eight to nine base-pairs downstream of the cleavage site, thus covering 10-13 base-pairs of duplex DNA. The footprint data are consistent with a model in which the substrate is bound in the DNA-binding groove such that the downstream duplex interacts with the helix-hairpin-helix motif of the enzyme. MD simulations to identify the substrate orientation in this model are consistent with the results of the enzyme activity assays on the methylphosphonate and 2'-O-methyl-modified substrates. To further refine the model, 5' flap endonuclease variants with alanine point substitutions at amino acid residues expected to contact phosphates in the substrate and one deletion mutant were tested in enzyme activity assays on the methylphosphonate-modified substrates. Changes in the enzyme footprint observed for two point mutants, R64A and R94A, and for the deletion mutant in the enzyme's beta(A)/beta(B) region, were interpreted as being the result of specific interactions in the enzyme/DNA complex and were used as distance restraints in MD simulations. The final structure suggests that the substrate's 5' flap interacts with the enzyme's helical arch and that the helix-hairpin-helix motif interacts with the template strand in the downstream duplex eight base-pairs from the cleavage site. This model suggests specific interactions between the 3' end of the upstream oligonucleotide and the enzyme. The proposed structure presents the first detailed description of substrate recognition by structure-specific 5' nucleases.

Base Sequence↗

PubMed2006

Gemini surfactants at the air/water interface: a fully atomistic molecular dynamics study.

Gemini surfactants typically consist of two single-chain surfactants chemically linked by a spacer molecule. We report herein the results of fully atomistic molecular dynamics (MD) simulations of a series of Gemini surfactants: CsH2s-alpha,omega-bis(C12H25N+(CH3)2Cl-), at the air/water interface with s = 3, 4, 6, 12, 14, and 16, at values of the initial surface area per surfactant AS = 70 A2, 77 A2, 95 A2, 151 A2, 133 A2, and 103 A2, respectively. The AS values employed were obtained from surface tension and neutron reflection experiments at the respective cmc of each surfactant. The Gemini surfactant corresponding to s = 3 was also simulated at AS = 105 A2, which is the experimentally derived value of surface area per surfactant at 1/10th of cmc. Only the surfactants with s = 12 and 14 and the surfactant with s = 3 at AS = 105 A2 gave a stable monolayer at the air/water interface. In other cases, we observe movement of some surfactant molecules from the air/water interface into the aqueous phase, resulting in a stable primary monolayer of surfactants at the air/water interface and a small concentration of surfactant molecules below it. The latter form aggregates, with their hydrophobic chains in the core. The density profiles along the normal to the interface are compared with the ones obtained from neutron reflection experiments. The MD simulations confirm the bending of the spacer toward the hydrophobic chains as the spacer length is increased and the spacer becomes more hydrophobic. The simulations have helped to shed light on the low-resolution picture which emerges from experimental analyses.

Air↗

PubMed2005

M.TaqI facilitates the base flipping via an unusual DNA backbone conformation.

MD simulations have been carried out to understand the dynamical behavior of the DNA substrate of the Thermus aquaticus DNA methyltransferase (M.TaqI) in the methylation process at N6 of adenine. As starting structures, an x-ray structure of M.TaqI in complex with DNA and cofactor analogue (PDB code: 1G 38) and free decamer d(GTTCGATGTC)(2) were taken. The x-ray structure shows two consecutive BII substates that are not observed in the free decamer. These consecutive BII substates are also observed during our simulation. Additionally, their facing backbones adopt the same conformations. These double facing BII substates are stable during the last 9 ns of the trajectories and result in a stretched DNA structure. On the other hand, protein-DNA contacts on 5' and 3' phosphodiester groups of the partner thymine of flipped adenine have changed. The sugar and phosphate parts of thymine have moved further into the empty space left by the flipping base without the influence of protein. Furthermore, readily high populated BII substates at the GpA step of palindromic tetrad TCGA rather than CpG step are observed in the free decamer. On the contrary, the BI substate at the GpA step is observed on the flipped adenine strand. A restrained MD simulation, reproducing the BI/BII pattern in the complex, demonstrated the influence of the unusual backbone conformation on the dynamical behavior of the target base. This finding along with the increased nearby interstrand phosphate distance is supportive to the N6-methylation mechanism.

Adenine↗

PubMed2003

The structure and stability of Si60 and Ge60 cages: a computational study.

Structural studies of fullerene-like Si(60) and Ge(60) cages using ab initio methods were augmented by density functional tight-binding molecular dynamics (DFTB-MD) simulations of finite temperature effects. Neither the perfect I(h) symmetry nor the distorted T(h) structures are true minima. The energies of both are high relative to distorted, lower symmetry minima, C(i) and T, respectively, which still preserve C(60)-type connectivity. Both Si(60) and Ge(60) favor C(i) symmetry cages in which Si and Ge vertexes exhibit either near-trigonal or pyramidal geometries. These structural variations imply significant reactivity differences between different positions. The small magnetic shielding effects (NICS) indicate that aromaticity is not important in these systems. The inorganic fullerene cages have lower stabilities compared with their carbon analogs. Si(60) is stable towards spontaneous disintegration up to 700 K according to DFTB-MD simulations, and thus has potential for experimental observation. In contrast, Ge(60) preserves its cage structure only up to 200 K.

Journal Article↗