Role of Conformational Dynamics of Sulfotransferases SULT1A1 and SULT1A3 in Substrate Specificity.

Sulfotransferases (SULTs) are phase II metabolizing enzymes catalyzing the sulfoconjugation from the co-factor 3'-Phosphoadenosine 5'-Phosphosulfate (PAPS) to a wide variety of endogenous compounds, drugs and natural products. Although SULT1A1 and SULT1A3 share 93% identity, SULT1A1, the most abundant SULT isoform in humans, exhibits a broad substrate range with specificity for small phenolic compounds, while SULT1A3 displays a high affinity toward monoamine neurotransmitters like dopamine. To elucidate the factors determining the substrate specificity of the SULT1 isoenzymes, we studied the dynamic behavior and structural specificities of SULT1A1 and SULT1A3 by using molecular dynamics (MD) simulations and ensemble docking of common and specific substrates of the two isoforms. Our results demonstrated that while SULT1A1 exhibits a relatively rigid structure by showing lower conformational flexibility except for the lip (loop L1), the loop L2 and the cap (L3) of SULT1A3 are extremely flexible. We identified protein residues strongly involved in the recognition of different substrates for the two isoforms. Our analyses indicated that being more specific and highly flexible, the structure of SULT1A3 has particularities in the binding site, which are crucial for its substrate selectivity.

Small GTPase Ran: Depicting the nucleotide-specific conformational landscape of the functionally important C-terminus.

The small GTPase Ran is the main regulator of the nucleo-cytoplasmic import and export through the nuclear pore complex. It functions as a molecular switch cycling between the GDP-bound inactive and GTP-bound active state. It consists of a globular (G) domain and a C-terminal region, which is bound to the G-domain in the inactive, GDP-bound states. Crystal structures of the GTP-bound active form complexed with Ran binding proteins (RanBP) show that the C-terminus undergoes a large conformational change, embracing Ran binding domains (RanBD). Whereas in the crystal structures of macromolecular complexes not containing RanBDs the structure of the C-terminal segment remains unresolved, indicating its large conformational flexibility. This movement could not have been followed either by experimental or simulation methods. Here, starting from the crystal structure of Ran in both GDP- and GTP-bound forms we show how rigid the C-terminal region in the inactive structure is during molecular dynamics (MD) simulations. Furthermore, we show how MD simulations of the active form are incapable of mapping the open conformations of the C-terminus. By using the MDeNM (Molecular Dynamics with excited Normal Modes) method, we were able to widely map the conformational surface of the C-terminus of Ran in the active GTP-bound form, which allows us to envisage how it can embrace RanBDs.

Sampling of Protein Conformational Space Using Hybrid Simulations: A Critical Assessment of Recent Methods.

Recent years have seen several hybrid simulation methods for exploring the conformational space of proteins and their complexes or assemblies. These methods often combine fast analytical approaches with computationally expensive full atomic molecular dynamics (MD) simulations with the goal of rapidly sampling large and cooperative conformational changes at full atomic resolution. We present here a systematic comparison of the utility and limits of four such hybrid methods that have been introduced in recent years: MD with excited normal modes (MDeNM), collective modes-driven MD (CoMD), and elastic network model (ENM)-based generation, clustering, and relaxation of conformations (ClustENM) as well as its updated version integrated with MD simulations (ClustENMD). We analyzed the predicted conformational spaces using each of these four hybrid methods, applied to four well-studied proteins, triosephosphate isomerase (TIM), 3-phosphoglycerate kinase (PGK), HIV-1 protease (PR) and HIV-1 reverse transcriptase (RT), which provide extensive ensembles of experimental structures for benchmarking and comparing the methods. We show that a rigorous multi-faceted comparison and multiple metrics are necessary to properly assess the differences between conformational ensembles and provide an optimal protocol for achieving good agreement with experimental data. While all four hybrid methods perform well in general, being especially useful as computationally efficient methods that retain atomic resolution, the systematic analysis of the same systems by these four hybrid methods highlights the strengths and limitations of the methods and provides guidance for parameters and protocols to be adopted in future studies.

Insights into the substrate binding mechanism of SULT1A1 through molecular dynamics with excited normal modes simulations.

Sulfotransferases (SULTs) are phase II drug-metabolizing enzymes catalyzing the sulfoconjugation from the co-factor 3'-phosphoadenosine 5'-phosphosulfate (PAPS) to a substrate. It has been previously suggested that a considerable shift of SULT structure caused by PAPS binding could control the capability of SULT to bind large substrates. We employed molecular dynamics (MD) simulations and the recently developed approach of MD with excited normal modes (MDeNM) to elucidate molecular mechanisms guiding the recognition of diverse substrates and inhibitors by SULT1A1. MDeNM allowed exploring an extended conformational space of PAPS-bound SULT1A1, which has not been achieved up to now by using classical MD. The generated ensembles combined with docking of 132 SULT1A1 ligands shed new light on substrate and inhibitor binding mechanisms. Unexpectedly, our simulations and analyses on binding of the substrates estradiol and fulvestrant demonstrated that large conformational changes of the PAPS-bound SULT1A1 could occur independently of the co-factor movements that could be sufficient to accommodate large substrates as fulvestrant. Such structural displacements detected by the MDeNM simulations in the presence of the co-factor suggest that a wider range of drugs could be recognized by PAPS-bound SULT1A1 and highlight the utility of including MDeNM in protein-ligand interactions studies where major rearrangements are expected.

Revealing the activation mechanism of autoinhibited RalF by integrated simulation and experimental approaches.

RalF is an Arf GEF from Legionella pneumophilia, the bacterium that causes severe pneumonia. In its crystal structure, RalF is in the autoinhibited form. A large-scale domain motion is expected to lift the autoinhibition, the mechanism of which is still unknown. Since RalF is activated in the presence of the membrane, its active structure and the structure of the RalF-Arf1 complex could not have been determined experimentally. On the simulation side, it has been proven that classical Molecular Dynamics (MD) alone is not efficient enough to map motions of such amplitude and determine the active conformation of RalF. In this article, using Molecular Dynamics with excited Normal Modes (MDeNM) combined with previous experimental findings we were able to determine the active RalF structure and the structure of the RalF-Arf1 complex in the presence of the membrane, bridging the gap between experiments and simulation.

Nucleotide-Specific Autoinhibition of Full-Length K-Ras4B Identified by Extensive Conformational Sampling.

K-Ras is one of the most frequently mutated oncogenes in human tumor cells. It consists of a well-conserved globular catalytic domain and a flexible tail-like hypervariable region (HVR) at its C-terminal end. It plays a key role in signaling networks in proliferation, differentiation, and survival, undergoing a conformational switch between the active and inactive states. It is regulated through the GDP-GTP cycle of the inactive GDP-bound and active GTP-bound states. Here, without imposing any prior constraints, we mapped the interaction pattern between the catalytic domain and the HVR using Molecular Dynamics with excited Normal Modes (MDeNM) starting from an initially extended HVR conformation for both states. Our sampling captured similar interaction patterns in both GDP- and GTP-bound states with shifted populations depending on the bound nucleotide. In the GDP-bound state, the conformations where the HVR interacts with the effector lobe are more populated than in the GTP-bound state, forming a buried thus autoinhibited catalytic site; in the GTP-bound state conformations where the HVR interacts with the allosteric lobe are more populated, overlapping the α3/α4 dimerization interface. The interaction of the GTP with Switch I and Switch II is stronger than that of the GDP in line with a decrease in the fluctuation upon GTP binding.

An allosteric signaling pathway of human 3-phosphoglycerate kinase from force distribution analysis.

3-Phosphogycerate kinase (PGK) is a two domain enzyme, which transfers a phosphate group between its two substrates, 1,3-bisphosphoglycerate bound to the N-domain and ADP bound to the C-domain. Indispensable for the phosphoryl transfer reaction is a large conformational change from an inactive open to an active closed conformation via a hinge motion that should bring substrates into close proximity. The allosteric pathway resulting in the active closed conformation has only been partially uncovered. Using Molecular Dynamics simulations combined with Force Distribution Analysis (FDA), we describe an allosteric pathway, which connects the substrate binding sites to the interdomain hinge region. Glu192 of alpha-helix 7 and Gly394 of loop L14 act as hinge points, at which these two secondary structure elements straighten, thereby moving the substrate-binding domains towards each other. The long-range allosteric pathway regulating hPGK catalytic activity, which is partially validated and can be further tested by mutagenesis, highlights the virtue of monitoring internal forces to reveal signal propagation, even if only minor conformational distortions, such as helix bending, initiate the large functional rearrangement of the macromolecule.

Importance of aspartate residues in balancing the flexibility and fine-tuning the catalysis of human 3-phosphoglycerate kinase.

The exact role of the metal ion, usually Mg(2+), in the catalysis of human 3-phosphoglycerate kinase, a well-studied two-domain enzyme, has not been clarified. Here we have prepared single and double alanine mutants of the potential metal-binding residues, D374 and D218. While all mutations weaken the catalytic interactions with Mg(2+), they surprisingly strengthen binding of both MgADP and MgATP, and the effects are even more pronounced for ADP and ATP. Thermodynamic parameters of binding indicate an increase in the binding entropy as a reason for the strengthening. In agreement with the experimental results, computer-simulated annealing calculations for the complexes of these mutants have supported the mobility of the nucleotide phosphates and, as a consequence, formation of their new interaction(s) within the active site. A similar type of mobility is suggested to be a characteristic feature of the nucleotide site of the wild-type enzyme, too, both in its inactive open conformation and in the active closed conformation. This mobility of the nucleotide phosphates that is regulated by the aspartate side chains of D218 and D374 through the complexing Mg(2+) is suggested to be essential in enzyme function.

Vibrational softening of a protein on ligand binding.

Neutron scattering experiments have demonstrated that binding of the cancer drug methotrexate softens the low-frequency vibrations of its target protein, dihydrofolate reductase (DHFR). Here, this softening is fully reproduced using atomic detail normal-mode analysis. Decomposition of the vibrational density of states demonstrates that the largest contributions arise from structural elements of DHFR critical to stability and function. Mode-projection analysis reveals an increase of the breathing-like character of the affected vibrational modes consistent with the experimentally observed increased adiabatic compressibility of the protein on complexation.

Ligand chirality effects on the dynamics of human 3-phosphoglycerate kinase: comparison between D- and L-nucleotides.

l-nucleoside analogues are now largely used as antiviral drugs for the treatment of viral infections like HBV, HCV and HIV. However, in order to be fully active, they need to be phosphorylated by cellular or viral kinases. Human 3-phosphogycerate kinase (hPGK) was shown to catalyze the last step of activation of l-enantiomers and thus constitutes an attractive target for theoretical predictions of its phosphorylation efficiency. Molecular dynamics simulations were carried out with four different nucleotides (d-/l-ADP and d-/l-CDP) in complex with hPGK and 1,3-bisphospho-d-glycerate (bPG). The binding affinities of CDPs (both enantiomers) for hPGK were found very weak while d- and l-ADP were better substrates. Interestingly, the binding affinity of the bPG substrate was found to be lower in presence of d-ADP than l-ADP which indicates a potential antagonistic effect on one substrate to the other. A detailed analysis of the simulations unravels important dynamic conditions for efficient phosphorylation. Indeed, as previously described for the natural substrate, the hinge bending motion of the domains upon substrates binding should be more correlated and directional. Interestingly, the unforeseen finding was the larger dynamics freedoms observed for the substrates that was favored by the protein atoms flexibility around the nucleobase binding site.

Substrate binding modifies the hinge bending characteristics of human 3-phosphoglycerate kinase: a molecular dynamics study.

3-Phosphogycerate kinase (PGK) is a two domain enzyme, with a binding site of the 1,3-bisphosphoglycerate on the N-domain and of the ADP on the C-domain. To transfer a phosphate group the enzyme has to undergo a hinge bending motion from open to closed conformation to bring the substrates to close proximity. Molecular dynamics simulation was used to elucidate the effect of ligand binding onto the domain motions of this enzyme. The simulation results of the apo form indicate a hinge bending motion in the ns timescale while the time period of the hinge bending motion of the complex form is clearly over the 20 ns simulation time. The apo form exhibits several hinge points that contribute to the hinge bending motion while upon binding the ligands, the hinge bending becomes strictly restrained with one dominant hinge point in the vicinity of the substrates. At the same time, ligand binding results in an enhanced correlation of internal domain motions.

The influence of interdomain interactions on the intradomain motions in yeast phosphoglycerate kinase: a molecular dynamics study.

A 3-ns molecular dynamics simulation in explicit solvent was performed to examine the inter- and intradomain motions of the two-domain enzyme yeast phosphoglycerate kinase without the presence of substrates. To elucidate contributions from individual domains, simulations were carried out on the complete enzyme as well as on each isolated domain. The enzyme is known to undergo a hinge-bending type of motion as it cycles from an open to a closed conformation to allow the phosphoryl transfer occur. Analysis of the correlation of atomic movements during the simulations confirms hinge bending in the nanosecond timescale: the two domains of the complete enzyme exhibit rigid body motions anticorrelated with respect to each other. The correlation of the intradomain motions of both domains converges, yielding a distinct correlation map in the enzyme. In the isolated domain simulations-in which interdomain interactions cannot occur-the correlation of domain motions no longer converges and shows a very small correlation during the same simulation time. This result points to the importance of interdomain contacts in the overall dynamics of the protein. The secondary structure elements responsible for interdomain contacts are also discussed.

Conformational heterogeneity and low-frequency vibrational modes of proteins.

Molecular dynamics simulation and normal mode analysis are used to calculate the vibrational density of states of dihydrofolate reductase complexed with nicotinamide adenine dinucleotide phosphate at 120 K and the results are compared with the experimental spectrum derived from inelastic neutron scattering. The simulation results indicate that the experimental spectrum arises from an average over proteins trapped in different conformations with structural differences mainly in the loop regions, and that these conformations have significantly different low-frequency (<20 cm(-1)) spectra. Thus, the experimentally measured spectrum is an average over the vibrational modes of different protein conformations and is thus inhomogeneously broadened. The implications of this broadening for future neutron scattering experiments and ligand binding calculations are discussed.

Direct determination of vibrational density of states change on ligand binding to a protein.

The change in the vibrational density of states of a protein (dihydrofolate reductase) on binding a ligand (methotrexate) is determined using inelastic neutron scattering. The vibrations of the complex soften significantly relative to the unbound protein. The resulting free-energy change, which is directly determined by the density of states change, is found to contribute significantly to the binding equilibrium.