Search

The role of ribonucleic acid (RNA) in biology continues to grow, but insight into important aspects of RNA behavior is lacking, such as dynamic structural ensembles in different environments, how flexibility is coupled to function, and how function might be modulated by small molecule binding. In the case of proteins, much progress in these areas has been made by complementing experiments with atomistic simulations, but RNA simulation methods and force fields are less mature. It remains challenging to generate stable RNA simulations, even for small systems where well-defined, thermostable structures have been established by experiments. Many different aspects of RNA energetics have been adjusted in force fields, seeking improvements that are transferable across a variety of RNA structural motifs. In this work, the role of weak CH···O interactions is explored, which are ubiquitous in RNA structure but have received less attention in RNA force field development. By comparing data extracted from high-resolution RNA crystal structures to energy profiles from quantum mechanics and force field calculations, it is shown that CH···O interactions are overly repulsive in the widely used Amber RNA force fields. A simple, targeted adjustment of CH···O repulsion that leaves the remainder of the force field unchanged was developed. Then, the standard and modified force fields were tested using molecular dynamics (MD) simulations with explicit water and salt, amassing over 300 μs of data for multiple RNA systems containing important features such as the presence of loops, base stacking interactions as well as canonical and noncanonical base pairing. In this work and others, standard force fields lead to reproducible unfolding of the NMR-based structures. Including a targeted CH···O adjustment in an otherwise identical protocol dramatically improves the outcome, leading to stable simulations for all RNA systems tested.

Ribonucleic acid (RNA) molecules can adopt a variety of secondary and tertiary structures in solution, with stem-loops being one of the more common motifs. Here, we present a systematic analysis of 15 RNA stem-loop sequences simulated with molecular dynamics simulations in an implicit solvent environment. Analysis of RNA cluster ensembles showed that the stem-loop structures can generally adopt the A-form RNA in the stem region. Loop structures are more sensitive, and experimental structures could only be reproduced with modification of CH···O interactions in the force field, combined with an implicit solvent nonpolar correction to better model base stacking interactions. Accurately modeling RNA with current atomistic physics-based models remains challenging, but the RNA systems studied herein may provide a useful benchmark set for testing other RNA modeling methods in the future.

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.

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.

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 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.

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.

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.

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.

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.

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.

Transfer of alanine from Escherichia coli alanyl-tRNA synthetase (AlaRS) to RNA minihelices that mimic the amino acid acceptor stem of tRNA(Ala) has been shown, by analysis of variant minihelix aminoacylation activities, to involve a transition state sensitive to changes in the 'discriminator' base at position 73. Solution NMR has indicated that this single-stranded nucleotide is predominantly stacked onto G1 of the first base pair of the alanine acceptor stem helix. We report the activity of a new variant with the adenine at position 73 substituted by its non-polar isostere 4-methylindole (M). Despite lacking N7, this analog is well tolerated by AlaRS. Molecular dynamics (MD) simulations show that the M substitution improves position 73 base stacking over G1, as measured by a stacking lifetime analysis. Additional MD simulations of wild-type microhelix(Ala) and six variants reveal a positive correlation between N73 base stacking propensity over G1 and aminoacylation activity. For the two DeltaN7 variants simulated we found that the propensity to stack over G1 was similar to the analogous variants that contain N7 and we conclude that the decrease in aminoacylation efficiency observed upon deletion of N7 is likely due to loss of a direct stabilizing interaction with the synthetase.