Molecular Dynamics of Apo-Adenylate Kinase: A Principal Component Analysis

Identification of Potent ADCK3 Inhibitors through Structure-Based Virtual Screening

ADCK3 is a member of the UbiB family of atypical protein kinases in humans, with homologues in archaea, bacteria, and eukaryotes. In lieu of protein kinase activity, ADCK3 plays a role in the biosynthesis of coenzyme Q10 (CoQ10), and inactivating mutations can cause a CoQ10 deficiency and ataxia. However, the exact functions of ADCK3 are still unclear, and small-molecule inhibitors could be useful as chemical probes to elucidate its molecular mechanisms. In this study, we applied structure-based virtual screening (VS) to discover a novel chemical series of ADCK3 inhibitors. Through extensive structural analysis of the active-site residues, we developed a pharmacophore model and applied it to a large-scale VS. Out of ∼170,000 compounds virtually screened, 800 top-ranking candidate compounds were selected and tested in both ADCK3 and p38 biochemical assays for hit validation. In total, 129 compounds were confirmed as ADCK3 inhibitors, and among them, 114 compounds are selective against p38, which was used as a counter-target. Molecular dynamics (MD) simulations were then conducted to predict the binding modes of the most potent compounds within the ADCK3 active site. Through metadynamics analysis, we successfully detected the key amino acid residues that govern intermolecular interactions. The findings provided in this study can serve as a promising starting point for drug development.

Predicting binding events in very flexible, allosteric, multi-domain proteins

Knowledge of the structures formed by proteins and small ligands is of fundamental importance for understanding molecular principles of chemotherapy and for designing new and more effective drugs. Due to the still high costs and to the several limitations of experimental techniques, it is most often desirable to predict these ligand-protein complexes in silico, particularly when screening for new putative drugs from databases of millions of compounds. While virtual screening based on molecular docking is widely used for this purpose, it generally fails in mimicking binding events associated with large conformational changes in the protein, particularly when the latter involve multiple domains. In this work, we describe a new methodology aimed at generating bound-like conformations of very flexible and allosteric proteins bearing multiple binding sites. Validation was performed on the enzyme adenylate kinase (ADK), a paradigmatic example of proteins that undergo very large conformational changes upon ligand binding. By only exploiting the unbound structure and the putative binding sites of the protein, we generated a significant fraction of bound-like structures, which employed in ensemble-docking calculations allowed to find native-like poses of substrates, inhibitors, and catalytically incompetent binders. Our protocol provides a general framework for the generation of bound-like conformations of flexible proteins that are suitable to host different ligands, demonstrating high sensitivity to the fine chemical details that regulate protein's activity. We foresee applications in virtual screening for difficult targets, prediction of the impact of amino acid mutations on structure and dynamics, and protein engineering.

Molecular insights into the inhibition mechanism of harringtonine against essential proteins associated with SARS-CoV-2 entry

Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) has recently posed a serious threat to global public health. Harringtonine (HT), as a small-molecule antagonist, has antiviral activity against a variety of viruses. There is evidence that HT can inhibit the SARS-CoV-2 entry into host cells by blocking the Spike protein and transmembrane protease serine 2 (TMPRSS2). However, the molecular mechanism underlying the inhibition effect of HT is largely elusive. Here, docking and all-atom molecular dynamics simulations were used to investigate the mechanism of HT against the receptor binding domain (RBD) of Spike, TMPRSS2, as well as the complex of RBD and angiotensin-converting enzyme 2 complex (RBD-ACE2). The results reveal that HT binds to all proteins primarily through hydrogen bond and hydrophobic interactions. Binding with HT influences the structural stability and dynamic motility processes of each protein. The interactions of HT with residues N33, H34 and K353 of ACE2, and residue K417 and Y453 of RBD contribute to disrupting the binding affinity between RBD and ACE2, which may hinder the virus entry into host cells. Our research provides molecular insights into the inhibition mechanism of HT against SARS-CoV-2 associated proteins, which will help for the novel antiviral drugs development.

A comparative study of 5- fluorouracil, doxorubicin, methotrexate, paclitaxel for their inhibition ability for Mpro of nCoV: Molecular docking and molecular dynamics simulations

A new corona virus (nCoV) is aetiological agent responsible for the viral pneumonia epidemic. Three is no specific therapeutic medicines available for the treatment of this condition and also effective treatment choices are few. In this work author tried to investigate some repurposing drug such as 5- fluorouracil, doxorubicin, methotrexate and paclitaxel against the main protease (Mpro) of nCoV by the computational model. Molecular docking was performed to screen out the best compound and doxorubicin was found to have minimum binding energy −121.89 kcal/mol. To further study, MD simulations were performed at 300 K and the result successfully corroborate the energy obtained by molecular docking. Temperature dependent MD simulation of the best molecule that is doxorubicin obtained from docking result was performed to check the variation in structural changes in Mpro of nCoV at 290 K, 310 K, 320 K and 325 K. It is sound that doxorubicin binds effectively with Mpro of nCoV at 290 K. Further ADME properties of the 5- fluorouracil, doxorubicin, methotrexate and paclitaxel were also evaluated to understand the bioavailability.

Importance of Three-Body Problems and Protein-Protein Interactions in Proteolysis-Targeting Chimera Modeling: Insights from Molecular Dynamics Simulations

Proteolysis-targeting chimeras (PROTACs) are a class of bifunctional molecules that can induce the ubiquitin degradation of its target protein by hijacking the E3 ligase to form a target protein-PROTAC-E3 ligase ternary complex. Its underlying principle has inspired the development of a wide range of protein degraders that are similar to or beyond PROTACs in recent years. The formation of the ternary complexes is the key to the success of PROTAC-induced protein degradation. Nevertheless, the lack of effective ternary complex modeling techniques has limited the application of computer-aided drug discovery tools to this emerging and fast developing new land in drug industry. Thus, in this study, we explored the application of the more physically sound molecular dynamics simulation and the molecular mechanics combined with the generalized Born and surface area continuum solvation (MM/GBSA) method to solve the underlying three-body problem in PROTAC modeling. We first verified the accuracy of our approach using a series of known Brd4 BD2 degraders. The calculated binding energy showed a good correlation with the experimental K_d values. The modeling of a unique property, namely, the α value, for PROTACs was also first and accurately performed to our best knowledge. The results also demonstrated the importance of PROTAC-induced protein-protein interactions in its modeling, either qualitatively or quantitatively. Finally, by standing on the success of earlier docking-based approaches, our protocol was also applied as a rescoring function in pose prediction. The results showed a notable improvement in reranking the initial poses generated from a modified Rosetta method, which was reportedly one of the best among a handful of PROTAC modeling approaches available in this field. We hope this work could provide a practical protocol and more insights to study the binding and the design of PROTACs and other protein degraders.

Tracking the ATP-binding response in adenylate kinase in real time

The biological function of proteins is critically dependent on dynamics inherent to the native structure. Such structural dynamics obey a predefined order and temporal timing to execute the specific reaction. Determination of the cooperativity of key structural rearrangements requires monitoring protein reactions in real time. In this work, we used time-resolved x-ray solution scattering (TR-XSS) to visualize structural changes in the Escherichia coli adenylate kinase (AdK) enzyme upon laser-induced activation of a protected ATP substrate. A 4.3-ms transient intermediate showed partial closing of both the ATP- and AMP-binding domains, which indicates a cooperative closing mechanism. The ATP-binding domain also showed local unfolding and breaking of an Arg¹³¹-Asp¹⁴⁶ salt bridge. Nuclear magnetic resonance spectroscopy data identified similar unfolding in an Arg¹³¹Ala AdK mutant, which refolded in a closed, substrate-binding conformation. The observed structural dynamics agree with a “cracking mechanism” proposed to underlie global structural transformation, such as allostery, in proteins.

Comparison between slow anisotropic LE4PD fluctuations and the principal component analysis modes of ubiquitin

The biological function and folding mechanisms of proteins are often guided by large-scale slow motions, which involve crossing high energy barriers. In a simulation trajectory, these slow fluctuations are commonly identified using a principal component analysis (PCA). Despite the popularity of this method, a complete analysis of its predictions based on the physics of protein motion has been so far limited. This study formally connects the PCA to a Langevin model of protein dynamics and analyzes the contributions of energy barriers and hydrodynamic interactions to the slow PCA modes of motion. To do so, we introduce an anisotropic extension of the Langevin equation for protein dynamics, called the LE4PD-XYZ, which formally connects to the PCA "essential dynamics." The LE4PD-XYZ is an accurate coarse-grained diffusive method to model protein motion, which describes anisotropic fluctuations in the alpha carbons of the protein. The LE4PD accounts for hydrodynamic effects and mode-dependent free-energy barriers. This study compares large-scale anisotropic fluctuations identified by the LE4PD-XYZ to the mode-dependent PCA predictions, starting from a microsecond-long alpha carbon molecular dynamics atomistic trajectory of the protein ubiquitin. We observe that the inclusion of free-energy barriers and hydrodynamic interactions has important effects on the identification and timescales of ubiquitin's slow modes.

Exploring Conformational Change of Adenylate Kinase by Replica Exchange Molecular Dynamic Simulation

Replica exchange molecular dynamics (REMD) simulation is a popular enhanced sampling method that is widely used for exploring the atomic mechanism of protein conformational change. However, the requirement of huge computational resources for REMD, especially with the explicit solvent model, largely limits its application. In this study, the availability and efficiency of a variant of velocity-scaling REMD (vsREMD) was assessed with adenylate kinase as an example. Although vsREMD achieved results consistent with those from conventional REMD and experimental studies, the number of replicas required for vsREMD (30) was much less than that for conventional REMD (80) to achieve a similar acceptance rate (∼0.2), demonstrating high efficiency of vsREMD to characterize the protein conformational change and associated free-energy profile. Thus, vsREMD is a highly efficient approach for studying the large-scale conformational change of protein systems.

Structural and Dynamic Characterizations Highlight the Deleterious Role of SULT1A1 R213H Polymorphism in Substrate Binding

Sulfotransferase 1A1 (SULT1A1) is responsible for catalyzing various types of endogenous and exogenous compounds. Accumulating data indicates that the polymorphism rs9282861 (R213H) is responsible for inefficient enzymatic activity and associated with cancer progression. To characterize the detailed functional consequences of this mutation behind the loss-of-function of SULT1A1, the present study deployed molecular dynamics simulation to get insights into changes in the conformation and binding energy. The dynamics scenario of SULT1A1 in both wild and mutated types as well as with and without ligand showed that R213H induced local conformational changes, especially in the substrate-binding loop rather than impairing overall stability of the protein structure. The higher conformational changes were observed in the loop3 (residues, 235-263), turning loop conformation to A-helix and B-bridge, which ultimately disrupted the plasticity of the active site. This alteration reduced the binding site volume and hydrophobicity to decrease the binding affinity of the enzyme to substrates, which was highlighted by the MM-PBSA binding energy analysis. These findings highlight the key insights of structural consequences caused by R213H mutation, which would enrich the understanding regarding the role of SULT1A1 mutation in cancer development and also xenobiotics management to individuals in the different treatment stages.

Namdinator - automatic molecular dynamics flexible fitting of structural models into cryo-EM and crystallography experimental maps

Model building into experimental maps is a key element of structural biology, but can be both time consuming and error prone for low-resolution maps. Here we present Namdinator, an easy-to-use tool that enables the user to run a molecular dynamics flexible fitting simulation followed by real-space refinement in an automated manner through a pipeline system. Namdinator will modify an atomic model to fit within cryo-EM or crystallography density maps, and can be used advantageously for both the initial fitting of models, and for a geometrical optimization step to correct outliers, clashes and other model problems. We have benchmarked Namdinator against 39 deposited cryo-EM models and maps, and observe model improvements in 34 of these cases (87%). Clashes between atoms were reduced, and the model-to-map fit and overall model geometry were improved, in several cases substantially. We show that Namdinator is able to model large-scale conformational changes compared to the starting model. Namdinator is a fast and easy tool for structural model builders at all skill levels. Namdinator is available as a web service (https://namdinator.au.dk), or it can be run locally as a command-line tool.

Electrostatic interactions determine entrance/release order of substrates in the catalytic cycle of adenylate kinase

Adenylate kinase is a monomeric phosphotransferase with important biological function in regulating concentration of adenosine triphosphate (ATP) in cells, by transferring the terminal phosphate group from ATP to adenosine monophosphate (AMP) and forming two adenosine diphosphate (ADP) molecules. During this reaction, the kinase may undergo a large conformational transition, forming different states with its substrates. Although many structures of the protein are available, atomic details of the whole process remain unclear. In this article, we use both conventional molecular dynamics (MD) simulation and an enhanced sampling technique called parallel cascade selection MD simulation to explore different conformational states of the Escherichia coli adenylate kinase. Based on the simulation results, we propose a possible entrance/release order of substrates during the catalytic cycle. The substrate-free protein prefers an open conformation, but changes to a closed state once ATP·Mg enters into its binding pocket first and then AMP does. After the reaction of ATP transferring the terminal phosphate group to AMP, ADP·Mg and ADP are released sequentially, and finally the whole catalyze cycle is completed. Detailed contact and distance analysis reveals that the entrance/release order of substrates may be largely controlled by electrostatic interactions between the protein and the substrates.

BIM Binding Remotely Regulates BAX Activation: Insights from the Free Energy Landscapes

Activation of the pro-apoptotic BAX protein, a BCL-2 family member, is known to trigger apoptosis by forming pores in the mitochondrial outer membrane (MOM). While in the cytosol, release of its transmembrane C-terminal helix (called α9 helix) from a well-characterized binding pocket (BC groove) and subsequent permeabilization of the MOM are understood to be the initiating events of the activation. Concerning what initiates BAX activation, so far one plausible suggestion has been that the transient attachment of BH3-only peptide at a distal site from the BC groove triggers the activation process. Yet how this pivotal step displaces α9 from the BC groove has remained poorly understood. Using a combination of standard molecular dynamics and enhanced sampling methods, the energy landscape of BIM (BH3-only peptide) induced BAX activation has been computed, and the molecular origin of those events is hereby reported in atomistic detail. The simulated transition pathway of α9 release reveals that BIM subdues the energetic cost of the process by reducing the activation energy barrier to some extent but mostly by minimizing the free energy difference between the active (α9-released) and inactive (α9-bound) states. Interestingly, the flexibility of the α9 helix itself plays a decisive role in this mechanism. The impact of BIM encounter at the distal site is found to propagate to the α9 (BC groove bound) mostly through conserved pathways of residue level interactions. Overall, the thermodynamic basis of the "hit-and-run" mechanism for activation of the BCL-2 family is presented reconciling the available biochemical observations.

Polymodal allosteric regulation of Type 1 Serine/Threonine Kinase Receptors via a conserved electrostatic lock

Type 1 Serine/Threonine Kinase Receptors (STKR1) transduce a wide spectrum of biological signals mediated by TGF-β superfamily members. The STKR1 activity is tightly controlled by their regulatory glycine-serine rich (GS) domain adjacent to the kinase domain. Despite decades of studies, it remains unknown how physiological or pathological GS domain modifications are coupled to STKR1 kinase activity. Here, by performing molecular dynamics simulations and free energy calculation of Activin-Like Kinase 2 (ALK2), we found that GS domain phosphorylation, FKBP12 dissociation, and disease mutations all destabilize a D354-R375 salt-bridge, which normally acts as an electrostatic lock to prevent coordination of adenosine triphosphate (ATP) to the catalytic site. We developed a WAFEX-guided principal analysis and unraveled how phosphorylation destabilizes this highly conserved salt-bridge in temporal and physical space. Using current-flow betweenness scores, we identified an allosteric network of residue-residue contacts between the GS domain and the catalytic site that controls the formation and disruption of this salt bridge. Importantly, our novel network analysis approach revealed how certain disease-causing mutations bypass FKBP12-mediated kinase inhibition to produce leaky signaling. We further provide experimental evidence that this salt-bridge lock exists in other STKR1s, and acts as a general safety mechanism in STKR1 to prevent pathological leaky signaling. In summary, our study provides a compelling and unifying allosteric activation mechanism in STKR1 kinases that reconciles a large number of experimental studies and sheds light on a novel therapeutic avenue to target disease-related STKR1 mutants.

Kinases play central role in essential physiological process and are attractive therapeutic drug targets. One of the important kinase families is Type 1 Serine/Threonine Kinase Receptors (STKR1), which control gene expression in response to extracellular growth factors. The activities of STKR1 are tightly controlled by their regulatory domain, which is distant from the kinase catalytic site. The underlying molecular mechanism is elucidated here. We identified that formation or disruption of a highly conserved charge-charge interaction located near the ATP binding site, mediates the physiological inhibition or activation of STKR1. We find that the stability of this charge-charge interaction is remotely controlled by interactions propagated from the distant regulatory domain. Several disease-causing mutations are located at the regulatory domain. We demonstrate how those mutations bypass these endogenous STKR1 inhibition mechanisms to produce pathological phenotypes. This study provides a general activation mechanism in STKR1 kinases, thus may benefit understanding the molecular mechanism of diseases and drug development.

JED: a Java Essential Dynamics Program for comparative analysis of protein trajectories

BACKGROUND

Essential Dynamics (ED) is a common application of principal component analysis (PCA) to extract biologically relevant motions from atomic trajectories of proteins. Covariance and correlation based PCA are two common approaches to determine PCA modes (eigenvectors) and their eigenvalues. Protein dynamics can be characterized in terms of Cartesian coordinates or internal distance pairs. In understanding protein dynamics, a comparison of trajectories taken from a set of proteins for similarity assessment provides insight into conserved mechanisms. Comprehensive software is needed to facilitate comparative-analysis with user-friendly features that are rooted in best practices from multivariate statistics.

RESULTS

We developed a Java based Essential Dynamics toolkit called JED to compare the ED from multiple protein trajectories. Trajectories from different simulations and different proteins can be pooled for comparative studies. JED implements Cartesian-based coordinates (cPCA) and internal distance pair coordinates (dpPCA) as options to construct covariance (Q) or correlation (R) matrices. Statistical methods are implemented for treating outliers, benchmarking sampling adequacy, characterizing the precision of Q and R, and reporting partial correlations. JED output results as text files that include transformed coordinates for aligned structures, several metrics that quantify protein mobility, PCA modes with their eigenvalues, and displacement vector (DV) projections onto the top principal modes. Pymol scripts together with PDB files allow movies of individual Q- and R-cPCA modes to be visualized, and the essential dynamics occurring within user-selected time scales. Subspaces defined by the top eigenvectors are compared using several statistical metrics to quantify similarity/overlap of high dimensional vector spaces. Free energy landscapes can be generated for both cPCA and dpPCA.

CONCLUSIONS

JED offers a convenient toolkit that encourages best practices in applying multivariate statistics methods to perform comparative studies of essential dynamics over multiple proteins. For each protein, Cartesian coordinates or internal distance pairs can be employed over the entire structure or user-selected parts to quantify similarity/differences in mobility and correlations in dynamics to develop insight into protein structure/function relationships.

Dynamic, structural and thermodynamic basis of insulin-like growth factor 1 kinase allostery mediated by activation loop phosphorylation

Despite the importance of kinases' catalytic activity regulation in cell signaling, detailed mechanisms underlying their activity regulation are poorly understood. Herein, using insulin-like growth factor 1 receptor kinase (IGF-1RK) as a model, the mechanisms of kinase regulation by its activation loop (A-loop) phosphorylation were investigated through molecular dynamics (MD) and alchemical free energy simulations. Analyses of the simulation results and free energy landscapes determined for the entire catalytic cycle of the kinase revealed that A-loop phosphorylation affects each step in the IGF-1RK catalytic cycle, including conformational change, substrate binding/product release and catalytic phosphoryl transfer. Specifically, the conformational equilibrium of the kinase is shifted by 13.2 kcal mol-1 to favor the active conformation after A-loop phosphorylation, which increases substrate binding affinity of the activated kinase. This free energy shift is achieved primarily via destabilization of the inactive conformation. The free energy of the catalytic reaction is also changed by 3.3 kcal mol-1 after the phosphorylation and in the end, facilitates product release. Analyses of MD simulations showed that A-loop phosphorylation produces these energetic effects by perturbing the side chain interactions around each A-loop tyrosine. These interaction changes are propagated to the remainder of the kinase to modify the orientations and dynamics of the αC-helix and A-loop, and together yield the observed free energy changes. Since many protein kinases share similar interactions identified in this work, the mechanisms of kinase allostery and catalysis unraveled here can be applicable to them.

Application of conventional molecular dynamics simulation in evaluating the stability of apomyoglobin in urea solution

In this study, we had exploited the advancement in computer technology to determine the stability of four apomyoglobin variants namely wild type, E109A, E109G and G65A/G73A by conducting conventional molecular dynamics simulations in explicit urea solution. Variations in RMSD, native contacts and solvent accessible surface area of the apomyoglobin variants during the simulation were calculated to probe the effect of mutation on the overall conformation of the protein. Subsequently, the mechanism leading to the destabilization of the apoMb variants was studied through the calculation of correlation matrix, principal component analyses, hydrogen bond analyses and RMSF. The results obtained here correlate well with the study conducted by Baldwin and Luo which showed improved stability of apomyoglobin with E109A mutation and contrariwise for E109G and G65A/G73A mutation. These positive observations showcase the feasibility of exploiting MD simulation in determining protein stability prior to protein expression.

Dynamics of fluoroquinolones induced resistance in DNA gyrase of Mycobacterium tuberculosis

DNA gyrase is a validated target of fluoroquinolones which are key components of multidrug resistance tuberculosis (TB) treatment. Most frequent occurring mutations associated with high level of resistance to fluoroquinolone in clinical isolates of TB patients are A90V, D94G, and A90V-D94G (double mutant [DM]), present in the larger subunit of DNA Gyrase. In order to explicate the molecular mechanism of drug resistance corresponding to these mutations, molecular dynamics (MD) and mechanics approach was applied. Structure-based molecular docking of complex comprised of DNA bound with Gyrase A (large subunit) and Gyrase C (small subunit) with moxifloxacin (MFX) revealed high binding affinity to wild type with considerably high Glide XP docking score of -7.88 kcal/mol. MFX affinity decreases toward single mutants and was minimum toward the DM with a docking score of -3.82 kcal/mol. Docking studies were also performed against 8-Methyl-moxifloxacin which exhibited higher binding affinity against wild and mutants DNA gyrase when compared to MFX. Molecular Mechanics/Generalized Born Surface Area method predicted the binding free energy of the wild, A90V, D94G, and DM complexes to be -55.81, -25.87, -20.45, and -12.29 kcal/mol, respectively. These complexes were further subjected to 30 ns long MD simulations to examine significant interactions and conformational flexibilities in terms of root mean square deviation, root mean square fluctuation, and strength of hydrogen bond formed. This comparative drug interaction analysis provides systematic insights into the mechanism behind drug resistance and also paves way toward identifying potent lead compounds that could combat drug resistance of DNA gyrase due to mutations.

Clarifying binding difference of ATP and ADP to extracellular signal-regulated kinase 2 by using molecular dynamics simulations

Extracellular signal-regulated kinase 2 is a promising target for designs and development of anticancer drugs. Molecular dynamics simulations and molecular mechanics Poisson-Boltzmann method were applied to study binding difference of ADP and ATP to extracellular signal-regulated kinase 2. The results prove that the binding ability of ATP to extracellular signal-regulated kinase 2 is stronger than that of ADP. Principal component analysis performed by using molecular dynamics trajectories suggests that binding of ADP and ATP to extracellular signal-regulated kinase 2 change motion directions of two helices α1 and α2. Residue-based free energy decomposition method was adopted to calculate contributions of separate residues to associations of ADP and ATP with extracellular signal-regulated kinase 2. The results show that ADP and ATP produce strong CH-π interactions with five residues Ile29, Val37, Ala50, Leu105, and Leu154. In addition, five hydrogen bonding interactions of ADP and ATP with residues Lys52, Gln103, Asp104, and Met106 also stabilize bindings of ADP and ATP to extracellular signal-regulated kinase 2. Overall, the CH-π interactions of ATP with five residues Ile29, Val37, Ala50, Leu105, and Leu154 are stronger than ADP. This study is expected to contribute a significant theoretical hint for designs of anticancer drugs targeting extracellular signal-regulated kinase 2.

Substrate Binding Specifically Modulates Domain Arrangements in Adenylate Kinase

The enzyme adenylate kinase (ADK) features two substrate binding domains that undergo large-scale motions during catalysis. In the apo state, the enzyme preferentially adopts a globally open state with accessible binding sites. Binding of two substrate molecules (AMP + ATP or ADP + ADP) results in a closed domain conformation, allowing efficient phosphoryl-transfer catalysis. We employed molecular dynamics simulations to systematically investigate how the individual domain motions are modulated by the binding of substrates. Two-dimensional free-energy landscapes were calculated along the opening of the two flexible lid domains for apo and holo ADK as well as for all single natural substrates bound to one of the two binding sites of ADK. The simulations reveal a strong dependence of the conformational ensembles on type and binding position of the bound substrates and a nonsymmetric behavior of the lid domains. Altogether, the ensembles suggest that, upon initial substrate binding to the corresponding lid site, the opposing lid is maintained open and accessible for subsequent substrate binding. In contrast, ATP binding to the AMP-lid induces global domain closing, preventing further substrate binding to the ATP-lid site. This might constitute a mechanism by which the enzyme avoids the formation of a stable but enzymatically unproductive state.

A principal component analysis of the dynamics of subdomains and binding sites in human serum albumin

The conformational dynamics of human serum albumin (HSA) was investigated by principal component analysis (PCA) applied to three molecular dynamics trajectories of 200 ns each. The overlap of the essential subspaces spanned by the first 10 principal components (PC) of different trajectories was about 0.3 showing that the PCA based on a trajectory length of 200 ns is not completely convergent for this protein. The contributions of the relative motion of subdomains and of the subdomains (internal) distortion to the first 10 PCs were found to be comparable. Based on the distribution of the first 3 PC, 10 protein conformers are identified showing relative root mean square deviations (RMSD) between 2.3 and 4.6 Å. The main PCs are found to be delocalized over the whole protein structure indicating that the motions of different protein subdomains are coupled. This coupling is considered as being related to the allosteric effects observed upon ligand binding to HSA. On the other hand, the first PC of one of the three trajectories describes a conformational transition of the protein domain I that is close to that experimentally observed upon myristate binding. This is a theoretical support for the older hypothesis stating that changes of the protein onformation favorable to binding can precede the ligand complexation. A detailed all atoms PCA performed on the primary Sites 1 and 2 confirms the multiconformational character of the HSA binding sites as well as the significant coupling of their motions.

Factor Analysis of Conformations and NMR Signals of Rotaxanes: AIMD and Polarizable MD Simulations

The interlocked ⟨rod | ring⟩ structures of pseudorotaxanes and [2]rotaxanes are usually maintained by the complex hydrogen-bonding (H-bonding) network between the rod and ring. Ab initio molecular dynamics (AIMD) using generalized energy-based fragmentation approach and polarizable force field (polar FF)-based molecular dynamics (MD) simulations were performed to investigate the conformational changes of mechanically interlocked systems and to obtain the ensemble-averaged NMR chemical shifts. Factor analysis (FA) demonstrates that the ring H-donor (2,6 pyridinedicarboxamide group) plays an important role in the ring-rod recognition. In comparison to the conventional fixed-charge force field, the polarization effect is crucial to account for the H-bonding interactions in supramolecular systems. In the hybrid scheme, the polar FF-based MD simulations are used to generate different initial states for the AIMD simulations, which are able to give better prediction of ensemble-averaged NMR signals for chemically equivalent amide protons. The magnitude of the deshielding shift of NMR signal is correlated with the length of hydrogen bond. The polar FF model with variable charges shows that the dipole-dipole interactions between the flexible diethylene glycol chain of ring and polar solvents induce the upfield shifts of NMR signals of rod H-donors and the directional distribution of the neighboring CH3CN solvents.

Revealing an outward-facing open conformational state in a CLC Cl(-)/H(+) exchange transporter

CLC secondary active transporters exchange Cl^- for H⁺. Crystal structures have suggested that the conformational change from occluded to outward-facing states is unusually simple, involving only the rotation of a conserved glutamate (Glu_ex) upon its protonation. Using ¹⁹F NMR, we show that as [H⁺] is increased to protonate Glu_ex and enrich the outward-facing state, a residue ~20 Å away from Glu_ex, near the subunit interface, moves from buried to solvent-exposed. Consistent with functional relevance of this motion, constriction via inter-subunit cross-linking reduces transport. Molecular dynamics simulations indicate that the cross-link dampens extracellular gate-opening motions. In support of this model, mutations that decrease steric contact between Helix N (part of the extracellular gate) and Helix P (at the subunit interface) remove the inhibitory effect of the cross-link. Together, these results demonstrate the formation of a previously uncharacterized 'outward-facing open' state, and highlight the relevance of global structural changes in CLC function.

DOI:http://dx.doi.org/10.7554/eLife.11189.001

Cells have transporter proteins on their surface to carry molecules in and out of the cell. For example, the CLC family of transporters move two chloride ions in one direction at the same time as moving one hydrogen ion in the opposite direction.

To be able to move these ions in opposite directions, transporters have to cycle through a series of shapes in which the ions can only access alternate sides of the membrane. First, the transporter adopts an 'outward-facing' shape when the ions first bind to the transporter, then it switches into the 'occluded' shape to move the ions through the membrane. Finally, the transporter takes on the 'inward-facing' shape to release the ions on the other side of the membrane. However, structural studies of CLCs suggest that the structures of these proteins do not change much while they are moving ions, which suggests that they might work in a different way.

Khantwal, Abraham et al. have now used techniques called “nuclear magnetic resonance” and "double electron-electron resonance" to investigate how a CLC from a bacterium moves ions. The experiments suggest that when the transporter adopts the outward-facing shape, points on the protein known as Y419 and D417 shift their positions. Chemically linking two regions of the CLC prevented this movement and inhibited the transport of chloride ions across the membrane.

Khantwal, Abraham et al. then used a computer simulation to model how the protein changes shape in more detail. This model predicts that two regions of the transporter undergo major rearrangements resulting in a gate-opening motion that widens a passage to allow the chloride ions to bind to the protein. Khantwal, Abraham et al.’s findings will prompt future studies to reveal the other shapes and how CLCs transition between them.

DOI:http://dx.doi.org/10.7554/eLife.11189.002

Structural and mechanistic insight into substrate binding from the conformational dynamics in apo and substrate-bound DapE enzyme

Conformational dynamics induced by substrate binding in DapE enzyme.

Extrapolation of Inter Domain Communications and Substrate Binding Cavity of Camel HSP70 1A: A Molecular Modeling and Dynamics Simulation Study

Heat shock protein 70 (HSP70) is an important chaperone, involved in protein folding, refolding, translocation and complex remodeling reactions under normal as well as stress conditions. However, expression of HSPA1A gene in heat and cold stress conditions associates with other chaperons and perform its function. Experimental structure for Camel HSP70 protein (cHSP70) has not been reported so far. Hence, we constructed 3D models of cHSP70 through multi- template comparative modeling with HSP110 protein of S. cerevisiae (open state) and with HSP70 protein of E. coli 70kDa DnaK (close state) and relaxed them for 100 nanoseconds (ns) using all-atom Molecular Dynamics (MD) Simulation. Two stable conformations of cHSP70 with Substrate Binding Domain (SBD) in open and close states were obtained. The collective mode analysis of different transitions of open state to close state and vice versa was examined via Principal Component Analysis (PCA) and Minimum Distance Matrix (MDM). The results provide mechanistic representation of the communication between Nucleotide Binding Domain (NBD) and SBD to identify the role of sub domains in conformational change mechanism, which leads the chaperone cycle of cHSP70. Further, residues present in the chaperon functioning site were also identified through protein-peptide docking. This study provides an overall insight into the inter domain communication mechanism and identification of the chaperon binding cavity, which explains the underlying mechanism involved during heat and cold stress conditions in camel.

The Atomistic Mechanism of Conformational Transition of Adenylate Kinase Investigated by Lorentzian Structure-Based Potential

We present a new all-atom structure-based method to study protein conformational transitions using Lorentzian attractive interactions based on native structures. The variability of each native contact is estimated based on evolutionary information using a machine learning method. To test the validity of this approach, we have investigated the conformational transition of adenylate kinase (ADK). The intrinsic boundedness of the Lorentzian attractive interactions facilitated frequent conformational transitions, and consequently we were able to observe more than 1000 structural interconversions between the open and closed states of ADK out of a total of 6 μs MD simulations. ADK has three domains: the nucleoside monophosphate (NMP) binding domain, the LID-domain, and the CORE domain, which catalyze the interconversion between ATP and ADP. We identified two transition states: a more frequent LID-closed-NMP-open (TS1) state and a less frequent LID-open-NMP-closed (TS2) state. The transition was found to be symmetric in both directions via TS1. We also obtained an off-pathway metastable state that was previously observed with physics-based all-atom simulations but not with coarse-grained models. In the metastable state, the LID domain was slightly twisted and formed contacts with the NMP domain. Our model correctly identified a total of 14 out of the top 16 residues with highest fluctuation by NMR experiment, thus showing excellent agreement with experimental NMR relaxation data and overwhelmingly better results than existing models.

Features of large hinge-bending conformational transitions. Prediction of closed structure from open state

We performed a detailed analysis of conformational transition pathways for a set of 10 proteins, which undergo large hinge-bending-type motions with 4-12 Å RMSD (root mean-square distance) between open and closed crystal structures. Anisotropic network model-Monte Carlo (ANM-MC) algorithm generates a targeted pathway between two conformations, where the collective modes from the ANM are used for deformation at each iteration and the conformational energy of the deformed structure is minimized via an MC algorithm. The target structure was approached successfully with an RMSD of 0.9-4.1 Å when a relatively low cutoff radius of 10 Å was used in ANM. Even though one predominant mode (first or second) directed the open-to-closed conformational transition, changes in the dominant mode character were observed for most cases along the transition. By imposing radius of gyration constraint during mode selection, it was possible to predict the closed structure for eight out of 10 proteins (with initial 4.1-7.1 Å and final 1.7-2.9 Å RMSD to target). Deforming along a single mode leads to most successful predictions. Based on the previously reported free energy surface of adenylate kinase, deformations along the first mode produced an energetically favorable path, which was interestingly facilitated by a change in mode shape (resembling second and third modes) at key points. Pathway intermediates are provided in our database of conformational transitions (http://safir.prc.boun.edu.tr/anmmc/method/1).

Computational study on the intramolecular self-organization of the macrorings of some ‘giant’ cyclodextrins (CDn, n = 40, 70, 85, 100)

A limited number of modes determine the overall deformations of the macrorings, which may have more than one cavity. Accordingly, they have the potential to accommodate more than one substrate molecule.

Understanding the recognition mechanisms of Zα domain of human editing enzyme ADAR1 (hZα(ADAR1)) and various Z-DNAs from molecular dynamics simulation

The Z-DNA-binding domain of human double-stranded RNA adenosine deaminase I (hZαADAR1) can specifically recognize the left-handed Z-DNA which preferentially occurs at alternating purine-pyrimidine repeats, especially the CG-repeats. The interactions of hZαADAR1 and Z-DNAs in different sequence contexts can affect many important biological functions including gene regulation and chromatin remodeling. Therefore it is of great necessity to fully understand their recognition mechanisms. However, most existing studies are aimed at the standard CG-repeat Z-DNA rather than the non-CG-repeats, and whether the molecular basis of hZαADAR1 binding to various Z-DNAs are identical or not is still unclear on the atomic level. Here, based on the recently determined crystal structures of three representative non-CG-repeat Z-DNAs (d(CACGTG)2, d(CGTACG)2 and d(CGGCCG)2) in complex with hZαADAR1, 40 ns molecular dynamics simulation together with binding free energy calculation were performed for each system. For comparison, the standard CG-repeat Z-DNA (d(CGCGCG)2) complexed with hZαADAR1 was also simulated. The consistent results demonstrate that nonpolar interaction is the driving force during the protein-DNA binding process, and that polar interaction mainly from helix α3 also provides important contributions. Five common hot-spot residues were identified, namely Lys169, Lys170, Asn173, Arg174 and Tyr177. Hydrogen bond analysis coupled with surface charge distribution further reveal the interfacial information between hZαADAR1 and Z-DNA in detail. All of the analysis illustrate that four complexes share the common key features and the similar binding modes irrespective of Z-DNA sequences, suggesting that Z-DNA recognition by hZαADAR1 is conformation-specific rather than sequence-specific. Additionally, by analyzing the conformational changes of hZαADAR1, we found that the binding of Z-DNA could effectively stabilize hZαADAR1 protein. Our study can provide some valuable information for better understanding the binding mechanism between hZαADAR1 or even other Z-DNA-binding protein and Z-DNA.

Molecular dynamics studies on the conformational transitions of adenylate kinase: a computational evidence for the conformational selection mechanism

Escherichia coli adenylate kinase (ADK) is a monomeric phosphotransferase enzyme that catalyzes reversible transfer of phosphoryl group from ATP to AMP with a large-scale domain motion. The detailed mechanism for this conformational transition remains unknown. In the current study, we performed long time-scale molecular dynamics simulations on both open and closed states of ADK. Based on the structural analyses of the simulation trajectories, we detected over 20 times conformational transitions between the open and closed states of ADK and identified two novel conformations as intermediate states in the catalytic processes. With these findings, we proposed a possible mechanism for the large-scale domain motion of Escherichia coli ADK and its catalytic process: (1) the substrate free ADK adopted an open conformation; (2) ATP bound with LID domain closure; (3) AMP bound with NMP domain closure; (4) phosphoryl transfer occurred with ATP, and AMP converted into two ADPs, and no conformational transition was detected in the enzyme; (5) LID domain opened with one ADP released; (6) another ADP released with NMP domain open. As both open and closed states sampled a wide range of conformation transitions, our simulation strongly supported the conformational selection mechanism for Escherichia coli ADK.

Exploring the Dynamic Functional Landscape of Adenylate Kinase Modulated by Substrates

Adenylate kinase (ADK) has been explored widely, through both experimental and theoretical studies. However, still less is known about how the functional dynamics of ADK is modulated explicitly by its natural substrates. Here, we report a quantitative study of the dynamic energy landscape for ADK responding to the substrate binding by integrating both experimental investigations and theoretical modeling. We make theoretical predictions which are in remarkable agreement with the single molecule experiments on the substrate-bound complex. With our combined models of ADK in its apo form, in the presence of AMP or ATP, and in complex with both substrates, we specifically address the following key questions: (1) Are there intermediate state(s) during their catalytic cycle and if so how many? (2) How many pathways are there along the open-to-closed transitions and what are their corresponding weights? (3) How do substrates influence the pathway weights and the stability of the intermediates? (4) Which lid's motion is rate-limiting along the turnover cycle, the NMP or the LID domain? Our models predict two major parallel stepwise pathways and two on-pathway intermediates which are denoted as IN (NMP domain open while LID domain closed) and IL (LID domain open and NMP domain closed), respectively. Further investigation of temperature effects suggests that the IN pathway is dominant at room temperature, but the IL pathway is dominant at the optimal temperature. This leads us to propose that the IL pathway is more dominant by entropy and IN pathway by enthalpy. Remarkably, our results show that even with maximum concentrations of natural substrates, ADK still fluctuates between multiple functional states, reflecting an intrinsic capability of large-scale conformational fluctuations which may be essential to its biological function. The results based on the dual-ligands model provide the theoretical validation of random bisubstrate biproducts (Bi-Bi) mechanism for the enzymatic reaction of ADK. Additionally, the pathway flux analysis strongly suggests that the motion of the NMP domain is the rate-determining step for the conformational cycle (opening and closing).

Minimum free energy path of ligand-induced transition in adenylate kinase

Large-scale conformational changes in proteins involve barrier-crossing transitions on the complex free energy surfaces of high-dimensional space. Such rare events cannot be efficiently captured by conventional molecular dynamics simulations. Here we show that, by combining the on-the-fly string method and the multi-state Bennett acceptance ratio (MBAR) method, the free energy profile of a conformational transition pathway in Escherichia coli adenylate kinase can be characterized in a high-dimensional space. The minimum free energy paths of the conformational transitions in adenylate kinase were explored by the on-the-fly string method in 20-dimensional space spanned by the 20 largest-amplitude principal modes, and the free energy and various kinds of average physical quantities along the pathways were successfully evaluated by the MBAR method. The influence of ligand binding on the pathways was characterized in terms of rigid-body motions of the lid-shaped ATP-binding domain (LID) and the AMP-binding (AMPbd) domains. It was found that the LID domain was able to partially close without the ligand, while the closure of the AMPbd domain required the ligand binding. The transition state ensemble of the ligand bound form was identified as those structures characterized by highly specific binding of the ligand to the AMPbd domain, and was validated by unrestrained MD simulations. It was also found that complete closure of the LID domain required the dehydration of solvents around the P-loop. These findings suggest that the interplay of the two different types of domain motion is an essential feature in the conformational transition of the enzyme.

Predicting flexible loop regions that interact with ligands: the challenge of accurate scoring

Flexible loop regions play a critical role in the biological function of many proteins and have been shown to be involved in ligand binding. In the context of structure-based drug design, using or predicting an incorrect loop configuration can be detrimental to the study if the loop is capable of interacting with the ligand. Three protein systems, each with at least one flexible loop region in close proximity to the known binding site, were selected for loop prediction using the CorLps program; a six residue loop region from phosphoribosylglycinamide formyltransferase (GART), two nine residue loop regions from cytochrome P450 (CYP) 119, and an 11 residue loop region from enolase were selected for loop prediction. The results of this study indicate that the statistically based DFIRE scoring function implemented in the CorLps program did not accurately rank native-like predicted loop configurations in any protein system. In an attempt to improve the ranking of the native-like predicted loop configurations, the MM/GBSA and the optimized MM/GBSA-dsr scoring functions were used to re-rank the predicted loops with and without bound ligand. In general, single snapshot MM/GBSA scoring provided the best ranking of native-like loop configurations. Based on the scoring function analyses presented, the optimal ranking of native-like loop configurations is still a difficult challenge and the choice of the "best" scoring function appears to be system dependent.

Computational studies on the conformations of some large-ring cyclodextrins (CDn, n = 20, 21, 22, 23)

Computational studies were carried out on the conformations of large-ring cyclodextrins with degree of polymerization from 20 to 23. Principal component analysis (PCA) was applied for postprocessing of trajectories from conformational search, based on 100.0 ns molecular dynamics simulations. The dominant PCA modes for concerted motions of the macroring atoms were monitored in a lower-dimensions subspace. The first six lowest indexed principal components contribute more than 90% of the total atomic motions in all cases, with about 70% (CD21) to 83% (CD22) contribution coming from the three highest-eigenvalue principal components. Representative average geometries of the cyclodextrin macrorings were also obtained for the whole simulation and for the ten 10.0 ns time intervals of the simulation. We concluded that resemblance exists of the representative conformations of these four cyclodextrins with the circularized three-turn single helical structure proposed for CD21 from small-angle X-ray scattering, as well as with the representative conformations of CD26.

Efficient incorporation of protein flexibility and dynamics into molecular docking simulations

Flexibility and dynamics are protein characteristics that are essential for the process of molecular recognition. Conformational changes in the protein that are coupled to ligand binding are described by the biophysical models of induced fit and conformational selection. Different concepts that incorporate protein flexibility into protein-ligand docking within the context of these two models are reviewed. Several computational studies that discuss the validity and possible limitations of such approaches will be presented. Finally, different approaches that incorporate protein dynamics, e.g., configurational entropy, and solvation effects into docking will be highlighted.