Search PubMed⌕ Search

SEARCH · Search PubMed

Results for “Molecular Dynamic simulation”

Search indexed PubMed citations on genomics, clinical trials, systematic reviews and public health. Explore titles, authors and supplied subject terms, then open the PubMed record.

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

At least 1,279 records · Page 71Linked to original sources
PubMed1998

The solution conformations of amino acids from molecular dynamics simulations of Gly-X-Gly peptides: comparison with NMR parameters.

The conformations that amino acids can adopt in the random coil state are of fundamental interest in the context of protein folding research and studies of protein-peptide interactions. To date, no detailed quantitative data from experimental studies have been reported; only nuclear magnetic resonance parameters such as chemical shifts and J coupling constants have been reported. These experimental nuclear magnetic resonance data represent averages over multiple conformations, and hence they do not provide unique structural information. I have performed relatively long (2.5 ns) molecular dynamics simulations of Gly-X-Gly tripeptides, surrounded by explicit water molecules, where X represents eight different amino acids with long side chains. From the trajectories one can calculate time averaged backbone chemical shifts and 3J(NH alpha) coupling constants and compare these with experimental data. These calculated quantities are quite close to the experimental values for most amino acids, suggesting that these simulations are a good model for the random coil state of the tripeptides. On the basis of my simulations I predict 3J(alphabeta) coupling constants and present dihedral distributions for the phi, psi, as well as chi1 and chi2 angles. Finally, I present correlation plots for these dihedral angles.

Amino Acids↗
PubMed1996

Molecular-dynamics simulation of domain movements in aspartate aminotransferase.

Mitochondrial aspartate aminotransferase is a homodimeric protein with 2 x 402 amino acid residues. The enzyme in solution undergoes ligand-induced and syncatalytic conformational changes which appear to correspond to shifts in the equilibrium between the crystallographically defined open and closed conformation. In the closed conformation, the small domain of each subunit has rotated as a rigid body by 13 degrees and 14 degrees towards the large coenzyme-binding domain and has closed the active-site pocket. Molecular dynamics simulations at 300 K of 120-ps duration were started from the crystal structures of the unliganded pyridoxal form (open form) of the dimeric enzyme and the 2-methylaspartate-liganded closed form in which the 2-methyl group had been removed. Both structures contained the crystal water molecules and were placed in a 5-A shell of water. The rms fluctuations of the individual C alpha atoms during the simulations agreed well with the corresponding B factors of the crystal structures. Superposition of the initial structures and the average structures of the last 20 ps showed in both simulations extensive C alpha deviations in the case of the whole subunit but much smaller changes in the individual large and small domains, indicating a movement of the two domains relative to each other. In the simulation of the open form, superposition of the large domains made evident a displacement of the small domain towards its position in the closed crystal structure, which can be described by a rotation of the small domain by about 13 degrees around the twofold symmetry (z) axis. A significantly less extensive rearrangement of parts of the small domain, i.e. a rotation of about 5 degrees around the z axis, was observed in the simulation of the substrate-liganded enzyme (closed form) which, in contrast to the open form, showed only small conformational changes around the active site. In both simulations an additional rotation of the small domain by 9 degrees around the x axis occurred. The actual domain movement is estimated to occur in a time range at least two orders of magnitude larger than the simulation time of 120 ps. Apparently, the surface tension of the unrestrained nonspherical shell of water accelerates the simulated conformational change which, however, quite closely imitates the geometric features of the extensive movement of the small domains (each approximately 130 residues).

Animals↗
PubMed2003

Understanding the acylation mechanisms of active-site serine penicillin-recognizing proteins: a molecular dynamics simulation study.

Beta-lactam antibiotics inhibit enzymes involved in the last step of peptidoglycan synthesis. These enzymes, also identified as penicillin-binding proteins (PBPs), form a long-lived acyl-enzyme complex with beta-lactams. Antibiotic resistance is mainly due to the production of beta-lactamases, which are enzymes that hydrolyze the antibiotics and so prevent them reaching and inactivating their targets, and to mutations of the PBPs that decrease their affinity for the antibiotics. In this study, we present a theoretical study of several penicillin-recognizing proteins complexed with various beta-lactam antibiotics. Hybrid quantum mechanical/molecular mechanical potentials in conjunction with molecular dynamics simulations have been performed to understand the role of several residues, and pK(a) calculations have also been done to determine their protonation state. We analyze the differences between the beta-lactamase TEM-1, the membrane-bound PBP2x of Streptococcus pneumoniae, and the soluble DD-transpeptidase of Streptomyces K15.

Acylation↗
PubMed1993

A comparison of 15N NMR relaxation measurements with a molecular dynamics simulation: backbone dynamics of the glucocorticoid receptor DNA-binding domain.

The rapid motions of the backbone of the DNA-binding domain of the glucocorticoid receptor (GR DBD) have been investigated using proton-detected heteronuclear NMR experiments on 15N-labeled protein at pH 6.0 and with a 200 psec molecular dynamics simulation of hydrated GR DBD. The experimental data were interpreted in terms of a generalized order parameter (S2) and an effective correlation time (tau e) for the internal motion of each amide bond. A back calculation, using the same model, yielded the [1H]-14N nuclear Overhauser effects (NOEs) and the 15N spin-lattice relaxation times (T1) from the simulated data. The rapid motions of the backbone turned out to be rather limited and uniform throughout the protein, with a somewhat reduced mobility in the two major alpha-helical regions and a slightly enhanced flexibility for some residues in the first zinc coordinating region. The agreement between the experimental and simulated S2-values was as good as quantitative for most of the residues, except for some residues that were subject to a more large-scale, and in the simulation thus poorly sampled, motion. Examples of such motions that were found in the simulation include jumps of the amide bond of Ile-487 between the charged oxygens of the side chain of Asp-485 and less distinct large scale motions for some of the residues in the extended regions, that were shown to give rise to noisy and/or fast decaying internal reorientational correlation functions. For these residues large differences in the simulated and experimental tau e-values were found, indicating that motions on different time scales were dominating in the experimental and simulated values. The lower (< 0.7) experimental NOEs for these residues could not be reproduced in the simulation and were shown to be a consequence of the lower tau e-values estimated in the simulation. By combining information from the simulation and the experiment a more complete picture of the motions for these residues can be obtained as is illustrated with an estimation of the jump angle and jump frequency for the amide bond of Ile-487.

Amino Acid Sequence↗
PubMed2001

Can non-mechanical proteins withstand force? Stretching barnase by atomic force microscopy and molecular dynamics simulation.

Atomic force microscopy (AFM) experiments have provided intriguing insights into the mechanical unfolding of proteins such as titin I27 from muscle, but will the same be possible for proteins that are not physiologically required to resist force? We report the results of AFM experiments on the forced unfolding of barnase in a chimeric construct with I27. Both modules are independently folded and stable in this construct and have the same thermodynamic and kinetic properties as the isolated proteins. I27 can be identified in the AFM traces based on its previous characterization, and distinct, irregular low-force peaks are observed for barnase. Molecular dynamics simulations of barnase unfolding also show that it unfolds at lower forces than proteins with mechanical function. The unfolding pathway involves the unraveling of the protein from the termini, with much more native-like secondary and tertiary structure being retained in the transition state than is observed in simulations of thermal unfolding or experimentally, using chemical denaturant. Our results suggest that proteins that are not selected for tensile strength may not resist force in the same way as those that are, and that proteins with similar unfolding rates in solution need not have comparable unfolding properties under force.

Animals↗
PubMed2005

Ultrasonic velocities, densities, viscosities, electrical conductivities, Raman spectra, and molecular dynamics simulations of aqueous solutions of Mg(OAc)2 and Mg(NO3)2: Hofmeister effects and ion pair formation.

The ultrasonic velocities, densities, viscosities, and electrical conductivities of aqueous solutions of magnesium nitrate and magnesium acetate have been measured from dilute to saturation concentrations at 0 < or = t/degrees C < or = 50. The temperature derivative of the isentropic compressibility, kappa(s), became zero at 2.28 and 2.90 mol kg(-1) for Mg(OAc)2 and Mg(NO3)2 solutions, respectively, at 25 degrees C. The total hydration numbers of the dissolved ions were estimated to be, respectively, 24.3 and 19.2 at these concentrations. Differences in kappa(s) for various M2+ salts, using the present and literature data, correlated with reported M2+-OH2 bond lengths and to a lesser extent with cationic charge densities (ionic radii). The influence of anions on kappa(s) appears to follow the Hofmeister series and also correlates approximately with the anionic charge density. Substantial differences between Mg(OAc)2(aq) and Mg(NO3)2(aq) occur with respect to their structural relaxation times (derived from compressibility and viscosity data) and their electrical conductivities. These differences were attributed to a much greater ion association in Mg(OAc)2 solutions. Raman spectra recorded at 28 degrees C confirmed the presence of various types of contact ion pairs including mono- and bidentate complexes in Mg(OAc)2(aq). In Mg(NO3)2(aq), only noncontact ion pairs appear to be formed even at high concentrations. The experimental results are supported by molecular dynamics simulations, which also reveal the much stronger tendency of OAc- compared to NO3- to associate with Mg2+ in aqueous solutions. The simulations also allow an evaluation of the ion-ion and ion-water radial distribution functions and cumulative sums and provide a molecular picture of ion hydration in Mg(OAc)2(aq) and Mg(NO3)2(aq) at varying concentrations.

Journal Article↗
PubMed2005

Molecular Dynamics Simulation of WSK-3, a Computationally Designed, Water-Soluble Variant of the Integral Membrane Protein KcsA.

Poor solubility and low expression levels often make membrane proteins difficult to study. An alternative to the use of detergents to solubilize these aggregation-prone proteins is the partial redesign of the sequence so as to confer water solubility. Recently, computationally assisted membrane protein solubilization (CAMPS) has been reported, where exposed hydrophobic residues on a protein's surface are computationally redesigned. Herein, the structure and fluctuations of a designed, water-soluble variant of KcsA (WSK-3) were studied using molecular dynamics simulations. The root mean square deviation of the protein from its starting structure, where the backbone coordinates are those of KcsA, was 1.8 angstroms. The structure of salt bridges involved in structural specificity and solubility were examined. The preferred configuration of ions and water in the selectivity filter of WSK-3 was consistent with the reported preferences for KcsA. The structure of the selectivity filter was maintained, which is consistent with WSK-3 having an affinity for agitoxin2 comparable to that of wild-type KcsA. In contrast to KcsA, the central cavity's side chains were observed to reorient, allowing water diffusion through the side of the cavity wall. These simulations provide an atomistic analysis of the CAMPS strategy and its implications for further investigations of membrane proteins.

Bacterial Proteins↗
PubMed2003

Probing cofactor specificity in phenylalanine hydroxylase by molecular dynamics simulations.

Phenylalanine hydroxylase (PAH) is a tetrahydrobiopterin-dependent enzyme that catalyzes the hydroxylation of L-phenylalanine (L-Phe) to L-tyrosine using dioxygen as an additional substrate. The requirement of PAH for a cofactor is absolute, but several cofactor analogs are able to substitute the natural cofactor in catalysis. However, it is only the natural cofactor 6R-tetrahydrobiopterin (6R-BH(4)) that induces a negative regulatory effect on the enzyme. In order to get further insights on the molecular basis for this specificity, we studied the structure of the cofactor-enzyme complex and the conformational changes induced by cofactor binding by molecular dynamics simulations. Simulations were carried out on the enzyme alone and complexed with 6R-BH(4) and with two cofactor analogs, 6S-BH(4) and 6-methyl-tetrahydropterin (6M-PH(4)). In the resting unbound enzyme Tyr377 in the catalytic domain is hydrogen bonded to both Ser23 and Glu21 of the autoregulatory N-terminal sequence. This hydrogen bonding network is disturbed by the binding of BH(4), which interacts with Ser23. By doing so, 6R-BH(4) facilitates an interaction between Glu21 and the active site iron, further pulling the N-terminal into the active site of PAH and blocking the L-Phe binding site. Thus, in the 6R-BH(4) complexed enzyme, the N-terminal functions as an intrinsic amino acid regulatory sequence (IARS). Neither 6M-PH(4) nor 6S-BH(4) can interact favorably with Ser23, and do not induce an inhibitory effect on PAH. These simulations thus explain the previous findings that the two hydroxyl groups in the side chain of the 6R epimer of BH(4) are essential for the inhibitory regulatory effect on PAH.

Binding Sites↗
PubMed2005

Solvation of tetraalkylammonium chlorides in acetonitrile-water mixtures: mass spectrometry and molecular dynamics simulations.

The solvation of tetramethylammonium chloride (Me4NCl) and tetra-n-butylammonium chloride (Bu4NCl) in water-acetonitrile mixtures was investigated by mass spectrometry of clusters isolated from the solution. As far as the positive ions are concerned, clusters composed of alkylammonium ions and acetonitrile molecules only were observed, even for mixtures with high water content. In contrast, for the negative ions, clusters composed of chloride with both water and/or acetonitrile molecules were observed. For the smaller system (Me4NCl) we performed quantum chemical calculations and molecular dynamics simulations. It was found that even though water is present in the solvation shell of Me4N+, only acetonitrile has a strong electrostatic interaction with the cation. Water molecules around Me4N+ form hydrogen bonds with other water molecules, and they interact with Me4N+ mainly via dispersive interactions. These results indicate that Me4N+ behaves like a hydrophobic solute. On the other hand, the interaction of Cl- with water and acetonitrile is of comparable strength and, in both cases, the electrostatic interaction dominates. Herein we demonstrate experimentally and theoretically that positive and negative ions give rise to characteristic solvation structures in mixed solvents: even a relatively small organic cation, such as Me4N+, exhibits a hydrophobic-like solvation shell.

Acetonitriles↗
PubMed2004

Molecular dynamics simulation links conformation of a pore-flanking region to hyperekplexia-related dysfunction of the inhibitory glycine receptor.

Inhibitory glycine receptors mediate rapid synaptic inhibition in mammalian spinal cord and brainstem. The previously identified hyperekplexia mutation GLRA1(P250T), located within the intracellular TM1-2 loop of the GlyR alpha1 subunit, results in altered receptor activation and desensitization. Here, elementary steps of ion channel function of alpha1(250) mutants were resolved and shown to correlate with hydropathy and molar volume of residue alpha1(250). Single-channel recordings and rapid activation kinetic studies using laser pulse photolysis showed reduced conductance but similar open probability of alpha1(P250T) mutant channels. Molecular dynamics simulation of a helix-turn-helix motif representing the intracellular TM1-2 domain revealed alterations in backbone conformation, indicating an increased flexibility in these mutants that paralleled changes in elementary steps of channel function. Thus, the architecture of the TM1-2 loop is a critical determinant of ion channel conductance and receptor desensitization.

Cell Line↗
PubMed2005

Molecular dynamics simulations of the adipocyte lipid binding protein reveal a novel entry site for the ligand.

The adipocyte lipid binding protein (ALBP) binds fatty acids (FA) in a cavity that is inaccessible from the bulk. Therefore, the penetration of the FA necessitates conformational changes whose nature is still unknown. It was suggested that the lipid first enters through a "portal region" which consists of the alphaII helix and the adjacent tight turns. The initial event in the ligand binding must be the interaction of the lipid with the protein surface. To analyze this interaction, we have carried out three molecular dynamics simulations of the apo-ALBP, with a palmitate ion initially located at different positions near the protein surface. The simulation indicated that the ligand could adsorb to the protein in more than one location. Yet, in one case, the ligand managed to penetrate the protein by entering a newly formed cavity some 10 A deep. The entry site is located near the N-terminus, at the junction between the loops connecting the beta-strands. Further progression of the penetration seems to be arrested by the need for desolvation of the COOH end of the headgroup. Evolutionary analysis showed that amino acids in this entry site are well conserved. On the basis of these observations, we suggest that the ligand may enter the protein from a site other than the portal region. Furthermore, the rate-limiting step is proposed to be the desolvation of the FA polar headgroup.

Adipocytes↗
PubMed2002

Molecular dynamics simulation of the RNA complex of a double-stranded RNA-binding domain reveals dynamic features of the intermolecular interface and its hydration.

The interaction between double-stranded RNA (dsRNA) and the third double-stranded domain (dsRBD) from Drosophila Staufen protein represents a paradigm to understand how the dsRBD protein family, one of the most common RNA-binding protein units, binds dsRNA. The nuclear magnetic resonance (NMR) structure of this complex and the x-ray structure of another family member revealed the stereochemical basis for recognition, but also raised new questions. Although the crystallographic studies revealed a highly ordered interface containing numerous water-mediated contacts, NMR suggested extensive residual motion at the interface. To address how interfacial motion contributes to molecular recognition in the dsRBD-dsRNA system, we conducted a 2-ns molecular dynamics simulation of the complex derived from Staufen protein and of the separate protein and RNA components. The results support the observation that a high degree of conformational flexibility is retained upon complex formation and that this involves interfacial residues that are critical for dsRBD-dsRNA binding. The structural origin of this residual flexibility is revealed by the analysis of the trajectory of motion. Individual basic side chains switch continuously from one RNA polar group to another with a residence time seldom exceeding 100 ps, while retaining favorable interaction with RNA throughout much of the simulation. Short-lived water molecules mediate some of these interactions for a large fraction of the trajectory studied here. This result indicates that water molecules are not statically associated with the interface, but continuously exchange with the bulk solvent on a 1-10-ps time scale. This work provides new insight into dsRBD-dsRNA recognition and builds upon a growing body of evidence, suggesting that short-lived dynamic interactions play important roles in protein-nucleic acid interactions.

Binding Sites↗
PubMed1992

Conformational states of a TT mismatch from molecular dynamics simulation of duplex d (CGCGATTCGCG).

The TT mismatch region in duplex d (CGCGATTCGCG) was studied using a 500-ps molecular dynamics (MD) simulation in water, and a series of 1-ps MD simulations and energy minimizations in vacuum. The DNA maintained its duplex structure, although the mismatch region showed significantly higher flexibility than the GC regions. The predominant conformation in the 500-ps MD simulation involved an average -42 degrees propeller twist between T6 and T'6, and a -22 degree buckle between A5 and T'7. One hydrogen bond was formed between T6 and T'6, and another between T6 and the O2 of T'7, with both Watson-Crick hydrogen bonds between A5 and T'7 remaining intact. The minimizations resulted in conformations with the equivalent hydrogen-bonding pattern, as well as ones with "wobble pair" hydrogen bonds between T6 and T'6. However, the wobble pair conformation was found to be unstable in the water simulation.

Base Sequence↗
PubMed2005

DASH: a novel analysis method for molecular dynamics simulation data. Analysis of ligands of PPAR-gamma.

A novel molecular dynamics (MD) analysis algorithm, DASH, is introduced in this paper. DASH has been developed to utilize the sequential nature of MD simulation data. By adjusting a set of parameters, the sensitivity of DASH can be controlled, allowing molecular motions of varying magnitudes to be detected or ignored as desired, with no knowledge of the number of conformations required being prerequisite. MD simulations of three synthetic ligands of the orphan nuclear receptor PPARgamma were generated in vacuo using Tripos's SYBYL and used as the training set for DASH. Two X-ray crystal structures of PPARgamma complexed with Rosiglitazone were compared to gain knowledge of the pharmacophoric conformation; this showed that the conformation of the ligand is significantly different between the two structures, indicating that there is no distinct conformation in which rosiglitazone binds to PPARgamma but multiple binding modes. An investigation into simulation length was carried out. A simulation of 5 ns was found to give highly variable results, whereas a simulation of 25 ns gave a representative window of motion for molecules of this size. DASH was compared with Ward's hierarchical cluster analysis method. The results show that DASH analysis is as good as Ward analysis in some areas (e.g. conformation identification) and is superior in others (e.g. speed and input size).

Algorithms↗
PubMed1996

Molecular dynamics simulations of the bis-intercalated complexes of ditercalinium and Flexi-Di with the hexanucleotide d(GCGCGC)2: theoretical analysis of the interaction and rationale for the sequence binding specificity.

The X-ray crystal structures of the complexes of ditercalinium and Flexi-Di with d(CGCG)2 have been studied by computational chemistry methods in an attempt to rationalize their distinct structural features. In addition, the complexes of these two bisintercalating drugs with d(GCGCGC)2 have been modeled and subjected to 0.5 ns of molecular dynamics simulations in explicit solvent with the aim of evaluating the relative importance of hydrogen bonding and stacking interactions in the sequence binding specificity of these compounds. According to our calculations, the electrostatic term is attractive for the stacking interactions between the pyridocarbazole chromophores of these drugs and the base pairs that make up the sandwiched GpC step. On the contrary, this energy term is repulsive for the base pairs that make up the boundaries of the bisintercalation site. This differential electrostatic binding energy component, which is shown to have a strong orientational dependence, could lie at the origin of the observed binding preferences of these drugs. In addition, both the Lennard-Jones and the electrostatic energy terms contribute to stabilizing the underwound central GpC step. The attractive electrostatic interactions between the linkers and the major groove are in concert with the stacking specificities for the sandwiched GpC step, which is thus very effectively stapled by the drugs. The hydrogen-bonding potential of the linkers, however, appears to be reduced in an aqueous medium due to competing interactions with water. Binding of either ditercalinium or Flexi-Di to d(GCGCGC)2 appears to favor the A-type conformation that this DNA molecule most likely adopts in the free state. The possible relevance of these findings to the process of bis-intercalation and to the pharmacological action of these compounds is discussed.

Carbazoles↗
PubMed1996

Molecular dynamics simulations of water within models of ion channels.

The transbilayer pores formed by ion channel proteins contain extended columns of water molecules. The dynamic properties of such waters have been suggested to differ from those of water in its bulk state. Molecular dynamics simulations of ion channel models solvated within and at the mouths of their pores are used to investigate the dynamics and structure of intra-pore water. Three classes of channel model are investigated: a) parallel bundles of hydrophobic (Ala20) alpha-helices; b) eight-stranded hydrophobic (Ala10) antiparallel beta-barrels; and c) parallel bundles of amphipathic alpha-helices (namely, delta-toxin, alamethicin, and nicotinic acetylcholine receptor M2 helix). The self-diffusion coefficients of water molecules within the pores are reduced significantly relative to bulk water in all of the models. Water rotational reorientation rates are also reduced within the pores, particularly in those pores formed by alpha-helix bundles. In the narrowest pore (that of the Ala20 pentameric helix bundle) self-diffusion coefficients and reorientation rates of intra-pore waters are reduced by approximately an order of magnitude relative to bulk solvent. In Ala20 helix bundles the water dipoles orient antiparallel to the helix dipoles. Such dipole/dipole interaction between water and pore may explain how water-filled ion channels may be formed by hydrophobic helices. In the bundles of amphipathic helices the orientation of water dipoles is modulated by the presence of charged side chains. No preferential orientation of water dipoles relative to the pore axis is observed in the hydrophobic beta-barrel models.

Alamethicin↗
PubMed2000

Homology modeling, molecular dynamics simulations, and analysis of CYP119, a P450 enzyme from extreme acidothermophilic archaeon Sulfolobus solfataricus.

The recent characterization of a thermophilic and barophilic CYP119 from Sulfolobus solfataricus offers a new opportunity to identify the origin of its stability by comparing it with mesophilic P450s with known structures. Since the three-dimensional structure of CYP119 is not yet available, homology modeling techniques were used to build model structures for this enzyme. The overall quality and stability of the models were assessed using three protein analysis programs and by monitoring structural stability during 1 ns of molecular dynamics simulations at 300 and 390 K. The results show the CYP119 models to be of good quality. Possible origins of the thermo- and barostability of CYP119 were then investigated by examining the amino acid compositions and the three-dimensional structure of CYP119 compared with the five mesophilic templates. Three possible factors were identified that could contribute to the enhanced stability of CYP119. The first was the higher relative population of salt bridges and the presence of a few unique salt bridges found in CYP119 that were absent in all five template CYP450s. The second factor was a decreased population of Ala and an increased population of Ile found in the interior of CYP119, which are likely to improve packing in CYP119. The third factor was a more extensive aromatic cluster seen in CYP119 which was not found in all five template P450s. In addition, the model CYP119 three-dimensional structures were also used to determine key properties related to its function. Specifically, binding site residues and surface residues for redox partner interactions were identified. These residues identified together with those residues found that might contribute to the increased stability are suggested for future mutagenesis studies. The results obtained from these experimental studies shall then provide further validation of the suggested origins of stability and the structure-function relationships derived from the model structures of this enzyme.

Amino Acid Sequence↗
PubMed2006

A molecular dynamics simulation study of an aminoglycoside/A-site RNA complex: conformational and hydration patterns.

Aminoglycoside antibiotics interfere with the translation mechanism by binding to the tRNA decoding site of the 16S ribosomal RNA. Crystallographic structures of aminoglycosides bound to A-site systems clarified many static aspects of RNA-ligand interactions. To gain some insight on the dynamic aspects of recognition phenomena, we conducted molecular dynamics simulations of the aminoglycoside paromomycin bound to a eubacterial ribosomal decoding A-site oligonucleotide. Results from 25 ns of simulation time revealed that: (i) the neamine part of the antibiotic represents the main anchor for binding, (ii) additional sugar rings provide limited and fragile contacts, (iii) long-resident water molecules present at the drug/RNA interface are involved in the recognition phenomena. The combination of MD simulations together with systematic structural information offers striking insights into the molecular recognition processes underlying RNA/aminoglycoside binding. Important methodological considerations related to the use of medium resolution starting structures and associated sampling problems are thoroughly discussed.

Aminoglycosides↗