Molecular dynamics and thermodynamics of protein-RNA interactions: mutation of a conserved aromatic residue modifies stacking interactions and structural adaptation in the U1A-stem loop 2 RNA complex

Design, Synthesis, and Characterization of a Novel 2'-5'-Linked Amikacin-Binding Aptamer: An Experimental and MD Simulation Study

To extend the approach of using RNA aptamers as transient protective groups for the synthesis of novel small-molecule drug derivatives from the existing aminoglycosides, we incorporated 2'-5' phosphodiester backbone modification in a structurally known neomycin RNA aptamer and studied the binding of a series of aminoglycosides using isothermal calorimetry (ITC) and molecular dynamics (MD) simulation. Experimental characterization of amikacin, a commercially available and widely used aminoglycoside for treating bacterial infections, shows that the aptamer A1 with a 2'-5' linkage between G15 and U16 exhibits a sevenfold increase in binding affinity with a lower binding energy compared to the native aptamer. Molecular dynamics (MD) simulation studies rationalize that this noncanonical linkage generates a narrower binding pocket by creating a superspiral RNA helical structure, which improves the ligand's fit in the binding pocket. These results provide new insights into applying 2'-5' linkages to diversify functional RNA aptamers as noncovalent protective groups in the synthesis of aminoglycoside derivatives, which can be further extended to other current drug molecules and complex natural compounds to make new pools of drug candidates more efficiently.

Pharmacological evaluation of newly synthesized organotin IV complex for antiulcer potential

The present study aims to investigate the newly synthesized organotin (IV) complex (2E, 2′E) dibutylstannanediyl bis (4-(4-nitrophenyl) amino)-4-oxobut-2-enoate (DTN) for its anti-ulcer potential. Characterization performed by carbon nuclear magnetic resonance spectroscopy proved that all values are in the expected ranges of the new compound. Gastroprotective activity of DTN was evaluated through in-silico, anti-H. pylori, in-vitro, in-vivo, and ex-vivo proteomic analysis. In-silico analysis shows that DTN possess stable binding with protein targets involved in gastric ulcer pathophysiology. DTN exhibited an inhibitory effect against 2,2-diphenyl-1-picrylhydrazyl, H. pylori and hydrogen potassium ATPase (H⁺/K⁺-ATPase). The antiulcer activity was performed using an ethanol-induced gastric ulcer model in rats. Anti-oxidant profile of DTN showed a significant increase in glutathione-S-transferase, glutathione and catalase levels whereas lipid peroxidation levels were reduced. Histopathological findings confirmed that DTN protected the gastric mucosa of rats. Inflammatory markers tumor necrosis factor-alpha, nuclear factor kappa B, cyclooxygenase-2, interleukin 6 and interleukin-1β were reduced and prostaglandin-E₂ restored expression of these cytokines in DTN pretreated animals when analyzed by using immunohistochemistry, enzyme-linked immunosorbent assay and western blot techniques. In real-time polymerase chain reaction technique, the expression of H⁺/K⁺-ATPase was downregulated in DTN pretreated group. DTN did not cause any mortality up to 400 mg/Kg. This study indicates that the newly synthesized compound DTN, possess stable binding against selected targets. DTN exhibits a gastro-protective effect, mediated via anti-H. pylori, H⁺/K⁺-ATPase inhibition, anti-oxidant and anti-inflammatory pathways, exploring its therapeutic potential in gastric ulcer management.

In Vitro Evaluation of the Anti-Diabetic Potential of Aqueous Acetone Helichrysum petiolare Extract (AAHPE) with Molecular Docking Relevance in Diabetes Mellitus

Diabetes mellitus (DM) is a chronic metabolic condition that can lead to significant complications and a high fatality rate worldwide. Efforts are ramping up to find and develop novel α-glucosidase and α-amylase inhibitors that are both effective and potentially safe. Traditional methodologies are being replaced with new techniques that are less complicated and less time demanding; yet, both the experimental and computational strategies are viable and complementary in drug discovery and development. As a result, this study was conducted to investigate the in vitro anti-diabetic potential of aqueous acetone Helichrysum petiolare and B.L Burtt extract (AAHPE) using a 2-NBDG, 2-(N-(7-Nitrobenz-2-oxa-1,3-diazol-4-yl) amino)-2-deoxy-d-glucose uptake assay. In addition, we performed molecular docking of the flavonoid constituents identified and quantified by liquid chromatography-mass spectrometry (LC-MS) from AAHPE with the potential to serve as effective and safe α-amylase and α-glucosidase inhibitors, which are important in drug discovery and development. The results showed that AAHPE is a potential inhibitor of both α-amylase and α-glucosidase, with IC₅₀ values of 46.50 ± 6.17 (µg/mL) and 37.81 ± 5.15 (µg/mL), respectively. This is demonstrated by a significant increase in the glucose uptake activity percentage in a concentration-dependent manner compared to the control, with the highest AAHPE concentration of 75 µg/mL of glucose uptake activity being higher than metformin, a standard anti-diabetic drug, in the insulin-resistant HepG2 cell line. The molecular docking results displayed that the constituents strongly bind α-amylase and α-glucosidase while achieving better binding affinities that ranged from ΔG = -7.2 to -9.6 kcal/mol (compared with acarbose ΔG = -6.1 kcal/mol) for α-amylase, and ΔG = -7.3 to -9.0 kcal/mol (compared with acarbose ΔG = -6.3 kcal/mol) for α-glucosidase. This study revealed the potential use of the H. petiolare plant extract and its phytochemicals, which could be explored to develop potent and safe α-amylase and α-glucosidase inhibitors to treat postprandial glycemic levels in diabetic patients.

Role of salt-bridging interactions in recognition of viral RNA by arginine-rich peptides

Interactions between RNA molecules and proteins are critical to many cellular processes and are implicated in various diseases. The RNA-peptide complexes are good model systems to probe the recognition mechanism of RNA by proteins. In this work, we report studies on the binding-unbinding process of a helical peptide from a viral RNA element using nonequilibrium molecular dynamics simulations. We explored the existence of various dissociation pathways with distinct free-energy profiles that reveal metastable states and distinct barriers to peptide dissociation. We also report the free-energy differences for each of the four pathways to be 96.47 ± 12.63, 96.1 ± 10.95, 91.83 ± 9.81, and 92 ± 11.32 kcal/mol. Based on the free-energy analysis, we further propose the preferred pathway and the mechanism of peptide dissociation. The preferred pathway is characterized by the formation of sequential hydrogen-bonding and salt-bridging interactions between several key arginine amino acids and the viral RNA nucleotides. Specifically, we identified one arginine amino acid (R8) of the peptide to play a significant role in the recognition mechanism of the peptide by the viral RNA molecule.

Estrogenic activity of multicyclic aromatic hydrocarbons in rainbow trout (Oncorhynchus mykiss) in vitro assays

A representative group of multicyclic aromatic hydrocarbons (MAHC) which can be further classified as bridged-ring (bridged-MAHC) or fused-ring (fused-MAHC) were examined for their ability to interact with the estrogen receptor of rainbow trout (rtER) in a hepatic cytosolic estrogen receptor competitive binding assay (cyto rtERαβ) and the vitellogenin (Vtg) mRNA gene activation liver slice assay. All five fused-MAHCs; naphthalene (NAFT), fluorene (FE), Fluoranthene (FAT), pyrene (PY), and 9,10-dihydroanthracene (DAC) had no estrogenic activity in the in vitro assays used. Five of the eight bridged-MAHCs; triphenylethylene (3PE), o-terphenyl (OTP), triphenylmethane (TPM), 1,1-diphenylethylene (DPE), and cis-stilbene (CSB) were positive in the rtER-binding assay. The additional three bridged-MAHC's; trans-stilbene (TSB), tetraphenylethylene (4PE), and 4,4-di-tertbutylphenyl (DtBB) were determined to be non-binders due to isomeric configuration, solubility limitation, and possible steric hinderance. It is possible that the bridged-MAHCs bind to the rtER through a proposed aromatic-aromatic stacking (π-π interaction) facilitated by perpendicular ring orientation achieved through free rotation of the bridged rings. The fused-ring structures are locked in a planar configuration which doesn't allow for rotation of rings perpendicular to one another. This first report of the rtER-binding of bridged-MAHCs in fish demonstrates binding for a class of chemicals normally not thought of as having an affinity for the estrogen receptor and further supports the versatility or promiscuity of ER ligand selectivity.

Molecular basis for the increased affinity of an RNA recognition motif with re-engineered specificity: A molecular dynamics and enhanced sampling simulations study

The RNA recognition motif (RRM) is the most common RNA binding domain across eukaryotic proteins. It is therefore of great value to engineer its specificity to target RNAs of arbitrary sequence. This was recently achieved for the RRM in Rbfox protein, where four mutations R118D, E147R, N151S, and E152T were designed to target the precursor to the oncogenic miRNA 21. Here, we used a variety of molecular dynamics-based approaches to predict specific interactions at the binding interface. Overall, we have run approximately 50 microseconds of enhanced sampling and plain molecular dynamics simulations on the engineered complex as well as on the wild-type Rbfox·pre-miRNA 20b from which the mutated systems were designed. Comparison with the available NMR data on the wild type molecules (protein, RNA, and their complex) served to establish the accuracy of the calculations.

Free energy calculations suggest that further improvements in affinity and selectivity are achieved by the S151T replacement.

RNA is an outstanding target for oncological intervention. Engineering the most common RNA binding motif in human proteins (called RRM) so as to bind to a specific RNA has an enormous pharmacological potential. Yet, it is highly non trivial to design RRM-bearing protein variants with RNA selectivity and affinity sufficiently high for clinical applications. Here we present an extensive molecular simulation study which shed light on the exquisite molecular recognition of the empirically-engineered complex between the RRM-bearing protein Rbfox and its RNA target pre-miR21. The simulations allow predicting a variant, the S151T, which may lead to further enhancement of selectivity and affinity for pre-miR21.

RNA Structural Dynamics As Captured by Molecular Simulations: A Comprehensive Overview

With both catalytic and genetic functions, ribonucleic acid (RNA) is perhaps the most pluripotent chemical species in molecular biology, and its functions are intimately linked to its structure and dynamics. Computer simulations, and in particular atomistic molecular dynamics (MD), allow structural dynamics of biomolecular systems to be investigated with unprecedented temporal and spatial resolution. We here provide a comprehensive overview of the fast-developing field of MD simulations of RNA molecules. We begin with an in-depth, evaluatory coverage of the most fundamental methodological challenges that set the basis for the future development of the field, in particular, the current developments and inherent physical limitations of the atomistic force fields and the recent advances in a broad spectrum of enhanced sampling methods. We also survey the closely related field of coarse-grained modeling of RNA systems. After dealing with the methodological aspects, we provide an exhaustive overview of the available RNA simulation literature, ranging from studies of the smallest RNA oligonucleotides to investigations of the entire ribosome. Our review encompasses tetranucleotides, tetraloops, a number of small RNA motifs, A-helix RNA, kissing-loop complexes, the TAR RNA element, the decoding center and other important regions of the ribosome, as well as assorted others systems. Extended sections are devoted to RNA-ion interactions, ribozymes, riboswitches, and protein/RNA complexes. Our overview is written for as broad of an audience as possible, aiming to provide a much-needed interdisciplinary bridge between computation and experiment, together with a perspective on the future of the field.

Computational and Pharmacological Evaluation of Ferrocene-Based Acyl Ureas and Homoleptic Cadmium Carboxylate Derivatives for Anti-diabetic Potential

We investigated possible anti-diabetic effect of ferrocene-based acyl ureas: 4-ferrocenyl aniline (PFA), 1-(4-chlorobenzoyl)-3-(4-ferrocenylphenyl) urea (DPC1), 1-(3-chlorobenzoyl)-3-(4-ferrocenylphenyl) urea (DMC1), 1-(2-chlorobenzoyl)-3-(4-ferrocenylphenyl) urea (DOC1) and homoleptic cadmium carboxylates: bis (diphenylacetato) cadmium (II) (DPAA), bis (4-chlorophenylacetato) cadmium (II) (CPAA), using in silico and in vivo techniques. PFA, DPC1, DMC1, DOC1, DPAA and CPAA exhibited high binding affinities (ACE ≥ −350 Kcal/mol) against targets: aldose reductase, peroxisome proliferator-activated receptor γ, 11β-hydroxysteroid dehydrogenase-1, C-alpha glucosidase and glucokinase, while showed moderate affinities (ACE ≥ −250 Kcal/mol) against N-alpha glucosidase, dipeptidyl peptidase-IV, phosphorylated-Akt, glycogen synthase kinase-3β, fructose-1,6-bisphosphatase and phosphoenolpyruvate carboxykinase, whereas revealed lower affinities (ACE < −250 Kcal/mol) vs. alpha amylase, protein tyrosine phosphatases 1B, glycogen phosphorylase and phosphatidylinositol 3 kinase. In alloxan (300 mg/Kg)-induced diabetic mice, DPAA and DPC1 (1–10 mg/Kg) at day 1, 5, 10, 15, and 20th decreased blood glucose levels, compared to diabetic control group and improved the treated animals body weight. DPAA (10 mg/Kg) and DPC1 (5 mg/Kg) in time-dependent manner (30–120 min.) enhanced tolerance of oral glucose overload in mice. DPAA and DPCI dose-dependently at 1, 5, and 10 mg/Kg decreased glycosylated hemoglobin levels in diabetic animals, as caused by metformin. These results indicate that aforementioned derivatives of ferrocene and cadmium possess anti-diabetic potential.

Molecular Dynamics Simulations for Deciphering the Structural Basis of Recognition of Pre-let-7 miRNAs by LIN28

LIN28 protein inhibits biogenesis of miRNAs belonging to the let-7 family by binding to precursor forms of miRNAs. Overexpression of LIN28 and low levels of let-7 miRNAs are associated with several forms of cancer cells. We have performed multiple explicit solvent molecular dynamics simulations ranging from 200 to 500 ns in length on different isoforms of preE-let-7 in complex with LIN28 and also in isolation to identify structural features and key specificity-determining residues (SDRs) that are important for the inhibitory role of LIN28. Our simulations suggest that a conserved structural feature of the loop regions of preE-let-7 miRNAs is more important for LIN28 recognition than sequence conservation among members of the let-7 family or the presence of the GGAG motif in the 3' region. The loop region consisting of a minimum of five nucleotides helps pre-miRNAs to acquire a conformation ideal for binding to LIN28, but pre-let-7c-2 prefers a conformation with a three-nucleotide loop. Thus, our simulations provide a theoretical rationale for the recent experimental observation of the escape of LIN28-mediated repression by pre-let-7c-2. The essential structural and sequence features highlighted in this study might aid in designing synthetic small molecule inhibitors for modulating LIN28-let-7 interaction in malignant cells. We have also identified crucial SDRs of the LIN28-preE-let-7 complex involving 13 residues of LIN28 and 10 residues of the pre-miRNA. On the basis of the conservation profile of these 13 SDRs, we have identified 10 novel proteins that are not annotated as LIN28 like but are similar in sequence, domain, or fold level to LIN28.

Synergy between NMR measurements and MD simulations of protein/RNA complexes: application to the RRMs, the most common RNA recognition motifs

RNA recognition motif (RRM) proteins represent an abundant class of proteins playing key roles in RNA biology. We present a joint atomistic molecular dynamics (MD) and experimental study of two RRM-containing proteins bound with their single-stranded target RNAs, namely the Fox-1 and SRSF1 complexes. The simulations are used in conjunction with NMR spectroscopy to interpret and expand the available structural data. We accumulate more than 50 μs of simulations and show that the MD method is robust enough to reliably describe the structural dynamics of the RRM-RNA complexes. The simulations predict unanticipated specific participation of Arg142 at the protein-RNA interface of the SRFS1 complex, which is subsequently confirmed by NMR and ITC measurements. Several segments of the protein-RNA interface may involve competition between dynamical local substates rather than firmly formed interactions, which is indirectly consistent with the primary NMR data. We demonstrate that the simulations can be used to interpret the NMR atomistic models and can provide qualified predictions. Finally, we propose a protocol for 'MD-adapted structure ensemble' as a way to integrate the simulation predictions and expand upon the deposited NMR structures. Unbiased μs-scale atomistic MD could become a technique routinely complementing the NMR measurements of protein-RNA complexes.

Can We Execute Stable Microsecond-Scale Atomistic Simulations of Protein-RNA Complexes?

We report over 30 μs of unrestrained molecular dynamics simulations of six protein-RNA complexes in explicit solvent. We utilize the AMBER ff99bsc0χ(OL3) RNA force field combined with the ff99SB protein force field and its more recent ff12SB version with reparametrized side-chain dihedrals. The simulations show variable behavior, ranging from systems that are essentially stable to systems with progressive deviations from the experimental structure, which we could not stabilize anywhere close to the starting experimental structure. For some systems, microsecond-scale simulations are necessary to achieve stabilization after initial sizable structural perturbations. The results show that simulations of protein-RNA complexes are challenging and every system should be treated individually. The simulations are affected by numerous factors, including properties of the starting structures (the initially high force field potential energy, resolution limits, conformational averaging, crystal packing, etc.), force field imbalances, and real flexibility of the studied systems. These factors, and thus the simulation behavior, differ from system to system. The structural stability of simulated systems does not correlate with the size of buried interaction surface or experimentally determined binding affinities but reflects the type of protein-RNA recognition. Protein-RNA interfaces involving shape-specific recognition of RNA are more stable than those relying on sequence-specific RNA recognition. The differences between the protein force fields are considerably smaller than the uncertainties caused by sampling and starting structures. The ff12SB improves description of the tyrosine side-chain group, which eliminates some problems associated with tyrosine dynamics.

Molecular dynamic simulations of protein/RNA complexes: CRISPR/Csy4 endoribonuclease

BACKGROUND

Many prokaryotic genomes comprise Clustered Regularly Interspaced Short Palindromic Repeats (CRISPRs) offering defense against foreign nucleic acids. These immune systems are conditioned by the production of small CRISPR-derived RNAs matured from long RNA precursors. This often requires a Csy4 endoribonuclease cleaving the RNA 3'-end.

METHODS

We report extended explicit solvent molecular dynamic (MD) simulations of Csy4/RNA complex in precursor and product states, based on X-ray structures of product and inactivated precursor (55 simulations; ~3.7μs in total).

RESULTS

The simulations identify double-protonated His29 and deprotonated terminal phosphate as the likely dominant protonation states consistent with the product structure. We revealed potential substates consistent with Ser148 and His29 acting as the general base and acid, respectively. The Ser148 could be straightforwardly deprotonated through solvent and could without further structural rearrangements deprotonate the nucleophile, contrasting similar studies investigating the general base role of nucleobases in ribozymes. We could not locate geometries consistent with His29 acting as general base. However, we caution that the X-ray structures do not always capture the catalytically active geometries and then the reactive structures may be unreachable by the simulation technique.

CONCLUSIONS

We identified potential catalytic arrangement of the Csy4/RNA complex but we also report limitations of the simulation technique. Even for the dominant protonation state we could not achieve full agreement between the simulations and the structural data.

GENERAL SIGNIFICANCE

Potential catalytic arrangement of the Csy4/RNA complex is found. Further, we provide unique insights into limitations of simulations of protein/RNA complexes, namely, the influence of the starting experimental structures and force field limitations. This article is part of a Special Issue entitled Recent developments of molecular dynamics.

Interaction investigations of HipA binding to HipB dimer and HipB dimer + DNA complex: a molecular dynamics simulation study

We carried out molecular dynamics simulations and free energy calculations for a series of ternary and diplex models for the HipA protein, HipB dimer, and DNA molecule to address the mechanism of HipA sequestration and the binding order of events from apo HipB/HipA to 2HipA + HipB dimer + DNA complex. The results revealed that the combination of DNA with the HipB dimer is energetically favorable for the combination of HipB dimer with HipA protein. The binding of DNA to HipB dimer induces a long-range allosteric communication from the HipB2 -DNA interface to the HipA-HipB2 interface, which involves the closeness of α1 helices of HipB dimer to HipA protein and formations of extra hydrogen bonds in the HipA-HipB2 interface through the extension of α2/3 helices in the HipB dimer. These simulated results suggested that the DNA molecule, as a regulative media, modulates the HipB dimer conformation, consequently increasing the interactions of HipB dimer with the HipA proteins, which explains the mechanism of HipA sequestration reported by the previous experiment. Simultaneously, these simulations also explored that the thermodynamic binding order in a simulated physiological environment, that is, the HipB dimer first bind to DNA to form HipB dimer + DNA complex, then capturing strongly the HipA proteins to form a ternary complex, 2HipA + HipB dimer + DNA, for sequestrating HipA in the nucleoid.

Analysis of protein-RNA complexes involving a RNA recognition motif engineered to bind hairpins with seven- and eight-nucleotide loops

U1A binds U1hpII, a hairpin RNA with a 10-nucleotide loop. A U1A mutant (ΔK50ΔM51) binds U1hpII-derived hairpins with shorter loops, making it an interesting scaffold for engineering or evolving proteins that bind similarly sized disease-related hairpin RNAs. However, a more detailed understanding of complexes involving ΔK50ΔM51 is likely a prerequisite to generating such proteins. Toward this end, we measured mutational effects for complexes involving U1A ΔK50ΔM51 and U1hpII-derived hairpin RNAs with seven- or eight-nucleotide loops and identified contacts that are critical to the stabilization of these complexes. Our data provide valuable insight into sequence-selective recognition of seven- or eight-nucleotide loop hairpins by an engineered RNA binding protein.

Distinct mechanisms of endocrine disruption of DDT-related pesticides toward estrogen receptor α and estrogen-related receptor γ

Dichlorodiphenyltrichloroethane (DDT) is ubiquitous in the environment, and the exposure to DDT and its related pesticides has long been linked to endocrine disruption. The mechanism of endocrine disruption toward targeted receptors, however, remains unclear. Probing the molecular recognition of DDT analogs by targeted receptors at the atomic level is critical for deciphering this mechanism. Molecular dynamics (MD) simulations were applied to probe the molecular recognition process of DDT and its five analogs, including dichlordiphenyldichloroethylene (DDE), dichlorodiphenyldichloroethane (DDD), methoxychlor (MXC), p,p'-hydroxy-DDT (HPTE), and dicofol by human estrogen receptor (ER) α and human ER-related receptor (ERR) γ. Van der Waals interactions mainly drive the interactions of DDT analogs with ERα ligand-binding domain (LBD) and ERRγ LBD. Minor structural changes of DDT analogs in the number and position of chlorine and phenolic hydroxyl moiety cause differences in binding modes through aromatic stacking and hydrogen bonding and thus affect differently conformational changes of ERα LBD and ERRγ LBD. The binding of DDT analogs affects the helix 12 orientation of ERα LBD but causes no rearrangement of helix 12 of ERRγ LBD. These results extend our understanding of how DDT analogs exert their estrogen-disrupting effects toward different receptors via multiple mechanisms.

Structural analysis of zinc-finger (TTK) + [Cu(BPA)]2+ /[Cu(IDB)]2+ + DNA complexes: an investigation by molecular dynamics simulation

In the present study, the molecular dynamics simulation technique is employed to investigate the hydrogen abstraction possibility from sugar of DNA in two designed complexes of copper-based chemical nuclease [Cu(BPA)](2+) bis(2-pyridylmethyl) amine (BPA) or [Cu(IDB)](2+) N,N-bis(2-benzimidazolylmethyl) amine (IDB) bound to the zinc finger protein Tramtrack (TTK). The simulated results show that each of the designed complexes can form a stable conformation within 30 ns of simulation time with the substrate OOH(-) and an 18-base pair (bp) DNA segment and is located in the major groove of the DNA segment. The active terminal O atom of the OOH(-) substrate is found in close proximity to the target C2'H, C3'H, C4'H or C5'H proton of the DNA in TTK + [Cu(BPA)OOH](+) + DNA or TTK + [Cu(IDB)OOH](+) + DNA complex, which is crucial to propose the hydrogen abstraction possibility that is responsible for the DNA cleavage. The positions of copper-based chemical nucleases bound to TTK may substantially influence the hydrogen abstraction possibility. The structures and sizes of ligands in copper-based nucleases are also found to have influence on the order of difficulty of the hydrogen abstraction from the sugars of DNA.

U1A protein-stem loop 2 RNA recognition: prediction of structural differences from protein mutations

Molecular dynamics (MD) simulations were carried out to compare the free and bound structures of wild type U1A protein with several Phe56 mutant U1A proteins that bind the target stem loop 2 (SL2) RNA with a range of affinities. The simulations indicate the free U1A protein is more flexible than the U1A-RNA complex for both wild type and Phe56 mutant systems. A complete analysis of the hydrogen-bonding (HB) and non-bonded (VDW) interactions over the course of the MD simulations suggested that changes in the interactions in the free U1A protein caused by the Phe56Ala and Phe56Leu mutations may stabilize the helical character in loop 3, and contribute to the weak binding of these proteins to SL2 RNA. Compared with wild type, changes in HB and VDW interactions in Phe56 mutants of the free U1A protein are global, and include differences in β-sheet, loop 1 and loop 3 interactions. In the U1A-RNA complex, the Phe56Ala mutation leads to a series of differences in interactions that resonate through the complex, while the Phe56Leu and Phe56Trp mutations cause local differences around the site of mutation. The long-range networks of interactions identified in the simulations suggest that direct interactions and dynamic processes in both the free and bound forms contribute to complex stability.

Aromatic stacking between nucleobase and enzyme promotes phosphate ester hydrolysis in dUTPase

Aromatic interactions are well-known players in molecular recognition but their catalytic role in biological systems is less documented. Here, we report that a conserved aromatic stacking interaction between dUTPase and its nucleotide substrate largely contributes to the stabilization of the associative type transition state of the nucleotide hydrolysis reaction. The effect of the aromatic stacking on catalysis is peculiar in that uracil, the aromatic moiety influenced by the aromatic interaction is relatively distant from the site of hydrolysis at the alpha-phosphate group. Using crystallographic, kinetics, optical spectroscopy and thermodynamics calculation approaches we delineate a possible mechanism by which rate acceleration is achieved through the remote π–π interaction. The abundance of similarly positioned aromatic interactions in various nucleotide hydrolyzing enzymes (e.g. most families of ATPases) raises the possibility of the reported phenomenon being a general component of the enzymatic catalysis of phosphate ester hydrolysis.

Induced fit or conformational selection for RNA/U1A folding

The hairpin II of U1 snRNA can bind U1A protein with high affinity and specificity. NMR spectra suggest that the loop region of apo-RNA is largely unstructured and undergoes a transition from unstructured to well-folded upon U1Abinding. However, the mechanism that RNA folding coupled protein binding is poorly understood. To get an insight into the mechanism, we have performed explicit-solvent molecular dynamics (MD) to study the folding kinetics of bound RNA and apo-RNA. Room-temperature MD simulations suggest that the conformation of bound RNA has significant adjustment and becomes more stable upon U1A binding. Kinetic analysis of high-temperature MD simulations shows that bound RNA and apo-RNA unfold via a two-state process, respectively. Both kinetics and free energy landscape analyses indicate that bound RNA folds in the order of RNA contracting, U1A binding, and tertiary folding. The predicted Phi-values suggest that A8, C10, A11, and G16 are key bases for bound RNA folding. Mutant Arg52Gln analysis shows that electrostatic interaction and hydrogen bonds between RNA and U1A (Arg52Gln) decrease. These results are in qualitative agreement with experiments. Furthermore, this method could be used in other studies about biomolecule folding upon receptor binding.

Adjustable ellipsoid nanoparticles assembled from re-engineered connectors of the bacteriophage phi29 DNA packaging motor

A 24 x 30 nm ellipsoid nanoparticle containing 84 subunits or 7 dodecamers of the re-engineered core protein of the bacteriophage phi29 DNA packaging motor was constructed. Homogeneous nanoparticles were obtained with simple one-step purification. Electron microscopy and analytical ultracentrifugation were employed to elucidate the structure, shape, size, and mechanism of assembly. The formation of this structure was mediated and stabilized by N-terminal peptide extensions. Reversal of the 84-subunit ellipsoid nanoparticle to its dodecamer subunit was controlled by the cleavage of the extended N-terminal peptide with a protease. The 84 outward-oriented C-termini were conjugated with a streptavidin binding peptide which can be used for the incorporation of markers. This further extends the application of this nanoparticle to pathogen detection and disease diagnosis by signal enhancement.

Prediction of salt and mutational effects on the association rate of U1A protein and U1 small nuclear RNA stem/loop II

We have developed a computational approach for predicting protein-protein association rates (Alsallaq and Zhou, Structure 2007, 15, 215). Here we expand the range of applicability of this approach to protein-RNA binding and report the first results for protein-RNA binding rates predicted from atomistic modeling. The system studied is the U1A protein and stem/loop II of the U1 small nuclear RNA. Experimentally it was observed that the binding rate is significantly reduced by increasing salt concentration while the dissociation changes little with salt concentration, and charges distant from the binding site make marginal contribution to the binding rate. These observations are rationalized. Moreover, predicted effects of salt and charge mutations are found to be in quantitative agreement with experimental results.

Steps towards the formation of a protocell: the possible role of short peptides

The paper deals with molecular self-organization leading to formation of a protocell. Plausible steps towards a protocell include: polymerization of peptides and oligonucleotides on mineral surfaces; coevolution of peptides and oligonucleotides with formation of collectively autocatalytic sets; self-organization of short peptides into vesicles; entrapment of the peptide/oligonucleotide systems in mixed peptide and simple amphiphile membranes; and formation of functioning protocells with metabolism and cell division. The established propensity of short peptides to self-ordering and to formation of vesicles makes this sequence plausible. We further suggest that evolution of a protocell produced cellular ancestors of viruses as well as ancestors of cellular organisms.

Do collective atomic fluctuations account for cooperative effects? Molecular dynamics studies of the U1A-RNA complex

A complete understanding of gene expression relies on a comprehensive understanding of the protein-RNA recognition process. However, the study of protein-RNA recognition is complicated by many factors that contribute to both binding affinity and specificity, including structure, energetics, dynamical motions, and cooperative interactions. Several recent studies have suggested that energetic coupling between residues contributes to formation of the complex between the U1A protein and stem loop 2 of U1 snRNA as a consequence of a cooperative network of interactions. We have performed molecular dynamics simulations on the U1A-RNA complex, including explicit water and counterions, and analyzed the results based on the calculated positional cross-correlations of atomic fluctuations. The results indicate that cross-correlations calculated on a per residue basis agree well with the observed inter-residue cooperativity and predict that the networks identified to date may also be coupled into an extensive hyper-network that reflects the intrinsic rigidity of the RNA recognition motif. In addition, we report a comparison of the MD calculated correlations with the results of a positional covariance analysis based on the sequences of 330 RNA recognition motifs, including U1A. The calculated inter-residue cross-correlations agree very well with the results of the sites exhibiting positional covariance. Collectively, these results strongly support the hypothesis that collective fluctuations contribute to cooperativity and the corresponding observed thermodynamic coupling. Predictions of additional sites in U1A that may be involved in cooperative networks are advanced.

Role of aromatic amino acids in protein-nucleic acid recognition

Statistical analysis of structures from the PBD has been used to examine the role that the aromatic amino acids play in protein-nucleic acid recognition. In protein-DNA complexes, the residues Phe and His are found to bind selectively to the DNA chain--Phe to A and T, and His to T and G. The preferred binding modes are identified, and the interactions involving Phe are shown to be important in the transcription process. In protein-RNA complexes, Phe is found to occur far less often and is instead replaced by Trp, which binds selectively to C and G, offering a possible mechanism for differentiation between the two nucleic acids. SASA analysis of the two sets of complexes suggests that all of the aromatic amino acids are more heavily involved in binding than would be expected on the balance of probability. Phe and Tyr occur approximately equal in both sets of data, whereas the proportions of His and Trp vary considerably, supporting the idea that these residues may be involved in differentiating between the two nucleic acids.

Affinity and specificity of protein U1A-RNA complex formation based on an additive component free energy model

An MM-GBSA computational protocol was used to investigate wild-type U1A-RNA and F56 U1A mutant experimental binding free energies. The trend in mutant binding free energies compared to wild-type is well-reproduced. Following application of a linear-response-like equation to scale the various energy components, the binding free energies agree quantitatively with observed experimental values. Conformational adaptation contributes to the binding free energy for both the protein and the RNA in these systems. Small differences in DeltaGs are the result of different and sometimes quite large relative contributions from various energetic components. Residual free energy decomposition indicates differences not only at the site of mutation, but throughout the entire protein. MM-GBSA and ab initio calculations performed on model systems suggest that stacking interactions may nearly, but not completely, account for observed differences in mutant binding affinities. This study indicates that there may be different underlying causes of ostensibly similar experimentally observed binding affinities of different mutants, and thus recommends caution in the interpretation of binding affinities and specificities purely by inspection.

Template-Controlled Face-to-Face Stacking of Olefinic and Aromatic Carboxylic Acids in the Solid State

A rigid 1,8-dipyridylnaphthalene template for supramolecular organization of unsaturated dicarboxylic acids in the solid state has been developed. Self-assembly of the template and aromatic dicarboxylic acids generates nondistorted paracyclophane and pyridinophane architectures with perfectly superimposed pi-systems. The dipyridyl template can be used to promote stereoselective solid-state dimerization of olefinic dicarboxylic acids.

Molecular dynamics simulations of RNA: an in silico single molecule approach

RNA molecules are now known to be involved in the processing of genetic information at all levels, taking on a wide variety of central roles in the cell. Understanding how RNA molecules carry out their biological functions will require an understanding of structure and dynamics at the atomistic level, which can be significantly improved by combining computational simulation with experiment. This review provides a critical survey of the state of molecular dynamics (MD) simulations of RNA, including a discussion of important current limitations of the technique and examples of its successful application. Several types of simulations are discussed in detail, including those of structured RNA molecules and their interactions with the surrounding solvent and ions, catalytic RNAs, and RNA-small molecule and RNA-protein complexes. Increased cooperation between theorists and experimentalists will allow expanded judicious use of MD simulations to complement conceptually related single molecule experiments. Such cooperation will open the door to a fundamental understanding of the structure-function relationships in diverse and complex RNA molecules. .

Ab initio calculations on π-stacked thiophene dimer, trimer, and tetramer: Structure, interaction energy, cooperative effects, and intermolecular electronic parameters

Pi-stacked complexes formed by two, three, and four thiophene rings have been investigated using abinitio quantum mechanical calculations. The relative orientation between the rings was investigated for each complex by exploring the corresponding potential energy surface at the MP2/6-31+G(d,p) level, the inter-ring distance, and the degree of tilting being examined in each case. Interaction energies were calculated at the MP2, MP3, MP4, and CCSD, levels of theory. Negligible or even slightly positive n-body effects have been predicted for the stacked thiophene arrangements studied in this work. This is consequence of the cancellation of favorable induction contribution by the destabilizing dispersion component. On the other hand analysis of the optimized geometries obtained for the trimer and tetramer revealed that the orientation of the rings presents a preferred degree of periodicity. Finally, we found that the lowest transition energy decreases when the size of the complex increases, this feature being attributed to desestabilization of the HOMO and stabilization of the LUMO that occur simultaneously.

A study of collective atomic fluctuations and cooperativity in the U1A-RNA complex based on molecular dynamics simulations

Cooperative interactions play an important role in recognition and binding in macromolecular systems. In this study, we find that cross-correlated atomic fluctuations can be used to identify cooperative networks in a protein-RNA system. The dynamics of the RRM-containing protein U1A-stem loop 2 RNA complex have been calculated theoretically from a 10 ns molecular dynamics (MD) simulation. The simulation was analyzed by calculating the covariance matrix of all atomic fluctuations. These matrix elements are then presented in the form of a two-dimensional grid, which displays fluctuations on a per residue basis. The results indicate the presence of strong, selective cross-correlated fluctuations throughout the RRM in U1A-RNA. The atomic fluctuations correspond well with previous biophysical studies in which a multiplicity of cooperative networks have been reported and indicate that the various networks identified in separate individual experiments are fluctuationally correlated into a hyper-network encompassing most of the RRM. The calculated results also correspond well with independent results from a statistical covariance analysis of 330 aligned RRM sequences. This method has significant implications as a predictive tool regarding cooperativity in the protein-nucleic acid recognition process.

The NMR structure of an internal loop from 23S ribosomal RNA differs from its structure in crystals of 50s ribosomal subunits

Internal loops play an important role in structure and folding of RNA and in recognition of RNA by other molecules such as proteins and ligands. An understanding of internal loops with propensities to form a particular structure will help predict RNA structure, recognition, and function. The structures of internal loops 5' 1009CUAAG1013 3'/3' 1168GAAGC1164 5' and 5' 998CUAAG1002 3'/3' 1157GAAGC1153 5' from helix 40 of the large subunit rRNA in Deinococcus radiodurans and Escherichia coli, respectively, are phylogenetically conserved, suggesting functional relevance. The energetics and NMR solution structure of the loop were determined in the duplex 5' 1GGCUAAGAC9 3'/3' 18CCGAAGCUG10 5'. The internal loop forms a different structure in solution and in the crystal structures of the ribosomal subunits. In particular, the crystal structures have a bulged out adenine at the equivalent of position A15 and a reverse Hoogsteen UA pair (trans Watson-Crick/Hoogsteen UA) at the equivalent of U4 and A14, whereas the solution structure has a single hydrogen bond UA pair (cis Watson-Crick/sugar edge A15U4) between U4 and A15 and a sheared AA pair (trans Hoogsteen/sugar edge A14A5) between A5 and A14. There is cross-strand stacking between A6 and A14 (A6/A14/A15 stacking pattern) in the NMR structure. All three structures have a sheared GA pair (trans Hoogsteen/sugar edge A6G13) at the equivalent of A6 and G13. The internal loop has contacts with ribosomal protein L20 and other parts of the RNA in the crystal structures. These contacts presumably provide the free energy to rearrange the base pairing in the loop. Evidently, molecular recognition of this internal loop involves induced fit binding, which could confer several advantages. The predicted thermodynamic stability of the loop agrees with the experimental value, even though the thermodynamic model assumes a Watson-Crick UA pair.

The role of RNA structure in the interaction of U1A protein with U1 hairpin II RNA

The N-terminal RNA Recognition Motif (RRM1) of the spliceosomal protein U1A interacting with its target U1 hairpin II (U1hpII) has been used as a paradigm for RRM-containing proteins interacting with their RNA targets. U1A binds to U1hpII via direct interactions with a 7-nucleotide (nt) consensus binding sequence at the 5' end of a 10-nt loop, and via hydrogen bonds with the closing C-G base pair at the top of the RNA stem. Using surface plasmon resonance (Biacore), we have examined the role of structural features of U1hpII in binding to U1A RRM1. Mutational analysis of the closing base pair suggests it plays a minor role in binding and mainly prevents "breathing" of the loop. Lengthening the stem and nontarget part of the loop suggests that the increased negative charge of the RNA might slightly aid association. However, this is offset by an increase in dissociation, which may be caused by attraction of the RRM to nontarget parts of the RNA. Studies of a single stranded target and RNAs with untethered loops indicate that structure is not very relevant for association but is important for complex stability. In particular, breaking the link between the stem and the 5' side of the loop greatly increases complex dissociation, presumably by hindering simultaneous contacts between the RRM and stem and loop nucleotides. While binding of U1A to a single stranded target is much weaker than to U1hpII, it occurs with nanomolar affinity, supporting recent evidence that binding of unstructured RNA by U1A has physiological significance.

A model for the role of short self-assembled peptides in the very early stages of the origin of life

The molecular basis of the origin of life is one of the most fundamental questions in modern biology. While the "RNA world" hypothesis offers a very sensible model for the evolvement of the current biochemical networks, there is a lack of knowledge about the early steps that led to the formation of the first RNA molecules. This issue is essential as it is practically impossible that complex molecules as functional RNA oligonucleotides had evolved spontaneously. It was recently demonstrated that peptide molecules as simple as dipeptides can self-assemble into well-ordered tubular, fibrilar, and closed-cage structures. Other studies have confirmed the ability of dipeptides to act as catalysts and the capability of other peptides, as short as tripeptides, to serve as a template for nucleotide binding and orientation. Unlike complex RNA molecules, the spontaneous formation of functional short peptides in the primordial earth conditions is very likely. We suggest a novel mechanism for the origin of life that is based on the ability of short peptides to form encapsulated structures, catalyst chemical reaction, and serve as highly ordered template for the assembly of nucleotides. This model may explain the early events that led to the formation of the current biochemical machinery that combines the elaborated and coordinated interaction between nucleic acids and proteins to allow the function of living systems.

Structural dynamics of precursor and product of the RNA enzyme from the hepatitis delta virus as revealed by molecular dynamics simulations

The hepatitis delta virus (HDV) ribozyme is a self-cleaving RNA enzyme involved in the replication of a human pathogen, the hepatitis delta virus. Recent crystal structures of the precursor and product of self-cleavage, together with detailed kinetic analyses, have led to hypotheses on the catalytic strategies employed by the HDV ribozyme. We report molecular dynamics (MD) simulations (approximately 120 ns total simulation time) to test the plausibility that specific conformational rearrangements are involved in catalysis. Site-specific self-cleavage requires cytidine in position 75 (C75). A precursor simulation with unprotonated C75 reveals a rather weak dynamic binding of C75 in the catalytic pocket with spontaneous, transient formation of a H-bond between U-1(O2') and C75(N3). This H-bond would be required for C75 to act as the general base. Upon protonation in the precursor, C75H+ has a tendency to move towards its product location and establish a firm H-bonding network within the catalytic pocket. However, a C75H+(N3)-G1(O5') H-bond, which would be expected if C75 acted as a general acid catalyst, is not observed on the present simulation timescale. The adjacent loop L3 is relatively dynamic and may serve as a flexible structural element, possibly gated by the closing U20.G25 base-pair, to facilitate a conformational switch induced by a protonated C75H+. L3 also controls the electrostatic environment of the catalytic core, which in turn may modulate C75 base strength and metal ion binding. We find that a distant RNA tertiary interaction involving a protonated cytidine (C41) becomes unstable when left unprotonated, leading to disruptive conformational rearrangements adjacent to the catalytic core. A Na ion temporarily compensates for the loss of the protonated hydrogen bond, which is strikingly consistent with the experimentally observed synergy between low pH and high Na+ concentrations in mediating residual self-cleavage of the HDV ribozyme in the absence of divalents.

Highly congested nondistorted diheteroarylnaphthalenes: model compounds for the investigation of intramolecular pi-stacking interactions

[structure: see text] The rigid, highly congested structure of 1,8-diacridylnaphthalenes has been studied in solution and in the solid state. The unique geometry of these compounds forces the acridyl rings to undergo face-to-face interactions while rendering T-shaped orientations and face-to-edge interactions impossible. Crystallographic analysis shows that splaying between the heteroaryl rings decreases while twisting between the cofacial rings increases as the acridyl nitrogens of the 1,8-diacridylnaphthalene framework are subsequently oxidized. The peri-acridyl rings are slightly splayed but remain perfectly planar in all cases. The significant decrease in splaying indicates enhanced pi-pi-attraction between the electron-rich acridyl N-oxide moieties, which is in agreement with recently reported symmetry-adapted perturbation theory calculations. The pi-stacking and the molecular geometry between the acridyl rings observed in the solid state have been confirmed through in-solution studies showing characteristic proton NMR upfield shifts and optical properties indicative of static intramolecular arene-arene interactions. Acridyl protons located directly above the adjacent aryl moiety as a consequence of twisting between the heteroaryl rings were identified by COSY NMR measurements and found to intrude into the pi-cloud and diamagnetic ring current of the neighboring acridine. Different shapes and strong red shifts of the fluorescence emission maxima of the diacridylnaphthalenes in comparison to parental acridyl monomers have been attributed to static excimer emission.

A single intermolecular contact mediates intramolecular stabilization of both RNA and protein

An arginine-rich peptide from the Jembrana disease virus (JDV) Tat protein is a structural "chameleon" that binds bovine immunodeficiency virus (BIV) or HIV TAR RNAs in two different binding modes, with an affinity for BIV TAR even higher than the cognate BIV peptide. We determined the NMR structure of the JDV Tat-BIV TAR high-affinity complex and found that the C-terminal tyrosine in JDV Tat forms a network of inter- and intramolecular hydrogen bonding and stacking interactions that simultaneously stabilize the beta-hairpin conformation of the peptide and a base triple in the RNA. A neighboring histidine also appears to help stabilize the peptide conformation. Induced fit binding is recurrent in protein-protein and protein-nucleic acid interactions, and the JDV Tat complex demonstrates how high affinity can be achieved not only by optimization of the binding interface but also by inducing new intramolecular contacts that stabilize each binding partner. Comparison to the cognate BIV Tat peptide-TAR complex shows how such a costabilization mechanism can evolve with only small changes to the peptide sequence. In addition, the bound structure of BIV TAR in the chameleon peptide complex is strikingly similar to the bound conformation of HIV TAR, suggesting new strategies for the development of HIV TAR binding molecules.

A binding mechanism in protein-nucleotide interactions: implication for U1A RNA binding

We present a close electronic view of the protein-base interface for the N-terminal domain of the human protein U1A. Combining accurate mixed quantum mechanics/molecular mechanics techniques and protein structure prediction methods, we provide a detailed electronic structure description of the protein-RNA stacking interactions. Our analysis indicates the evolution of the protein structure optimizing the interaction between Asp-92 and the RNA bases. The results show a direct coupling of the C-terminal tail and Asp-92, providing a direct rationalization of the experimentally determined role of the C-terminal domain in RNA binding. Here, we propose a mechanism where a protein side chain, with a delocalized electronic pi system, assists in the nucleotide binding. The binding mechanism involves a short-range interaction of the side chain with the nucleotide base and an electronic long-range interaction through a sandwich-stacking motif. The structural motif of the binding mechanism is observed in similar protein-RNA interactions and in various protein-ATP-binding sites.

The snRNP 15.5K protein folds its cognate K-turn RNA: a combined theoretical and biochemical study

The human 15.5K protein binds to the 5' stem-loop of U4 snRNA, promotes the assembly of the spliceosomal U4/U6 snRNP, and is required for the recruitment of the 61K protein and the 20/60/90K protein complex to the U4 snRNA. In the crystallographic structure of the 15.5K-U4 snRNA complex, the conformation of the RNA corresponds to the family of kink-turn (K-turn) structural motifs. We simulated the complex and the free RNA, showing how the protein binding and the intrinsic flexibility contribute to the RNA folding process. We found that the RNA is significantly more flexible in the absence of the 15.5K protein. Conformational transitions such as the interconversion between alternative purine stacking schemes, the loss of G-A base pairs, and the opening of the K-turn occur only in the free RNA. Furthermore, the stability of one canonical G-C base pair is influenced both by the binding of the 15.5K protein and the nature of the adjacent structural element in the RNA. We performed chemical RNA modification experiments and observed that the free RNA lacks secondary structure elements, a result in excellent agreement with the simulations. Based on these observations, we propose a protein-assisted RNA folding mechanism in which the RNA intrinsic flexibility functions as a catalyst.

Site-directed mutagenesis reveals the thermodynamic requirements for single-stranded DNA recognition by the telomere-binding protein Cdc13

The essential Saccharomyces cerevisiae protein Cdc13 binds the conserved single-stranded overhang at the end of telomeres and mediates access of protein complexes involved in both end-capping and telomerase activity. The single-stranded DNA-binding domain (ssDBD) of Cdc13 exhibits both high affinity (K(d) of 3 pM) and sequence specificity for the GT-rich sequences present at yeast telomeres. We have used the ssDBD of Cdc13 to understand the sequence-specific recognition of extended single-stranded DNA (ssDNA). The recent structure of the Cdc13 DNA-binding domain revealed that ssDNA is recognized by a large protein surface containing an oligonucleotide/oligosaccharide-binding fold (OB-fold) augmented by an extended 30-amino acid loop. Contacts to ssDNA occur via a contiguous surface of aromatic, hydrophobic, and basic residues. A complete alanine scan of the binding interface has been used to determine the contribution of each contacting side chain to binding affinity. Substitution of any aromatic or hydrophobic residue at the interface was deleterious to binding (20 to >700-fold decrease in binding affinity), while tolerance for replacement of basic residues was observed. The important aromatic and hydrophobic contacts are spread throughout the extended interface, indicating that the entire surface is both structurally and thermodynamically required for binding. While all of these contacts are important, several of the individual alanine substitutions that abolish binding cluster to one region of the protein surface. This region is vital for recognition of four bases at the 5' end of the DNA and constitutes a "hotspot" of binding affinity.

Design of an adenosine analogue that selectively improves the affinity of a mutant U1A protein for RNA

The RNA recognition motif (RRM), one of the most common RNA binding domains, contains three highly conserved aromatic amino acids that participate in stacking interactions with RNA bases. We have investigated the contribution of these highly conserved aromatic amino acids to the affinity of the complex formed between the N-terminal RRM of the U1A protein and stem loop 2 of U1 snRNA. Previously, we found that substitution of one of these conserved aromatic amino acids, Phe56, with Ala resulted in a large destabilization of the complex. Here, we have modified A6, the base in stem loop 2 RNA that stacks with Phe56, to compensate for a portion of the destabilization caused by the Phe56Ala mutation. We have designed two modified adenosines, A-3CPh and A-4CPh, in which a phenyl group is linked to the adenosine such that it may replace the phenyl group that is eliminated by the Phe56Ala mutation in the complex. We have found that incorporation of A-3CPh into stem loop 2 RNA stabilizes the complex formed with Phe56Ala by 0.6 kcal/mol, while incorporation of A-4CPh into stem loop 2 RNA stabilizes this complex by 1.8 kcal/mol. Either base modification destabilizes the wild-type complex by 0.8-0.9 kcal/mol. Experiments with other U1A mutant proteins suggest that the stabilization of the complex between the Phe56Ala U1A protein and stem loop 2 RNA is due to a specific interaction between the Phe56Ala U1A protein and A6-4CPh stem loop 2 RNA.

Theoretical methods for the simulation of nucleic acids

Different theoretical methods for the description of nucleic acid structures are reviewed. Firstly, we introduce the concept of classical force-field in the context of nucleic acid structures, discussing their accuracy. We then examine theoretical approaches to the description of nucleic acids based on: i) a rigid or quasi-rigid description of the molecule, ii) molecular mechanics optimization, and iii) molecular dynamics. Special emphasis is made ion current state of the art molecular dynamics simulations of nucleic acids structures.

Molecular dynamics simulation studies of induced fit and conformational capture in U1A-RNA binding: do molecular substates code for specificity?

Molecular dynamics (MD) simulations on stem loop 2 of U1 small nuclear RNA and a construct of the U1A protein were carried out to obtain predictions of the structures for the unbound forms in solution and to elucidate dynamical aspects of induced fit upon binding. A crystal structure of the complex between the U1A protein and stem loop 2 RNA and an NMR structure for the uncomplexed form of the U1A protein are available from Oubridge et al. (Nature, 1994, Vol. 372, pp. 432-438) and Avis et al. (Journal of Molecular Biology, 1996, Vol. 257, pp. 398-411), respectively. As a consequence, U1A-RNA binding is a particularly attractive case for investigations of induced fit in protein-nucleic acid complexation. When combined with the available structural data, the results from simulations indicate that structural adaptation of U1A protein and RNA define distinct mechanisms for induced fit. For the protein, the calculations indicate that induced fit upon binding involves a non-native thermodynamic substate in which the structure is preorganized for binding. In contrast, induced fit of the RNA involves a distortion of the native structure in solution to an unstable form. However, the RNA solution structures predicted from simulation show evidence that structures in which groups of bases are favorably oriented for binding the U1A protein are thermally accessible. These results, which quantify with computational modeling recent proposals on induced fit and conformational capture by Leuillot and Varani (Biochemistry, 2001, Vol. 40, pp. 7947-7956) and by Williamson (Nature Structural Biology, 2000, Vol. 7, pp. 834-837) suggest an important role for intrinsic molecular architecture and substates other than the native form in the specificity of protein-RNA interactions.

Substitution of an essential adenine in the U1A-RNA complex with a non-polar isostere

The RNA recognition motif (RRM) binds to single-stranded RNA target sites of diverse sequences and structures. A conserved mode of base recognition by the RRM involves the simultaneous formation of a network of hydrogen bonds with the base functional groups and a stacking interaction between the base and a highly conserved aromatic amino acid. We have investigated the energetic contribution of the functional groups involved in the recognition of an essential adenine, A6, in stem-loop 2 of U1 snRNA by the N-terminal RRM of the U1A protein. Previously, we found that elimination of individual hydrogen bond donors and acceptors on A6 destabilized the complex by 0.8-1.9 kcal/mol, while mutation of the aromatic amino acid (Phe56) that stacks with A6 to Ala destabilized the complex by 5.5 kcal/mol. Here we continue to probe the contribution of A6 to complex stability through mutation of both the RNA and protein. We have removed two hydrogen-bonding functional groups by introducing a U1A mutation, Ser91Ala, and replacing A6 with tubercidin, purine, or 1-deazaadenine. We find that the complex is destabilized an additional 1.2-2.6 kcal/mol by the elimination of the second hydrogen bond donor or acceptor. Surprisingly, deletion of all of the functional groups involved in hydrogen bonds with the U1A protein by substituting adenine with 4-methylindole reduced the binding free energy by only 2.0 kcal/mol. Experiments with U1A proteins containing mutations of Phe56 suggested that improved stacking interactions due to the greater hydrophobicity of 4-methylindole than adenine may be partly responsible for the small destabilization of the complex upon substitution of 4-methylindole for A6. The data imply that hydrophobic interactions can compensate energetically for the disruption of the complex hydrogen-bonding network between nucleotide and protein.

Analysis of the pi-pi stacking interactions between the aminoglycoside antibiotic kinase APH(3')-IIIa and its nucleotide ligands

A key contact in the active site of an aminoglycoside phosphotransferase enzyme (APH(3')-IIIa) is a pi-pi stacking interaction between Tyr42 and the adenine ring of bound nucleotides. We investigated the prevalence of similar Tyr-adenine contacts and found that many different protein systems employ Tyr residues in the recognition of the adenine ring. The geometry of these stacking interactions suggests that electrostatics play a role in the attraction between these aromatic systems. Kinetic and calorimetric experiments on wild-type and mutant forms of APH(3')-IIIa yielded further experimental evidence of the importance of electrostatics in the adenine binding region and suggested that the stacking interaction contributes approximately 2 kcal/mol of binding energy. This type of information concerning the forces that govern nucleotide binding in APH(3')-IIIa will facilitate inhibitor design strategies that target the nucleotide binding site of APH-type enzymes.

Complex role of the beta 2-beta 3 loop in the interaction of U1A with U1 hairpin II RNA

RNA recognition motifs (RRMs) are characterized by highly conserved regions located centrally on a beta-sheet, which forms the RNA binding surface. Variable flanking regions, such as the loop connecting beta-strands 2 and 3, are thought to be important in determining the RNA-binding specificities of individual RRMs. The N-terminal RRM of the spliceosomal U1A protein mediates binding to an RNA hairpin (U1hpII) in the U1 small nuclear RNA. In this complex, the beta(2)-beta(3) loop protrudes through the 10-nucleotide RNA loop. Shortening of the RNA loop strongly perturbs binding, suggesting that an optimal "fit" of the beta(2)-beta(3) loop into the RNA loop is an important factor in complexation. To understand this interaction further, we mutated or deleted loop residues Lys(50) and Met(51), which protrude centrally into the RNA loop but do not make any direct contacts to the bases. Using BIACORE, we analyzed the ability of these U1A mutants to bind to wild type RNAs, or RNAs with shortened loops. Alanine replacement mutations only modestly affected binding to wild type U1hpII. Interestingly, simultaneous replacement of Lys(50) and Met(51) with alanine appeared to alleviate the loss of binding caused by shortening of the RNA loop. Deletion of Lys(50) or Met(51) caused a dramatic loss in stability of the U1A.U1hpII complex. However, deletion of both residues simultaneously was much less deleterious. Simulated annealing molecular dynamics analyses suggest this is due to the ability of this mutant to rearrange flanking amino acids to substitute for the two deleted residues. The double deletion mutant also exhibited substantially reduced negative effects of RNA loop shortening, suggesting the rearranged loop is better able to accommodate a short RNA loop. Our results indicate that one of the roles of the beta(2)-beta(3) loop is to provide a steric fit into the RNA loop, thereby stabilizing the RNA.protein complex.

Free-energy component analysis of 40 protein-DNA complexes: a consensus view on the thermodynamics of binding at the molecular level

Noncovalent association of proteins to specific target sites on DNA--a process central to gene expression and regulation--has thus far proven to be idiosyncratic and elusive to generalizations on the nature of the driving forces. The spate of structural information on protein--DNA complexes sets the stage for theoretical investigations on the molecular thermodynamics of binding aimed at identifying forces responsible for specific macromolecular recognition. Computation of absolute binding free energies for systems of this complexity transiting from structural information is a stupendous task. Adopting some recent progresses in treating atomic level interactions in proteins and nucleic acids including solvent and salt effects, we have put together an energy component methodology cast in a phenomenological mode and amenable to systematic improvements and developed a computational first atlas of the free energy contributors to binding in approximately 40 protein-DNA complexes representing a variety of structural motifs and functions. Illustrating vividly the compensatory nature of the free energy components contributing to the energetics of recognition for attaining optimal binding, our results highlight unambiguously the roles played by packing, electrostatics including hydrogen bonds, ion and water release (cavitation) in protein-DNA binding. Cavitation and van der Waals contributions without exception favor complexation. The electrostatics is marginally unfavorable in a consensus view. Basic residues on the protein contribute favorably to binding despite the desolvation expense. The electrostatics arising from the acidic and neutral residues proves unfavorable to binding. An enveloping mode of binding to short stretches of DNA makes for a strong unfavorable net electrostatics but a highly favorable van der Waals and cavitation contribution. Thus, noncovalent protein-DNA association is a system-specific fine balancing act of these diverse competing forces. With the advances in computational methods as applied to macromolecular recognition, the challenge now seems to be to correlate the differential (initial vs. final) energetics to substituent effects in drug design and to move from affinity to specificity.

Investigation of a conserved stacking interaction in target site recognition by the U1A protein

Three highly conserved aromatic residues in RNA recognition motifs (RRM) participate in stacking interactions with RNA bases upon binding RNA. We have investigated the contribution of one of these aromatic residues, Phe56, to the complex formed between the N-terminal RRM of the spliceosomal protein U1A and stem-loop 2 of U1 snRNA. Previous work showed that the aromatic group is important for high affinity binding. Here we probe how mutation of Phe56 affects the kinetics of complex dissociation, the strength of the hydrogen bonds formed between U1A and the base that stacks with Phe56 (A6) and specific target site recognition. Substitution of Phe56 with Trp or Tyr increased the rate of dissociation of the complex, consistent with previously reported results. However, substitution of Phe56 with His decreased the rate of complex association, implying a change in the initial formation of the complex. Simultaneous modification of residue 56 and A6 revealed energetic coupling between the aromatic group and the functional groups of A6 that hydrogen bond to U1A. Finally, mutation of Phe56 to Leu reduced the ability of U1A to recognize stem-loop 2 correctly. Taken together, these experiments suggest that Phe56 contributes to binding affinity by stacking with A6 and participating in networks of energetically coupled interactions that enable this conserved aromatic amino acid to play a complex role in target site recognition.