分子动力学模拟研究A1408G突变对核糖体16S rRNA片段的影响

利用分子动力学模拟方法, 研究了原核生物核糖体小亚基中的16S rRNA片段与氨基糖苷类抗生素巴龙霉素复合物结构的柔性. 结果表明, 16S rRNA片段中的1408位点的腺嘌呤(A)突变为鸟嘌呤(G), 改变了与tRNA中反密码子环识别相关的2个腺嘌呤A1492和A1493的空间构象, 阻碍了氨基糖苷类抗生素与核糖体的结合, 从而影响原核生物蛋白转录过程. 模拟结果与实验测定的晶体结构相吻合, 可为基于核糖体16S rRNA的药物分子设计提供较可靠的结构信息.     

Correlation between the conformation of nucleoside antibiotics and their biological activity

Nucleoside antibiotics, which result as a consequence of minor modifications in pyrimidine and purine nucleosides, exhibit a wide variety of antiviral, antibacterial, antitumor, and cancerostatic properties. The conformational properties of a number of these antibiotics have been investigated by using the quantummechanical PCILO method, and the results indicate that the nucleoside antibiotics and their parent nucleosides have very similar conformational preferences. This similarity is strikingly marked in the situations which prevail in an aqueous medium. As a result, these antibiotics easily get incorporated in growing chains of RNA and DNA by mimicking their parent nucleosides and then bring about the inhibition of protein, RNA, or DNA syntheses. The experimental observations corroborate these deductions, and thus a correlation has been obtained between the conformation and the biological activity of nucleoside antibiotics; it is the striking conformational similarity between the nucleoside antibiotics and their parent nucleosides which gives rise to their biological activity. The PCILO investigations carried out on two 3-deazapurine nucleosides demonstrate that the converse of the above correlation also holds true.     

Reparameterizations of the χ Torsion and Lennard‐Jones σ Parameters Improve the Conformational Characteristics of Modified Uridines

The currently available force field parameters for modified RNA residues in AMBER show significant deviations in conformational properties from experimental observations. The examination of the transferability of the recently revised torsion parameters revealed that there was an overall improvement in the conformational properties for some of the modifications but the improvements were still insufficient in describing the sugar pucker preferences (J. Chem. Inf. Model. 2014, 54, 1129–1142). Here, we report an approach for the development and fine tuning of the AMBER force field parameters for 2‐thiouridine, 4‐thiouridine, and pseudouridine with diverse conformational preferences. The χ torsion parameters were reparameterized at the individual nucleoside level. The effect of combining the revised γ torsion parameter and modifying the Lennard‐Jones σ parameters were also tested by directly comparing the conformational preferences obtained from our extensive molecular dynamics simulations with those from experimental observations. © 2016 Wiley Periodicals, Inc.     

Molecular properties from conformational ensembles. 1. Dipole moments of molecules with multiple internal rotations

We present a novel method for constructing the stable conformational space of small molecules with many rotatable bonds that uses our iterative stochastic elimination (ISE) algorithm, a robust stochastic search method capable of finding ensembles of best solutions for large combinatorial problems. To validate the method, we show that ISE reproduces the best conformers found in a fully exhaustive search, as well as compare computed dipole moments to experimental values, based on molecular ensembles and their Boltzmann distributions. Results were also compared to the alternative molecular dynamics and simulated annealing methods. Our results clarify that many low energy conformations may be required to reproduce molecular properties, while single low energy conformers or ensembles of low energy conformers cannot account for the experimental properties of flexible molecules. Whereas ISE well reproduces conformations that are not separated by very large energy barriers, it has not been successful in reproducing conformations of strained molecules.     

Energetic Tuning by tRNA Modifications Ensures Correct Decoding of Isoleucine and Methionine on the Ribosome

Chemical modifications of tRNAs are critical for accurate translation of the genetic code on the ribosome. The discrimination between isoleucine (AUA) and methionine (AUG) codons depends on such modifications of the wobble position in isoleucine tRNA anticodon loops, in all kingdoms of life. Bacteria and archaea employ functionally similar lysine‐ and agmatine‐conjugated cytidine derivatives to ensure decoding fidelity, but the thermodynamics underlying codon discrimination remains unknown. Here, we report structure‐based computer simulations that quantitatively reveal the energetics of this decoding strategy in archaea. The results further show that the agmatidine modification confers tRNA specificity primarily by desolvation of the incorrect codon in the non‐cognate complex. Tautomerism is found to play no significant role in this decoding system as the usual amino form of the modified tRNA is by far the most stable.     

Additive CHARMM force field for naturally occurring modified ribonucleotides

More than 100 naturally occurring modified nucleotides have been found in RNA molecules, in particular in tRNAs. We have determined molecular mechanics force field parameters compatible with the CHARMM36 all‐atom additive force field for all these modifications using the CHARMM force field parametrization strategy. Emphasis was placed on fine tuning of the partial atomic charges and torsion angle parameters. Quantum mechanics calculations on model compounds provided the initial set of target data, and extensive molecular dynamics simulations of nucleotides and oligonucleotides in aqueous solutions were used for further refinement against experimental data. The presented parameters will allow for computational studies of a wide range of RNAs containing modified nucleotides, including the ribosome and transfer RNAs. © 2016 The Authors. Journal of Computational Chemistry Published by Wiley Periodicals, Inc.     

A new set of molecular mechanics parameters for hydroxyproline and its use in molecular dynamics simulations of collagen-like peptides

Recently, the importance of proline ring pucker conformations in collagen has been suggested in the context of hydroxylation of prolines. The previous molecular mechanics parameters for hydroxyproline, however, do not reproduce the correct pucker preference. We have developed a new set of parameters that reproduces the correct pucker preference. Our molecular dynamics simulations of proline and hydroxyproline monomers as well as collagen-like peptides, using the new parameters, support the theory that the role of hydroxylation in collagen is to stabilize the triple helix by adjusting to the right pucker conformation (and thus the right phi angle) in the Y position.     

Conformational aspects and biological functions of biomolecules

For studying the molecular mechanism of the biological function of transfer ribonucleic acid (tRNA, essential cell component for translation), NMR spin coupling constants, lanthanide-induced shifts and relaxation rates, and nuclear Overhauser effects have been analyzed of uridine phosphates and modified uridine derivatives in aqueous solution. On the basis of the conformational characteristics of uridine derivatives, the molecular mechanism of enhanced thermostability of extreme thermophile tRNA species and the biological significance of the two types of modifications of uridine in the first position of anticodon have been elucidated.     

Conformational sensitivity of polyether macrocycles to electrostatic potential: Partial atomic charges,molecular mechanics,and conformational prediction

Molecular recognition (whether by enzymes, the immune system, or chelating ligands) depends critically on molecular conformation. Molecular mechanics predicts energetically favorable molecular conformations by locating low energy conformations using an empirical fit of molecular potential energy as a function of internal coordinates. Molecular mechanics analysis of 18-crown-6 demonstrates that the nonbonded term (primarily the electrostatic part) is the largest contributor to the conformational energy. Nevertheless, common methods of treating the electrostatic interaction for 18-crown-6 yield inconsistent values for conformational energies partly because partial charges assigned to each atom can change with conformation due to through-space inductive effects which are not considered in most molecular mechanics programs. Similar findings from several other groups are reviewed to support our conclusions. We argue for care and caution in predicting conformational preferences of molecules with two or more highly polar atoms. We also discuss the desirability of using an empirical method of partial charge determination such as the charge equilibration algorithm of Rappé and Goddard (or a suitable generalization which includes polarization) as a method of including these effects in molecular mechanics and molecular dynamics calculations.     

A Method to Explore Protein Side Chain Conformational Variability Using Experimental Data

Experimentally measured values of molecular properties or observables of biomolecules such as proteins are generally averages over time and space, which do not contain su?cient information to determine the underlying conformational distribution of the molecules in solution. The relationship between experimentally measured NMR ³J‐coupling values and the corresponding dihedral angle values is a particularly complicated case due to its nonlinear, multiple‐valued nature. Molecular dynamics (MD) simulations at constant temperature can generate Boltzmann ensembles of molecular structures that are free from a priori assumptions about the nature of the underlying conformational distribution. They suffer, however, from limited sampling with respect to time and conformational space. Moreover, the quality of the obtained structures is dependent on the choice of force ?eld and solvation model. A recently proposed method that uses time‐averaging with local‐elevation (LE) biasing of the conformational search provides an elegant means of overcoming these three problems. Using a set of side chain ³J‐coupling values for the FK506 binding protein (FKBP), we ?rst investigate the uncertainty in the angle values predicted theoretically. We then propose a simple MD‐based technique to detect inconsistencies within an experimental data set and identify degrees of freedom for which conformational averaging takes place or for which force ?eld parameters may be de?cient. Finally, we show that LE MD is the best method for producing ensembles of structures that, on average, ?t the experimental data.     

On-the-Fly Ring-Polymer Molecular Dynamics Calculations of the Dissociative Photodetachment Process of the Oxalate Anion

The dissociative photodetachment dynamics of the oxalate anion, C₂O₄H⁻ + hν → CO₂ + HOCO + e⁻, were theoretically studied using the on-the-fly path-integral and ring-polymer molecular dynamics methods, which can account for nuclear quantum effects at the density-functional theory level in order to compare with the recent experimental study using photoelectron–photofragment coincidence spectroscopy. To reduce computational time, the force acting on each bead of ring-polymer was approximately calculated from the first and second derivatives of the potential energy at the centroid position of the nuclei beads. We find that the calculated photoelectron spectrum qualitatively reproduces the experimental spectrum and that nuclear quantum effects are playing a role in determining spectral widths. The calculated coincidence spectrum is found to reasonably reproduce the experimental spectrum, indicating that a relatively large energy is partitioned into the relative kinetic energy between the CO₂ and HOCO fragments. This is because photodetachment of the parent anion leads to Franck–Condon transition to the repulsive region of the neutral potential energy surface. We also find that the dissociation dynamics are slightly different between the two isomers of the C₂O₄H⁻ anion with closed- and open-form structures.     

Nonfluorescent quenchers to correlate single-molecule conformational and compositional dynamics

Single-molecule F?rster resonance energy transfer (smFRET) is a powerful method for studying the conformational dynamics of a biomolecule in real-time. However, studying how interacting ligands correlate with and regulate the conformational dynamics of the biomolecule is extremely challenging because of the availability of a limited number of fluorescent dyes with both high quantum yield and minimal spectral overlap. Here we report the use of a nonfluorescent quencher (Black Hole Quencher, BHQ) as an acceptor for smFRET. Using a Cy3/BHQ pair, we can accurately follow conformational changes of the ribosome during elongation in real time. We demonstrate the application of single-color FRET to correlate the conformational dynamics of the ribosome with the compositional dynamics of tRNA. We use the normal Cy5 FRET acceptor to observe arrival of a fluorescently labeled tRNA with a concomitant transition of the ribosome from the locked to the unlocked conformation. Our results illustrate the potential of nonfluorescent quenchers in single-molecule correlation studies.     

Conformation of nucleoside antibiotics

Minor modifications or substitutions in the sugar or in the base part of pyrimidine and purine nucleosides have a profound effect on their biological activity. These modified nucleosides usually become antiviral, antibacterial, or cancerostatic agents and they are collectively called nucleoside antibiotics. The conformational properties of some of these nucleoside antibiotics have been studied by the PCILO method. The results obtained from such study indicate that the conformational preferences of these nucleoside antibiotics are very similar to those of their parent nucleosides and especially so in the situations that occur in aqueous solutions. The important biological significance of these results is that these nucleoside antibiotics can easily get incorporated into growing chains of DNA and RNA by mimicking their parent nucleosides and can interfere with the protein synthesis of RNA or DNA synthesis.     

Force field parameters for the simulation of modified histone tails

We describe the development of force field parameters for methylated lysines and arginines, and acetylated lysine for the CHARMM all‐atom force field. We also describe a CHARMM united‐atom force field for modified sidechains suitable for use with fragment‐based docking methods. The development of these parameters is based on results of ab initio quantum mechanics calculations of model compounds with subsequent refinement and validation by molecular mechanics and molecular dynamics simulations. The united‐atom parameters are tested by fragment docking to target proteins using the MCSS procedure. The all‐atom force field is validated by molecular dynamics simulations of multiple experimental structures. In both sets of calculations, the computational predictions using the force field were compared to the corresponding experimental structures. We show that the parameters yield an accurate reproduction of experimental structures. Together with the existing CHARMM force field, these parameters will enable the general modeling of post‐translational modifications of histone tails. © 2010 Wiley Periodicals, Inc. J Comput Chem, 2010     

Molecular mechanics models for tetracycline analogs

Tetracyclines (Tcs) are an important family of antibiotics that bind to the ribosome and several proteins. To model Tc interactions with protein and RNA, we have developed a molecular mechanics force field for 12 tetracyclines, consistent with the CHARMM force field. We considered each Tc variant in its zwitterionic tautomer, with and without a bound Mg(2+). We used structures from the Cambridge Crystallographic Data Base to identify the conformations likely to be present in solution and in biomolecular complexes. A conformational search by simulated annealing was undertaken, using the MM3 force field, for tetracycline, anhydrotetracycline, doxycycline, and tigecycline. Resulting, low-energy structures were optimized with an ab initio method. We found that Tc and its analogs all adopt an extended conformation in the zwitterionic tautomer and a twisted one in the neutral tautomer, and the zwitterionic-extended state is the most stable in solution. Intermolecular force field parameters were derived from a standard supermolecule approach: we considered the ab initio energies and geometries of a water molecule interacting with each Tc analog at several different positions. The final, rms deviation between the ab initio and force field energies, averaged over all forms, was 0.35 kcal/mol. Intramolecular parameters were adopted from either the standard CHARMM force field, the ab initio structure, or the earlier, plain Tc force field. The model reproduces the ab initio geometry and flexibility of each Tc. As tests, we describe MD and free energy simulations of a solvated complex between three Tcs and the Tet repressor protein.     

PEPCAT—A new tool for conformational analysis of peptides

We report a new technique for the efficient analysis and visualization of peptide and protein conformations and conformational relationships, which we have implemented in a computer program called PEPCAT. PEPCAT (an abbreviation for Peptide Conformational Analysis Tool) provides a simple, graphical, and flexible framework that allows the user to define a specific structural feature or juxtaposition of amino acids and to follow the fate of the motif during a molecular dynamics simulation. Here we describe the PEPCAT analysis of the effects of environmental and chemical modifications on conformational preferences of a regulator of hemopoiesis, namely the pentapeptide pyro‐EEDCK, and of a conformational transition in the immunosuppressant drug cyclosporin A. PEPCAT, however, can be applied to the conformational analysis of peptides and proteins in general. © 2000 John Wiley & Sons, Inc. J Comput Chem 21: 446–461, 2000     

Aminoglycoside-induced reduction in nucleotide mobility at the ribosomal RNA A-site as a potentially key determinant of antibacterial activity

Steady-state and time-resolved fluorescence techniques have been used to characterize the energetics and dynamics associated with the interaction of an E. coli 16 S rRNA A-site model oligonucleotide and four aminoglycoside antibiotics that exhibit a broad range of antibacterial activity. The results of these characterizations suggest that aminoglycoside-induced reduction in the mobility of an adenine residue at position 1492 of the rRNA A-site is a more important determinant of antibacterial activity than drug affinity for the A-site. This observation is consistent with a recently proposed model for the mechanism of protein synthesis inhibition by aminoglycosides that invokes a drug-induced alteration in the conformational equilibrium of the rRNA A-site (centered around the conserved adenine residues at positions 1492 and 1493), which, in turn, promotes an enhanced interaction between the rRNA and the minihelix formed by the tRNA anticodon and the mRNA codon, even when the anticodon is noncognate. Regarded as a whole, the results reported here indicate that the rational design of antibiotics that target the 16 S rRNA A-site requires consideration of not only the structure and energetics of the drug-RNA complex but also the dynamics associated with that complex.     

HMGB1–Carbenoxolone Interactions: Dynamics Insights from Combined Nuclear Magnetic Resonance and Molecular Dynamics

The interplay of protein dynamics and molecular recognition is of fundamental importance in biological processes. Atomic‐resolution insights into these phenomena may provide new opportunities for drug discovery. Herein, we have combined NMR relaxation experiments and residual dipolar coupling (RDC) measurements with molecular dynamics (MD) simulations to study the effects of the anti‐inflammatory drug carbenoxolone (CBNX) on the conformational properties and on the internal dynamics of a subdomain (box A) of high‐mobility group B protein (HMGB1). ¹⁵N relaxation data show that CBNX binding enhances the fast pico‐ to nanosecond motions of a loop and partially removes the internal motional anisotropy of the first two helices of box A. Dipolar wave analysis of amide RDC data shows that ligand binding induces helical distortions. In parallel, increased mobility of the loop upon ligand binding is highlighted by the essential dynamics analysis (EDA) of MD simulations. Moreover, simulations detect two possible orientations for CBNX, which induces two possible conformations of helix H3, one being similar to the free form and the second one causing a partial helical distortion. Finally, we introduce a new approach for the analysis of the internal coordination of protein residues that is consistent with experimental data and allows us to pinpoint which substructures of box A are dynamically affected by CBNX. The observations reported here may be useful for understanding the role of protein dynamics in binding at atomic resolution.