Disorder-to-order transition of an intrinsically disordered region of sortase revealed by multiscale enhanced sampling

Molecular functions of intrinsically disordered proteins (IDPs) or intrinsically disordered regions (IDRs), such as molecular recognition and cellular signaling, are ascribed to dynamic changes in the conformational space in response to binding of target molecules. Sortase, a transpeptitase in Gram-positive bacteria, has an IDR in a loop which undergoes a disordered-to-ordered transition (called "disordered loop"), accompanying a tilt of another loop ("dynamic loop"), upon binding of a signal peptide and a calcium ion. In this study, all-atom conformational ensembles of sortase were calculated for the four different binding states (with/without the peptide and with/without a calcium ion) by the multiscale enhanced sampling (MSES) simulation to examine how the binding of the peptide and/or calcium influences the conformational ensemble. The MSES is a multiscale and multicopy simulation method that allows an enhanced sampling of the all-atom model of large proteins including explicit solvent. A 100 ns MSES simulation of the ligand-free sortase using 20 replicas (in total 2 μs) demonstrated large flexibility in both the disordered and dynamic loops; however, their distributions were not random but had a clear preference which populates the N-terminal part of the disordered loop near the bound form. The MSES simulations of the three binding states clarified the allosteric mechanism of sortase: the N- and C-terminal parts of the disordered loop undergo a disorder-to-order transition independently of each other upon binding of the peptide and a calcium ion, respectively; however, upon binding of both ligands, the two parts work cooperatively to stabilize the bound peptide.     

Incorporating dynamics in E. coli dihydrofolate reductase enhances structure-based drug discovery

Escherichia coli dihydrofolate reductase (DHFR) is a long-standing target for enzyme studies. The influence of protein motion on its catalytic cycle is significant, and the conformation of the M20 loop is of particular interest. We present receptor-based pharmacophore models-an equivalent of solvent-mapping of binding hotspots-based on ensembles of protein conformations from molecular dynamics simulations of DHFR.NADPH in both the closed and open conformation of the M20 loop. The optimal models identify DHFR inhibitors over druglike non-inhibitors; furthermore, high-affinity inhibitors of E. coli DHFR are preferentially identified over general DHFR inhibitors. As expected, models resulting from simulations with DHFR in the productive conformation with a closed M20 loop have better performance than those from the open-loop simulations. Model performance improves with increased dynamic sampling, indicating that including a greater degree of protein flexibility can enhance the quest for potent inhibitors. This was compared to the limited conformational sampling seen in crystal structures, which were suboptimal for this application.     

Relative stability of the open and closed conformations of the active site loop of streptavidin

The eight-residue surface loop, 45-52 (Ser, Ala, Val, Gly, Asn, Ala, Glu, Ser), of the homotetrameric protein streptavidin has a "closed" conformation in the streptavidin-biotin complex, where the corresponding binding affinity is one of the strongest found in nature (ΔG ～ -18 kcal∕mol). However, in most of the crystal structures of apo (unbound) streptavidin, the loop conformation is "open" and typically exhibits partial disorder and high B-factors. Thus, it is plausible to assume that the loop structure is changed from open to closed upon binding of biotin, and the corresponding difference in free energy, ΔF = F(open) - F(closed) in the unbound protein, should therefore be considered in the total absolute free energy of binding. ΔF (which has generally been neglected) is calculated here using our "hypothetical scanning molecular-dynamics" (HSMD) method. We use a protein model in which only the atoms closest to the loop are considered (the "template") and they are fixed in the x-ray coordinates of the free protein; the x-ray conformation of the closed loop is attached to the same (unbound) template and both systems are capped with the same sphere of TIP3P water. Using the force field of the assisted model building with energy refinement (AMBER), we carry out two separate MD simulations (at temperature T = 300 K), starting from the open and closed conformations, where only the atoms of the loop and water are allowed to move (the template-water and template-loop interactions are considered). The absolute F(open) and F(closed) (of loop + water) are calculated from these trajectories, where the loop and water contributions are obtained by HSMD and a thermodynamic integration (TI) process, respectively. The combined HSMD-TI procedure leads to total (loop + water) ΔF = -27.1 ± 2.0 kcal∕mol, where the entropy TΔS constitutes 34% of ΔF, meaning that the effect of S is significant and should not be ignored. Also, ΔS is positive, in accord with the high flexibility of the open loop observed in crystal structures, while the energy ΔE is unexpectedly negative, thus also adding to the stability of the open loop. The loop and the 250 capped water molecules are the largest system studied thus far, which constitutes a test for the efficiency of HSMD-TI; this efficiency and technical issues related to the implementation of the method are also discussed. Finally, the result for ΔF is a prediction that will be considered in the calculation of the absolute free energy of binding of biotin to streptavidin, which constitutes our next project.     

End-to-end vs interior loop formation kinetics in unfolded polypeptide chains

The conformational search for favorable intramolecular interactions during protein folding is limited by intrachain diffusion processes. Recent studies on the dynamics of loop formation in unfolded polypeptide chains have focused on loops involving residues near the chain ends. During protein folding, however, most contacts are formed between residues in the interior of the chain. We compared the kinetics of end-to-end loop formation (type I loops) to the formation of end-to-interior (type II loops) and interior-to-interior loops (type III loops) using triplet-triplet energy transfer from xanthone to naphthylalanine. The results show that formation of type II and type III loops is slower compared to type I loops of the same size and amino acid sequence. The rate constant for type II loop formation decreases with increasing overall chain dimensions up to a limiting value, at which loop formation is about 2.5-fold slower for type II loops compared to type I loops. Comparing type II loops of different loop size and amino acid sequence shows that the ratio of loop dimension over total chain dimension determines the rate constant for loop formation. Formation of type III loops is 1.7-fold slower than formation of type II loops, indicating that local chain motions are strongly coupled to motions of other chain segments which leads to faster dynamics toward the chain ends. Our results show that differences in the kinetics of formation of type I, type II, and type III loops are mainly caused by differences in internal flexibility at the different positions in the polypeptide chain. Interactions of the polypeptide chain with the solvent contribute to the kinetics of loop formation, which are strongly viscosity-dependent. However, the observed differences in the kinetics of formation of type I, type II, and type III loops are not due to the increased number of peptide-solvent interactions in type II and type III loops compared to type I loops as indicated by identical viscosity dependencies for the kinetics of formation of the different types of loops.     

(Thermo)dynamic Role of Receptor Flexibility,Entropy, and Motional Correlation in Protein–Ligand Binding

The binding of 2‐amino‐5‐methylthiazole to the W191G cavity mutant of cytochrome c peroxidase is an ideal test case to investigate the entropic contribution to the binding free energy due to changes in receptor flexibility. The dynamic and thermodynamic role of receptor flexibility are studied by 50 ns‐long explicit‐solvent molecular dynamics simulations of three separate receptor ensembles: W191G binding a K⁺ ion, W191G–2a5mt complex with a closed 190–195 gating loop, and apo with an open loop. We employ a method recently proposed to estimate accurate absolute single‐molecule configurational entropies and their differences for systems undergoing conformational transitions. We find that receptor flexibility plays a generally underestimated role in protein–ligand binding (thermo)dynamics and that changes of receptor motional correlation determine such large entropy contributions.     

A general strategy for creating "inactive-conformation" abl inhibitors

Kinase inhibitors that bind to the ATP cleft can be broadly classified into two groups: those that bind exclusively to the ATP site with the kinase assuming a conformation otherwise conducive to phosphotransfer (type I), and those that exploit a hydrophobic site immediately adjacent to the ATP pocket made accessible by a conformational rearrangement of the activation loop (type II). To date, all type II inhibitors were discovered by using structure-activity-guided optimization strategies. Here, we describe a general pharmacophore model of type II inhibition that enables a rational "hybrid-design" approach whereby a 3-trifluoromethylbenzamide functionality is appended to four distinct type I scaffolds in order to convert them into their corresponding type II counterparts. We demonstrate that the designed compounds function as type II inhibitors by using biochemical and cellular kinase assays and by cocrystallography with Abl.     

Assessing protein loop flexibility by hierarchical Monte Carlo sampling

Loop flexibility is often crucial to protein biological function in solution. We report a new Monte Carlo method for generating conformational ensembles for protein loops and cyclic peptides. The approach incorporates the triaxial loop closure method which addresses the inverse kinematic problem for generating backbone move sets that do not break the loop. Sidechains are sampled together with the backbone in a hierarchical way, making it possible to make large moves that cross energy barriers. As an initial application, we apply the method to the flexible loop in triosephosphate isomerase that caps the active site, and demonstrate that the resulting loop ensembles agree well with key observations from previous structural studies. We also demonstrate, with 3 other test cases, the ability to distinguish relatively flexible and rigid loops within the same protein.     

Insights into the dynamics of HIV-1 protease: a kinetic network model constructed from atomistic simulations

The conformational dynamics in the flaps of HIV-1 protease plays a crucial role in the mechanism of substrate binding. We develop a kinetic network model, constructed from detailed atomistic simulations, to determine the kinetic mechanisms of the conformational transitions in HIV-1 PR. To overcome the time scale limitation of conventional molecular dynamics (MD) simulations, our method combines replica exchange MD with transition path theory (TPT) to study the diversity and temperature dependence of the pathways connecting functionally important states of the protease. At low temperatures the large-scale flap opening is dominated by a small number of paths; at elevated temperatures the transition occurs through many structurally heterogeneous routes. The expanded conformation in the crystal structure 1TW7 is found to closely mimic a key intermediate in the flap-opening pathways at low temperature. We investigated the different transition mechanisms between the semi-open and closed forms. The calculated relaxation times reveal fast semi-open ? closed transitions, and infrequently the flaps fully open. The ligand binding rate predicted from this kinetic model increases by 38-fold from 285 to 309 K, which is in general agreement with experiments. To our knowledge, this is the first application of a network model constructed from atomistic simulations together with TPT to analyze conformational changes between different functional states of a natively folded protein.     

Modeling loop backbone flexibility in receptor‐ligand docking simulations

The relevance of receptor conformational change during ligand binding is well documented for many pharmaceutically relevant receptors, but is still not fully accounted for in in silico docking methods. While there has been significant progress in treatment of receptor side chain flexibility sampling of backbone flexibility remains challenging because the conformational space expands dramatically and the scoring function must balance protein–protein and protein–ligand contributions. Here, we investigate an efficient multistage backbone reconstruction algorithm for large loop regions in the receptor and demonstrate that treatment of backbone receptor flexibility significantly improves binding mode prediction starting from apo structures and in cross docking simulations. For three different kinase receptors in which large flexible loops reconstruct upon ligand binding, we demonstrate that treatment of backbone flexibility results in accurate models of the complexes in simulations starting from the apo structure. At the example of the DFG‐motif in the p38 kinase, we also show how loop reconstruction can be used to model allosteric binding. Our approach thus paves the way to treat the complex process of receptor reconstruction upon ligand binding in docking simulations and may help to design new ligands with high specificity by exploitation of allosteric mechanisms. © 2012 Wiley Periodicals, Inc.     

Simulation of the opening and closing of Hsp70 chaperones by coarse-grained molecular dynamics

Heat-shock proteins 70 (Hsp70s) are key molecular chaperones which assist in the folding and refolding/disaggregation of proteins. Hsp70s, which consist of a nucleotide-binding domain (NBD, consisting of NBD-I and NBD-II subdomains) and a substrate-binding domain [SBD, further split into the β-sheet (SBD-β) and α-helical (SBD-α) subdomains], occur in two major conformations having (a) a closed SBD, in which the SBD and NBD domains do not interact, (b) an open SBD, in which SBD-α interacts with NBD-I and SBD-β interacts with the top parts of NBD-I and NBD-II. In the SBD-closed conformation, SBD is bound to a substrate protein, with release occurring after transition to the open conformation. While the transition from the closed to the open conformation is triggered efficiently by binding of adenosine triphosphate (ATP) to the NBD, it also occurs, although less frequently, in the absence of ATP. The reverse transition occurs after ATP hydrolysis. Here, we report canonical and multiplexed replica exchange simulations of the conformational dynamics of Hsp70s using a coarse-grained molecular dynamics approach with the UNRES force field. The simulations were run in the following three modes: (i) with the two halves of the NBD unrestrained relative to each other, (ii) with the two halves of the NBD restrained in an "open" geometry as in the SBD-closed form of DnaK (2KHO), and (iii) the two halves of NBD restrained in a "closed" geometry as in known experimental structures of ATP-bound NBD forms of Hsp70. Open conformations, in which the SBD interacted strongly with the NBD, formed spontaneously during all simulations; the number of transitions was largest in simulations carried out with the "closed" NBD domain, and smallest in those carried out with the "open" NBD domain; this observation is in agreement with the experimentally-observed influence of ATP-binding on the transition of Hsp70's from the SBD-closed to the SBD-open form. Two kinds of open conformations were observed: one in which SBD-α interacts with NBD-I and SBD-β interacts with the top parts of NBD-I and NBD-II (as observed in the structures of nucleotide exchange factors), and another one in which this interaction pattern is swapped. A third type of motion, in which SBD-α binds to NBD without dissociating from SBD-β was also observed. It was found that the first stage of interdomain communication (approach of SBD-β, to NBD) is coupled with the rotation of the long axes of NBD-I and NBD-II towards each other. To the best of our knowledge, this is the first successful simulation of the full transition of an Hsp70 from the SBD-closed to the SBD-open conformation.     

Binding of low molecular weight inhibitors promotes large conformational changes in the dengue virus NS2B-NS3 protease: fold analysis by pseudocontact shifts

The two-component dengue virus NS2B-NS3 protease (DEN NS2B-NS3pro) is an established drug target, but inhibitor design is hampered by the lack of a crystal structure of the protease in its fully active form. In solution and without inhibitors, the functionally important C-terminal segment of the NS2B cofactor is dissociated from DEN NS3pro ("open state"), necessitating a large structural change to produce the "closed state" thought to underpin activity. We analyzed the fold of DEN NS2B-NS3pro in solution with and without bound inhibitor by nuclear magnetic resonance (NMR) spectroscopy. Multiple paramagnetic lanthanide tags were attached to different sites to generate pseudocontact shifts (PCS). In the face of severe spectral overlap and broadening of many signals by conformational exchange, methods for assignment of (15)N-HSQC cross-peaks included selective mutation, combinatorial isotope labeling, and comparison of experimental PCSs and PCSs back-calculated for a structural model of the closed conformation built by using the structure of the related West Nile virus (WNV) protease as a template. The PCSs show that, in the presence of a positively charged low-molecular weight inhibitor, the enzyme assumes a closed state that is very similar to the closed state previously observed for the WNV protease. Therefore, a model of the protease built on the closed conformation of the WNV protease is a better template for rational drug design than available crystal structures, at least for positively charged inhibitors. To assess the open state, we created a binding site for a Gd(3+) complex and measured paramagnetic relaxation enhancements. The results show that the specific open conformation displayed in the crystal of DEN NS2B-NS3pro is barely populated in solution. The techniques used open an avenue to the fold analysis of proteins that yield poor NMR spectra, as PCSs from multiple sites in combination with model building generate powerful information even from incompletely assigned (15)N-HSQC spectra.     

The large scale conformational change of the human DPP III-substrate prefers the "closed" form

Human dipeptidyl peptidase III (DPP III) is a two domain metallo-peptidase from the M49 family. The wide interdomain cleft and broad substrate specificity suggest that this enzyme could experience significant conformational change. Long (>100 ns) molecular dynamics (MD) simulations of DPP III revealed large range conformational changes of the protein, suggesting the pre-existing equilibrium model for a substrate binding. The binding free energy calculations revealed tighter binding of the preferred synthetic substrate Arg-Arg-2-naphtylamide to the "closed" than to the "open" DPP III conformation. Our assumption that Asp372 plays a crucial role in the large scale interdomain closure was proved by the MD simulations of the Asp372Ala variant. During the same simulation time, the variant remained more "open" than the wild type protein. Apparently, Ala was not as efficient as Asp in establishing the interdomain interactions. According to the MM-PBSA calculations, the electrostatic component of the free energy of solvation turned out to be higher for the "closed" protein than for its less compact form. However, the gain in entropy due to water released from the interdomain cleft nicely balanced this negative effect.     

Coarse-grained modeling of the structural states and transition underlying the powerstroke of dynein motor domain

This study aims to model a minimal dynein motor domain capable of motor function, which consists of the linker domain, six AAA+ modules (AAA1-AAA6), coiled coil stalk, and C-terminus domain. To this end, we have used the newly solved X-ray structures of dynein motor domain to perform a coarse-grained modeling of dynein's post- and pre-powerstroke conformation and the conformational transition between them. First, we have used normal mode analysis to identify a single normal mode that captures the coupled motions of AAA1-AAA2 closing and linker domain rotation, which enables the ATP-driven recovery stroke of dynein. Second, based on the post-powerstroke conformation solved crystallographically, we have modeled dynein's pre-powerstroke conformation by computationally inducing AAA1-AAA2 closing and sliding of coiled coil stalk, and the resulting model features a linker domain near the pre-powerstroke position and a slightly tilted stalk. Third, we have modeled the conformational transition from pre- to post-powerstroke conformation, which predicts a clear sequence of structural events that couple microtubule binding, powerstroke and product release, and supports a signaling path from stalk to AAA1 via AAA3 and AAA4. Finally, we have found that a closed AAA3-AAA4 interface (compatible with nucleotide binding) is essential to the mechano-chemical coupling in dynein. Our modeling not only offers unprecedented structural insights to the motor function of dynein as described by past single-molecule, fluorescence resonance energy transfer, and electron microscopy studies, but also provides new predictions for future experiments to test.     

The amino terminus of cGMP-dependent protein kinase Iβ increases the dynamics of the protein's cGMP-binding pockets

The type I cGMP-dependent protein kinases play critical roles in regulating vascular tone, platelet activation and synaptic plasticity. PKG I α and PKG Iβ differ in their first ~100 amino acids giving each isoform unique dimerization and autoinhibitory domains with identical cGMP-binding pockets and catalytic domains. The N-terminal leucine zipper and autoinhibitory domains have been shown to mediate isoform specific affinity for cGMP. PKG Iα has a >10 fold higher affinity for cGMP than PKG Iβ, and PKG Iβ that is missing its leucine zipper has a three-fold decreased affinity for cGMP. The exact mechanism through which the N-terminus of PKG alters cGMP-affinity is unknown. In the present study, we have used deuterium exchange mass spectrometry to study how PKG Iβ's N-terminus affects the conformation and dynamics of its cGMP-binding pockets. We found that the N-terminus increases the rate of deuterium exchange throughout the cGMP-binding domain. Our results suggest that the N-terminus shifts the conformational dynamics of the binding pockets, leading to an "open" conformation that has an increased affinity for cGMP.     

Discovering New Classes of <Emphasis Type="Italic">Brugia malayi</Emphasis> Asparaginyl-tRNA Synthetase Inhibitors and Relating Specificity to Conformational Change

SLIDE software, which models the flexibility of protein and ligand side chains while docking, was used to screen several large databases to identify inhibitors of Brugia malayi asparaginyl-tRNA synthetase (AsnRS), a target for anti-parasitic drug design. Seven classes of compounds identified by SLIDE were confirmed as micromolar inhibitors of the enzyme. Analogs of one of these classes of inhibitors, the long side-chain variolins, cannot bind to the adenosyl pocket of the closed conformation of AsnRS due to steric clashes, though the short side-chain variolins identified by SLIDE␣apparently bind isosterically with adenosine. We hypothesized that an open conformation of the motif 2 loop also permits the long side-chain variolins to bind in the adenosine pocket and that their selectivity for Brugia relative to human AsnRS can be explained by differences in the sequence and conformation of this loop. Loop flexibility sampling using Rigidity Optimized Conformational Kinetics (ROCK) confirms this possibility, while scoring of the relative affinities of the different ligands by SLIDE correlates well with the compounds’ ranks in inhibition assays. Combining ROCK and SLIDE provides a promising approach for exploiting conformational flexibility in structure-based screening and design of species selective inhibitors.     

Backbone motions of free and pheromone-bound major urinary protein I studied by molecular dynamics simulation

Molecular motions of free and pheromone-bound mouse major urinary protein I, previously investigated by NMR relaxation, were simulated in 30 ns molecular dynamics (MD) runs. The backbone flexibility was described in terms of order parameters and correlation times, commonly used in the NMR relaxation analysis. Special attention was paid to the effect of conformational changes on the nanosecond time scale. Time-dependent order parameters were determined in order to separate motions occurring on different time scales. As an alternative approach, slow conformational changes were identified from the backbone torsion angle variances, and "conformationally filtered" order parameters were calculated for well-defined conformation states. A comparison of the data obtained for the free and pheromone-bound protein showed that some residues are more rigid in the bound form, but a larger portion of the protein becomes more flexible upon the pheromone binding. This finding is in general agreement with the NMR results. The higher flexibility observed on the fast (fs-ps) time scale was typically observed for the residues exhibiting higher conformational freedom on the ns time scale. An inspection of the hydrogen bond network provided a structural explanation for the flexibility differences between the free and pheromone-bound proteins in the simulations.     

分子动力学模拟别构抑制剂Efavirenz对HIV-1逆转录酶的作用

为了理解非核苷类逆转录酶抑制剂(NNRTIs)与HIV-1逆转录酶(RT)的相互作用机制,利用新力场ff12SB对未结合和结合Efavirenz (EFV)逆转录酶的三种RT大分子体系分别进行了100 ns的长时间动力学模拟。通过分析EFV对RT结构的影响、不同残基柔性和不同体系构象的动力学行为等,发现EFV的结合会导致RT结构变化,从而影响RT的活性;证实了EFV的“分子楔”作用;还发现EFV的结合不但引起“拇指关节炎”,而且引起轻度“手指关节炎”;整个模拟过程中没有出现不同构象间的跃迁,但是无别构分子时的RT张开构象表现出明显的闭合倾向。这些结果有助于理解NNRTIs的抑制机制和RT构象变化的动力学性质。另外,还比较分析了模拟方法对计算结果的影响,对大分子体系的动力学模拟具有重要借鉴意义。     

Single molecule conformational dynamics of adenylate kinase: energy landscape, structural correlations, and transition state ensembles

We developed a coarse grained two-well model to study the single molecule protein conformational dynamics in microscopic detail at the residue level, overcoming the often encountered computational bottleneck. In particular, we explored the underlying conformational energy landscape of adenylate kinase, a crucial protein for signal transduction in the cell, and identified two major kinetic pathways for the conformational switch between open and closed states through either the intermediate state or the transient state. Based on the parameters fitted to the room-temperature experimental data, we predicted open and closed kinetic rates at the whole temperature ranges from 10 to 50 degrees C, which agree well with the experimental turnover numbers. After uncovering the underlying mechanism for conformational dynamics and exploring the structural correlations, we found the crucial dynamical interplay between the nucleoside monophosphate binding domain (NMP) and the ATP-binding domain (LID) in controlling the conformational switch. The key residues and contacts responsible for the conformational transitions are identified by following the time evolution of the two-dimensional spatial contact maps and characterizing the transition state as well as intermediate structure ensembles through phi value analysis. Our model provides a general framework to study the conformational dynamics of biomolecules and can be applied to many other systems.     

Characterization of the dynamics of an essential helix in the U1A protein by time-resolved fluorescence measurements

The RNA recognition motif (RRM), one of the most common RNA-binding domains, recognizes single-stranded RNA. A C-terminal helix that undergoes conformational changes upon binding is often an important contributor to RNA recognition. The N-terminal RRM of the U1A protein contains a C-terminal helix (helix C) that interacts with the RNA-binding surface of a beta-sheet in the free protein (closed conformation), but is directed away from this beta-sheet in the complex with RNA (open conformation). The dynamics of helix C in the free protein have been proposed to contribute to binding affinity and specificity. We report here a direct investigation of the dynamics of helix C in the free U1A protein on the nanosecond time scale using time-resolved fluorescence anisotropy. The results indicate that helix C is dynamic on a 2-3 ns time scale within a 20 degrees range of motion. Steady-state fluorescence experiments and molecular dynamics simulations suggest that the dynamical motion of helix C occurs within the closed conformation. Mutation of a residue on the beta-sheet that contacts helix C in the closed conformation dramatically destabilizes the complex (Phe56Ala) and alters the steady-state fluorescence, but not the time-resolved fluorescence anisotropy, of a Trp in helix C. Mutation of Asp90 in the hinge region between helix C and the remainder of the protein to Ala or Gly subtly alters the dynamics of the U1A protein and destabilizes the complex. Together these results show that helix C maintains a dynamic closed conformation that is stable to these targeted protein modifications and does not equilibrate with the open conformation on the nanosecond time scale.