Electrostatic potential of aminoacyl-tRNA synthetase navigates tRNA on its pathway to the binding site

Enhanced sampling without borders: on global biasing functions and how to reweight them

Molecular dynamics (MD) simulations are a powerful tool to follow the time evolution of biomolecular motions in atomistic resolution. However, the high computational demand of these simulations limits the timescales of motions that can be observed. To resolve this issue, so called enhanced sampling techniques are developed, which extend conventional MD algorithms to speed up the simulation process. Here, we focus on techniques that apply global biasing functions. We provide a broad overview of established enhanced sampling methods and promising new advances. As the ultimate goal is to retrieve unbiased information from biased ensembles, we also discuss benefits and limitations of common reweighting schemes. In addition to concisely summarizing critical assumptions and implications, we highlight the general application opportunities as well as uncertainties of global enhanced sampling.

Stimuli-Free Transcuticular Delivery of Zn Microelement Using Biopolymeric Nanovehicles: Experimental, Theoretical, and In Planta Studies

This paper reports one-step synthesis of polysaccharide-based nanovehicles, capable of transporting ionic zinc via plant cuticle without auxiliary stimulation. Delivery of highly hydrophilic nutritive microelements via the hydrophobic cuticle of plant foliage is one of the major challenges in modern agriculture. In traditional nutrition via roots, up to 80% of microelements permeate to soil and get wasted; therefore, foliar treatment is an environmentally and economically preferable alternative. Carboxymethyl cellulose (CMC) was modified to amphiphilic N-octylamide-derivative (CMC-8), which spontaneously self-assemble to nanovehicles. It was found that hydrophobic substituents endow a biopolymer with unexpected affinity toward a hydrophilic payload. CMC-8 nanovehicles effectively encapsulated ionic zinc (ZnSO₄) and delivered it upon foliar application to pepper (Capsicum annuum) and tomato (Solanum lycopersicum) plants. Zinc uptake and translocation in plants were monitored by SEM-EDS and fluorescence microscopic methods. In planta monitoring of the carrier was done by labeling nanovehicles with fluorescent carbon dots. Three-dimensional (3-D) structural modeling and conformational dynamics explained the CMC-8 self-assembly mechanism and zinc coordination phenomenon upon introduction of hydrophobic substituents.

Computational investigation of the binding of a designed peptide to λ light chain amyloid fibril

Systemic light chain amyloidosis (AL) causes a malignant pathology associated with the formation of amyloid fibrils that deposit in human organs and tissues, leading to dysfunction and severe morbidity. Amyloid fibril-reactive antibodies have been used to remove amyloid from organs and are effective in restoring organ function in patients with AL amyloidosis. Unfortunately, antibodies do not bind amyloid in all AL patients, nor do they efficiently bind many other forms of amyloid. Recently, a synthetic peptide P62 was developed, which binds many forms of systemic amyloidosis and can be further modified and fused to a high-affinity peptide epitope to expand its utility as a novel amyloid immunotherapeutic. However, the molecular-level details of P62-fibril binding mechanisms, critical for future peptide design, are unclear. Here, we combine protein docking, all-atom molecular dynamics simulation and umbrella sampling to study the dynamical interactions between peptide P62 and a structural model of the λ light chain in systemic amyloidosis. We found that P62 only binds to the canonical interface of the fibril where the peptide inserts into the fibril groove and its two termini are more mobile than the helix core. Our results also revealed an important role of the lysine residues of P62 in the binding process by forming initial contacts with aspartic acids on the fibril surface. Collectively, our computational study provided molecular-level insights into the binding mechanism between an amyloid fibril model and peptide P62, which could lay a foundation for rational design of peptides for improved amyloid diagnosis and immunotherapy.

Role of the transfer ribonucleic acid (tRNA) bound magnesium ions in the charging step of aminoacylation reaction in the glutamyl tRNA synthetase and the seryl tRNA synthetase bound with cognate tRNA

Aminoacylation reaction is the first step of protein biosynthesis. Transfer RNA (tRNA) is charged with an amino acid in this reaction and the reaction is catalyzed by aminoacyl tRNA synthetase enzyme (aaRS). In the present work, we use classical molecular dynamics simulation to show that the tRNA bound Mg²⁺ ions significantly influence the charging step of class I ^TtGluRS: Glu-AMP: tRNA^Glu and class II dimeric ^TtSerRS: Ser-AMP: tRNA^Ser. The CCA end of the acceptor terminal is disordered in the absence of coordinated Mg²⁺ ions and the CCA end can freely explore beyond the specific conformational space of the tRNA in its precharging state. A balance between the conformational disorder of the tRNA and the restriction imposed on the CCA terminal via coordination with the Mg²⁺ ions is needed for the placement of the CCA terminal in a precharging state organization. This result provides a molecular-level explanation of the experimental observation that the presence of Mg²⁺ ions is a necessary condition for a successful aminoacylation reaction.Communicated by Ramaswamy H. Sarma.

Shark Antibody Variable Domains Rigidify Upon Affinity Maturation-Understanding the Potential of Shark Immunoglobulins as Therapeutics

Sharks and other cartilaginous fish are the phylogenetically oldest living organisms that have antibodies as part of their adaptive immune system. As part of their humoral adaptive immune response, they produce an immunoglobulin, the so-called immunoglobulin new antigen receptor (IgNAR), a heavy-chain only antibody. The variable domain of an IgNAR, also known as V NAR , binds the antigen as an independent soluble domain. In this study, we structurally and dynamically characterized the affinity maturation mechanism of the germline and somatically matured (PBLA8) V NAR to better understand their function and their applicability as therapeutics. We observed a substantial rigidification upon affinity maturation, which is accompanied by a higher number of contacts, thereby contributing to the decrease in flexibility. Considering the static x-ray structures, the observed rigidification is not obvious, as especially the mutated residues undergo conformational changes during the simulation, resulting in an even stronger network of stabilizing interactions. Additionally, the simulations of the V NAR in complex with the hen egg-white lysozyme show that the V NAR antibodies evidently follow the concept of conformational selection, as the binding-competent state already preexisted even without the presence of the antigen. To have a more detailed description of antibody-antigen recognition, we also present here the binding/unbinding mechanism between the hen egg-white lysozyme and both the germline and matured V NAR s. Upon maturation, we observed a substantial increase in the resulting dissociation-free energy barrier. Furthermore, we were able to kinetically and thermodynamically describe the binding process and did not only identify a two-step binding mechanism, but we also found a strong population shift upon affinity maturation toward the native binding pose.

Constructive Neutral Evolution 20 Years Later

Evolution has led to a great diversity that ranges from elegant simplicity to ornate complexity. Many complex features are often assumed to be more functional or adaptive than their simpler alternatives. However, in 1999, Arlin Stolzfus published a paper in the Journal of Molecular Evolution that outlined a framework in which complexity can arise through a series of non-adaptive steps. He called this framework Constructive Neutral Evolution (CNE). Despite its two-decade-old roots, many evolutionary biologists still appear to be unaware of this explanatory framework for the origins of complexity. In this perspective piece, we explain the theory of CNE and how it changes the order of events in narratives that describe the evolution of complexity. We also provide an extensive list of cellular features that may have become more complex through CNE. We end by discussing strategies to determine whether complexity arose through neutral or adaptive processes.

Aminoacyl-tRNA synthetases

The aminoacyl-tRNA synthetases are an essential and universally distributed family of enzymes that plays a critical role in protein synthesis, pairing tRNAs with their cognate amino acids for decoding mRNAs according to the genetic code. Synthetases help to ensure accurate translation of the genetic code by using both highly accurate cognate substrate recognition and stringent proofreading of noncognate products. While alterations in the quality control mechanisms of synthetases are generally detrimental to cellular viability, recent studies suggest that in some instances such changes facilitate adaption to stress conditions. Beyond their central role in translation, synthetases are also emerging as key players in an increasing number of other cellular processes, with far-reaching consequences in health and disease. The biochemical versatility of the synthetases has also proven pivotal in efforts to expand the genetic code, further emphasizing the wide-ranging roles of the aminoacyl-tRNA synthetase family in synthetic and natural biology.

Protein-Protein Binding as a Two-Step Mechanism: Preselection of Encounter Poses during the Binding of BPTI and Trypsin

Biomolecular recognition between proteins follows complex mechanisms, the understanding of which can substantially advance drug discovery efforts. Here, we track each step of the binding process in atomistic detail with molecular dynamics simulations using trypsin and its inhibitor bovine pancreatic trypsin inhibitor (BPTI) as a model system. We use umbrella sampling to cover a range of unbinding pathways. Starting from these simulations, we subsequently seed classical simulations at different stages of the process and combine them to a Markov state model. We clearly identify three kinetically separated states (an unbound state, an encounter state, and the final complex) and describe the mechanisms that dominate the binding process. From our model, we propose the following sequence of events. The initial formation of the encounter complex is driven by long-range interactions because opposite charges in trypsin and BPTI draw them together. The encounter complex features the prealigned binding partners with binding sites still partially surrounded by solvation shells. Further approaching leads to desolvation and increases the importance of van der Waals interactions. The native binding pose is adopted by maximizing short-range interactions. Thereby side-chain rearrangements ensure optimal shape complementarity. In particular, BPTI’s P1 residue adapts to the S1 pocket and prime site residues reorient to optimize interactions. After the paradigm of conformation selection, binding-competent conformations of BPTI and trypsin are already present in the apo ensembles and their probabilities increase during this proposed two-step association process. This detailed characterization of the molecular forces driving the binding process includes numerous aspects that have been discussed as central to the binding of trypsin and BPTI and protein complex formation in general. In this study, we combine all these aspects into one comprehensive model of protein recognition. We thereby contribute to enhance our general understanding of this fundamental mechanism, which is particularly critical as the development of biopharmaceuticals continuously gains significance.

Breaking a single hydrogen bond in the mitochondrial tRNAPhe -PheRS complex leads to phenotypic pleiotropy of human disease

Various pathogenic variants in both mitochondrial tRNA^Phe and Phenylalanyl‐tRNA synthetase mitochondrial protein coding gene (FARS2) gene encoding for the human mitochondrial PheRS have been identified and associated with neurological and/or muscle‐related pathologies. An important Guanine‐34 (G34)A anticodon mutation associated with myoclonic epilepsy with ragged red fibers (MERRF) syndrome has been reported in hmit‐tRNA^Phe. The majority of G34 contacts in available aaRSs‐tRNAs complexes specifically use that base as an important tRNA identity element. The network of intermolecular interactions providing its specific recognition also largely conserved. However, their conservation depends also on the invariance of the residues in the anticodon binding domain (ABD) of human mitochondrial Phenylalanyl‐tRNA synthetase (hmit‐PheRS). A defect in recognition of the anticodon of tRNA^Phe may happen not only because of G34A mutation, but also due to mutations in the ABD. Indeed, a pathogenic mutation in FARS2 has been recently reported in a 9‐year‐old female patient harboring a p.Asp364Gly mutation. Asp364 is hydrogen bonded (HB) to G34 in WT hmit‐PheRS. Thus, there are two pathogenic variants disrupting HB between G34 and Asp364: one is associated with G34A mutation, and the other with Asp364Gly mutation. We have measured the rates of tRNA^Phe aminoacylation catalyzed by WT hmit‐PheRS and mutant enzymes. These data ranked the residues making a HB with G34 according to their contribution to activity and the signal transduction pathway in the hmit‐PheRS‐tRNA^Phe complex. Furthermore, we carried out extensive MD simulations to reveal the interdomain contact topology on the dynamic trajectories of the complex, and gaining insight into the structural and dynamic integrity effects of hmit‐PheRS complexed with tRNA^Phe.

Database

Structural data are available in PDB database under the accession number(s): 3CMQ, 3TUP, 5MGH, 5MGV.

Solvation Free Energy as a Measure of Hydrophobicity: Application to Serine Protease Binding Interfaces

Solvation and hydrophobicity play a key role in a variety of biological mechanisms. In substrate binding, but also in structure-based drug design, the thermodynamic properties of water molecules surrounding a given protein are of high interest. One of the main algorithms devised in recent years to quantify thermodynamic properties of water is the grid inhomogeneous solvation theory (GIST), which calculates these features on a grid surrounding the protein. Despite the inherent advantages of GIST, the computational demand is a major drawback, as calculations for larger systems can take days or even weeks. Here, we present a GPU accelerated version of the GIST algorithm, which facilitates efficient estimates of solvation free energy even of large biomolecular interfaces. Furthermore, we show that GIST can be used as a reliable tool to evaluate protein surface hydrophobicity. We apply the approach on a set of nine different proteases calculating localized solvation free energies on the surface of the binding interfaces as a measure of their hydrophobicity. We find a compelling agreement with the hydrophobicity of their substrates, i.e., peptides, binding into the binding cleft, and thus our approach provides a reliable description of hydrophobicity characteristics of these biological interfaces.

A Mechanism of Modulating the Direction of Flagellar Rotation in Bacteria by Fumarate and Fumarate Reductase

Fumarate, an electron acceptor in anaerobic respiration of Escherichia coli, has an additional function of assisting the flagellar motor to shift from counterclockwise to clockwise rotation, with a consequent modulation of the bacterial swimming behavior. Fumarate transmits its effect to the motor via the fumarate reductase complex (FrdABCD), shown to bind to FliG-one of the motor's switch proteins. How binding of the FrdABCD respiratory enzyme to FliG enhances clockwise rotation and how fumarate is involved in this activity have remained puzzling. Here we show that the FrdA subunit in the presence of fumarate is sufficient for binding to FliG and for clockwise enhancement. We further demonstrate by in vitro binding assays and super-resolution microscopy in vivo that the mechanism by which fumarate-occupied FrdA enhances clockwise rotation involves its preferential binding to the clockwise state of FliG (FliG_cw). Continuum electrostatics combined with docking analysis and conformational sampling endorsed the experimental conclusions and suggested that the FrdA-FliG_cw interaction is driven by the positive electrostatic potential generated by FrdA and the negatively charged areas of FliG. They further demonstrated that fumarate changes FrdA's conformation to one that can bind to FliG_cw. These findings also show that the reason for the failure of the succinate dehydrogenase flavoprotein SdhA (an almost-identical analog of FrdA shown to bind to FliG equally well) to enhance clockwise rotation is that it has no binding preference for FliG_cw. We suggest that this mechanism is physiologically important as it can modulate the magnitude of ΔG⁰ between the clockwise and counterclockwise states of the motor to tune the motor to the growth conditions of the bacteria.

Electrostatic recognition in substrate binding to serine proteases

Serine proteases of the Chymotrypsin family are structurally very similar but have very different substrate preferences. This study investigates a set of 9 different proteases of this family comprising proteases that prefer substrates containing positively charged amino acids, negatively charged amino acids, and uncharged amino acids with varying degree of specificity. Here, we show that differences in electrostatic substrate preferences can be predicted reliably by electrostatic molecular interaction fields employing customized GRID probes. Thus, we are able to directly link protease structures to their electrostatic substrate preferences. Additionally, we present a new metric that measures similarities in substrate preferences focusing only on electrostatics. It efficiently compares these electrostatic substrate preferences between different proteases. This new metric can be interpreted as the electrostatic part of our previously developed substrate similarity metric. Consequently, we suggest, that substrate recognition in terms of electrostatics and shape complementarity are rather orthogonal aspects of substrate recognition. This is in line with a 2‐step mechanism of protein‐protein recognition suggested in the literature.

Kinetic and structural changes in HsmtPheRS, induced by pathogenic mutations in human FARS2

Mutations in the mitochondrial aminoacyl-tRNA synthetases (mtaaRSs) can cause profound clinical presentations, and have manifested as diseases with very selective tissue specificity. To date most of the mtaaRS mutations could be phenotypically recognized, such that clinicians could identify the affected mtaaRS from the symptoms alone. Among the recently reported pathogenic variants are point mutations in FARS2 gene, encoding the human mitochondrial PheRS. Patient symptoms range from spastic paraplegia to fatal infantile Alpers encephalopathy. How clinical manifestations of these mutations relate to the changes in three-dimensional structures and kinetic characteristics remains unclear, although impaired aminoacylation has been proposed as possible etiology of diseases. Here, we report four crystal structures of HsmtPheRS mutants, and extensive MD simulations for wild-type and nine mutants to reveal the structural changes on dynamic trajectories of HsmtPheRS. Using steady-state kinetic measurements of phenylalanine activation and tRNA^Phe aminoacylation, we gained insight into the structural and kinetic effects of mitochondrial disease-related mutations in FARS2 gene.

BindUP: a web server for non-homology-based prediction of DNA and RNA binding proteins

Gene expression is a multi-step process involving many layers of regulation. The main regulators of the pathway are DNA and RNA binding proteins. While over the years, a large number of DNA and RNA binding proteins have been identified and extensively studied, it is still expected that many other proteins, some with yet another known function, are awaiting to be discovered. Here we present a new web server, BindUP, freely accessible through the website http://bindup.technion.ac.il/, for predicting DNA and RNA binding proteins using a non-homology-based approach. Our method is based on the electrostatic features of the protein surface and other general properties of the protein. BindUP predicts nucleic acid binding function given the proteins three-dimensional structure or a structural model. Additionally, BindUP provides information on the largest electrostatic surface patches, visualized on the server. The server was tested on several datasets of DNA and RNA binding proteins, including proteins which do not possess DNA or RNA binding domains and have no similarity to known nucleic acid binding proteins, achieving very high accuracy. BindUP is applicable in either single or batch modes and can be applied for testing hundreds of proteins simultaneously in a highly efficient manner.

Chimeric human mitochondrial PheRS exhibits editing activity to discriminate nonprotein amino acids

Mitochondria are considered as the primary source of reactive oxygen species (ROS) in nearly all eukaryotic cells during respiration. The harmful effects of these compounds range from direct neurotoxicity to incorporation into proteins producing aberrant molecules with multiple physiological problems. Phenylalanine exposure to ROS produces multiple oxidized isomers: tyrosine, Levodopa, ortho-Tyr, meta-Tyr (m-Tyr), and so on. Cytosolic phenylalanyl-tRNA synthetase (PheRS) exerts control over the translation accuracy, hydrolyzing misacylated products, while monomeric mitochondrial PheRS lacks the editing activity. Recently we showed that "teamwork" of cytosolic and mitochondrial PheRSs cannot prevent incorporation of m-Tyr and l-Dopa into proteins. Here, we present human mitochondrial chimeric PheRS with implanted editing module taken from EcPheRS. The monomeric mitochondrial chimera possesses editing activity, while in bacterial and cytosolic PheRSs this type of activity was detected for the (αβ)2 architecture only. The fusion protein catalyzes aminoacylation of tRNA(Phe) with cognate phenylalanine and effectively hydrolyzes the noncognate aminoacyl-tRNAs: Tyr-tRNA(Phe) and m-Tyr-tRNA(Phe) .

Method for the analysis of contribution of sliding and hopping to a facilitated diffusion of DNA-binding protein: Application to in vivo data

DNA-binding protein searches for its target, a specific site on DNA, by means of diffusion. The search process consists of many recurrent steps of one-dimensional diffusion (sliding) along the DNA chain and three-dimensional diffusion (hopping) after dissociation of a protein from the DNA chain. Here we propose a computational method that allows extracting the contribution of sliding and hopping to the search process in vivo from the measurements of the kinetics of the target search by the lac repressor in Escherichia coli [P. Hammar et al., Science 336, 1595 (2012)]. The method combines lattice Monte Carlo simulations with the Brownian excursion theory and includes explicitly steric constraints for hopping due to the helical structure of DNA. The simulation results including all experimental data reveal that the in vivo target search is dominated by sliding. The short-range hopping to the same base pair interrupts one-dimensional sliding while long-range hopping does not contribute significantly to the kinetics of the search of the target in vivo.

Discovering RNA-protein interactome by using chemical context profiling of the RNA-protein interface

RNA-protein (RNP) interactions generally are required for RNA function. At least 5% of human genes code for RNA-binding proteins. Whereas many approaches can identify the RNA partners for a specific protein, finding the protein partners for a specific RNA is difficult. We present a machine-learning method that scores a protein's binding potential for an RNA structure by utilizing the chemical context profiles of the interface from known RNP structures. Our approach is applicable even when only a single RNP structure is available. We examined 801 mammalian proteins and find that 37 (4.6%) potentially bind transfer RNA (tRNA). Most are enzymes involved in cellular processes unrelated to translation and were not known to interact with RNA. We experimentally tested six positive and three negative predictions for tRNA binding in vivo, and all nine predictions were correct. Our computational approach provides a powerful complement to experiments in discovering new RNPs.

Emergence and evolution

The aminoacyl-tRNA synthetases (aaRSs) are essential components of the protein synthesis machinery responsible for defining the genetic code by pairing the correct amino acids to their cognate tRNAs. The aaRSs are an ancient enzyme family believed to have origins that may predate the last common ancestor and as such they provide insights into the evolution and development of the extant genetic code. Although the aaRSs have long been viewed as a highly conserved group of enzymes, findings within the last couple of decades have started to demonstrate how diverse and versatile these enzymes really are. Beyond their central role in translation, aaRSs and their numerous homologs have evolved a wide array of alternative functions both inside and outside translation. Current understanding of the emergence of the aaRSs, and their subsequent evolution into a functionally diverse enzyme family, are discussed in this chapter.

Thermodynamic properties distinguish human mitochondrial aspartyl-tRNA synthetase from bacterial homolog with same 3D architecture

In the mammalian mitochondrial translation apparatus, the proteins and their partner RNAs are coded by two genomes. The proteins are nuclear-encoded and resemble their homologs, whereas the RNAs coming from the rapidly evolving mitochondrial genome have lost critical structural information. This raises the question of molecular adaptation of these proteins to their peculiar partner RNAs. The crystal structure of the homodimeric bacterial-type human mitochondrial aspartyl-tRNA synthetase (DRS) confirmed a 3D architecture close to that of Escherichia coli DRS. However, the mitochondrial enzyme distinguishes by an enlarged catalytic groove, a more electropositive surface potential and an alternate interaction network at the subunits interface. It also presented a thermal stability reduced by as much as 12°C. Isothermal titration calorimetry analyses revealed that the affinity of the mitochondrial enzyme for cognate and non-cognate tRNAs is one order of magnitude higher, but with different enthalpy and entropy contributions. They further indicated that both enzymes bind an adenylate analog by a cooperative allosteric mechanism with different thermodynamic contributions. The larger flexibility of the mitochondrial synthetase with respect to the bacterial enzyme, in combination with a preserved architecture, may represent an evolutionary process, allowing nuclear-encoded proteins to cooperate with degenerated organelle RNAs.

Comprehensive comparative analysis and identification of RNA-binding protein domains: multi-class classification and feature selection

RNA-protein interaction plays an important role in various cellular processes, such as protein synthesis, gene regulation, post-transcriptional gene regulation, alternative splicing, and infections by RNA viruses. In this study, using Gene Ontology Annotated (GOA) and Structural Classification of Proteins (SCOP) databases an automatic procedure was designed to capture structurally solved RNA-binding protein domains in different subclasses. Subsequently, we applied tuned multi-class SVM (TMCSVM), Random Forest (RF), and multi-class ℓ1/ℓq-regularized logistic regression (MCRLR) for analysis and classifying RNA-binding protein domains based on a comprehensive set of sequence and structural features. In this study, we compared prediction accuracy of three different state-of-the-art predictor methods. From our results, TMCSVM outperforms the other methods and suggests the potential of TMCSVM as a useful tool for facilitating the multi-class prediction of RNA-binding protein domains. On the other hand, MCRLR by elucidating importance of features for their contribution in predictive accuracy of RNA-binding protein domains subclasses, helps us to provide some biological insights into the roles of sequences and structures in protein-RNA interactions.

Crystal structure of human mitochondrial PheRS complexed with tRNA(Phe) in the active "open" state

Monomeric human mitochondrial phenylalanyl-tRNA synthetase (PheRS), or hmPheRS, is the smallest known enzyme exhibiting aminoacylation activity. HmPheRS consists of only two structural domains and differs markedly from heterodimeric eukaryotic cytosolic and bacterial analogs both in the domain organization and in the mode of tRNA binding. Here, we describe the first crystal structure of mitochondrial aminoacyl-tRNA synthetase (aaRS) complexed with tRNA at a resolution of 3.0 Å. Unlike bacterial PheRSs, the hmPheRS recognizes C74, the G1-C72 base pair, and the "discriminator" base A73, proposed to contribute to tRNA(Phe) identity in the yeast mitochondrial enzyme. An interaction of the tRNA acceptor stem with the signature motif 2 residues of hmPheRS is of critical importance for the stabilization of the CCA-extended conformation and its correct placement in the synthetic site of the enzyme. The crystal structure of hmPheRS-tRNA(Phe) provides direct evidence that the formation of the complex with tRNA requires a significant rearrangement of the anticodon-binding domain from the "closed" to the productive "open" state. Global repositioning of the domain is tRNA modulated and governed by long-range electrostatic interactions.

Multistep kinetics of the U1A-SL2 RNA complex dissociation

The U1A-SL2 RNA complex is a model system for studying interactions between RNA and the RNA recognition motif (RRM), which is one of the most common RNA binding domains. We report here kinetic studies of dissociation of the U1A-SL2 RNA complex, using laser temperature jump and stopped-flow fluorescence methods with U1A proteins labeled with the intrinsic chromophore tryptophan. An analysis of the kinetic data suggests three phases of dissociation with time scales of ∼100 μs, ∼50 ms, and ∼2 s. We propose that the first step of dissociation is a fast rearrangement of the complex to form a loosely bound complex. The intermediate step is assigned to be the dissociation of the U1A-SL2 RNA complex, and the final step is assigned to a reorganization of the U1A protein structure into the conformation of the free protein. These assignments are consistent with previous proposals based on thermodynamic, NMR, and surface plasmon resonance experiments and molecular dynamics simulations. Together, these results begin to build a comprehensive model of the complex dynamic processes involved in the formation and dissociation of an RRM-RNA complex.

Analysis of electric moments of RNA-binding proteins: implications for mechanism and prediction

Background

Protein-RNA interactions play important role in many biological processes such as gene regulation, replication, protein synthesis and virus assembly. Although many structures of various types of protein-RNA complexes have been determined, the mechanism of protein-RNA recognition remains elusive. We have earlier shown that the simplest electrostatic properties viz. charge, dipole and quadrupole moments, calculated from backbone atomic coordinates of proteins are biased relative to other proteins, and these quantities can be used to identify DNA-binding proteins. Closely related, RNA-binding proteins are investigated in this study. In particular, discrimination between various types of RNA-binding proteins, evolutionary conservation of these bulk electrostatic features and effect of conformational changes by complex formation are investigated. Basic binding mechanism of a putative RNA-binding protein (HI1333 from Haemophilus influenza) is suggested as a potential application of this study.

Results

We found that similar to DNA-binding proteins (DBPs), RNA-binding proteins (RBPs) also show significantly higher values of electric moments. However, higher moments in RBPs are found to strongly depend on their functional class: proteins binding to ribosomal RNA (rRNA) constitute the only class with all three of the properties (charge, dipole and quadrupole moments) being higher than control proteins. Neural networks were trained using leave-one-out cross-validation to predict RBPs from control data as well as pair-wise classification capacity between proteins binding to various RNA types. RBPs and control proteins reached up to 78% accuracy measured by the area under the ROC curve. Proteins binding to rRNA are found to be best distinguished (AUC = 79%). Changes in dipole and quadrupole moments between unbound and bound structures were small and these properties are found to be robust under complex formation.

Conclusions

Bulk electric moments of proteins considered here provide insights into target recognition by RNA-binding proteins, as well as ability to recognize one type of RBP from others. These results help in understanding the mechanism of protein-RNA recognition, and identifying RNA-binding proteins.

Structure of human cytosolic phenylalanyl-tRNA synthetase: evidence for kingdom-specific design of the active sites and tRNA binding patterns

The existence of three types of phenylalanyl-tRNA synthetase (PheRS), bacterial (alphabeta)(2), eukaryotic/archaeal cytosolic (alphabeta)(2), and mitochondrial alpha, is a prominent example of structural diversity within the aaRS family. PheRSs have considerably diverged in primary sequences, domain compositions, and subunit organizations. Loss of the anticodon-binding domain B8 in human cytosolic PheRS (hcPheRS) is indicative of variations in the tRNA(Phe) binding and recognition as compared to bacterial PheRSs. We report herein the crystal structure of hcPheRS in complex with phenylalanine at 3.3 A resolution. A novel structural module has been revealed at the N terminus of the alpha subunit. It stretches out into the solvent of approximately 80 A and is made up of three structural domains (DBDs) possessing DNA-binding fold. The dramatic reduction of aminoacylation activity for truncated N terminus variants coupled with structural data and tRNA-docking model testify that DBDs play crucial role in hcPheRS activity.

Orientation and distance dependent chiral discrimination in the first step of the aminoacylation reaction: integrated molecular orbital and semi-empirical method (ONIOM) based calculation

Aminoacylation is a vital step in natural biosynthesis process of peptide and is the key step in correlating the realm of protein with the RNA world. Incorrect aminoacylation might lead to misacylation of d-amino acid in the tRNA which might cause synthesis of a hetero-peptide rather than natural homopeptide leading to the altered functionality of the peptide. However, the accuracy of this process is remarkable and leads to the attachment of the correct enantiomer of the amino acid with their cognate tRNA. Thus, the chiral discrimination is stringent. In the present work, we presented a combined ONIOM (ab initio/semi-empirical) study of the chiral discrimination in the first step of aminoacylation reaction based on a model of crystal structure of the oligomeric complex of histidyl-tRNA synthetase (HisRS) from Escherichia coli complexed with ATP and histidinol and histidyl-adenylate. The study reveals that the molecular mechanism of the chiral discrimination involves the amino acid, ATP as well as surrounding residues of the synthetase. Several factors are noted to be responsible for discrimination and explain the high level of stereospecificity of the process. The chirality of the amino acid of the substrate and its (principally) electrostatic interaction with the ATP is important for discrimination. The distance and orientational changes involved in the approach of the d-His towards the ATP is energetically unfavorable. The charge distributions on the His and ATP are important for the discrimination. Removal of the charges in the model drastically reduces the discrimination. Restricted nature of the mutual orientation within the cavity of the active site where the His and ATP are located during the change in orientation for the approach to form the adenylate makes the resultant interaction profile as different for l-His and d-His also influences chiral discrimination. The analysis of the transition state structure revealed that alteration of the chirality of the His destabilize the transition state by removing the favorable electrostatic interaction between the Glu-83 and NH(3)(+) group of the His substrate. The proximity of the surrounding residues as present in the active site of the synthetase with the His and ATP (the separation is of nanometer range) has influence of discrimination. The study provides a molecular mechanism of the retention of biological homochirality.

Intra-protein compensatory mutations analysis highlights the tRNA recognition regions in aminoacyl-tRNA synthetases

The aminoacyl-tRNA synthetases (aaRSs) covalently attach amino acids to their corresponding nucleic acid adapter molecules, tRNAs. The interactions in the tRNA-aaRSs complexes are mostly non-specific, and largely electrostatic. Tracing a way of aaRS-tRNA mutual adaptation throughout evolution offers a clearer view of understanding how aaRS-tRNA systems preserve patterns of tRNA recognition and binding. In this study, we used the compensatory mutations analysis to explore adaptation of aaRSs in respond to random mutations that can occur in the tRNA-recognition area. We showed that the frequency of compensatory mutations among residues that belong to the recognition region is 1.75-fold higher than that of the exposed residues. The highest frequencies of compensatory mutations are observed for pairs of charged residues, wherein one residue is located within the tRNA-recognition area, while the second is placed outside of the area, and contributes to the formation of the aaRS electrostatic landscape. Given charged residues are compensated by buried charge residues in more than 60% of the analyzed mutations. The cytoplasmatic and mitochondrial aaRSs preserve similar patterns of compensatory mutations in the tRNA recognition areas. Moreover, we found that mitochondrial aaRSs demonstrate a significant increase in the frequency of compensatory mutations in the area. Our findings shed light on the physical nature of compensatory mutations in aaRSs, thereby keeping unchanged tRNA-recognition patterns.

Idiosyncratic helix-turn-helix motif in Methanosarcina barkeri seryl-tRNA synthetase has a critical architectural role

All seryl-tRNA synthetases (SerRSs) are functional homodimers with a C-terminal active site domain typical for class II aminoacyl-tRNA synthetases and an N-terminal domain involved in tRNA binding. The recently solved three-dimensional structure of Methanosarcina barkeri SerRS revealed the idiosyncratic features of methanogenic-type SerRSs; that is, an active site zinc ion, a unique tRNA binding domain, and an insertion of approximately 30 residues in the catalytic domain, which adopt a helix-turn-helix (HTH) fold. Here, we present biochemical evidence for multiple roles of the HTH motif; it is important for dimerization of the enzyme, contributes to the overall stability, and is critical for the proper positioning of the tRNA binding domain relative to the catalytic domain. The changes in intrinsic fluorescence during denaturation of the wild-type M. barkeri SerRS and of the mutated variant lacking the HTH motif combined with cross-linking and gel analysis of protein subunits during various stages of the unfolding process revealed significantly reduced stability of the mutant dimers. In vitro kinetic analysis of enzymes, mutated in one of the N-terminal helices and the HTH motif, shows impaired tRNA binding and aminoacylation and emphasizes the importance of this domain for the overall architecture of the enzyme. The role of the idiosyncratic HTH motif in dimer stabilization and association between the catalytic and tRNA binding domain has been additionally confirmed by a yeast two-hybrid approach. Furthermore, we provide experimental evidence that tRNA binds across the dimer.

Toward a more complete view of tRNA biology

Transfer RNAs are ancient molecules present in all domains of life. In addition to translating the genetic code into protein and defining the second genetic code together with aminoacyl-tRNA synthetases, tRNAs act in many other cellular functions. Robust phenomenological observations on the role of tRNAs in translation, together with massive sequence and crystallographic data, have led to a deeper physicochemical understanding of tRNA architecture, dynamics and identity. In vitro studies complemented by cell biology data already indicate how tRNA behaves in cellular environments, in particular in higher Eukarya. From an opposite approach, reverse evolution considerations suggest how tRNAs emerged as simplified structures from the RNA world. This perspective discusses what basic questions remain unanswered, how these answers can be obtained and how a more rational understanding of the function and dysfunction of tRNA can have applications in medicine and biotechnology.

The tRNA-induced conformational activation of human mitochondrial phenylalanyl-tRNA synthetase

All class II aminoacyl-tRNA synthetases (aaRSs) are known to be active as functional homodimers, homotetramers, or heterotetramers. However, multimeric organization is not a prerequisite for phenylalanylation activity, as monomeric mitochondrial phenylalanyl-tRNA synthetase (PheRS) is also active. We herein report the structure, at 2.2 A resolution, of a human monomeric mitPheRS complexed with Phe-AMP. The smallest known aaRS, which is, in fact, 1/5 of a cytoplasmic analog, is a chimera of the catalytic module of the alpha and anticodon binding domain (ABD) of the bacterial beta subunit of (alphabeta)2 PheRS. We demonstrate that the ABD located at the C terminus of mitPheRS overlaps with the acceptor stem of phenylalanine transfer RNA (tRNAPhe) if the substrate is positioned in a manner similar to that seen in the binary Thermus thermophilus complex. Thus, formation of the PheRS-tRNAPhe complex in human mitochondria must be accompanied by considerable rearrangement (hinge-type rotation through approximately 160 degrees) of the ABD upon tRNA binding.

Classifying RNA-binding proteins based on electrostatic properties

Protein structure can provide new insight into the biological function of a protein and can enable the design of better experiments to learn its biological roles. Moreover, deciphering the interactions of a protein with other molecules can contribute to the understanding of the protein's function within cellular processes. In this study, we apply a machine learning approach for classifying RNA-binding proteins based on their three-dimensional structures. The method is based on characterizing unique properties of electrostatic patches on the protein surface. Using an ensemble of general protein features and specific properties extracted from the electrostatic patches, we have trained a support vector machine (SVM) to distinguish RNA-binding proteins from other positively charged proteins that do not bind nucleic acids. Specifically, the method was applied on proteins possessing the RNA recognition motif (RRM) and successfully classified RNA-binding proteins from RRM domains involved in protein–protein interactions. Overall the method achieves 88% accuracy in classifying RNA-binding proteins, yet it cannot distinguish RNA from DNA binding proteins. Nevertheless, by applying a multiclass SVM approach we were able to classify the RNA-binding proteins based on their RNA targets, specifically, whether they bind a ribosomal RNA (rRNA), a transfer RNA (tRNA), or messenger RNA (mRNA). Finally, we present here an innovative approach that does not rely on sequence or structural homology and could be applied to identify novel RNA-binding proteins with unique folds and/or binding motifs.

Gene expression in all living organisms is regulated by a complex set of events at both transcriptional and posttranscriptional levels. RNA-binding proteins play a key role in posttranscriptional events including splicing, stability, transport, and translation. Nowadays, there is increasing evidence that many other cellular processes may be mediated by RNA. Identifying new proteins involved in interaction with RNA is thus essential to unraveling the cellular processes in which these interactions are involved. In the current study we present a successful computational approach for classifying RNA-binding proteins and distinguishing them from other proteins based on structural and electrostatic properties. We test the method on a unique protein domain, the RNA recognition motif (RRM), which mediates both RNA and protein interactions. We show that we can discriminate RNA-binding RRMs from protein-binding RRMs. Further, we demonstrate that we can classify known RNA-binding proteins based on their RNA target (mRNA, rRNA, or tRNA). Our method does not rely on any kind of evolutionary information and thus can be applied to identify RNA-binding proteins with novel modes of RNA recognition.

Radial scan of the molecular electrostatic potential of RNA double helices: an application to the enzyme-tRNA recognition

We introduced a method to characterize quantitatively the molecular electrostatic potential (MEP) of the minor and major grooves of base pairs located at nucleic acid double helices. By means of a radial MEP scan, we obtained a n-tuple of potential values corresponding to each groove, which can be analyzed by plotting the MEP values as a function of the angle in the radial scan. We studied base pairs of two different tRNAs, relevant in the recognition process with their cognate aminoacyl tRNA synthetases (aaRSs), in order to correlate their electrostatic behavior with the corresponding aminoacylation activity. We analyzed the first three base pairs of the Escherichia coli tRNA(Ala) acceptor stem, finding several cases where the MEP profiles obtained from the plots are in agreement with the reported aminoacylation activities. Additionally, a non-hierarchical clustering performed over the MEP n-tuples resulted in meaningful classifications that correlate with the activity and with the predicted stereochemistry of the reaction. We also studied the first two base pairs of the E. coli tRNA(Thr) acceptor stem but constraining the analysis to the angle intervals that seem relevant for the binding sites of the enzyme. These intervals were deduced from the ThrRS-tRNA(Thr) complex crystal structure. In this case, we also found a good agreement between the MEP profiles and the activity, supporting the idea that the tRNA identity elements function is to allow an optimal electrostatic complementarity between the aminoacyl-tRNA synthetase and the tRNA.

Interaction of aminoacyl-tRNA synthetases with tRNA: general principles and distinguishing characteristics of the high-molecular-weight substrate recognition

This review summarizes results of numerous (mainly functional) studies that have been accumulated over recent years on the problem of tRNA recognition by aminoacyl-tRNA synthetases. Development and employment of approaches that use synthetic mutant and chimeric tRNAs have demonstrated general principles underlying highly specific interaction in different systems. The specificity of interaction is determined by a certain number of nucleotides and structural elements of tRNA (constituting the set of recognition elements or specificity determinants), which are characteristic of each pair. Crystallographic structures available for many systems provide the details of the molecular basis of selective interaction. Diversity and identity of biochemical functions of the recognition elements make substantial contribution to the specificity of such interactions.

Why are proteins charged? Networks of charge-charge interactions in proteins measured by charge ladders and capillary electrophoresis

Almost all proteins contain charged amino acids. While the function in catalysis or binding of individual charges in the active site can often be identified, it is less clear how to assign function to charges beyond this region. Are they necessary for solubility? For reasons other than solubility? Can manipulating these charges change the properties of proteins? A combination of capillary electrophoresis (CE) and protein charge ladders makes it possible to study the roles of charged residues on the surface of proteins outside the active site. This method involves chemical modification of those residues to generate a large number of derivatives of the protein that differ in charge. CE separates those derivatives into groups with the same number of modified charged groups. By studying the influence of charge on the properties of proteins using charge ladders, it is possible to estimate the net charge and hydrodynamic radius and to infer the role of charged residues in ligand binding and protein folding.

The role of positively charged amino acids and electrostatic interactions in the complex of U1A protein and U1 hairpin II RNA

Previous kinetic investigations of the N-terminal RNA recognition motif (RRM) domain of spliceosomal protein U1A, interacting with its RNA target U1 hairpin II, provided experimental evidence for a ‘lure and lock’ model of binding in which electrostatic interactions first guide the RNA to the protein, and close range interactions then lock the two molecules together. To further investigate the ‘lure’ step, here we examined the electrostatic roles of two sets of positively charged amino acids in U1A that do not make hydrogen bonds to the RNA: Lys20, Lys22 and Lys23 close to the RNA-binding site, and Arg7, Lys60 and Arg70, located on ‘top’ of the RRM domain, away from the RNA. Surface plasmon resonance-based kinetic studies, supplemented with salt dependence experiments and molecular dynamics simulation, indicate that Lys20 predominantly plays a role in association, while nearby residues Lys22 and Lys23 appear to be at least as important for complex stability. In contrast, kinetic analyses of residues away from the RNA indicate that they have a minimal effect on association and stability. Thus, well-positioned positively charged residues can be important for both initial complex formation and complex maintenance, illustrating the multiple roles of electrostatic interactions in protein–RNA complexes.