Molecular Mechanics Calculations of Proteins. Comparison of Different Energy Minimization Strategies

Abstract

A general strategy for performing energy minimization of proteins using the SYBYL molecular modelling program has been developed. The influence of several variables including energy minimization procedure, solvation, dielectric function and dielectric constant have been investigated in order to develop a general method, which is capable of producing high quality protein structures. Avian pancreatic polypeptide (APP) and bovine pancreatic phospholipase A2 (BP PLA2) were selected for the calculations, because high quality X-ray structures exist and because all classes of secondary structure are represented in the structures. The energy minimized structures were evaluated relative to the corresponding X-ray structures. The overall similarity was checked by calculating RMS distances for all atom positions. Backbone conformation was checked by Ramachandran plots and secondary structure elements evaluated by the length on hydrogen bonds. The dimensions of active site in BP PLA2 is very dependent on electrostatic interactions, due to the presence of the positively charged calcium ion. Thus, the distances between calcium and the calcium-coordinating groups were used as a quality index for this protein. Energy minimized structures of the trimeric PLA2 from Indian cobra (N.n.n. PLA2) were used for assessing the impact of protein-protein interactions. Based on the above mentioned criteria, it could be concluded that using the following conditions: Dielectric constant ? = 4 or 20; a distance dependent dielectric function and stepwise energy minimization, it is possible to reproduce X-ray structures very accurately without including explicit solvent molecules.     

Lipid II-independent trans Editing of Mischarged tRNAs by the Penicillin Resistance Factor MurM

Streptococcus pneumoniae is a causative agent of nosocomial infections such as pneumonia, meningitis, and septicemia. Penicillin resistance in S. pneumoniae depends in part upon MurM, an aminoacyl-tRNA ligase that attaches l-serine or l-alanine to the stem peptide lysine of Lipid II in cell wall peptidoglycan. To investigate the exact substrates the translation machinery provides MurM, quality control by alanyl-tRNA synthetase (AlaRS) was investigated. AlaRS mischarged serine and glycine to tRNA^Ala, as observed in other bacteria, and also transferred alanine, serine, and glycine to tRNA^Phe. S. pneumoniae tRNA^Phe has an unusual U4:C69 mismatch in its acceptor stem that prevents editing by phenylalanyl-tRNA synthetase (PheRS), leading to the accumulation of misaminoacylated tRNAs that could serve as substrates for translation or for MurM. Although the peptidoglycan layer of S. pneumoniae tolerates a combination of both branched and linear muropeptides, deletion of MurM results in a reversion to penicillin sensitivity in strains that were previously resistant. However, because MurM is not required for cell viability, the reason for its functional conservation across all strains of S. pneumoniae has remained elusive. We now show that MurM can directly function in translation quality control by acting as a broad specificity lipid-independent trans editing factor that deacylates tRNA. This activity of MurM does not require the presence of its second substrate, Lipid II, and can functionally substitute for the activity of widely conserved editing domain homologues of AlaRS, termed AlaXPs proteins, which are themselves absent from S. pneumoniae.     

Distinct tRNA modifications in the thermo-acidophilic archaeon, Thermoplasma acidophilum

Thermoplasma acidophilum is a thermo-acidophilic archaeon. We purified tRNA^Leu (UAG) from T. acidophilum using a solid-phase DNA probe method and determined the RNA sequence after determining via nucleoside analysis and m⁷G-specific aniline cleavage because it has been reported that T. acidophilum tRNA contains m⁷G, which is generally not found in archaeal tRNAs. RNA sequencing and liquid chromatography–mass spectrometry revealed that the m⁷G modification exists at a novel position 49. Furthermore, we found several distinct modifications, which have not previously been found in archaeal tRNA, such as 4-thiouridine9, archaeosine13 and 5-carbamoylmethyuridine34. The related tRNA modification enzymes and their genes are discussed.     

A voltage-gated pore for translocation of tRNA

Very little is known about how nucleic acids are translocated across membranes. The multi-subunit RNA Import Complex (RIC) from mitochondria of the kinetoplastid protozoon Leishmania tropica induces translocation of tRNAs across artificial or natural membranes, but the nature of the translocation pore remains unknown. We show that subunits RIC6 and RIC9 assemble on the membrane in presence of subunit RIC4A to form complex R3. Atomic Force Microscopy of R3 revealed particles with an asymmetric surface groove of ∼20 nm rim diameter and ∼1 nm depth. R3 induced translocation of tRNA into liposomes when the pH of the medium was lowered to ∼6 in the absence of ATP. R3-mediated tRNA translocation could also be induced at neutral pH by a K⁺ diffusion potential with an optimum of 60–70 mV. Point mutations in the Cys₂–His₂ Fe-binding motif of RIC6, which is homologous to the respiratory Complex III Fe–S protein, abrogated import induced by low pH but not by K⁺ diffusion potential. These results indicate that the R3 complex forms a pore that is gated by a proton-generated membrane potential and that the Fe–S binding region of RIC6 has a role in proton translocation. The tRNA import complex of L. tropica thus contains a novel macromolecular channel distinct from the mitochondrial protein import pore that is apparently involved in tRNA import in some species.     

Asymmetric structure and domain binding interfaces of human tyrosyl‐tRNA synthetase studied by molecular dynamics simulations

Human tyrosyl‐tRNA synthetase (HsTyrRS) is composed of two structural modules: N‐terminal catalytic core and an EMAP II‐like C‐terminal domain. The structures of these modules are known, but no crystal structure of the full‐length HsTyrRS is currently available. An all‐atom model of the full‐length HsTyrRS was developed in this work. The structure, dynamics, and domain binding interfaces of HsTyrRS were investigated by extensive molecular dynamics (MD) simulations. Our data suggest that HsTyrRS in solution consists of a number of compact asymmetric conformations, which differ significantly by their rigidity, internal mobility, orientation of C‐terminal modules, and the strength of interdomain binding. Interfaces of domain binding obtained in MD simulations are in perfect agreement with our previous coarse‐grained hierarchical rotations technique simulations. Formation of the hydrogen bonds between R93 residue of the ELR cytokine motif and the residues A340 and E479 in the C‐module was observed. This observation supports the idea that the lack of cytokine activity in the full‐length HsTyrRS is explained by interactions between N‐modules and C‐modules, which block the ELR motif. Copyright © 2013 John Wiley & Sons, Ltd.     

Mechanisms of the tRNA wobble cytidine modification essential for AUA codon decoding in prokaryotes

Bacteria and archaea have 2-lysylcytidine (L or lysidine) and 2-agmatinylcytidine (agm²C or agmatidine), respectively, at the first (wobble) position of the anticodon of the AUA codon-specific tRNA^Ile. These lysine- or agmatine-conjugated cytidine derivatives are crucial for the precise decoding of the genetic code. L is synthesized by tRNA^Ile-lysidine synthetase (TilS), which uses l-lysine and ATP as substrates. Agm²C formation is catalyzed by tRNA^Ile-agm²C synthetase (TiaS), which uses agmatine and ATP for the reaction. Despite the fact that TilS and TiaS synthesize structurally similar cytidine derivatives, these enzymes belong to non-related protein families. Therefore, these enzymes modify the wobble cytidine by distinct catalytic mechanisms, in which TilS activates the C2 carbon of the wobble cytidine by adenylation, while TiaS activates it by phosphorylation. In contrast, TilS and TiaS share similar tRNA recognition mechanisms, in which the enzymes recognize the tRNA acceptor stem to discriminate tRNA^Ile and tRNA^Met.     

Escherichia coli tRNAs Are Resistant to the Hyperprocessing Reaction of Homologous E. coli Ribonuclease P Ribozyme

Bacterial ribonuclease P RNA ribozyme can do the hyperprocessing reaction, the internal cleavage reaction of some floppy eukaryotic tRNAs. The hyperprocessing reaction can be used as a detection tool to examine the stability of the cloverleaf shape of tRNA. Until now, the hyperprocessing reaction has been observed in the heterologous combination of eukaryotic tRNAs and bacterial RNase P enzymes. In this paper, we examined the hyperprocessing reaction of Escherichia coli tRNAs by homologous E. coli RNase P, to find that these homologous tRNAs were resistant to the toxic hyperprocessing reaction. Our results display the evidence for molecular co-evolution between homologous tRNAs and RNase P in the bacterium E. coli.