RESUMO
Community network analysis derived from molecular dynamics simulations is used to identify and compare the signaling pathways in a bacterial glutamyl-tRNA synthetase (GluRS):tRNA(Glu) and an archaeal leucyl-tRNA synthetase (LeuRS):tRNA(Leu) complex. Although the 2 class I synthetases have remarkably different interactions with their cognate tRNAs, the allosteric networks for charging tRNA with the correct amino acid display considerable similarities. A dynamic contact map defines the edges connecting nodes (amino acids and nucleotides) in the physical network whose overall topology is presented as a network of communities, local substructures that are highly intraconnected, but loosely interconnected. Whereas nodes within a single community can communicate through many alternate pathways, the communication between monomers in different communities has to take place through a smaller number of critical edges or interactions. Consistent with this analysis, there are a large number of suboptimal paths that can be used for communication between the identity elements on the tRNAs and the catalytic site in the aaRS:tRNA complexes. Residues and nucleotides in the majority of pathways for intercommunity signal transmission are evolutionarily conserved and are predicted to be important for allosteric signaling. The same monomers are also found in a majority of the suboptimal paths. Modifying these residues or nucleotides has a large effect on the communication pathways in the protein:RNA complex consistent with kinetic data.
Assuntos
Archaea/enzimologia , Bactérias/enzimologia , Glutamato-tRNA Ligase/metabolismo , Leucina-tRNA Ligase/metabolismo , Aminoacil-RNA de Transferência/metabolismo , Simulação por Computador , Estrutura Secundária de Proteína , Aminoacilação de RNA de TransferênciaRESUMO
Elongation factor Tu (EF-Tu) binds to all standard aminoacyl transfer RNAs (aa-tRNAs) and transports them to the ribosome while protecting the ester linkage between the tRNA and its cognate amino acid. We use molecular dynamics simulations to investigate the dynamics of the EF-Tu.guanosine 5'-triphosphate.aa-tRNA(Cys) complex and the roles played by Mg2+ ions and modified nucleosides on the free energy of protein.RNA binding. Individual modified nucleosides have pronounced effects on the structural dynamics of tRNA and the EF-Tu.Cys-tRNA(Cys) interface. Combined energetic and evolutionary analyses identify the coevolution of residues in EF-Tu and aa-tRNAs at the binding interface. Highly conserved EF-Tu residues are responsible for both attracting aa-tRNAs as well as providing nearby nonbonded repulsive energies that help fine-tune molecular attraction at the binding interface. In addition to the 3' CCA end, highly conserved tRNA nucleotides G1, G52, G53, and U54 contribute significantly to EF-Tu binding energies. Modification of U54 to thymine affects the structure of the tRNA common loop resulting in a change in binding interface contacts. In addition, other nucleotides, conserved within certain tRNA specificities, may be responsible for tuning aa-tRNA binding to EF-Tu. The trend in EF-Tu.Cys-tRNA(Cys) binding energies observed as the result of mutating the tRNA agrees with experimental observation. We also predict variations in binding free energies upon misacylation of tRNA(Cys) with d-cysteine or O-phosphoserine and upon changing the protonation state of l-cysteine. Principal components analysis in each case reveals changes in the communication network across the protein.tRNA interface and is the basis for the entropy calculations.