Functional asymmetry in the lysyl-tRNA synthetase explored by molecular dynamics,free energy calculations and experiment

Background  

Charging of transfer-RNA with cognate amino acid is accomplished by the aminoacyl-tRNA synthetases, and proceeds through an aminoacyl adenylate intermediate. The lysyl-tRNA synthetase has evolved an active site that specifically binds lysine and ATP. Previous molecular dynamics simulations of the heat-inducible Escherichia coli lysyl-tRNA synthetase, LysU, have revealed differences in the binding of ATP and aspects of asymmetry between the nominally equivalent active sites of this dimeric enzyme. The possibility that this asymmetry results in different binding affinities for the ligands is addressed here by a parallel computational and biochemical study.     

Asymmetry in membrane responses to electric shocks: insights from bidomain simulations

Models of myocardial membrane dynamics have not been able to reproduce the experimentally observed negative bias in the asymmetry of transmembrane potential changes (DeltaVm) induced by strong electric shocks delivered during the action potential plateau. The goal of this study is to determine what membrane model modifications can bridge this gap between simulation and experiment. We conducted simulations of shocks in bidomain fibers and sheets with membrane dynamics represented by the LRd'2000 model. We found that in the fiber, the negative bias in DeltaVm asymmetry could not be reproduced by addition of electroporation only, but by further addition of hypothetical outward current, Ia, activated upon strong shock-induced depolarization. Furthermore, the experimentally observed rectangularly shaped positive DeltaVm, negative-to-positive DeltaVm ratio (asymmetry ratio) = approximately 2, electroporation occurring at the anode only, and the increase in positive DeltaVm caused by L-type Ca2+-channel blockade were reproduced in the strand only if Ia was assumed to be a part of K+ flow through the L-type Ca2+-channel. In the sheet, Ia not only contributed to the negative bias in DeltaVm asymmetry at sites polarized by physical and virtual electrodes, but also restricted positive DeltaVm. Inclusion of Ia and electroporation is thus the bridge between experiment and simulation.     

The duality of LysU, a catalyst for both Ap4A and Ap3A formation

Heat shock inducible lysyl-tRNA synthetase of Escherichia coli (LysU) is known to be a highly efficient diadenosine 5',5'-P1,P4-tetraphosphate (Ap4A) synthase. However, we use an ion-exchange HPLC technique to demonstrate that active LysU mixtures actually have a dual catalytic activity, initially producing Ap4A from ATP, before converting that tetraphosphate to a triphosphate. LysU appears to be an effective diadenosine 5',5'-P1,P3-triphosphate (Ap3A) synthase. Mechanistic investigations reveal that Ap3A formation requires: (a) that the second step of Ap4A formation is slightly reversible, thereby leading to a modest reappearance of adenylate intermediate; and (b) that phosphate is present to trap the intermediate (either as inorganic phosphate, as added ADP, or as ADP generated in situ from inorganic phosphate). Ap3A forms readily from Ap4A in the presence of such phosphate-based adenylate traps (via a 'reverse-trap' mechanism). LysU is also clearly demonstrated to exist in a phosphorylated state that is more physically robust as a catalyst of Ap4A formation than the nonphosphorylated state. However, phosphorylated LysU shows only marginally improved catalytic efficiency. We note that Ap3A effects have barely been studied in prokaryotic organisms. By contrast, there is a body of literature that describes Ap3A and Ap4A having substantially different functions in eukaryotic cells. Our data suggest that Ap3A and Ap4A biosynthesis could be linked together through a single prokaryotic dual 'synthase' enzyme. Therefore, in our view there is a need for new research into the effects and impact of Ap3A alone and the intracellular [Ap3A]/[Ap4A] ratio on prokaryotic organisms.     

The lrp gene product regulates expression of lysU in Escherichia coli K-12.

In Escherichia coli K-12, expression of the lysU gene is regulated by the lrp gene product, as indicated by an increase in the level of lysyl-tRNA synthetase activity and LysU protein in an lrp mutant. Comparison of the patterns of protein expression visualized by two-dimensional gel electrophoresis indicated that LysU is present at higher levels in an lrp strain than in its isogenic lrp+ parent. The purified lrp gene product was shown to bind to sites upstream of the lysU gene and to protect several sites against DNase I digestion. A region extending over 100 nucleotides, between 60 and 160 nucleotides upstream from the start of the lysU coding sequence, showed altered sensitivity to DNase I digestion in the presence of the Lrp protein. The extent of protected DNA suggests a complex interaction of Lrp protein and upstream lysU DNA.     

Simulated annealing coupled replica exchange molecular dynamics—An efficient conformational sampling method

Molecular dynamics simulated annealing (SA-MD) simulations are frequently used for refinement and optimization of peptide and protein structures. Depending on the simulation conditions and simulation length SA-MD simulations can be trapped in locally stable conformations far from the global optimum. As an alternative replica exchange molecular dynamics (RexMD) simulations can be used which allow exchanges between high and low simulation temperatures at all stages of the simulation. A significant drawback of RexMD simulations is, however, the rapid increase of the replica number with increasing system size to cover a desired temperature range. A combined SA-MD and RexMD approach termed SA-RexMD is suggested that employs a small number of replicas (4) and starts initially with a set of high simulation temperatures followed by gradual cooling of the set of temperatures until a target temperature has been reached. The protocol has been applied for the folding of several peptide systems and for the refinement of protein model structures. In all the cases, the SA-RexMD method turned out to be significantly more efficient in reaching low energy structures and also structures close to experiment compared to continuous MD simulations at the target temperature and to SA-MD simulations at the same computational demand. The approach is well suited for applications in structure refinement and for systematic force field improvement.     

Homology modeling and molecular dynamics simulations of the glycine receptor ligand binding domain

We present a homology based model of the ligand binding domain (LBD) of the homopentameric alpha1 glycine receptor (GlyR). The model is based on multiple sequence alignment with other members of the nicotinicoid ligand gated ion channel superfamily and two homologous acetylcholine binding proteins (AChBP) from the freshwater (Lymnaea stagnalis) and saltwater (Aplysia californica) snails with known high resolution structure. Using two template proteins with known structure to model three dimensional structure of a target protein is especially advantageous for sequences with low homology as in the case presented in this paper. The final model was cross-validated by critical evaluation of experimental and published mutagenesis, functional and other biochemical studies. In addition, a complex structure with strychnine antagonist in the putative binding site is proposed based on docking simulation using Autodock program. Molecular dynamics (MD) simulations with simulated annealing protocol are reported on the proposed LBD of GlyR, which is stable in 5 ns simulation in water, as well as for a deformed LBD structure modeled on the corresponding domain determined in low-resolution cryomicroscopy structure of the alpha subunit of the full-length acetylcholine receptor (AChR). Our simulations demonstrate that the beta-sandwich central core of the protein monomer is fairly rigid in the simulations and resistant to deformations in water.     

Partial characterization of a lysU mutant of Escherichia coli K-12.

The Escherichia coli K-12 strain GNB10181 shows no inducible lysyl-tRNA synthetase (LysRS) activity. Two-dimensional gel electrophoretic analysis of the polypeptides synthesized by this strain indicates that the normal lysU gene product, LysU, is absent. When both GNB10181 and its parent, MC4100, were grown at elevated temperatures (42 to 45 degrees C) no significant difference between their growth rates was observed. The lysU mutation was transferred to other E. coli K-12 backgrounds by using P1 transduction. The lysU transductants behaved comparably to their lysU+ parents at different growth temperatures. Therefore, the LysU proteins does not appear to be essential for growth at high temperatures, at least under the conditions examined here. In addition, lysU transductants were found to be defective for inducible lysine decarboxylase, (LDC), inducible arginine decarboxylase (ADI), and melibiose utilization (Mel), which are all missing in GNB10181. Complementation of the above missing functions was achieved by using the Clarke-Carbon plasmids pLC4-5 (LysU LDC) and pLC17-38 (LysU Mel ADI). From these experiments, it appears that GNB10181 has suffered a chromosomal deletion between 93.4 and 93.7 min, which includes the lysU gene. By using plasmid pLC17-38, the position of ADI on two-dimensional gels was identified. Finally, lysS delta lysU double mutants were constructed which can potentially be used as positive selection agents for the isolation of LysRS genes from other sources.     

Three-dimensional model of a substrate-bound SARS chymotrypsin-like cysteine proteinase predicted by multiple molecular dynamics simulations: catalytic efficiency regulated by substrate binding

Severe acute respiratory syndrome (SARS) is a contagious and deadly disease caused by a new coronavirus. The protein sequence of the chymotrypsin-like cysteine proteinase (CCP) responsible for SARS viral replication has been identified as a target for developing anti-SARS drugs. Here, I report the ATVRLQ(p1)A(p1')-bound CCP 3D model predicted by 420 different molecular dynamics simulations (2.0 ns for each simulation with a 1.0-fs time step). This theoretical model was released at the Protein Data Bank (PDB; code: 1P76) before the release of the first X-ray structure of CCP (PDB code: 1Q2W). In contrast to the catalytic dyad observed in X-ray structures of CCP and other coronavirus cysteine proteinases, a catalytic triad comprising Asp187, His41, and Cys145 is found in the theoretical model of the substrate-bound CCP. The simulations of the CCP complex suggest that substrate binding leads to the displacement of a water molecule entrapped by Asp187 and His41, thus converting the dyad to a more efficient catalytic triad. The CCP complex structure has an expanded active-site pocket that is useful for anti-SARS drug design. In addition, this work demonstrates that multiple molecular dynamics simulations are effective in correcting errors that result from low-sequence-identity homology modeling.     

Molecular dynamics simulations of peptides and proteins with amplified collective motions

We present a novel method that uses the collective modes obtained with a coarse-grained model/anisotropic network model to guide the atomic-level simulations. Based on this model, local collective modes can be calculated according to a single configuration in the conformational space of the protein. In the molecular dynamics simulations, the motions along the slowest few modes are coupled to a higher temperature by the weak coupling method to amplify the collective motions. This amplified-collective-motion (ACM) method is applied to two test systems. One is an S-peptide analog. We realized the refolding of the denatured peptide in eight simulations out of 10 using the method. The other system is bacteriophage T4 lysozyme. Much more extensive domain motions between the N-terminal and C-terminal domain of T4 lysozyme are observed in the ACM simulation compared to a conventional simulation. The ACM method allows for extensive sampling in conformational space while still restricting the sampled configurations within low energy areas. The method can be applied in both explicit and implicit solvent simulations, and may be further applied to important biological problems, such as long timescale functional motions, protein folding/unfolding, and structure prediction.     

Detecting and adjusting for small-study effects in meta-analysis

Publication bias and related types of small-study effects threaten the validity of systematic reviews. The existence of small-study effects has been demonstrated in empirical studies. Small-study effects are graphically diagnosed by inspection of the funnel plot. Though observed funnel plot asymmetry cannot be easily linked to a specific reason, tests based on funnel plot asymmetry have been proposed. Beyond a vast range of funnel plot tests, there exist several methods for adjusting treatment effect estimates for these biases. In this article, we consider the trim-and-fill method, the Copas selection model, and more recent regression-based approaches. The methods are exemplified using a meta-analysis from the literature and compared in a simulation study, based on binary response data. They are also applied to a large set of meta-analyses. Some fundamental differences between the approaches are discussed. An assumption common to the trim-and-fill method and the Copas selection model is that the small-study effect is caused by selection. The trim-and-fill method corresponds to an unknown implicit model generated by the symmetry assumption, whereas the Copas selection model is a parametric statistical model. However, it requires a sensitivity analysis. Regression-based approaches are easier to implement and not based on a specific selection model. Both simulations and applications suggest that in the presence of strong selection both the trim-and-fill method and the Copas selection model may not fully eliminate bias, while regression-based approaches seem to be a promising alternative.     

Validation of an Allosteric Binding Site of Src Kinase Identified by Unbiased Ligand Binding Simulations

Allostery plays a primary role in regulating protein activity, making it an important mechanism in human disease and drug discovery. Identifying allosteric regulatory sites to explore their biological significance and therapeutic potential is invaluable to drug discovery; however, identification remains a challenge. Allosteric sites are often “cryptic” without clear geometric or chemical features. Since allosteric regulatory sites are often less conserved in protein kinases than the orthosteric ATP binding site, allosteric ligands are commonly more specific than ATP competitive inhibitors. We present a generalizable computational protocol to predict allosteric ligand binding sites based on unbiased ligand binding simulation trajectories. We demonstrate the feasibility of this protocol by revisiting our previously published ligand binding simulations using the first identified viral proto-oncogene, Src kinase, as a model system. The binding paths for kinase inhibitor PP1 uncovered three metastable intermediate states before binding the high-affinity ATP-binding pocket, revealing two previously known allosteric sites and one novel site. Herein, we validate the novel site using a combination of virtual screening and experimental assays to identify a V-type allosteric small-molecule inhibitor that targets this novel site with specificity for Src over closely related kinases. This study provides a proof-of-concept for employing unbiased ligand binding simulations to identify cryptic allosteric binding sites and is widely applicable to other protein–ligand systems.     

Molecular dynamics simulation and experimental validation of the effect of pH on protein desorption in hydrophobic charge induction chromatography

The ligands in hydrophobic charge induction chromatography (HCIC) are hydrophobic and ionisable. Thus, the pH is crucial for the separation performance in HCIC, especially for elution. However, it is difficult to obtain the microscopic information in HCIC through experimental means. In this work, molecular dynamics simulations are performed to examine the effect of pH on elution and protein conformational transition in HCIC, using a 46-bead β-barrel coarse-grained model protein and an HCIC adsorbent pore model constructed in an earlier work. Corresponding experiments are carried out for the validation of simulation results, using lysozyme and MEP Hypercel. Both the activities and fluorescence of lysozyme are examined to evaluate the conformational transition. The simulations indicate that the elution efficiency of protein increases with decreasing pH value in a non-linear manner. This is qualitatively consistent with the experimental results. MD simulations indicate that protein unfolding occurs in elution at all pH values. However, the experimental data show that the activity and conformation of lysozyme is independent of pH of the elution buffer. The microscopic information from simulation shows that protein unfolding is mainly observed on the adsorbent surface, but it cannot be detected in the experiments that only probe the proteins in the bulk liquid.     

Fluid-gel coexistence in lipid membranes under differential stress

A widely conserved property of many biological lipid bilayers is their asymmetry. In addition to having distinct compositions on its two sides, a membrane can also exhibit different tensions in its two leaflets, a state known as differential stress. Here, we examine how this stress can influence the phase behavior of the constituent lipid monolayers of a single-component membrane. For temperatures sufficiently close to, but still above, the main transition, molecular dynamics simulations show the emergence of finite gel domains within the compressed leaflet. We describe the thermodynamics of this phenomenon by adding two empirical single-leaflet free energies for the fluid-gel transition, each evaluated at its respective asymmetry-dependent lipid density. Finite size effects arising in simulation are included in the theory through a geometry-dependent interfacial term. Our model reproduces the phase coexistence observed in simulation. It could therefore be used to connect the “hidden variable” of differential stress to experimentally observable properties of the main phase transition. These ideas could be generalized to any first-order bilayer phase transition in the presence of asymmetry, including liquid-ordered/liquid-disordered phase separation.     

Phosphatidylserine Asymmetry Promotes the Membrane Insertion of a Transmembrane Helix

The plasma membrane (PM) contains an asymmetric distribution of lipids between the inner and outer bilayer leaflets. A lipid of special interest in eukaryotic membranes is the negatively charged phosphatidylserine (PS). In healthy cells, PS is actively sequestered to the inner leaflet of the PM, but PS redistributes to the outer leaflet when the cell is damaged or at the onset of apoptosis. However, the influence of PS asymmetry on membrane protein structure and folding are poorly understood. The pH low insertion peptide (pHLIP) adsorbs to the membrane surface at a neutral pH, but it inserts into the membrane at an acidic pH. We have previously observed that in symmetric vesicles, PS affects the membrane insertion of pHLIP by lowering the pH midpoint of insertion. Here, we studied the effect of PS asymmetry on the membrane interaction of pHLIP. We developed a modified protocol to create asymmetric vesicles containing PS and employed Annexin V labeled with an Alexa Fluor 568 fluorophore as a new probe to quantify PS asymmetry. We observed that the membrane insertion of pHLIP was promoted by the asymmetric distribution of negatively charged PS, which causes a surface charge difference between bilayer leaflets. Our results indicate that lipid asymmetry can modulate the formation of an α-helix on the membrane. A corollary is that model studies using symmetric bilayers to mimic the PM may fail to capture important aspects of protein-membrane interactions.     

FF12MC: A revised AMBER forcefield and new protein simulation protocol

Specialized to simulate proteins in molecular dynamics (MD) simulations with explicit solvation, FF12MC is a combination of a new protein simulation protocol employing uniformly reduced atomic masses by tenfold and a revised AMBER forcefield FF99 with (i) shortened C? H bonds, (ii) removal of torsions involving a nonperipheral sp³ atom, and (iii) reduced 1–4 interaction scaling factors of torsions ? and ψ. This article reports that in multiple, distinct, independent, unrestricted, unbiased, isobaric–isothermal, and classical MD simulations FF12MC can (i) simulate the experimentally observed flipping between left‐ and right‐handed configurations for C14–C38 of BPTI in solution, (ii) autonomously fold chignolin, CLN025, and Trp‐cage with folding times that agree with the experimental values, (iii) simulate subsequent unfolding and refolding of these miniproteins, and (iv) achieve a robust Z score of 1.33 for refining protein models TMR01, TMR04, and TMR07. By comparison, the latest general‐purpose AMBER forcefield FF14SB locks the C14–C38 bond to the right‐handed configuration in solution under the same protein simulation conditions. Statistical survival analysis shows that FF12MC folds chignolin and CLN025 in isobaric–isothermal MD simulations 2–4 times faster than FF14SB under the same protein simulation conditions. These results suggest that FF12MC may be used for protein simulations to study kinetics and thermodynamics of miniprotein folding as well as protein structure and dynamics. Proteins 2016; 84:1490–1516. © 2016 The Authors Proteins: Structure, Function, and Bioinformatics Published by Wiley Periodicals, Inc.     

Effect of Ceramide Tail Length on the Structure of Model Stratum Corneum Lipid Bilayers

Lipid bilayers composed of non-hydroxy sphingosine ceramide (CER NS), cholesterol (CHOL), and free fatty acids (FFAs), which are components of the human skin barrier, are studied via molecular dynamics simulations. Since mixtures of these lipids exist in dense gel phases with little molecular mobility at physiological conditions, care must be taken to ensure that the simulations become decorrelated from the initial conditions. Thus, we propose and validate an equilibration protocol based on simulated tempering, in which the simulation takes a random walk through temperature space, allowing the system to break out of metastable configurations and hence become decorrelated from its initial configuration. After validating the equilibration protocol, which we refer to as random-walk molecular dynamics, the effects of the lipid composition and ceramide tail length on bilayer properties are studied. Systems containing pure CER NS, CER NS + CHOL, and CER NS + CHOL + FFA, with the CER NS fatty acid tail length varied within each CER NS-CHOL-FFA composition, are simulated. The bilayer thickness is found to depend on the structure of the center of the bilayer, which arises as a result of the tail-length asymmetry between the lipids studied. The hydrogen bonding between the lipid headgroups and with water is found to change with the overall lipid composition, but is mostly independent of the CER fatty acid tail length. Subtle differences in the lateral packing of the lipid tails are also found as a function of CER tail length. Overall, these results provide insight into the experimentally observed trend of altered barrier properties in skin systems where there are more CERs with shorter tails present.     

In silico chaperonin-like cycle helps folding of proteins for structure prediction

Currently, one of the most serious problems in protein-folding simulations for de novo structure prediction is conformational sampling of medium-to-large proteins. In vivo, folding of these proteins is mediated by molecular chaperones. Inspired by the functions of chaperonins, we designed a simple chaperonin-like simulation protocol within the framework of the standard fragment assembly method: in our protocol, the strength of the hydrophobic interaction is periodically modulated to help the protein escape from misfolded structures. We tested this protocol for 38 proteins and found that, using a certain defined criterion of success, our method could successfully predict the native structures of 14 targets, whereas only those of 10 targets were successfully predicted using the standard protocol. In particular, for non-α-helical proteins, our method yielded significantly better predictions than the standard approach. This chaperonin-inspired protocol that enhanced de novo structure prediction using folding simulations may, in turn, provide new insights into the working principles underlying the chaperonin system.     

Assembly of a tetrameric alpha-helical bundle: computer simulations on an intermediate-resolution protein model.

Discontinuous molecular dynamics (DMD) simulation on an intermediate-resolution protein model is used to study the folding of an isolated, small model peptide to an amphipathic alpha-helix and the assembly of four of these model peptides into a four-helix bundle. A total of 129 simulations were performed on the isolated peptide, and 50 simulations were performed on the four-peptide system. Simulations efficiently sample conformational space allowing complete folding trajectories from random initial configurations to be observed within 15 min for the one-peptide system and within 15 h for the four-peptide system on a 500-MHz workstation. The native structures of both the alpha-helix and the four-helix bundle are consistent with experimental characterization studies and with results from previous simulations on these model peptides. In both the one- and four-peptide systems, the native state is achieved during simulations within an optimal temperature range, a phenomenon also observed experimentally. The ease with which our simulations yield reasonable estimates of folded structures demonstrates the power of the intermediate-resolution model developed for this work and the DMD algorithm and suggests that simulations of very long times and of multiprotein systems may be possible with this model.     

Control of membrane fusion mechanism by lipid composition: predictions from ensemble molecular dynamics

Membrane fusion is critical to biological processes such as viral infection, endocrine hormone secretion, and neurotransmission, yet the precise mechanistic details of the fusion process remain unknown. Current experimental and computational model systems approximate the complex physiological membrane environment for fusion using one or a few protein and lipid species. Here, we report results of a computational model system for fusion in which the ratio of lipid components was systematically varied, using thousands of simulations of up to a microsecond in length to predict the effects of lipid composition on both fusion kinetics and mechanism. In our simulations, increased phosphatidylcholine content in vesicles causes increased activation energies for formation of the initial stalk-like intermediate for fusion and of hemifusion intermediates, in accordance with previous continuum-mechanics theoretical treatments. We also use our large simulation dataset to quantitatively compare the mechanism by which vesicles fuse at different lipid compositions, showing a significant difference in fusion kinetics and mechanism at different compositions simulated. As physiological membranes have different compositions in the inner and outer leaflets, we examine the effect of such asymmetry, as well as the effect of membrane curvature on fusion. These predicted effects of lipid composition on fusion mechanism both underscore the way in which experimental model system construction may affect the observed mechanism of fusion and illustrate a potential mechanism for cellular regulation of the fusion process by altering membrane composition.     

Conformational sampling by NMR solution structures calculated with the program DIANA evaluated by comparison with long-time molecular dynamics calculations in explicit water

The NMR solution structure of bovine pancreatic trypsin inhibitor (BPTI) obtained by distance geometry calculations with the program DIANA is compared with groups of conformers generated by molecular dynamics (MD) simulations in explicit water at ambient temperature and pressure. The MD simulations started from a single conformer and were free or restrained either by the experimental NOE distance restraints or by time-averaged restraints; the groups of conformers were collected either in 10 ps intervals during 200 ps periods of simulation, or in 50 ps intervals during a 1 ns period of simulation. Overall, these comparisons show that the standard protein structure determination protocol with the program DIANA provides a picture of the protein structure that is in agreement with MD simulations using “realistic” potential functions over a nanosecond timescale. For well-constrained molecular regions there is a trend in the free MD simulation of duration 1 ns that the sampling of the conformation space is slightly increased relative to the DIANA calculations. In contrast, for surface-exposed side-chains that are less extensively constrained by the NMR data, the DIANA conformers tend to sample larger regions of conformational space than conformers selected from any of the MD trajectories. Additional insights into the behavior of surface side-chains come from comparison of the MD runs of 200 ps or 1 ns duration. In this time range the sampling of conformation space by the protein surface depends strongly on the length of the simulation, which indicates that significant side-chain transitions occur on the nanosecond timescale and that much longer simulations will be needed to obtain statistically significant data on side-chain dynamics.