AmberTools.

AmberTools is a free and open-source collection of programs used to set up, run, and analyze molecular simulations. The newer features contained within AmberTools23 are briefly described in this Application note.

Phosphomimetic Mutation at Ser165 of α-Tubulin Promotes the Persistence of GTP Caps in Microtubules.

Protein kinase C (PKC)-mediated phosphorylation of α-tubulin at Ser165 or expression of phosphomimetic (S165D)-α-tubulin stimulates microtubule (MT) polymerization (Cytoskeleton 2014, 71, 257-272). Ser165 lies near the interface between adjacent αß-tubulin heterodimers and helix H8, which contains Glu254, the catalytic residue in α-tubulin that hydrolyzes the exchangeable GTP in ß-tubulin (ß:GTP) and triggers MT depolymerization. It was hypothesized that S165D, a phosphomimetic variant of α-tubulin, perturbs the alignment of α:Glu254 with respect to ß:GTP, thereby impairing its hydrolysis. Molecular simulations were performed with cryoEM structures of MTs (PDB ID: 3J6E) in which phosphomimetic S165D or control S165N had been substituted. Unlike native and S165N structures, the distance between S165D and α:Glu254 increased by 0.6 Å, while the distance between α:Glu254 and ß:GTP decreased by 0.4 Å. Rotation of ß:GTP by 4 Å occurred in the S165D variant, whereas ß:GTP in the S165N control was unchanged from the native structure. Additionally, the S165D variant exhibited an altered pattern of H-bonding to ß:GTP, including the loss of three H-bonds. The significance of these findings to ß:GTP hydrolysis was analyzed in MCF-10A human breast cells treated with an antibody that detects GTP-bound tubulin. Compared with controls, GTP-tubulin signals were at higher levels in cells that ectopically expressed S165D-α-tubulin (TUBA1C) or had been treated with PKC activator DAG-lactone to induce phosphorylation of Ser165 in native α-tubulin. These findings support a model whereby conformational changes induced by Ser165 phosphorylation alter the spatial relationship between ß:GTP and α:Glu254, thereby slowing GTP hydrolysis and promoting GTP caps.

The Impacts of the Molecular Education and Research Consortium in Undergraduate Computational Chemistry on the Careers of Women in Computational Chemistry.

The Molecular Education and Research Consortium in Undergraduate Computational Chemistry (MERCURY) has supported a diverse group of faculty and students for over 20 years by providing computational resources as well as networking opportunities and professional support. The consortium comprises 38 faculty (42% women) at 34 different institutions, who have trained nearly 900 undergraduate students, more than two-thirds of whom identify as women and one-quarter identify as students of color. MERCURY provides a model for the support necessary for faculty to achieve professional advancement and career satisfaction. The range of experiences and expertise of the consortium members provides excellent networking opportunities that allow MERCURY faculty to support each other's teaching, research, and service needs, including generating meaningful scientific advancements and outcomes with undergraduate researchers as well as being leaders at the departmental, institutional, and national levels. While all MERCURY faculty benefit from these supports, the disproportionate number of women in the consortium, relative to their representation in computational sciences generally, produces a sizable impact on advancing women in the computational sciences. In this report, the women of MERCURY share how the consortium has benefited their careers and the careers of their students.

Women Make COMP: Mentoring the Next Generation of Women in Computational Chemistry.

The Women Make COMP symposium (258th American Chemical Society Meeting) aims at inspiring, motivating, and supporting young women in computational and theoretical chemistry. As a role model of the event, Ada Lovelace (1815-1852) was an English mathematician and writer, known for having founded computing science.

Molecular dynamics simulations of human tRNA Lys,3 UUU: the role of modified bases in mRNA recognition.

Accuracy in translation of the genetic code into proteins depends upon correct tRNA-mRNA recognition in the context of the ribosome. In human tRNA(Lys,3)UUU three modified bases are present in the anticodon stem-loop--2-methylthio-N6-threonylcarbamoyladenosine at position 37 (ms2t6A37), 5-methoxycarbonylmethyl-2-thiouridine at position 34 (mcm5s2U34) and pseudouridine (psi) at position 39--two of which, ms2t6A37 and mcm5s2U34, are required to achieve wild-type binding activity of wild-type human tRNA(Lys,3)UUU [C. Yarian, M. Marszalek, E. Sochacka, A. Malkiewicz, R. Guenther, A. Miskiewicz and P. F. Agris (2000) Biochemistry, 39, 13390-13395]. Molecular dynamics simulations of nine tRNA anticodon stem-loops with different combinations of nonstandard bases were performed. The wild-type simulation exhibited a canonical anticodon stair-stepped conformation. The ms2t6 modification at position 37 is required for maintenance of this structure and reduces solvent accessibility of U36. Ms2t6A37 generally hydrogen bonds across the loop and may prevent U36 from rotating into solution. A water molecule does coordinate to psi39 most of the simulation time but weakly, as most of the residence lifetimes are <40 ps.

Backbone-base interactions critical to quantum stabilization of transfer RNA anticodon structure.

Transfer RNA (tRNA) anticodons adopt a highly ordered 3'-stack without significant base overlap. Density functional theory at the M06-2X/6-31+G(d,p) level in combination with natural bond orbital analysis was utilized to calculate the intramolecular interactions within the tRNA anticodon that are responsible for stabilizing the stair-stepped conformation. Ten tRNA X-ray crystal structures were obtained from the PDB databank and were trimmed to include only the anticodon bases. Hydrogenic positions were added and optimized for the structures in the stair-stepped conformation. The sugar-phosphate backbone has been retained for these calculations, revealing the role it plays in RNA structural stability. It was found that electrostatic interactions between the sugar-phosphate backbone and the base provide the most stability, rather than the traditionally studied interbase stacking. Base-stacking interactions, though present, were weak and inconsistent. Aqueous solvation was found to have little effect on the intramolecular interactions.

Water, shape recognition, salt bridges, and cation-pi interactions differentiate peptide recognition of the HIV rev-responsive element.

Recognition of the human immunodeficiency virus Rev-responsive element (RRE) RNA by the Rev protein is an essential step in the viral life cycle. Formation of the Rev-RRE complex signals nucleocytoplasmic export of unspliced and partially spliced viral RNA. Essential components of the complex have been localized to a minimal arginine-rich Rev peptide and stem IIB of RRE. In vitro selection studies have identified a synthetic peptide known as RSG 1.2 that binds with better specificity and affinity to RRE than the Rev peptide. NMR structures of both peptide-RNA complexes of Rev and RSG 1.2 bound to RRE stem IIB have been solved and reveal gross structural differences between the two bound complexes. Molecular dynamics simulations of the Rev and RSG 1.2 peptides in complex with RRE stem IIB have been simulated to better understand on an atomic level how two arginine-rich peptides of similar length recognize the same sequence of RNA with such different structural motifs. While the Rev peptide employs some base-specific hydrogen bonding for recognition of RRE, shape recognition, through contact with the sugar-phosphate backbone, and cation-pi interactions are also important. Molecular dynamics simulations suggest that RSG 1.2 binds more tightly to the RRE sequence than Rev by forming more base-specific contacts, using water to mediate peptide-RNA contacts, and is held in place by a strong salt bridge network spanning the major groove of the RNA.

Pre-transfer editing by class II prolyl-tRNA synthetase: role of aminoacylation active site in "selective release" of noncognate amino acids.

Aminoacyl-tRNA synthetases catalyze the attachment of cognate amino acids to specific tRNA molecules. To prevent potential errors in protein synthesis caused by misactivation of noncognate amino acids, some synthetases have evolved editing mechanisms to hydrolyze misactivated amino acids (pre-transfer editing) or misacylated tRNAs (post-transfer editing). In the case of post-transfer editing, synthetases employ a separate editing domain that is distinct from the site of amino acid activation, and the mechanism is believed to involve shuttling of the flexible CCA-3' end of the tRNA from the synthetic active site to the site of hydrolysis. The mechanism of pre-transfer editing is less well understood, and in most cases, the exact site of pre-transfer editing has not been conclusively identified. Here, we probe the pre-transfer editing activity of class II prolyl-tRNA synthetases from five species representing all three kingdoms of life. To locate the site of pre-transfer editing, truncation mutants were constructed by deleting the insertion domain characteristic of bacterial prolyl-tRNA synthetase species, which is the site of post-transfer editing, or the N- or C-terminal extension domains of eukaryotic and archaeal enzymes. In addition, the pre-transfer editing mechanism of Escherichia coli prolyl-tRNA synthetase was probed in detail. These studies show that a separate editing domain is not required for pre-transfer editing by prolyl-tRNA synthetase. The aminoacylation active site plays a significant role in preserving the fidelity of translation by acting as a filter that selectively releases non-cognate adenylates into solution, while protecting the cognate adenylate from hydrolysis.

Efficient aminoacylation of the tRNA(Ala) acceptor stem: dependence on the 2:71 base pair.

Specific aminoacylation by aminoacyl-tRNA synthetases requires accurate recognition of cognate tRNA substrates. In the case of alanyl-tRNA synthetase (AlaRS), RNA duplexes that mimic the acceptor stem of the tRNA are efficient substrates for aminoacylation in vitro. It was previously shown that recognition by AlaRS is severely affected by a simple base pair transversion of the G2:C71 pair at the second position in the RNA helix. In this study, we determined the aminoacylation efficiencies of 50 variants of the tRNA(Ala) acceptor stem containing substitutions at the 2:71 position. We find that there is not a single functional group of the wild-type G2:C71 base pair that is critical for positive recognition. Rather, we observed that base-pair orientation plays an important role in recognition. In particular, pyrimidine2:purine71 combinations generally resulted in decreased aminoacylation efficiency compared to the corresponding purine:pyrimidine pair. Moreover, the activity of a pyrimidine:purine variant could be partially restored by the presence of a major groove amino group at position 71. In an attempt to understand this result further, dielectric continuum electrostatic calculations were carried out, in some cases with additional inclusion of van der Waals interaction energies, to determine interaction potentials of the wild-type duplexAla and seven 2:71 variants. This analysis revealed a positive correlation between major groove negative electrostatic potential in the vicinity of the 3:70 base pair and measured aminoacylation efficiency.