Structural modelling of the complex of leucyl-tRNA synthetase and mis-aminoacylated tRNA

To assure fidelity of translation, class Ia aminoacyl-tRNA synthetases (aaRSs) edit mis-aminoacylated tRNAs. Mis-attached amino acids and structural water molecules are not included simultaneously in the current crystal structures of the aaRS•tRNA complexes, where the 3′-ends (adenine 76; A76) are bound to the editing sites. A structural model of the completely solvated leucyl-tRNA synthetase complexed with valyl-tRNA^Leu was constructed by exploiting molecular dynamics simulations modified for the present modelling. The results showed that the ribose conformation of A76 is distinct from those observed in the above-mentioned crystal structures, which could be derived from structural constraints in a sandwiched manner induced by the mis-attached valine and tRNA^Leu.     

Exit Strategies for Charged tRNA from GluRS

For several class I aminoacyl-tRNA synthetases (aaRSs), the rate-determining step in aminoacylation is the dissociation of charged tRNA from the enzyme. In this study, the following factors affecting the release of the charged tRNA from aaRSs are computationally explored: the protonation states of amino acids and substrates present in the active site, and the presence and the absence of AMP and elongation factor Tu.Through molecular modeling, internal pK_a calculations, and molecular dynamics simulations, distinct, mechanistically relevant post-transfer states with charged tRNA bound to glutamyl-tRNA synthetase from Thermus thermophilus (Glu-tRNA^Glu) are considered. The behavior of these nonequilibrium states is characterized as a function of time using dynamical network analysis, local energetics, and changes in free energies to estimate transitions that occur during the release of the tRNA. The hundreds of nanoseconds of simulation time reveal system characteristics that are consistent with recent experimental studies.Energetic and network results support the previously proposed mechanism in which the transfer of amino acid to tRNA is accompanied by the protonation of AMP to H-AMP. Subsequent migration of proton to water reduces the stability of the complex and loosens the interface both in the presence and in the absence of AMP. The subsequent undocking of AMP or tRNA then proceeds along thermodynamically competitive pathways. Release of the tRNA acceptor stem is further accelerated by the deprotonation of the α-ammonium group on the charging amino acid. The proposed general base is Glu41, a residue binding the α-ammonium group that is conserved in both structure and sequence across nearly all class I aaRSs. This universal handle is predicted through pK_a calculations to be part of a proton relay system for destabilizing the bound charging amino acid following aminoacylation. Addition of elongation factor Tu to the aaRS·tRNA complex stimulates the dissociation of the tRNA core and the tRNA acceptor stem.     

The Role of L1 Stalk-tRNA Interaction in the Ribosome Elongation Cycle

The ribosomal L1 stalk is a mobile structure implicated in directing tRNA movement during translocation through the ribosome. This article investigates three aspects of L1 stalk-tRNA interaction. First, by combining data from cryo electron microscopy, X-ray crystallography, and molecular dynamics simulations through the molecular dynamics flexible fitting method, we obtained atomic models of different tRNAs occupying the hybrid P/E state interacting with the L1 stalk. These models confirm the assignment of fluorescence resonance energy transfer states from previous single-molecule investigations of L1 stalk dynamics. Second, the models reconcile how initiator tRNA^fMet interacts less strongly with the L1 stalk compared to elongator tRNA^Phe, as seen in previous single-molecule experiments. Third, results from a simulation of the entire ribosome in which the L1 stalk is moved from a half-closed conformation to its open conformation are found to support the hypothesis that L1 stalk opening is involved in tRNA release from the ribosome.     

Classical molecular dynamics simulation of seryl tRNA synthetase and threonyl tRNA synthetase bound with tRNA and aminoacyl adenylate

Lacunae of understanding exist concerning the active site organization during the charging step of the aminoacylation reaction. We present here a molecular dynamics simulation study of the dynamics of the active site organization during charging step of subclass IIa dimeric SerRS from Thermus thermophilus (^ttSerRS) bound with ^tttRNA^Ser and dimeric ThrRS from Escherichia coli (^ecThrRS) bound with ^ectRNA^Thr. The interactions between the catalytically important loops and tRNA contribute to the change in dynamics of tRNA in free and bound states, respectively. These interactions help in the development of catalytically effective organization of the active site. The A76 end of the ^tttRNA^Ser exhibits fast dynamics in free State, which is significantly slowed down within the active site bound with adenylate. The loops change their conformation via multimodal dynamics (a slow diffusive mode of nanosecond time scale and fast librational mode of dynamics in picosecond time scale). The active site residues of the motif 2 loop approach the proximal bases of tRNA and adenylate by slow diffusive motion (in nanosecond time scale) and make conformational changes of the respective side chains via ultrafast librational motion to develop precise hydrogen bond geometry. Presence of bound Mg²⁺ ions around tRNA and dynamically slow bound water are other common features of both aaRSs. The presence of dynamically rigid Zinc ion coordination sphere and bipartite mode of recognition of ^ectRNA^Thr are observed.     

Open interface and large quaternary structure movements in 3D domain swapped proteins: insights from molecular dynamics simulations of the C-terminal swapped dimer of ribonuclease A

Bovine pancreatic ribonuclease (RNase A) forms two three-dimensional (3D) domain swapped dimers. Crystallographic investigations have revealed that these dimers display completely different quaternary structures: one dimer (N-dimer), which presents the swapping of the N-terminal helix, is characterized by a compact structure, whereas the other (C-dimer), which is stabilized by the exchange of the C-terminal end, shows a rather loose assembly of the two subunits. The dynamic properties of monomeric RNase A and of the N-dimer have been extensively characterized. Here, we report a molecular dynamics investigation carried out on the C-dimer. This computational experiment indicates that the quaternary structure of the C-dimer undergoes large fluctuations. These motions do not perturb the proper folding of the two subunits, which retain the dynamic properties of RNase A and the N-dimer. Indeed, the individual subunits of the C-dimer display the breathing motion of the beta-sheet structure, which is important for the enzymatic activity of pancreatic-like ribonucleases. In contrast to what has been observed for the N-dimer, the breathing motion of the two subunits of the C-dimer is not coupled. This finding suggests that the intersubunit communications in a 3D domain swapped dimer strongly rely on the extent of the interchain interface. Furthermore, the observation that the C-dimer is endowed with a high intrinsic flexibility holds interesting implications for the specific properties of 3D domain swapped dimers. Indeed, a survey of the quaternary structures of the other 3D domain swapped dimers shows that large variations are often observed when the structural determinations are conducted in different experimental conditions. The 3D domain swapping phenomenon coupled with the high flexibility of the quaternary structure may be relevant for protein-protein recognition, and in particular for the pathological aggregations.     

Molecular dynamics simulations on SDF-1alpha: binding with CXCR4 receptor

Insights into the interacting mode of CXCR4 with SDF-1alpha are crucial in understanding the structural and functional characteristics of CXCR4 receptor. In this paper a computational pipeline, integrating protein structure prediction, molecular dynamics simulations, automated molecular docking, and Brownian dynamics simulations were employed to investigate the dynamic and energetic aspects of CXCR4 associating with SDF-1alpha. The entire simulation revealed the surface distribution feature of electrostatic potentials and conformational "open-close" process of the receptor. The possible binding conformation of CXCR4 was identified, and the CXCR4-SDF-1alpha binding complex was generated. Arg188-Glu277 salt bridge plays an important role for both the extracellular domain conformational change and SDF-1alpha binding. Two binding sites were mapped at the extracellular domain (Site 1) and inside the transmembrane domain (Site 2), which are composed of conserved residues. Sites 1 and 2 contribute approximately 60% and 40% to the binding affinity with SDF-1alpha, respectively. The binding model is in agreement with most of the experimental data. Transmembrane VI has more significant motion in the harmonious conformational transition of CXCR4 during SDF-1alpha binding, which may be possibly associated with signal transduction. Based on the modeling and simulation, a binding mechanism hypothesis between CXCR4 and SDF-1alpha and its relationship to the signal transduction has been proposed.     

Low-mass molecular dynamics simulation: A simple and generic technique to enhance configurational sampling

CLN025 is one of the smallest fast-folding proteins. Until now it has not been reported that CLN025 can autonomously fold to its native conformation in a classical, all-atom, and isothermal–isobaric molecular dynamics (MD) simulation. This article reports the autonomous and repeated folding of CLN025 from a fully extended backbone conformation to its native conformation in explicit solvent in multiple 500-ns MD simulations at 277 K and 1 atm with the first folding event occurring as early as 66.1 ns. These simulations were accomplished by using AMBER forcefield derivatives with atomic masses reduced by 10-fold on Apple Mac Pros. By contrast, no folding event was observed when the simulations were repeated using the original AMBER forcefields of FF12SB and FF14SB. The results demonstrate that low-mass MD simulation is a simple and generic technique to enhance configurational sampling. This technique may propel autonomous folding of a wide range of miniature proteins in classical, all-atom, and isothermal–isobaric MD simulations performed on commodity computers—an important step forward in quantitative biology.     

Molecular recognition of Cullin3 by KCTDs: Insights from experimental and computational investigations

Recent investigations have highlighted a key role of the proteins of the KCTD (K-potassium channel tetramerization domain containing proteins) family in several fundamental biological processes. Despite the growing importance of KCTDs, our current understanding of their biophysical and structural properties is very limited. Biochemical characterizations of these proteins have shown that most of them act as substrate adaptor in E3 ligases during protein ubiquitination. Here we present a characterization of the KCTD5-Cullin3 interactions which are mediated by the KCTD5 BTB domain. Isothermal titration calorimetry experiments reveal that KCTD5 avidly binds the Cullin3 (Cul3). The complex presents a 5:5 stoichiometry and a dissociation constant of 59 nM. Molecular modeling and molecular dynamics simulations clearly indicate that the two proteins form a stable (KCTD5–Cul3)₅ pinwheel-shaped heterodecamer in which two distinct KCTD5 subunits cooperate in the binding of each cullin chain. Molecular dynamics simulations indicate that different types of interactions contribute to the stability of the assembly. Interestingly, residues involved in Cul3 recognitions are conserved in the KCTD5 orthologs and paralogs implicated in important biological processes. These residues are also rather well preserved in most of the other KCTD proteins. By using molecular modeling techniques, the entire ubiquitination system including the E3 ligase, the E2 conjugating enzyme and ubiquitin was generated. The analysis of the molecular architecture of this complex machinery provides insights into the ubiquitination processes which involve E3 ligases with a high structural complexity.     

Structure, dynamics, and stability of assemblies of the human prion fragment SNQNNF

Misfolding of the prion protein (PrP) is associated with the development of Transmissible Spongiform Encephalopathies. The recent crystal structure of ‘steric zipper’ aggregates of the peptide SNQNNF (human PrP fragment 170-175) has highlighted its potential involvement in the misfolding process. A detailed molecular dynamics investigation on SNQNNF aggregates has been performed to analyze the behavior of the assemblies in a non-crystalline context. Stability, dynamics, and structural features suggest that SNQNNF assemblies are very good candidates to be involved in the structure of PrP fibrils. In addition, the analysis of small aggregates shows that steric zipper interfaces are able to stabilize assemblies composed of four strands per sheet. Altogether, the present findings indicate that steric zipper may play a key role in prion diseases. This suggestion is also corroborated by MD analyses of point mutations within the region 170-175.     

Direct detection of the myosin super-relaxed state and interacting-heads motif in solution

The interacting-heads motif (IHM) is a structure of myosin that has been proposed to modulate cardiac output by occluding myosin molecules from undergoing the force-generating cycle. It is hypothesized to be the structural basis for the super-relaxed state (SRX), a low-ATPase kinetic state thought to be cardioprotective. The goal of the present study was to test this hypothesis by determining directly and quantitatively the fractions of myosin in the IHM and SRX under the same conditions in solution. To detect the structural IHM, we used time-resolved fluorescence resonance energy transfer to quantitate two distinct populations. One population was observed at a center distance of 2.0 nm, whereas the other was not detectable by fluorescence resonance energy transfer, implying a distance greater than 4 nm. We confirmed the IHM assignment to the 2.0-nm population by applying the same cross-linking protocol used previously to image the IHM by electron microscopy. Under the same conditions, we also measured the fraction of myosin in the SRX using stopped-flow kinetics. Our results show that the populations of SRX and IHM myosin were similar, unless treated with mavacamten, a drug that recently completed phase III clinical trials to treat hypertrophic cardiomyopathy and is proposed to act by stabilizing both the SRX and IHM. However, we found that mavacamten had a much greater effect on the SRX (55% increase) than on the IHM (4% increase). We conclude that the IHM structure is sufficient but not necessary to produce the SRX kinetic state.     

Human trans-editing enzyme displays tRNA acceptor-stem specificity and relaxed amino acid selectivity

Accurate translation of genetic information into proteins is vital for cell sustainability. ProXp-ala prevents proteome-wide Pro-to-Ala mutations by hydrolyzing misacylated Ala-tRNA^Pro, which is synthesized by prolyl-tRNA synthetase. Bacterial ProXp-ala was previously shown to combine a size-based exclusion mechanism with conformational and chemical selection for the recognition of the alanyl moiety, whereas tRNA^Pro is selected via recognition of tRNA acceptor-stem elements G72 and A73. The identity of these critical bases changed during evolution with eukaryotic cytosolic tRNA^Pro possessing a cytosine at the corresponding positions. The mechanism by which eukaryotic ProXp-ala adapted to these changes remains unknown. In this work, recognition of the aminoacyl moiety and tRNA acceptor stem by human (Homo sapiens, or Hs) ProXp-ala was examined. Enzymatic assays revealed that Hs ProXp-ala requires C72 and C73 in the context of Hs cytosolic tRNA^Pro for efficient deacylation of mischarged Ala-tRNA^Pro. The strong dependence on these bases prevents cross-species deacylation of bacterial Ala-tRNA^Pro or of Hs mitochondrial Ala-tRNA^Pro by the human enzyme. Similar to the bacterial enzyme, Hs ProXp-ala showed strong tRNA acceptor-stem recognition but differed in its amino acid specificity profile relative to bacterial ProXp-ala. Changes at conserved residues in both the Hs and bacterial ProXp-ala substrate-binding pockets modulated this specificity. These results illustrate how the mechanism of substrate selection diverged during the evolution of the ProXp-ala family, providing the first example of a trans-editing domain whose specificity evolved to adapt to changes in its tRNA substrate.     

Structural Events in a Bacterial Uniporter Leading to Translocation of Glucose to the Cytosol

Among classes of sugar transporters, there exists a comparatively new family of transporters named SWEET transporters (semi-SWEETS in bacteria) that are uniport transmembrane proteins. It is hypothesized that sugar is transported from the extracellular side (via outward-open state) to intracellular side (inward-open state) through intermediate occluded state (both extracellular and intracellular gates closed). In our study, extensive unbiased all-atom molecular dynamics simulations were carried out with the outward-open and inward-open conformations to study this transition mechanism. We find that after 100?ns, the outward-open structure without sugar bound starts changing to the occluded form leading to closure of extracellular gates stabilized by electrostatic and hydrophobic interactions. Further simulations (up to 7?μs) have led to a transition toward the inward-open form and suggest that there exists more than one intermediate occluded conformation. We have also performed 5-μs simulations on the glucose-docked structure to identify the putative substrate-bound translocation pathway. Glucose binds to semi-SWEET with strong hydrogen bonds to Asn66 and Trp50. Comparative simulations of substrate bound, and unbound forms suggested that glucose, the putative substrate, facilitates relatively rapid conformational transitions. For the first time, we captured the release of glucose to the cytosol, in this family of transporters. We find that prior to release of glucose, the glucose forms interactions with polar residues near the intracellular gate which may facilitate its release. The distance between the residues Asn31 and Gly34 of the other protomer was found to play a decisive role in the transport of glucose to the cytoplasmic side.     

Coupling and dynamics of subunits in the hexameric AAA+ chaperone ClpB

The bacterial AAA+ protein ClpB and its eukaryotic homologue Hsp104 ensure thermotolerance of their respective organisms by reactivating aggregated proteins in cooperation with the Hsp70/Hsp40 chaperone system. Like many members of the AAA+ superfamily, the ClpB protomers form ringlike homohexameric complexes. The mechanical energy necessary to disentangle protein aggregates is provided by ATP hydrolysis at the two nucleotide-binding domains of each monomer. Previous studies on ClpB and Hsp104 show a complex interplay of domains and subunits resulting in homotypic and heterotypic cooperativity. Using mutations in the Walker A and Walker B nucleotide-binding motifs in combination with mixing experiments we investigated the degree of inter-subunit coupling with respect to different aspects of the ClpB working cycle. We find that subunits are tightly coupled with regard to ATPase and chaperone activity, but no coupling can be observed for ADP binding. Comparison of the data with statistical calculations suggests that for double Walker mutants, approximately two in six subunits are sufficient to abolish chaperone and ATPase activity completely. In further experiments, we determined the dynamics of subunit reshuffling. Our results show that ClpB forms a very dynamic complex, reshuffling subunits on a timescale comparable to steady-state ATP hydrolysis. We propose that this could be a protection mechanism to prevent very stable aggregates from becoming suicide inhibitors for ClpB.     

Ensemble docking and molecular dynamics identify knoevenagel curcumin derivatives with potent anti-EGFR activity

Epidermal growth factor receptor tyrosine kinase (EGFR-TK) is an attractive target for cancer therapy. Despite a number of effective EGFR inhibitors that are constantly expanding and different methods being employed to obtain novel compounds, the search for newer EGFR inhibitors is still a major scientific challenge. In the present study, a molecular docking and molecular dynamics investigation has been carried out with an ensemble of EGFR-TK structures against a synthetically feasible library of curcumin analogs to discover potent EGFR inhibitors. To resolve protein flexibility issue we have utilized 5 EGFR wild type crystal structures during docking as this gives improved possibility of identifying an active compound as compared to using a single crystal structure. We then identified five curcumin analogs representing different scaffolds that can serve as lead molecules. Finally, the 5 ns molecular dynamics simulation shows that knoevenagel condensate of curcumin specifically C29 and C30 can be used as starting blocks for developing effective leads capable of inhibiting EGFR.     

Propagation of dynamic changes in barnase upon binding of barstar: an NMR and computational study

NMR spectroscopy and computer simulations were used to examine changes in chemical shifts and in dynamics of the ribonuclease barnase that result upon binding to its natural inhibitor barstar. Although the spatial structures of free and bound barnase are very similar, binding results in changes of the dynamics of both fast side-chains, as revealed by (2)H relaxation measurements, and NMR chemical shifts in an extended beta-sheet that is located far from the binding interface. Both side-chain dynamics and chemical shifts are sensitive to variations in the ensemble populations of the inter-converting molecular states, which can escape direct structural observation. Molecular dynamics simulations of free barnase and barnase in complex with barstar, as well as a normal mode analysis of barnase using a Gaussian network model, reveal relatively rigid domains that are separated by the extended beta-sheet mentioned above. The observed changes in NMR parameters upon ligation can thus be rationalized in terms of changes in inter-domain dynamics and in populations of exchanging states, without measurable structural changes. This provides an alternative model for the propagation of a molecular response to ligand binding across a protein that is based exclusively on changes in dynamics.     

在活细胞中应用FRET技术研究TGFβ/TβRI信号传导

转移生长因子β(TGFβ)信号传导通路参与调节细胞的增殖、分化、凋亡、细胞迁移等一系列细胞过程,与骨代谢疾病的发病机制密切相关.本研究根据荧光共振能量转移(FRET)技术原理,构建包含CFP-TβRI-YFP融合蛋白的TβRI生物传感器,转染293T细胞,观察转染效率.以TGFβ1为诱导剂,激活TGFβ/TβRI信号传导通路,在活细胞生理条件下,动态监测TβRI生物传感器的FRET效应.结果表明,成功构建了TβRI生物传感器,转染细胞效率达50%,在TGFβ1诱导刺激6min后,FRET效率增加并达到最大值,该过程经历9 min后,随时间的延长,FRET效率下降.研究结果表明:在活细胞生理条件下,TGFβ1/TβRI信号转导过程存在一定的时间特异性.     

Single-Molecule Studies of the Im7 Folding Landscape

Under appropriate conditions, the four-helical Im7 (immunity protein 7) folds from an ensemble of unfolded conformers to a highly compact native state via an on-pathway intermediate. Here, we investigate the unfolded, intermediate, and native states populated during folding using diffusion single-pair fluorescence resonance energy transfer by measuring the efficiency of energy transfer (or proximity or P ratio) between pairs of fluorophores introduced into the side chains of cysteine residues placed in the center of helices 1 and 4, 1 and 3, or 2 and 4. We show that while the native states of each variant give rise to a single narrow distribution with high P values, the distributions of the intermediates trapped at equilibrium (denoted I^eqm) are fitted by two Gaussian distributions. Modulation of the folding conditions from those that stabilize the intermediate to those that destabilize the intermediate enabled the distribution of lower P value to be assigned to the population of the unfolded ensemble in equilibrium with the intermediate state. The reduced stability of the I^eqm variants allowed analysis of the effect of denaturant concentration on the compaction and breadth of the unfolded state ensemble to be quantified from 0 to 6 M urea. Significant compaction is observed as the concentration of urea is decreased in both the presence and absence of sodium sulfate, as previously reported for a variety of proteins. In the presence of Na₂SO₄ in 0 M urea, the P value of the unfolded state ensemble approaches that of the native state. Concurrent with compaction, the ensemble displays increased peak width of P values, possibly reflecting a reduction in the rate of conformational exchange among iso-energetic unfolded, but compact conformations. The results provide new insights into the initial stages of folding of Im7 and suggest that the unfolded state is highly conformationally constrained at the outset of folding.     

The importance of G.A hydrogen bonding in the metal ion- and protein-induced folding of a kink turn RNA

The kink turn (K-turn) is a common motif in RNA structure, found in many RNA species important in translation, RNA modification and splicing, and the control of gene expression. In general the K-turn comprises a three nucleotide bulge followed by trans sugar-Hoogsteen G·A pairs. The RNA adopts a tightly kinked conformation, and is a common target for binding proteins, exemplified by the L7Ae family. We have measured the rates of association and dissociation for the binding of L7Ae to the Kt-7 kink turn, from which we calculate an affinity of K_D = 10 pM. This high affinity is consistent with the role of this binding as the first stage in the assembly of key functional nucleoproteins such as box C/D snoRNP. Kink-turn RNA undergoes a two-state transition between the kinked conformation, and a more extended structure, and folding into the kinked form is induced by divalent metal ions, or by binding of proteins of the L7Ae class. The K-turn provides an excellent, simple model for RNA folding, which can be dissected at the atomic level. We have analyzed the contributions of the hydrogen bonds that form the G·A pairs to the ion- and protein-induced folding of the K-turn. We find that all four hydrogen bonds are important to the stability of the kinked form of the RNA, and we can now define all the important hydrogen bonding interactions that stabilize the K-turn. The high affinity of L7Ae binding is coupled to the induced folding of the K-turn, allowing some sub-optimal variants to adopt the kinked geometry. However, in all such cases the affinity is lowered, and the results underline the importance of both G·A pairs to the stability of the K-turn.     

Protein-protein interactions as a tool for site-specific labeling of proteins

Probing structures and dynamics within biomolecules using ensemble and single-molecule fluorescence resonance energy transfer requires the conjugation of fluorophores to proteins in a site-specific and thermodynamically nonperturbative fashion. Using single-molecule fluorescence-aided molecular sorting and the chymotrypsin inhibitor 2-subtilisin BPN' complex as an example, we demonstrate that protein-protein interactions can be exploited to afford site-specific labeling of a recombinant double-cysteine variant of CI2 without the need for extensive and time-consuming chromatography. The use of protein-protein interactions for site-specific labeling of proteins is compatible with and complementary to existing chemistries for selective labeling of N-terminal cysteines, and could be extended to label multiple positions within a given polypeptide chain.     

Identification of the nucleophilic factors and the productive complex for the editing reaction by leucyl-tRNA synthetase

To ensure fidelity of translation, several aminoacyl-tRNA synthetases (aaRSs) possess editing capability to hydrolyse mis-aminoacylated tRNAs. In this report, based on our previously-modelled structure of leucyl-tRNA synthetase (LeuRS) complexed with valyl-tRNA^Leu, further structural modelling has been performed along with molecular dynamics simulations. This enabled the identification of the nucleophile, which is different from that suggested by the crystal structure of the LeuRS • Nva2AA complex. Our results revealed that the 3′ hydroxyl group of A76 acts as a “gate” to regulate the accessibility of the nucleophile; thus, the opening of the gate leads to the productive complex for the reaction.