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 559 records · Page 31Linked to original sources
PubMed2005

The role of electrostatic interaction in triggering the unraveling of stable helix 1 in normal prion protein. A molecular dynamics simulation investigation.

The conversion of normal prion protein (PrPC) into scrapie isoform (PrPSc) is a key event in the pathogenesis of prion diseases. However, the conversion mechanism has given rise to much controversy. For instance, there is much debate on the behavior of helix 1 (H1) in the conversion. A series of experiments demonstrated that H1 in isolated state was very stable under a variety of conditions. But, other experiments indicated that helices 2 and 3 rather than H1 were retained in PrPSc. In this paper, molecular dynamics (MD) simulation is employed to investigate the dynamic behavior of H1. It is revealed that although the helix 1 of Human PrPC (HuPrPC) is very stable in the isolated state, it becomes unstable when incorporated into native HuPrPC, which likely results from the long-range electrostatic interaction between Asp147 and Arg208 located in the helices 1 and 3, respectively. This explanation is supported by experimental evaluation and MD simulation on D147N mutant of HuPrPC that the mutant becomes a little more stable than the wild type HuPrPC. This finding not only help to reconcile the existing debate on the role of helix 1 in the PrPC-->PrPSc transition, but also reveals a possible mechanism for triggering the PrPC-->PrPSc conversion.

Aspartic Acid↗
PubMed2004

Influence of a fluorobenzene nucleobase analogue on the conformational flexibility of RNA studied by molecular dynamics simulations.

Chemically modified bases are frequently used to stabilize nucleic acids, to study the driving forces for nucleic acid structure formation and to tune DNA and RNA hybridization conditions. In particular, fluorobenzene and fluorobenzimidazole base analogues can act as universal bases able to pair with any natural base and to stabilize RNA duplex formation. Although these base analogues are compatible with an A-form RNA geometry, little is known about the influence on the fine structure and conformational dynamics of RNA. In the present study, nano-second molecular dynamics (MD) simulations have been performed to characterize the dynamics of RNA duplexes containing a central 1'-deoxy-1'-(2,4-difluorophenyl)-beta-D-ribofuranose base pair or opposite to an adenine base. For comparison, RNA with a central uridine:adenine pair and a 1'-deoxy-1'-(phenyl)-beta-D-ribofuranose opposite to an adenine was also investigated. The MD simulations indicate a stable overall A-form geometry for the RNAs with base analogues. However, the presence of the base analogues caused a locally enhanced mobility of the central bases inducing mainly base pair shear and opening motions. No stable 'base-paired' geometry was found for the base analogue pair or the base analogue:adenine pairs, which explains in part the universal base character of these analogues. Instead, the conformational fluctuations of the base analogues lead to an enhanced accessibility of the bases in the major and minor grooves of the helix compared with a regular base pair.

Adenine↗
PubMed2004

Site-specific variations of carbonyl chemical shift anisotropies in proteins.

By combining molecular dynamics (MD) simulation with a novel extended density functional theory method, we calculate site-specific carbonyl chemical shift tensors in the SMN Tudor domain. We formulate a simple model for the C' chemical shift anisotropy (CSA) based solely on the isotropic chemical shift. Using this simple chemical shift tensor model and the MD simulation, an accurate prediction of transverse C'/N-H cross-correlated relaxation rates can be obtained.

Anisotropy↗
PubMed2005

Free-energy perturbation methods to study structure and energetics of DNA adducts: results for the major N2-dG adduct of benzo[a]pyrene in two conformations and different sequence contexts.

The potent mutagen/carcinogen benzo[a]pyrene (B[a]P) is activated to (+)-anti-B[a]PDE, which induces a variety of mutations (e.g., G --> T, G --> A, etc.) via its major adduct [+ta]-B[a]P-N2-dG. One hypothesis is that adducts (such as [+ta]-B[a]P-N2-dG) induce different mutations via different conformations, probably when replicated by different lesion-bypass DNA polymerases (DNAPs). We showed that Escherichia coli DNAP V was responsible for G --> T mutations with [+ta]-B[a]P-N2-dG in a 5'-TGT sequence (Yin et al., (2004) DNA Repair 3, 323), so we wish to study conformations of this adduct/sequence context by molecular modeling. The development of a CHARMM-based molecular dynamics (MD) simulations protocol with free-energy calculations in the presence of solvent and counterions is described. A representative base-pairing and base-displaced conformation of [+ta]-B[a]P-N2-dG in the 5'-TGT sequence are used: (1) BPmi5, which has the B[a]P moiety in the minor groove pointing toward the base on the 5'-side of the adduct, and (2) Gma5, which has the B[a]P moiety stacked with the surrounding base pairs and the dG moiety displaced into the major groove. The MD output structures are reasonable when compared to known NMR structures. Changes in DNA sequence context dramatically affect the biological consequences (e.g., mutagenesis) of [+ta]-B[a]P-N2-dG. Consequently, we also developed a MD-based free-energy perturbation (FEP) protocol to study DNA sequence changes. FEP involves the gradual "fading-out" of atoms in a starting structure (A) and "fading-in" of atoms in a final structure (B), which allows a realistic assessment of the energetic and structural changes when two structures A and B are closely related. Two DNA sequence changes are described: (1) 5'-TGT --> 5'-TGG, which involves two steps [T:A --> T:C --> G:C], and (2) 5'-TGT --> 5'-TGC, which involves three steps [T:A --> T:2AP --> C:2AP --> C:G], where 2AP (2-aminopurine) is included, because T:2AP and C:2AP retain more-or-less normal pairing orientations between complementary bases. FEP is also used to evaluate the impact that a 5'-TGT to 5'-UGT sequence change might have on mutagenesis with [+ta]-B[a]P-N2-dG. In summary, we developed (1) a CHARMM-based molecular dynamics (MD) simulations protocol with free-energy calculations in the presence of solvent and counterions to study B[a]P-N2-dG adducts in DNA duplexes, and (2) a MD-based free-energy perturbation (FEP) protocol to study DNA sequence context changes around B[a]P-N2-dG adducts.

Base Sequence↗
PubMed2005

Single-strand conformation polymorphism (SSCP) of oligodeoxyribonucleotides: an insight into solution structural dynamics of DNAs provided by gel electrophoresis and molecular dynamics simulations.

Studies on the solution structure dynamics of RNA/DNA are becoming crucially important. The phenomena of SSCP (single-strand conformation polymorphism), small RNA dynamics in a cell, and others can be related to the conformational changes of single-stranded (ss) RNAs/DNAs in solution. However, little is known about those dynamics. Only the intra-structural transition of ssDNAs in solution has been reported based on Watson-Crick (W-C) base-pairing. Here, we found a general feature of the SSCP phenomenon by studying the simpler molecules of ss-oligodeoxyribonucleotides. A single base substitution or a positional exchange of nucleotide in a highly homologous series of ss-dodecanucleotides led to a change in the mobility-in-gel. This was unexpected, since most of these nucleotides [such as d(A(11)G) or d(A(11)C)] have no possibility of forming W-C base-pairing. MD (molecular dynamics) experiments revealed differences in shape and size between the dynamic structures of these molecules which could affect their mobility-in-gel. In addition, a high correlation was observed between the electrophoretic mobility and the size-related parameters such as end-to-end distance obtained from MD simulations. Because the simulation was considerably shorter (nanosecond) than the experimental time-scale (second), the result must be considered conservatively; but it is nevertheless encouraging for utilizing MD simulation for structural analysis of oligonucleotides.

Base Sequence↗
PubMed2006

Influence of the disulfide bond configuration on the dynamics of the spin label attached to cytochrome c.

A series of multi-nanosecond molecular dynamics (MD) simulations of wild-type cytochrome c and its spin-labeled variants with the methanethiosulfonate moiety attached at position C102 were performed (1) to elucidate the effect of the spin probe presence on the protein structure and (2) to describe the structure and dynamics of the spin-label moiety. Comparisons with the reference crystal structure of cytochrome c (PDB entry: 1YCC) indicate that the protein secondary structure is well preserved during simulations of the wild-type cytochrome c but slightly changed in simulations of the cytochrome c labeled at position C102. At the time scale covered in our simulations, the spin label exhibits highly dynamical behavior. The number of observed distinct conformations of the spin label moiety is between 3 and 13. The spin probe was found to form short-lived hydrogen bonds with the protein. Temporary hydrophobic interactions between the probe and the protein were also detected. The MD simulations directly show that the disulfide bond in the tether linking a spin probe with a protein strongly influence the behavior of the nitroxide group. The conformational flexibility and interaction with the protein are different for each of the two low energy conformations of the disulfide bond.

Computer Simulation↗
PubMed2005

Effects of molecular association on polarizability relaxation in liquid mixtures of benzene and hexafluorobenzene.

In this work we have studied the relaxation dynamics of the many-body polarizability anisotropy in liquid mixtures of benzene (Bz) and hexafluorobenzene (Hf) at room temperature by femtosecond optical heterodyne-detected Raman-induced Kerr effect spectroscopy (OHD-RIKES) experiments and molecular dynamics (MD) simulations. The computed polarizability response arising from intermolecular interactions was included using the first-order dipole-induced-dipole model with the molecular polarizability distributed over the carbon sites of each molecule. We found good qualitative agreement between experiments and simulations in the features exhibited by the nuclear response function R(t) for pure liquids and mixtures. The long-time diffusive decay of R(t) was observed to vary substantially with composition, slowing down noticeably with dilution of each of the species as compared with that in the corresponding pure liquids. MD simulation shows that the effect on R(t) is due to the formation of strong and localized intermolecular association between Bz and Hf species that hinder the rotational diffusive dynamics. The formation of these Bz-Hf complexes in the liquid mixtures also modifies the rotational diffusive dynamics of the component species in such a way that cannot be explained solely in terms of a viscosity effect. Even though the computed orientational diffusive relaxation times associated with Bz and Hf are larger by a factor of approximately 2 than those from experiments, we found similar trends in experiments and simulations for these characteristic times as a function of composition. Namely, the collective and single-molecule orientational correlation times associated with Bz are observed to grow monotonically with the dilution of Bz, while those corresponding to Hf species exhibit a maximum at the equimolar composition. We attribute the quantitative discrepancy between experiments and simulations to the use of the Williams potential, which seems to overestimate the intermolecular interactions and thus predicts not only a slower translational dynamics but also a slower rotational diffusion dynamics than in real fluids.

Journal Article↗
PubMed2006

Nanopore-facilitated, voltage-driven phosphatidylserine translocation in lipid bilayers--in cells and in silico.

Nanosecond, megavolt-per-meter pulses--higher power but lower total energy than the electroporative pulses used to introduce normally excluded material into biological cells--produce large intracellular electric fields without destructively charging the plasma membrane. Nanoelectropulse perturbation of mammalian cells causes translocation of phosphatidylserine (PS) to the outer face of the cell, intracellular calcium release, and in some cell types a subsequent progression to apoptosis. Experimental observations and molecular dynamics (MD) simulations of membranes in pulsed electric fields presented here support the hypothesis that nanoelectropulse-induced PS externalization is driven by the electric potential that appears across the lipid bilayer during a pulse and is facilitated by the poration of the membrane that occurs even during pulses as brief as 3 ns. MD simulations of phospholipid bilayers in supraphysiological electric fields show a tight association between PS externalization and membrane pore formation on a nanosecond time scale that is consistent with experimental evidence for electropermeabilization and anode-directed PS translocation after nanosecond electric pulse exposure, suggesting a molecular mechanism for nanoelectroporation and nanosecond PS externalization: electrophoretic migration of the negatively charged PS head group along the surface of nanometer-diameter electropores initiated by field-driven alignment of water dipoles at the membrane interface.

Biological Transport↗
PubMed2006

Catalytic mechanism of cyclophilin as observed in molecular dynamics simulations: pathway prediction and reconciliation of X-ray crystallographic and NMR solution data.

Cyclophilins are proteins that catalyze X-proline cis-trans interconversion, where X represents any amino acid. Its mechanism of action has been investigated over the past years but still generates discussion, especially because until recently structures of the ligand in the cis and trans conformations for the same system were lacking. X-ray crystallographic structures for the complex cyclophilin A and HIV-1 capsid mutants with ligands in the cis and trans conformations suggest a mechanism where the N-terminal portion of the ligand rotates during the cis-trans isomerization. However, a few years before, a C-terminal rotating ligand was proposed to explain NMR solution data. In the present study we use molecular dynamics (MD) simulations to generate a trans structure starting from the cis structure. From simulations starting from the cis and trans structures obtained through the rotational pathways, the seeming contradiction between the two sets of experimental data could be resolved. The simulated N-terminal rotated trans structure shows good agreement with the equivalent crystal structure and, moreover, is consistent with the NMR data. These results illustrate the use of MD simulation at atomic resolution to model structural transitions and to interpret experimental data.

Catalysis↗
PubMed2004

Molecular modeling of green fluorescent protein: structural effects of chromophore deprotonation.

Molecular dynamics (MD) simulations were carried out to study the conformational rearrangement induced by deprotonation of the fluorescent chromophore in GFP, as well as the associated changes in the hydrogen-bonding network. For both the structures with either a neutral or an anionic chromophore, it was found that the beta-barrel was stable and rigid, and the conformation of the chromophore was consistent with the available x-ray structure. The conformational change in Thr203 due to deprotonation was also found to be consistent with the three-state isomerization model. Although GFP is highly fluorescent, denatured-GFP is nonfluorescent, indicating that the environment of the protein plays an important role in its fluorescence behavior. Our MD simulations, which explore the effect of the protein shell on the conformation of the chromophore, find the flexibility of the central chromophore to be significantly restricted due to the rigid nature of the protein shell. The hydrogen-bonding between the chromophore and neighboring residues was also shown to contribute to the chromophore rigidity. In addition to the MD studies, quantum mechanics/molecular mechanics (QM/MM) ONIOM calculations were carried out to investigate the effect of the beta-barrel on the internal rotation in the chromophore. Along with providing quantitative values for torsional rotation barriers about the bridging bond in the chromophore, the ONIOM calculations also validate our MD force field parameters.

Computer Simulation↗
PubMed2006

X-ray and neutron scattering studies of the hydration structure of alkali ions in concentrated aqueous solutions.

The presence of ions in water provides a rich and varied environment in which many natural processes occur with important consequences in biology, geology and chemistry. This article will focus on the structural properties of ions in water and it will be shown how the 'difference' methods of neutron diffraction with isotopic substitution (NDIS) and anomalous X-ray diffraction (AXD) can be used to obtain direct information regarding the radial pair distribution functions of many cations and anions in solution. This information can subsequently be used to calculate coordination numbers and to determine ion-water conformation in great detail. As well as enabling comparisons to be made amongst ions in particular groups in the periodic table, such information can also be contrasted with results provided by molecular dynamics (MD) simulation techniques. To illustrate the power of these 'difference' methods, reference will be made to the alkali group of ions, all of which have been successfully investigated by the above methods, with the exception of the radioactive element francium. Additional comments will be made on how NDIS measurements are currently being combined with MD simulations to determine the structure around complex ions and molecules, many of which are common in biological systems.

Ions↗
PubMed2004

How does Mg(2+) affect the binding of anhydrotetracycline in the TetR protein?

The binding of anhydrotetracycline (atc) in wild-type TetR(D), TetR(B), and four single tryptophan mutants of TetR(B) was investigated by UV/vis absorption, steady state, and time-resolved fluorescence spectroscopy. From absorption titration experiments with Mg2+, we conclude that binding of one [atc-Mg]+ complex in the homodimer causes changes in the protein conformation around the second binding pocket. In the presence of absence of Mg2+, several different groups of atc-protein arrangements must exist, each with a characteristic atc fluorescence decay time. Taking into account the results of molecular dynamics (MD) simulations, we propose as one possible origin for such a differentiation teh extent of hydrogen bonding between atc and the surrounding amino acids. Binding of Mg2+ should change the arrangement of the surrounding amino acids such that some of the excited atc molecules do not undergo the relaxation process typical for free atc. The MD simulations also show that the pattern of intra- and intermolecular hydrogen bonding in the two monomeric units is no correlated, thereby leading to different fluorescence kinetics for atc in the two monomeric units. Furthermore, it is suggested that hydrogen bonding between Arg104 and O10 of anhydrotetracycline could regulate the relaxation processes of excited anhydrotetracycline.

Bacterial Proteins↗
PubMed2005

Effect on DNA relaxation of the single Thr718Ala mutation in human topoisomerase I: a functional and molecular dynamics study.

The functional and dynamical properties of the human topoisomerase I Thr718Ala mutant have been compared to that of the wild-type enzyme using functional assays and molecular dynamics (MD) simulations. At physiological ionic strength, the cleavage and religation rates, evaluated on oligonucleotides containing the preferred topoisomerase I DNA sequence, are almost identical for the wild-type and the mutated enzymes, as is the cleavage/religation equilibrium. On the other hand, the Thr718Ala mutant shows a decreased efficiency in a DNA plasmid relaxation assay. The MD simulation, carried out on the enzyme complexed with its preferred DNA substrate, indicates that the mutant has a different dynamic behavior compared to the wild-type enzyme. Interestingly, no changes are observed in the proximity of the mutation site, whilst a different flexibility is detected in regions contacting the DNA scissile strand, such as the linker and the V-shaped alpha helices. Taken together, the functional and simulation results indicate a direct communication between the mutation site and regions located relatively far away, such as the linker domain, that with their altered flexibility confer a reduced DNA relaxation efficiency. These results provide evidence that the comprehension of the topoisomerase I dynamical properties are an important element in the understanding of its complex catalytic cycle.

Alanine↗
PubMed1998

Coarse-grained simulations of conformational dynamics of proteins: application to apomyoglobin.

A coarse-grained dynamic Monte Carlo method is proposed for investigating the conformational dynamics of proteins. Each residue is represented by two interaction sites, one at the alpha-carbon, and the other on the amino acid sidechain. Geometry and energy parameters extracted from databank structures are used. The method is applied to the crystal structure of apomyoglobin (apo-Mb). Equilibrium and dynamic properties of apo-Mb are characterized within computation times one order of magnitude shorter than conventional molecular dynamics (MD) simulations. The calculated rms fluctuations in alpha-carbons are in good agreement with crystallographic temperature factors. Regions exhibiting enhanced conformational mobilities are identified. Among the loops connecting the eight helices A to H, the loop CD undergoes the fastest motions, leading to partial unwinding of helix D. Helix G is the most stable helix on the basis of the kinetic stability of dihedral angles, followed by the respective helices A, E, H, and B. These results, in agreement with H/D exchange and two-dimensional NMR experiments, as well as with MD simulations, lend support to the use of the proposed approach as an efficient, yet physically plausible, means of characterizing protein conformational dynamics.

Animals↗
PubMed1996

An extended sampling of the configurational space of HPr from E. coli.

Recently, we developed a method (Amadei et al., J. Biomol. Str. Dyn. 13: 615-626; de Groot et al., J. Biomol. Str. Dyn. 13: 741-751, 1996) to obtain an extended sampling of the configurational space of proteins, using an adapted form of molecular dynamics (MD) simulations, based on the essential dynamics (ED) (Amadei et al., Proteins 17:412-425, 1993) method. In the present study, this ED sampling technique is applied to the histidine-containing phosphocarrier protein HPr from Escherichia coli. We find a cluster of conformations that is an order of magnitude larger than that found for a usual MD simulation of comparable length. The structures in this cluster are geometrically and energetically comparable to NMR structures. Moreover, on average, this large cluster satisfies nearly all NMR-derived distance restraints.

Bacterial Proteins↗
PubMed2006

NMR chemical shift powder pattern recoupling at high spinning speed and theoretical tensor evaluation applied to silk fibroin.

The NMR pulse sequence RAI (recoupling of anisotropy information) has been improved to obtain powder patterns at high MAS spinning speeds. The 2D iso-aniso experiment displays the static chemical shift spectra on the indirect dimension and the MAS spectra on the direct dimension; hence overlapping chemical shift tensor patterns can be well resolved. This efficient technique is applicable to compounds containing (13)C sp(3) (C(alpha), C(beta)) and sp(2) (C=O) sites with higher chemical shift (CS) anisotropy (CSA), and the reliability of the method was tested here on the (13)C chemical shift tensors of polycrystalline glycine, alanine, and serine. Subsequently, the same experiment was applied to the native silk protein fibroin from Bombyx mori, which consists mainly of these three amino acids. Molecular dynamics (MD) simulations of the silk II crystal structure of Takahashi et al. (Takahashi et al. Int. J. Biol. Macromol. 1999, 24, 127-138) were carried out to study the influence of motions on the chemical shift tensors. The (13)C chemical shift tensors were calculated using the bond polarization theory BPT on 200 structures created by an MD simulation. Very good agreement of the theoretical chemical shift anisotropy values with the experimental NMR results was obtained. The tensor orientations in the protein structure could thus be reliably derived.

Animals↗
PubMed2002

Oxygenated metabolites of anandamide and 2-arachidonoylglycerol: conformational analysis and interaction with cannabinoid receptors, membrane transporter, and fatty acid amide hydrolase.

This study was aimed at finding structural requirements for the interaction of the acyl chain of endocannabinoids with cannabinoid receptors, membrane transporter protein, and fatty acid amide hydrolase (FAAH). To this end, the flexibility of the acyl chain was restricted by introduction of an 1-hydroxy-2Z,4E-pentadiene system in anandamide (N-arachidonoylethanolamine, AEA) and 2-arachidonoylglycerol (2-AG) at various positions using different lipoxygenases. This brought about selectivity and attenuated the binding potency of AEA and 2-AG. Although the displacement constants were modest, 15(S)-hydroxy-eicosa-5Z,8Z,11Z,13E-tetraenoyl-N-(2-hydroxyethyl)amine was found to bind selectively to the CB(1) receptor, whereas its 1-arachidonoyl-sn-glycerol analogue and 13(S)-hydroxy-octadeca-9Z,11E-dienoyl-N-(2-hydroxyethyl)amine could selectively bind to the CB(2) receptor. 11(S)-Hydroxy-eicosa-5Z,8Z,12E,14Z-tetraenoyl-N-(2-hydroxyethyl)amine did not bind to either receptor, whereas 12(S)-hydroxy-eicosa-5Z,8Z,10E,14Z-tetraenoyl-N-(2-hydroxyethyl)amine did bind to both CB receptors with an affinity similar to that of AEA. All oxygenated anandamide derivatives were good inhibitors of FAAH (low micromolar K(i)) but were ineffective on the AEA transporter. 2-AG rapidly isomerizes into 1(3)-arachidonoyl-sn-glycerol. Both 1- and 3-arachidonoyl-sn-glycerol did not bind to either CB receptor and did not interfere with AEA transport. Thus, after it is isomerized, 2-AG is inactivated, thereby decreasing effective concentrations of 2-AG. Analysis of (1)H NMR spectra revealed that chloroform did not induce notably different conformations in the acyl chain of 15(S)-hydroxy-eicosa-5Z,8Z,11Z,13E-tetraenoic acid as compared with water. Molecular dynamics (MD) simulations of AEA and its analogues in the presence of explicit water molecules revealed that a tightly folded conformation of the acyl chain is not the only requirement for CB(1) binding. Structural details of the C(2)-C(15) loop, such as an sp(2) carbon at position 11, are necessary for receptor binding. The MD simulations may suggest that the average orientations of the pentyl tail of AEA and 12(S)-hydroxy-eicosa-5Z,8Z,10E,14Z-tetraenoyl-N-(2-hydroxyethyl)amine are different from that of the low-affinity, inactive ligands.

Amidohydrolases↗
PubMed2002

A structure-based design approach for the identification of novel inhibitors: application to an alanine racemase.

We report a new structure-based strategy for the identification of novel inhibitors. This approach has been applied to Bacillus stearothermophilus alanine racemase (AlaR), an enzyme implicated in the biosynthesis of the bacterial cell wall. The enzyme catalyzes the racemization of L- and D-alanine using pyridoxal 5'-phosphate (PLP) as a cofactor. The restriction of AlaR to bacteria and some fungi and the absolute requirement for D-alanine in peptidoglycan biosynthesis make alanine racemase a suitable target for drug design. Unfortunately, known inhibitors of alanine racemase are not specific and inhibit the activity of other PLP-dependent enzymes, leading to neurological and other side effects. This article describes the development of a receptor-based pharmacophore model for AllaR, taking into account receptor flexibility (i.e. a 'dynamic' pharmacophore model). In order to accomplish this, molecular dynamics (MD) simulations were performed on the full AlaR dimer from Bacillus stearothermophilus (PDB entry, 1 sft) with a D-alanine molecule in one active site and the non-covalent inhibitor, propionate, in the second active site of this homodimer. The basic strategy followed in this study was to utilize conformations of the protein obtained during MD simulations to generate a dynamic pharmacophore model using the property mapping capability of the LigBuilder program. Compounds from the Available Chemicals Directory that fit the pharmacophore model were identified and have been submitted for experimental testing. The approach described here can be used as a valuable tool for the design of novel inhibitors of other biomolecular targets.

Alanine Racemase↗