Conformational regulation of the EGFR kinase core by the juxtamembrane and C-terminal tail: a molecular dynamics study

The catalytic domain of epidermal growth factor receptor (EGFR) is activated by dimerization, which requires allosteric coupling between distal dimerization and catalytic sites. Although crystal structures of EGFR kinases, solved in various conformational states, have provided important insights into EGFR activation by dimerization, the atomic details of how dimerization signals are dynamically coupled to catalytic regions of the kinase core are not fully understood. In this study, we have performed unrestrained and targeted molecular dynamics simulations on the active and inactive states of EGFR, followed by principal component analysis on the simulated trajectories, to identify correlated motions in the EGFR kinase domain upon dimerization. Our analysis reveals that the conformational changes associated with the catalytic functions of the kinase core are highly correlated with motions in the juxtamembrane (JM) and C-terminal tail, two flexible structural elements that play an active role in EGFR kinase activation and dimerization. In particular, the opening and closing of the ATP binding lobe relative to the substrate binding lobe is highly correlated with motions in the JM and C-terminal tail, suggesting that ATP and substrate binding can be coordinated with dimerization through conformational changes in the JM and C-terminal tail. Our study pinpoints key residues involved in this conformational coupling, and provides new insights into the role of the JM and C-terminal tail segments in EGFR kinase functions.     

A comparison of the dynamics of pantothenate synthetase from M. tuberculosis and E. coli: computational studies

Pantothenate synthetase (PS) catalyzes the final step of the pantothenate pathway, in which pantothenate is formed from pantoate and β-alanine in an ATP-dependent reaction. Mycobacterium tuberculosis PS (MTB PS) is functionally a dimer and a potential target for novel antitubercular drugs. Molecular dynamics simulations show that the functional dynamics of the enzyme are dominated by motions of a flexible gate loop in the N-terminal domain and of the C-terminal domain. The gate loop motions dominate in MTB PS while the C-terminal domain motion dominates in Escherichia coli PS. Simulations also show that the correlated motions of the domains are severely compromised in the monomeric forms. Mutations that reduce the mobility of the gate loop in MTB PS and increased it in E. coli PS were designed and validated through simulations.     

A machine learning approach for the prediction of protein surface loop flexibility

Proteins often undergo conformational changes when binding to each other. A major fraction of backbone conformational changes involves motion on the protein surface, particularly in loops. Accounting for the motion of protein surface loops represents a challenge for protein-protein docking algorithms. A first step in addressing this challenge is to distinguish protein surface loops that are likely to undergo backbone conformational changes upon protein-protein binding (mobile loops) from those that are not (stationary loops). In this study, we developed a machine learning strategy based on support vector machines (SVMs). Our SVM uses three features of loop residues in the unbound protein structures-Ramachandran angles, crystallographic B-factors, and relative accessible surface area-to distinguish mobile loops from stationary ones. This method yields an average prediction accuracy of 75.3% compared with a random prediction accuracy of 50%, and an average of 0.79 area under the receiver operating characteristic (ROC) curve using cross-validation. Testing the method on an independent dataset, we obtained a prediction accuracy of 70.5%. Finally, we applied the method to 11 complexes that involve members from the Ras superfamily and achieved prediction accuracy of 92.8% for the Ras superfamily proteins and 74.4% for their binding partners.     

Evolved streptavidin mutants reveal key role of loop residue in high-affinity binding

We have performed a detailed analysis of streptavidin variants with altered specificity towards desthiobiotin. In addition to changes in key residues which widen the ligand binding pocket and accommodate the more structurally flexible desthiobiotin, the data revealed the role of a key, non-active site mutation at the base of the flexible loop (S52G) which slows dissociation of this ligand by approximately sevenfold. Our data suggest that this mutation results in the loss of a stabilizing contact which keeps this loop open and accessible in the absence of ligand. When this mutation was introduced into the wild-type protein, destabilization of the opened loop conferred a ～10-fold decrease in both the on-rate and off-rate for the ligand biotin-4-fluoroscein. A similar effect was observed when this mutation was added to a monomeric form of this protein. Our results provide key insight into the role of the streptavidin flexible loop in ligand binding and maintaining high affinity interactions.     

Strand 6B deformation and residues exposure towards N-terminal end of helix B during proteinase inhibition by Serpins

Serine Protease inhibitors (Serpins) like antithrombin, antitrypsin, neuroserpin, antichymotrypsin, protein C-inhibitor and plasminogen activator inhibitor is involved in important biological functions like blood coagulation, fibrinolysis, inflammation, cell migration and complement activation. Serpins native state is metastable, which undergoes transformation to a more stable state during the process of protease inhibition. Serpins are prone to conformation defects, however little is known about the factors and mechanisms which promote its conformational change and misfolding. Helix B region in serpins is with several point mutations which result in pathological conditions due to polymerization. Helix B analysis for residue burial and cavity was undertaken to understand its role in serpin structure function. A structural overlap and an accessible surface area analysis showed the deformation of strand 6B and exposure of helix B at N-terminal end in cleaved conformation but not in the native and latent conformation of various inhibitory serpins. A cleaved polymer like conformation of antitrypsin also showed deformation of s6B and helix B exposure. Cavity analysis showed that helix B residues were part of the largest cavity in most of the serpins in the native state which increase in size during the transformation to cleaved and latent states. These data for the first time show the importance of strand 6B deformation and exposure of helix B in smooth insertion of the reactive center loop during serpin inhibition and indicate that helix B exposure due to variants may increase its polymer propensity. ABBREVIATIONS: serpin -serine protease inhibitors RCL -reactive center loop ASA -accessible surface area.     

Substrate specificity and mutational analysis of Kluyveromyces lactis gamma-toxin, a eukaryal tRNA anticodon nuclease

tRNA anticodon damage inflicted by the Kluyveromyces lactis γ-toxin underlies an RNA-based innate immune system that distinguishes self from nonself species. γ-toxin arrests the growth of Saccharomyces cerevisiae by incising a single phosphodiester 3' of the wobble base of tRNA(Glu(UUC)) to generate a break with 2',3'-cyclic phosphate and 5'-OH ends. Recombinant γ-toxin cleaves oligonucleotide substrates in vitro that mimic the anticodon stem-loop of tRNA(Glu). A single 2'-deoxy sugar substitution at the wobble nucleoside abolishes anticodon nuclease activity. To gain further insights to γ-toxin's substrate specificity, we tested deoxynucleoside effects at positions other than the site of transesterification. The results attest to a stringent requirement for a ribonucleoside at the uridine 5' of the wobble base. In contrast, every other nonwobble ribonucleoside in the anticodon loop can be replaced by a deoxy without significantly affecting γ-toxin's cleavage activity. Whereas either the 5' half or the 3' half of the anticodon stem can be replaced en bloc with DNA without a major effect, simultaneously replacing both strands with DNA interfered strongly, signifying that γ-toxin requires an A-form helical conformation of the anticodon stem. We purified γ-toxin mutants identified previously as nontoxic in vivo and gauged their anticodon nuclease activities in vitro. The results highlight Glu9 and Arg151 as candidate catalytic residues, along with His209 implicated previously. By analogy to other endoribonucleases, we speculate that γ-toxin drives transesterification by general acid-base catalysis (via His209 and Glu9) and transition-state stabilization (via Arg151).     

The length of the reactive center loop modulates the latency transition of plasminogen activator inhibitor-1

Plasminogen activator inhibitor-1 (PAI-1) belongs to the serine protease inhibitor (serpin) protein family, which has a common tertiary structure consisting of three beta-sheets and several alpha-helices. Despite the similarity of its structure with those of other serpins, PAI-1 is unique in its conformational lability, which allows the conversion of the metastable active form to a more stable latent conformation under physiological conditions. For the conformational conversion to occur, the reactive center loop (RCL) of PAI-1 must be mobilized and inserted into the major beta-sheet, A sheet. In an effort to understand how the structural conversion is regulated in this conformationally labile serpin, we modulated the length of the RCL of PAI-1. We show that releasing the constraint on the RCL by extension of the loop facilitates a conformational transition of PAI-1 to a stable state. Biochemical data strongly suggest that the stabilization of the transformed conformation is owing to the insertion of the RCL into A beta-sheet, as in the known latent form. In contrast, reducing the loop length drastically retards the conformational change. The results clearly show that the constraint on the RCL is a factor that regulates the conformational transition of PAI-1.