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 289 records · Page 16Linked to original sources
PubMed2002

OmpA: a pore or not a pore? Simulation and modeling studies.

The bacterial outer membrane protein OmpA is composed of an N-terminal 171-residue beta-barrel domain (OmpA(171)) that spans the bilayer and a periplasmic, C-terminal domain of unknown structure. OmpA has been suggested to primarily serve a structural role, as no continuous pore through the center of the barrel can be discerned in the crystal structure of OmpA(171). However, several groups have recorded ionic conductances for bilayer-reconstituted OmpA(171). To resolve this apparent paradox we have used molecular dynamics (MD) simulations on OmpA(171) to explore the conformational dynamics of the protein, in particular the possibility of transient formation of a central pore. A total of 19 ns of MD simulations of OmpA(171) have been run, and the results were analyzed in terms of 1) comparative behavior of OmpA(171) in different bilayer and bilayer-mimetic environments, 2) solvation states of OmpA(171), and 3) pore characteristics in different MD simulations. Significant mobility was observed for residues and water molecules within the beta-barrel. A simulation in which putative gate region side chains of the barrel interior were held in a non-native conformation led to an open pore, with a predicted conductance similar to experimental measurements. The OmpA(171) pore has been shown to be somewhat more dynamic than suggested by the crystal structure. A gating mechanism is proposed to explain its documented channel properties, involving a flickering isomerization of Arg138, forming alternate salt bridges with Glu52 (closed state) and Glu128 (open state).

Arginine↗
PubMed2005

Molecular dynamics study of water penetration in staphylococcal nuclease.

The ionization properties of Lys and Glu residues buried in the hydrophobic core of staphylococcal nuclease (SN) suggest that the interior of this protein behaves as a highly polarizable medium with an apparent dielectric constant near 10. This has been rationalized previously in terms of localized conformational relaxation concomitant with the ionization of the internal residue, and with contributions by internal water molecules. Paradoxically, the crystal structure of the SN V66E variant shows internal water molecules and the structure of the V66K variant does not. To assess the structural and dynamical character of interior water molecules in SN, a series of 10-ns-long molecular dynamics (MD) simulations was performed with wild-type SN, and with the V66E and V66K variants with Glu66 and Lys66 in the neutral form. Internal water molecules were identified based on their coordination state and characterized in terms of their residence times, average location, dipole moment fluctuations, hydrogen bonding interactions, and interaction energies. The locations of the water molecules that have residence times of several nanoseconds and display small mean-square displacements agree well with the locations of crystallographically observed water molecules. Additional, relatively disordered water molecules that are not observed crystallographically were found in internal hydrophobic locations. All of the interior water molecules that were analyzed in detail displayed a distribution of interaction energies with higher mean value and narrower width than a bulk water molecule. This underscores the importance of protein dynamics for hydration of the protein interior. Further analysis of the MD trajectories revealed that the fluctuations in the protein structure (especially the loop elements) can strongly influence protein hydration by changing the patterns or strengths of hydrogen bonding interactions between water molecules and the protein. To investigate the dynamical response of the protein to burial of charged groups in the protein interior, MD simulations were performed with Glu66 and Lys66 in the charged state. Overall, the MD simulations suggest that a conformational change rather than internal water molecules is the dominant determinant of the high apparent polarizability of the protein interior.

Algorithms↗
PubMed2005

Stability and fluctuations of amide hydrogen bonds in a bacterial cytochrome c: a molecular dynamics study.

Molecular dynamics (MD) simulations on a bacterial cytochrome c were performed to investigate the lifetime and fluctuations of backbone hydrogen bonds and to correlate these data with protection factors for hydrogen exchange measured by NMR spectroscopy (Bartalesi et al. in Biochemistry, 42:10923-10930, 2003). The MD simulations provide a consistent pattern in that long lifetimes of hydrogen bonds go along with small amplitude fluctuations. In agreement with experiments, differences in stability were found with a rather flexible N-terminal segment as compared with a more rigid C-terminal part. Protection factors of backbone hydrogen exchange correlate strongly with the number of contacts but also with hydrogen-bond occupancy, hydrogen-bond survival times, as well as the inverse of fluctuations of backbone atoms and hydrogen-bond lengths derived from MD simulation data. We observed a conformational transition in the C-terminal loop, and significant motion in the N-terminal loop, which can be interpreted as being the structural units involved in the onset of the protein unfolding process in agreement with experimental evidence on mitochondrial cytochrome c.

Amides↗
PubMed2001

Conformational dynamics of an alanine dipeptide analog: an ab initio molecular dynamics study.

An ab initio molecular dynamics (MD) simulation technique employing the Born-Oppenheimer approach in the framework of a Gaussian implementation of Kohn-Sham density functional theory is used to study the gas-phase conformational dynamics of an alanine dipeptide analog. It is found that conformational transformation between C5 and C7(eq) occurs on the picosecond time scale. Classical MD simulations using most of the popular force fields do not yield a transition even after nanoseconds. An analysis is given of the difference, for this small gas-phase system, between ab initio MD and traditional MD simulation using force fields.

Alanine↗
PubMed1997

Prediction of titration properties of structures of a protein derived from molecular dynamics trajectories.

This paper explores the dependence of the molecular dynamics (MD) trajectory of a protein molecule on the titration state assigned to the molecule. Four 100-ps MD trajectories of bovine pancreatic trypsin inhibitor (BPTI) were generated, starting from two different structures, each of which was held in two different charge states. The two starting structures were the X-ray crystal structure and one of the solution structures determined by NMR, and the charge states differed only in the ionization state of N terminus. Although it is evident that the MD simulations were too short to sample fully the equilibrium distribution of structures in each case, standard Poisson-Boltzmann titration state analysis of the resulting configurations shows general agreement between the overall titration behavior of the protein and the charge state assumed during MD simulation: at pH 7, the total net charge of the protein resulting from the titration analysis is consistently lower for the protein with the N terminus assumed to be neutral than for the protein with the N terminus assumed to be charged. For most of the ionizable residues, the differences in the calculated pKaS among the four trajectories are statistically negligible and remain in good agreement with the data obtained by crystal structure titration and by experiment. The exceptions include the N terminus, which responds directly to the change of its imposed charge; the C terminus, which in the NMR structure interacts strongly with the former; and a few other residues (Arg 1, Glu 7, Tyr 35, and Arg 42) whose pKaS reflect the initial structure and the limited trajectory lengths. This study illustrates the importance of the careful assignment of protonation states at the start of MD simulations and points to the need for simulation methods that allow for the variation of the protonation state in the calculation of equilibrium properties.

Aprotinin↗
PubMed2006

Combining docking and molecular dynamic simulations in drug design.

A rational approach is needed to maximize the chances of finding new drugs, and to exploit the opportunities of potential new drug targets emerging from genomic and proteomic initiatives, and from the large libraries of small compounds now readily available through combinatorial chemistry. Despite a shaky early history, computer-aided drug design techniques can now be effective in reducing costs and speeding up drug discovery. This happy outcome results from development of more accurate and reliable algorithms, use of more thoughtfully planned strategies to apply them, and greatly increased computer power to allow studies with the necessary reliability to be performed. Our review focuses on applications and protocols, with the main emphasis on critical analysis of recent studies where docking calculations and molecular dynamics (MD) simulations were combined to dock small molecules into protein receptors. We highlight successes to demonstrate what is possible now, but also point out drawbacks and future directions. The review is structured to lead the reader from the simpler to more compute-intensive methods. Thus, while inexpensive and fast docking algorithms can be used to scan large compound libraries and reduce their size, more accurate but expensive MD simulations can be applied when a few selected ligand candidates remain. MD simulations can be used: during the preparation of the protein receptor before docking, to optimize its structure and account for protein flexibility; for the refinement of docked complexes, to include solvent effects and account for induced fit; to calculate binding free energies, to provide an accurate ranking of the potential ligands; and in the latest developments, during the docking process itself to find the binding site and correctly dock the ligand a priori.

Computer Simulation↗
PubMed2005

Molecular dynamics in the isothermal-isobaric ensemble: the requirement of a "shell" molecule. II. Simulation results.

The results of a series of constant pressure and temperature molecular-dynamics (MD) simulation studies based on the rigorous shell particle formulation of the isothermal-isobaric (NpT) ensemble are presented. These MD simulations validate the newly proposed constant pressure equations of motion in which a "shell" particle is used to define uniquely the volume of the system [M. J. Uline and D. S. Corti, J. Chem. Phys. (to be published), preceding paper]. Ensemble averages obtained with the new MD NpT algorithm match the ensemble averages obtained using the previously derived shell particle Monte Carlo NpT method [D. S. Corti, Mol. Phys. 100, 1887 (2002)]. In addition, we also verify that the Hoover NpT MD algorithm [W. G. Hoover, Phys. Rev. A 31, 1695 (1985); 34, 2499 (1986)] generates the correct ensemble averages, though only when periodic boundary conditions are employed. The extension of the shell particle MD algorithm to multicomponent systems is also discussed, in which we show for equilibrium properties that the identity of the shell particle is completely arbitrary when periodic boundary conditions are applied. Self-diffusion coefficients determined with the shell particle equations of motion are also identical to those obtained in other ensembles. Finally, since the mass of the shell particle is known, the system itself, and not a piston of arbitrary mass, controls the time scales for internal pressure and volume fluctuations. We therefore consider the effects of the shell particle on the dynamics of the system. Overall, the shell particle MD algorithm is an effective simulation method for studying systems exposed to a constant external pressure and may provide an advantage over other existing constant pressure approaches when developing nonequilibrium MD methods.

Journal Article↗
PubMed2000

Structure of the Sm binding site from human U4 snRNA derived from a 3 ns PME molecular dynamics simulation.

A molecular dynamics simulation of the Sm binding site from human U4 snRNA was undertaken to determine the conformational flexibility of this region and to identify RNA conformations that were important for binding of the Sm proteins. The RNA was fully-solvated (>9,000 water molecules) and charge neutralized by inclusion of potassium ions. A three nanosecond MD simulation was conducted using AMBER with long-range electrostatic forces considered using the particle mesh Ewald summation method. The initial model of the Sm binding site region had the central and 3' stem-loops that flanked the Sm site co-axial with one another, and with the single-stranded Sm binding site region ([I] conformation). During the course of the trajectory, the axes of the 3' stem-loop, and later the central stem-loop, became roughly orthogonal from their original anti-parallel orientation. As these conformational changes occurred, the snRNA adopted first an [L] conformation, and finally a [U] conformation. The [U] conformation was more stable than either the [I] or [L] conformations, and persisted for the final 1 ns of the trajectory. Analysis of the structure resulting from the MD simulations revealed the bulged nucleotide, U114, and the mismatched Ag91-G110 base pair provided distinctive structural features that may enhance Sm protein binding. Based on the results of the MD simulation and the available experimental data, we proposed a mechanism for the binding of the Sm protein sub-complexes to the snRNA. In this model, the D1/D2 and E/F/G Sm protein sub-complexes first bind the snRNA in the [U] conformation, followed by conformational re-arrangement to the [I] conformation and binding of the D3/B Sm protein sub-complex.

Base Sequence↗
PubMed1995

A simple protocol for identification of helical and mobile residues in membrane proteins.

A simple molecular dynamics (MD) simulations is protocol shown to predict whether a residue is in a structured alpha-helical segment or in a mobile loop or terminal segment of a membrane protein. The results are verified by comparisons with experimental NMR data. The protocol consists of performing several independent MD simulations on a polypeptide sequence of interest in a dielectric continuum with a relative permittivity epsilon = 2. The time histories of the individual angles between NH bond vectors at time 0 and time t later are calculated, and then Gaussian smoothing of typically 50 ps is applied. The smoothed data are subtracted from the original data to yield the short time fluctuations of the NH bond vectors, and then the rms deviations of the angles are calculated and compared to experimental NMR results. Pf1 coat protein and the c subunit of the E. coli F1F0 ATP synthase are used as examples of membrane proteins. The calculated NH bond rms fluctuations are in qualitative agreement with experimental NMR data in showing that each of these proteins has a mobile segment connecting two helices, as well as mobile N and C-terminal regions. This MD simulations protocol can demonstrate the presence of both the amphipathic and hydrophobic helices while hydropathy plots are able to detect only the hydrophobic helices present in membrane proteins.

Capsid↗
PubMed2004

A computational study of the binding of propidium to the peripheral anionic site of human acetylcholinesterase.

Combined docking and molecular dynamics (MD) simulations were carried out in order to investigate the binding mode of propidium at the human acetylcholinesterase (HuAChE) peripheral site. Two different docking protocols followed by cluster analyses were performed, allowing the identification of five high-populated and low-energy configuration families. To dynamically explore the behavior of the ligand at the peripheral HuAChE binding site, six complexes (five low-energy and one high-energy) were submitted to 2.5 ns of MD simulations. The representative propidium/HuAChE binding modes were chosen on the basis of both the docking energy score and the dynamic stability of the complexes throughout the MD simulations. The most stable poses of propidium at HuAChE PAS were similar to those experimentally determined with the murine enzyme. We therefore suggest that the present modeling protocol might be used in the dynamic investigation of the interactions of a small-molecule inhibitor with a surface-like binding site of a protein. Finally, because of the biological relevance of the target studied here, the present results can be of interest for the rational design of molecules potentially useful in the treatment of the Alzheimer's disease.

Acetylcholinesterase↗
PubMed1993

Molecular dynamics simulation of 1,2-dilauroyl-L-phosphatidylethanolamine binding to phospholipase A2: an attempt to explain the selective hydrolysis of substrate fatty acid ester at position 2.

To improve our understanding of why phospholipase A2 (PLA2) specifically catalyzes the hydrolysis of the fatty acid ester bond at position 2, not at position 1, of 1,2-diacyl-3-sn-phosphoglycerides, the binding of each fatty acid chain of 1,2-dilauroyl-L-phosphatidyl-ethanolamine (DLPE), a natural substrate, to bovine pancreas PLA2 was examined by molecular dynamics (MD) simulations. Two different binding modes were considered, i.e., the respective hydrocarbon chains of 1- and 2-lauroyl fatty acid esters were located at the PLA2 binding sites usually observed in the complex crystals (Form A2) and at the reverse sites (Form A1). Although the total energies of both forms fluctuated within nearly the same range during the 80 ps MD simulations, the binding mode of DLPE to the PLA2 catalytic site through the coordination to Ca2+ was much more advantageous in Form A2 than that in Form A1; significant deviation of the Ca2+ position from its starting structure was observed in the MD simulation of Form A1. The result suggests the importance of Ca2+ in the selective recognition and catalytic function of PLA2 toward the 2-positioned fatty acid ester of phosphoglyceride substrates.

Animals↗
PubMed2004

Molecular dynamics simulation of a decasaccharide fragment of heparin in aqueous solution.

Molecular dynamics (MD) simulations on heparin-water-sodium systems were carried out in order to establish a simulation protocol able to represent heparin solution conformation under physiological conditions. Atomic charges suitable for heparin oligosaccharides were obtained from ab initio quantum-mechanical computations, at the 6-31G(**) level. The GROMACS forcefield, the SPC, and SPC/E water models were employed. Also heparin was simulated with IdoA residues in 1C(4) or 2S(0) conformational states. The results of the performed MD simulations are in agreement with the available experimental data, suggesting that this approach can be applied for the study of heparin interactions with its target proteins and thus play a role in the development of new antithrombotic agents.

Carbohydrate Conformation↗
PubMed2001

Computer aided study of ligand binding with catalytic domain of Avian sarcoma virus integrase and its ligand binding loops.

We have carried out 1 nanosecond (ns) Molecular Dynamics (MD) simulations of the drug Y3 (4-acetylamino-5-hydroxynaphthalene- 2, 7-disulfonic acid) complexed with catalytic domain of Avian sarcoma virus Integrase (ASV-IN), both in vacuum and in the presence of explicit solvent. Starting models were obtained on the basis of PDB co- ordinates (1A5X) of ASV-IN-Y3 complex. Mn(2+) cation was present in the active site. To neutralize the positive charge in the presence of explicit solvent, eight Cl(-)anions were added. Energy Minimization (EM) and MD simulations, for both the systems, were carried out using Sander's module of AMBER5.0 with all atom force field. We also carried out 1 ns MD simulation on two flexible loops--L1 (Gly54-Gln62) and L2 (Trp138-Met155) playing crucial role in interaction of IN with the drug, under differing environmental conditions (vacuum, aqueous and organic solvent methanol). Comparison of the conformational changes in the loops, monomer and dimer is presented in the paper. Our results showed that the conformation of the loop region was closest to crystallographic data in case of monomer and constrained loops in aqueous environment. However, the dimer in vacuum was more stable than monomer. The beta sheet structure of the monomer in aqueous environment was unstable. Latter also took long time for equilibration. The box formed by loops L1 and L2 from two sub units (IINA and INB) of the dimer satisfies prerequisites for ligand recognition site and seems to be the functional biological unit.

Amino Acid Sequence↗
PubMed1996

Selection of peptides that bind to the HLA-A2.1 molecule by molecular modelling.

Cytotoxic T lymphocytes recognize antigenic peptides in association with major histocompatibility complex class I proteins. Although a large set of class I binding peptides has been described, it is not yet easy to search for potentially antigenic peptides without synthesis of a panel of peptides, and subsequent binding assays. In order to predict HLA-A2.1-restricted antigenic epitopes, a computer model of the HLA-A2.1 molecule was established using X-ray crystallography data. In this model nonameric peptide sequences were aligned. In a molecular dynamics (MD) simulation with two sets of peptides known to be presented by HLA-A2.1, it was important to know the anchor amino acid residue preference and the distance between the anchor residues. We show here that the peptides bound to the HLA-A2.1 model structure possess a side chain of C-terminal anchor residue oriented into the binding groove with different distances between the two anchor residues from 15 to 21A. We also synthesized a set of nonamer peptides containing amino acid sequences of Hepatitis B virus protein that were selected on the basis of previously described HLA-A2.1 specific motifs. When results obtained from the MD simulation were compared with functional binding assays using the TAP-deficient cell line T2, it was evident that the MD simulation method improves prediction of the HLA-A2.1 binding epitope sequence. These results suggest that this approach can provide a way to predict peptide epitopes and search for antigenic regions in sequences in a variety of antigens without screening a large number of synthetic peptides.

Amino Acid Sequence↗
PubMed2005

Molecular dynamics study on the interaction of a mithramycin dimer with a decanucleotide duplex.

The complex of a minor groove binding drug mithramycin (MTR) and the self-complementary d(TAGCTAGCTA) 10-mer duplex was investigated by molecular dynamics (MD) simulations using the AMBER 7.0 suite of programs. There is one disaccharide and trisaccharide segment projecting from opposite ends of an aglycone chromophore of MTR. A MTR dimer complex (MTR)2Mg2+ is formed in the presence of a coordinated ion Mg2+. A NMR solution structure of two (MTR)2Mg2+ complexes bound with one DNA duplex, namely, the 2:1 duplex complex, was taken as the starting structure for the MD simulation. The partial charge on each atom was calculated using the multiple-RESP fitting procedure, and all of the missing parameters in the Parm99 force field used were adapted comparably from the literature. The length of the MD simulation was 5 ns, and the binding free energy for the formation of a 1:1 or 2:1 duplex complex was determined from the last 4 ns of the simulation. The binding free energies were decomposed to components of the contributions from different energy types, and the changes in the helical parameters of the bound DNA duplex plus the glycosidic linkages between sugar residues of the bound MTR dimer were determined. It was found that binding of the first (MTR)2Mg2+ complex with the DNA duplex to form a 1:1 duplex complex does not cause stiffening of the duplex especially in the unoccupied site of the duplex. However, the overall flexibility of the DNA duplex is reduced substantially once the second (MTR)2Mg2+ complex is bound with the unoccupied site to form the 2:1 duplex complex. The van der Waals interactions were found to be dominant in the central part of the DNA duplex where sugar residues from each bound (MTR)2Mg2+ complex were inwardly pointing and the corresponding minor groove was widened.

Base Pairing↗
PubMed2005

Computational study of IAG-nucleoside hydrolase: determination of the preferred ground state conformation and the role of active site residues.

The mechanism of action of inosine-adenosine-guanosine nucleoside hydrolase (IAG-NH) has been investigated by long-term molecular dynamics (MD) simulation in TIP3P water using stochastic boundary conditions. Special attention has been given to the role of leaving group pocket residues and conformation of the bound substrate at the active site of IAG-NH. We also describe the positioning of the residues of an important flexible loop at the active site, which was previously unobservable by X-ray crystallography due to high B-factors. Five MD simulations have been performed with the Enzyme x Substrate complexes: Enzyme x anti-Adenosine with Asp40-COOH [E(40H) x Ade(a)], Enzyme x anti-Adenosine with Asp40-COO- [E(40) x Ade(a)], Enzyme x syn-Adenosine with Asp40-COOH [E(40H) x Ade(s)], Enzyme x syn-Adenosine with Asp40-COO- [E(40) x Ade(s)], and Enzyme x anti-Inosine with Asp40-COO- [E(40) x Ino(a)]. Overall, the structures generated from the MD simulation of E(40H) x Ade(s) preserve the catalytically important hydrogen bonds as well as electrostatic and hydrophobic interactions to provide a plausible catalytic structure. When deprotonated Asp40 (Asp4-COO-) is present, the active site is open to water solvent which interferes with the base stacking between Trp83 and nucleobase. A calculation using Poisson-Boltzmann equation module supports that Asp40 indeed has an elevated pK(app). Solvent accessible surface area (SASA) calculations on all the five MD structures shows that systems with protonated Asp40, namely, E(40H) x Ade(a) and E(40H) x Ade(s), have zero SASA. It is found that a water molecule is hydrogen-bonded to the N7 of the nucleobase and is probably the essential general acid to protonate the departing nucleobase anion. The N7-bonded water is in turn hydrogen-bonded to waters in a channel, held in place by the residues of the flexible loop, Tyr257, His247, and Cys245. Using normal-mode analysis with elastic network model, we find that the flexible loop explores a conformational space much larger than in the MD trajectory, leading to a "gating"-like motion with respect to the active site.

Adenosine↗
PubMed1998

Conformer selection and differential restriction of ligand mobility by a plant lectin--conformational behaviour of Galbeta1-3GlcNAcbeta1-R, Galbeta1-3GalNAcbeta1-R and Galbeta1-2Galbeta1-R' in the free state and complexed with galactoside-specific mistletoe lectin as revealed by random-walk and conformational-clustering molecular-mechanics.

To study conformational parameters of ligands before and after complex formation with the galactoside-binding agglutinin of Viscum album L. (VAA) in solution, combined computer-assisted random walk molecular mechanics (RAMM) calculations extended by conformational clustering analysis (CCA), molecular dynamics (MD) simulations as well as two-dimensional rotating-frame nuclear Overhauser effect (ROE) and two-dimensional nuclear Overhauser effect (NOE) spectroscopy NMR experiments were employed. Derivatives of the naturally occurring disaccharides Galbeta1-3GlcNAcbeta1-R and Galbeta1-3GalNAcbeta1-R as well as of a synthetic high-affinity binding partner, i.e. the disaccharide Galbeta1-2Galbeta1-R', were chosen as ligands in this study. The disaccharides displayed inherent flexibility in the valley of the global minimum between phi/psi combinations of (40 degrees/60 degrees) and (40 degrees/-60 degrees). Calculations of the de-N-acetylated sugars revealed that presence of this group did not markedly influence the distribution of low-energy conformers in the phi, psi, epsilon plot. Occupation of side minima at phi/psi (180 degrees/0 degrees) or (0 degrees/180 degrees) is either unlikely or low according to the results of MD simulations and RAMM calculations extended by CCA. Notably, these side minima define conformations which are not stable during a MD simulation. Transitions to other minima occur already a few picoseconds after the start of the simulation. NMR experiments of the free-state ligand confirmed the validity of the data sets obtained by the calculations. Following the description of the conformational space in the free-state NMR experiments were performed for these disaccharides complexed with VAA. They yielded two interresidual contacts for Galbeta1-3GlcNAcbeta1-R and Galbeta1-2Galbeta1-R'. The ligand conformations in the complex did not deviate markedly from those of a minimum conformation in the free state. One- and two-dimensional transferred nuclear Overhauser enhancement (TRNOE) experiments at different mixing times excluded the influence of spin-diffusion effects. When the NOE build-up curves in the three studied cases were compared, the residual mobility of the penultimate carbohydrate unit of Galbeta1-3GalNAcbeta1-R was observed to be higher than that of the respective hexopyranose unit of the other two bound ligands. Due to the availability of the conformational parameters of Galbeta1-2Galbeta1-R' in association with a galectin, namely the beta-galactoside-binding protein from chicken liver, it is remarkable to note that this ligand displays different conformations in the binding sites of either the plant or the animal lectin. They correspond to local energy-minimum conformations in the phi,psi, epsilon plot and substantiate differential conformer selection by these two lectins with identical nominal monosaccharide specificity.

Binding Sites↗
PubMed2000

Molecular dynamics simulations of the d(CCAACGTTGG)(2) decamer: influence of the crystal environment.

Molecular dynamics (MD) simulations of the DNA duplex d(CCAACGTTGG)(2) were used to study the relationship between DNA sequence and structure. Two crystal simulations were carried out; one consisted of one unit cell containing two duplexes, and the other of two unit cells containing four duplexes. Two solution simulations were also carried out, one starting from canonical B-DNA and the other starting from the crystal structure. For many helicoidal parameters, the results from the crystal and solution simulations were essentially identical. However, for other parameters, in particular, alpha, gamma, delta, (epsilon - zeta), phase, and helical twist, differences between crystal and solution simulations were apparent. Notably, during crystal simulations, values of helical twist remained comparable to those in the crystal structure, to include the sequence-dependent differences among base steps, in which values ranged from 20 degrees to 50 degrees per base step. However, in the solution simulations, not only did the average values of helical twist decrease to approximately 30 degrees per base step, but every base step was approximately 30 degrees, suggesting that the sequence-dependent information may be lost. This study reveals that MD simulations of the crystal environment complement solution simulations in validating the applicability of MD to the analysis of DNA structure.

Computer Simulation↗