The versatile functions of the heat shock protein 70 (Hsp70) family of molecular chaperones rely on allosteric interactions between their nucleotide-binding and substrate-binding domains, NBD and SBD. Understanding the mechanism of interdomain allostery is essential to rational design of Hsp70 modulators. Yet, despite significant progress in recent years, how the two Hsp70 domains regulate each other's activity remains elusive. Covariance data from experiments and computations emerged in recent years as valuable sources of information towards gaining insights into the molecular events that mediate allostery. In the present study, conservation and covariance properties derived from both sequence and structural dynamics data are integrated with results from Perturbation Response Scanning and in vivo functional assays, so as to establish the dynamical basis of interdomain signal transduction in Hsp70s. Our study highlights the critical roles of SBD residues D481 and T417 in mediating the coupled motions of the two domains, as well as that of G506 in enabling the movements of the α-helical lid with respect to the β-sandwich. It also draws attention to the distinctive role of the NBD subdomains: Subdomain IA acts as a key mediator of signal transduction between the ATP- and substrate-binding sites, this function being achieved by a cascade of interactions predominantly involving conserved residues such as V139, D148, R167 and K155. Subdomain IIA, on the other hand, is distinguished by strong coevolutionary signals (with the SBD) exhibited by a series of residues (D211, E217, L219, T383) implicated in DnaJ recognition. The occurrence of coevolving residues at the DnaJ recognition region parallels the behavior recently observed at the nucleotide-exchange-factor recognition region of subdomain IIB. These findings suggest that Hsp70 tends to adapt to co-chaperone recognition and activity via coevolving residues, whereas interdomain allostery, critical to chaperoning, is robustly enabled by conserved interactions.
The Hsp70 family of molecular chaperones assists in protein folding, degradation, assembly/disassembly of some complexes, and intracellular trafficking. These activities in the cell are accomplished by coupled conformational switches/signals between their nucleotide-binding and substrate-binding domains (NBD and SBD), assisted by cognate co-chaperones. Despite significant progress in the field, the molecular basis of Hsp70 machinery and the key interactions that regulate interdomain communication are not fully understood. Using a combination of experimental and computational methods, including in vivo functional assays, sequence- and structure-based analyses and perturbation response scanning, we identified a network of conserved interactions in subdomain IA of the NBD, which plays a key (effector) role in propagating signals between the ATP-binding and substrate-binding sites. Subdomain IIA, on the other hand, appears to adapt to J-domain co-chaperone binding by virtue of a broadly distributed cluster of co-evolving residues on the recognition surface.
Biomolecular conformational transitions are essential to biological functions. Most experimental methods report on the long-lived functional states of biomolecules, but information about the transition pathways between these stable states is generally scarce. Such transitions involve short-lived conformational states that are difficult to detect experimentally. For this reason, computational methods are needed to produce plausible hypothetical transition pathways that can then be probed experimentally. Here we propose a simple and computationally efficient method, called ANMPathway, for constructing a physically reasonable pathway between two endpoints of a conformational transition. We adopt a coarse-grained representation of the protein and construct a two-state potential by combining two elastic network models (ENMs) representative of the experimental structures resolved for the endpoints. The two-state potential has a cusp hypersurface in the configuration space where the energies from both the ENMs are equal. We first search for the minimum energy structure on the cusp hypersurface and then treat it as the transition state. The continuous pathway is subsequently constructed by following the steepest descent energy minimization trajectories starting from the transition state on each side of the cusp hypersurface. Application to several systems of broad biological interest such as adenylate kinase, ATP-driven calcium pump SERCA, leucine transporter and glutamate transporter shows that ANMPathway yields results in good agreement with those from other similar methods and with data obtained from all-atom molecular dynamics simulations, in support of the utility of this simple and efficient approach. Notably the method provides experimentally testable predictions, including the formation of non-native contacts during the transition which we were able to detect in two of the systems we studied. An open-access web server has been created to deliver ANMPathway results.
Many biomolecules are like tiny molecular machines that need to change their shapes and visit many states to perform their biological functions. For a complete molecular understanding of a biological process, one needs to have information on the relevant stable states of the system in question, as well as the pathways by which the system travels from one state to another. We report here an efficient computational method that uses the knowledge of experimental structures of a pair of stable states in order to construct an energetically favoravle pathway between them. We adopt a simple representation of the molecular system by replacing the atoms with beads connected by springs and constructing an energy function with two minima around the end-states. We searched for the structure with highest energy that the system is most likely to visit during the transition and created two paths starting from this structure and proceeding toward the end-states. The combined result of these two paths is the minimum energy pathway between the two stable states. We apply this method to study important structural changes in one enzyme and three large proteins that transport small molecules and ions across the cell membrane.
Recognition of injured mitochondria for degradation by macroautophagy is essential for cellular health, but the mechanisms remain poorly understood. Cardiolipin is an inner mitochondrial membrane phospholipid. We found that rotenone, staurosporine, 6-hydroxydopamine and other pro-mitophagy stimuli caused externalization of cardiolipin to the mitochondrial surface in primary cortical neurons and SH-SY5Y cells. RNAi knockdown of cardiolipin synthase or of phospholipid scramblase-3, which transports cardiolipin to the outer mitochondrial membrane, decreased mitochondrial delivery to autophagosomes. Furthermore, we found that the autophagy protein microtubule-associated-protein-1-light chain-3 (LC3), which mediates both autophagosome formation and cargo recognition, contains cardiolipin-binding sites important for the engulfment of mitochondria by the autophagic system. Mutation of LC3 residues predicted as cardiolipin-interaction sites by computational modeling inhibited its participation in mitophagy. These data indicate that redistribution of cardiolipin serves as an “eat-me” signal for the elimination of damaged mitochondria from neuronal cells.
Recent studies suggest that protein motions observed in molecular simulations are related to biochemical activities, although the computed time scales do not necessarily match those of the experimentally observed processes. The molecular origin of this conflicting observation is explored here for a test protein, cyanovirin-N (CV-N), through a series of molecular dynamics simulations that span a time range of three orders of magnitude up to 0.4 microseconds. Strikingly, increasing the simulation time leads to an approximately uniform amplification of the motional sizes, while maintaining the same conformational mechanics. Residue fluctuations exhibit amplitudes of 1–2 Å in the nanosecond simulations, while their average sizes increase by a factor of 4–5 in the microsecond regime. The mean-square displacements averaged over all residues (y) exhibit a power law dependence of the form y ∝ x0.26 on the simulation time (x). Essential dynamics analysis of the trajectories, on the other hand, demonstrates that CV-N has robust preferences to undergo specific types of motions that already can be detected at short simulation times, provided that multiple runs are performed and carefully analyzed.
structure-encoded dynamics; molecular dynamics simulations; power law; global motions; equilibrium fluctuations of cyanovirin-N
Domain swapping creates protein oligomers by exchange of structural units between identical monomers. At present, no unifying molecular mechanism of domain swapping has emerged. Here we used the protein Cyanovirin-N and 19F-NMR to investigate the process of domain swapping. CV-N is an HIV inactivating protein that can exist as a monomer or a domain-swapped dimer. We measured thermodynamic and kinetic parameters of the conversion process and determined the size of the energy barrier between the two species. The barrier is very large and of similar magnitude to that for equilibrium unfolding of the protein. Therefore, for CV-N, overall unfolding of the polypeptide is required for domain swapping.
Background: Mechanisms of sodium-coupled substrate binding/release are largely unknown in the family of glutamate transporters.
Results: Opening of helical hairpin HP2 leads to spontaneous substrate release from an aspartate transporter in the inward-facing state.
Conclusion: Substrate dissociation is preceded by the release of a sodium ion.
Significance: Exploring the dynamics of substrate gating is crucial for understanding the mechanisms of the full transport cycle.
Sodium-coupled neurotransmitter transporters play a key role in neuronal signaling by clearing excess transmitter from the synapse. Structural data on a trimeric archaeal aspartate transporter, GltPh, have provided valuable insights into structural features of human excitatory amino acid transporters. However, the time-resolved mechanisms of substrate binding and release, as well as that of coupling to sodium co-transport, remain largely unknown for this important family. We present here the results of the most extensive simulations performed to date for GltPh in both outward-facing and inward-facing states by taking advantage of significant advances made in recent years in molecular simulation technology. The generated multiple microsecond trajectories consistently show that the helical hairpin HP2, not HP1, serves as an intracellular gate (in addition to its extracellular gating role). In contrast to previous proposals, HP1 can neither initiate nor accommodate neurotransmitter release without prior opening of HP2 by at least 4.0 Å. Aspartate release invariably follows that of a sodium ion located near the HP2 gate entrance. Asp-394 on TM8 and Arg-276 on HP1 emerge as key residues that promote the reorientation and diffusion of substrate toward the cell interior. These findings underscore the significance of examining structural dynamics, as opposed to static structure(s), to make inferences on the mechanisms of transport and key interactions.
Aspartate; Membrane Transport; Neurotransmitter Transport; Sodium Transport; Transporters; Inward-facing; Sodium-coupled; Substrate Release
Ionotropic glutamate receptors (iGluRs) harbor two extracellular domains: the membrane-proximal ligand-binding domain (LBD) and the distal N-terminal domain (NTD). These are involved in signal sensing: the LBD binds L-glutamate, which activates the receptor channel. Ligand binding to the NTD modulates channel function in the NMDA receptor subfamily of iGluRs, which has not been observed for the AMPAR subfamily to date. Structural data suggest that AMPAR NTDs are packed into tight dimers and have lost their signaling potential. Here, we assess NTD dynamics from both subfamilies, using a variety of computational tools. We describe the conformational motions that underly NMDAR NTD allosteric signaling. Unexpectedly, AMPAR NTDs are capable of undergoing similar dynamics; although dimerization imposes restrictions, the two subfamilies sample similar, interconvertible conformational subspaces. Finally, we solve the crystal structure of AMPAR GluA4 NTD, and combined with molecular dynamics simulations, we characterize regions pivotal for an as-yet-unexplored dynamic spectrum of AMPAR NTDs.
► Crystal structure newly resolved for AMPAR GluA4 N-terminal domain (NTD) ► Dynamics underlying NMDA NR2B NTD allostery ► Similar global dynamics between AMPA- and NMDA-receptor NTDs ► Critical role of dimer interface contacts for GluA3 NTD allostery
Dutta et al. examine the dynamics and allostery of ionotropic glutamate receptors N-terminal domains (NTDs). Systematic comparison of the members of the AMPAR subfamily, and across the AMPAR and NMDAR subfamilies, reveals mechanisms of global and local motions favored by the NTDs and allosteric potential in AMPARs.
Catalytic loop motions facilitate substrate recognition and binding in many enzymes. While these motions appear to be highly flexible, their functional significance suggests that structure-encoded preferences may play a role in selecting particular mechanisms of motions. We performed an extensive study on a set of enzymes to assess whether the collective/global dynamics, as predicted by elastic network models (ENMs), facilitates or even defines the local motions undergone by functional loops. Our dataset includes a total of 117 crystal structures for ten enzymes of different sizes and oligomerization states. Each enzyme contains a specific functional/catalytic loop (10–21 residues long) that closes over the active site during catalysis. Principal component analysis (PCA) of the available crystal structures (including apo and ligand-bound forms) for each enzyme revealed the dominant conformational changes taking place in these loops upon substrate binding. These experimentally observed loop reconfigurations are shown to be predominantly driven by energetically favored modes of motion intrinsically accessible to the enzyme in the absence of its substrate. The analysis suggests that robust global modes cooperatively defined by the overall enzyme architecture also entail local components that assist in suitable opening/closure of the catalytic loop over the active site.
Protein loops have critical roles in ligand binding and catalysis. An unresolved issue in this context is the extent to which the intrinsic dynamics of proteins predispose loops to perform their molecular function. In this work, we (i) critically examine the structural changes undergone by functional/catalytic loops based on a set of enzyme crystal structures in the presence/absence of a ligand, and (ii) examine to what extent those motions are correlated with, or driven by, the global modes that are predictable using simplified, physics-based models. Using a dataset of 117 structures for ten enzymes of different sizes and oligomerization states, we show that the collective modes defined by the protein topology favor loop rearrangements in reasonable agreement with those experimentally observed upon activation. These results suggest that simple but robust motions encoded by the entire architecture, not the local binding site only, assist in binding of the ligand, positioning of the catalytic loop, and/or sequestration of the catalytic site, which in turn, enable efficient catalysis.
Background: Crystal structures of an inward-facing CaiT conformation with four substrate sites are reported.
Results: This conformation is substrate-releasing. Arg262 facilitates release from the primary site.
Conclusion: Unbinding pathway from the primary site is determined by a broken helical portion of a transmembrane domain.
Significance: Mechanism of substrate release may be relevant to structural homologues such as neurotransmitter transporters.
CaiT is a homotrimeric antiporter that exchanges l-carnitine (CRN) with γ-butyrobetaine (GBB) across the bacterial membrane. Three structures have been resolved to date for CaiT, all in the inward-facing state: CRN-bound (with four CRNs per subunit), GBB-bound (two GBBs per subunit), and apo. One of the reported binding sites is the counterpart of the primary site observed in structurally similar transporters. However, the mechanism and pathway(s) of CRN/GBB unbinding and translocation, or even the ability of the substrates to dislodge from the reported binding sites, are yet to be determined. To shed light on these issues, we performed a total of 1.3 μs of molecular dynamics simulations and examined the dynamics of substrate-bound CaiT structures under different conditions. We find that both CRN and GBB are able to dissociate completely from their primary site into the cytoplasm. Substrate molecules initially located at the secondary sites dissociate even faster (within tens of nanoseconds) into the extra- or intracellular regions. Interestingly, the unbinding pathway from the primary site appears to be dictated by the geometry of the unwound part of the transmembrane (TM) helix 3, mostly around Thr100 therein. Arg262 on TM7, which apparently mimics the role of Na+ in CaiT structural homologues, plays a key role in triggering the dissociation of the substrate away from the primary site and guiding its release to the cytoplasm provided that the unwound part of TM3 switches from a shielding to a yielding pose.
Biophysics; Carnitine; Membrane Proteins; Molecular Dynamics; Transport
Biochemical activity and core stability are essential properties of proteins, maintained usually by conserved amino acids. Structural dynamics emerged in recent years as another essential aspect of protein functionality. Structural dynamics enable the adaptation of the protein to binding substrates and to undergo allosteric transitions, while maintaining the native fold. Key residues that mediate structural dynamics would thus be expected to be conserved or exhibit coevolutionary patterns at least. Yet, the correlation between sequence evolution and structural dynamics is yet to be established. With recent advances in efficient characterization of structural dynamics, we are now in a position to perform a systematic analysis. In the present study, a set of 34 enzymes representing various folds and functional classes is analyzed using information theory and elastic network models. Our analysis shows that the structural regions distinguished by their coevolution propensity as well as high mobility are predisposed to serve as substrate recognition sites, whereas residues acting as global hinges during collective dynamics are often supported by conserved residues. We propose a mobility scale for different types of amino acids, which tends to vary inversely with amino acid conservation. Our findings suggest the balance between physical adaptability (enabled by structure-encoded motions) and chemical specificity (conferred by correlated amino acid substitutions) underlies the selection of a relatively small set of versatile folds by proteins.
protein dynamics; sequence evolution; Gaussian network model; specificity and adaptability
PUMA (p53 upregulated modulator of apoptosis) is a Bcl-2 homology 3 (BH3)-only Bcl-2 family member and a key mediator of apoptosis induced by a wide variety of stimuli. PUMA is particularly important in initiating radiation-induced apoptosis and damage in the gastrointestinal and hematopoietic systems. Unlike most BH3-only proteins, PUMA neutralizes all five known antiapoptotic Bcl-2 members through high affinity interactions with its BH3 domain to initiate mitochondria-dependent cell death. Using structural data on the conserved interactions of PUMA with Bcl-2-like proteins, we developed a pharmacophore model that mimics these interactions. In silico screening of the ZINC 8.0 database with this pharmacophore model yielded 142 compounds that could potentially disrupt these interactions. Thirteen structurally diverse compounds with favorable in silico ADME/Toxicity profiles have been retrieved from this set. Extensive testing of these compounds using cell-based and cell-free systems identified lead compounds that confer considerable protection against PUMA-dependent and radiation-induced apoptosis, and inhibit the interaction between PUMA and Bcl-xL.
Inhibition of PUMA-induced apoptosis; Bcl-2 protein family; BH3 domain; protein-protein interactions; pharmacophore modeling; druggability; virtual screening of libraries of small compounds
Glutamate transporters regulate excitatory amino acid neurotransmission across neuronal and glial cell membranes by coupling the translocation of their substrate (aspartate or glutamate) into the intracellular (IC) medium to the energetically favorable transport of sodium ions or other cations. The first crystallographically resolved structure of this family, the archaeal aspartate transporter, GltPh, has served as a structural paradigm for elucidating the mechanism of substrate translocation by these transporters. Two helical hairpins, HP2 and HP1, at the core domains of the three subunits that form this membrane protein have been proposed to act as the respective extracellular and IC gates for substrate intake and release. Molecular dynamics simulations using the outward-facing structure have confirmed that the HP2 loop acts as an EC gate. The mechanism of substrate release at atomic scale, however, remained unknown due to the lack of structural data until the recent determination of the inward-facing structure of GltPh. In the present study, we use this recently resolved structure to simulate the release of substrate to the cytoplasm and the roles of HP1 and HP2 in this process. The highly flexible HP2 loop is observed to serve as an activator (or initiator) prompting the release of a gatekeeper Na+ to the cytoplasm and promoting the influx of water molecules from the cytoplasm, which effectively disrupt substrate–protein interactions and drive the dislodging of the substrate from its binding site. The completion of substrate release and exit, however, entails the opening of the highly stable HP1 loop as well. Overall, the unique conformational flexibility of the HP2 loop, the dissociation of a Na+, the hydration of binding pocket, and final yielding of the HP1 loop 3-Ser motif emerge as the successive events controlling the release of the bound substrate to the cell interior by glutamate transporters.
Oligomerization is a functional requirement for many proteins. The interfacial interactions and the overall packing geometry of the individual monomers are viewed as important determinants of the thermodynamic stability and allosteric regulation of oligomers. The present study focuses on the role of the interfacial interactions and overall contact topology in the dynamic features acquired in the oligomeric state. To this aim, the collective dynamics of enzymes belonging to the amino acid kinase family both in dimeric and hexameric forms are examined by means of an elastic network model, and the softest collective motions (i.e., lowest frequency or global modes of motions) favored by the overall architecture are analyzed. Notably, the lowest-frequency modes accessible to the individual subunits in the absence of multimerization are conserved to a large extent in the oligomer, suggesting that the oligomer takes advantage of the intrinsic dynamics of the individual monomers. At the same time, oligomerization stiffens the interfacial regions of the monomers and confers new cooperative modes that exploit the rigid-body translational and rotational degrees of freedom of the intact monomers. The present study sheds light on the mechanism of cooperative inhibition of hexameric N-acetyl-L-glutamate kinase by arginine and on the allosteric regulation of UMP kinases. It also highlights the significance of the particular quaternary design in selectively determining the oligomer dynamics congruent with required ligand-binding and allosteric activities.
Protein function requires a three-dimensional structure with specific dynamic features for catalytic and binding events, and, in many cases, the structure results from the assembly of more than one polypeptide chain (also called monomer or subunit) to form an oligomer or multimer. Proteins such as hemoglobin or chaperonin GroEL are oligomers formed by 2 and 14 subunits, respectively, whereas virus capsids are multimers composed of hundreds of monomers. In these cases, the architecture of the interface between the subunits and the overall assembly geometry are essential in determining the functional motions that these sophisticated structures are able to perform under physiological conditions. Here we present results from our computational study of the large-amplitude motions of dimeric and hexameric proteins that belong to the Amino Acid Kinase family. Our study reveals that the monomers in these oligomeric proteins are arranged in such a way that the oligomer inherits the intrinsic dynamic features of its components. The packing geometry additionally confers the ability to perform highly cooperative conformational changes that involve all monomers and enable the biological activity of the multimer. The study highlights the significance of the quaternary design in favoring the oligomer dynamics that enables ligand-binding and allosteric regulation functions.
Many proteins require bound metals to achieve their function. We take advantage of increasing structural data on metal-binding proteins to elucidate three properties: the involvement of metal-binding sites in the global dynamics of the protein, predicted by elastic network models, their exposure/burial to solvent, and their signal-processing properties indicated by Markovian stochastics analysis. Systematic analysis of a dataset of 145 structures reveals that the residues that coordinate metal ions enjoy remarkably efficient and precise signal transduction properties. These properties are rationalized in terms of their physical properties: participation in hinge sites that control the softest modes collectively accessible to the protein and occupancy of central positions minimally exposed to solvent. Our observations suggest that metal-binding sites may have been evolutionary selected to achieve optimum allosteric communication. They also provide insights into basic principles for designing metal-binding sites, which are verified to be met by recently designed de novo metal-binding proteins.
metal-binding sites; elastic network models; GNM; signal propagation; allosteric communication
Summary: We developed a Python package, ProDy, for structure-based analysis of protein dynamics. ProDy allows for quantitative characterization of structural variations in heterogeneous datasets of structures experimentally resolved for a given biomolecular system, and for comparison of these variations with the theoretically predicted equilibrium dynamics. Datasets include structural ensembles for a given family or subfamily of proteins, their mutants and sequence homologues, in the presence/absence of their substrates, ligands or inhibitors. Numerous helper functions enable comparative analysis of experimental and theoretical data, and visualization of the principal changes in conformations that are accessible in different functional states. ProDy application programming interface (API) has been designed so that users can easily extend the software and implement new methods.
Availability: ProDy is open source and freely available under GNU General Public License from http://www.csb.pitt.edu/ProDy/.
Contact: firstname.lastname@example.org; email@example.com
A detailed analysis of high-resolution structural data and computationally predicted dynamics was carried out for a designed sugar binding protein. The mean-square-deviations in the positions of residues derived from Nuclear Magnetic Resonance (NMR) models, and those inferred from X-ray crystallographic B-factors for two different crystal forms were compared with the predictions based on the Gaussian Network Model (GNM), and the results from molecular dynamics (MD) simulations. GNM systematically yielded a higher correlation than MD, with experimental data, suggesting that the lack of atomistic details in the coarse-grained GNM is more than compensated for by the mathematically exact evaluation of fluctuations using the native contacts topology. Evidence is provided that particular loop motions are curtailed by intermolecular contacts in the crystal environment causing a discrepancy between theory and experiments. Interestingly, the information conveyed by X-ray crystallography becomes more consistent with NMR models and computational predictions when ensembles of X-ray models are considered. Less precise (broadly distributed) ensembles indeed appear to describe the accessible conformational space under native state conditions better than B-factors. Our results highlight the importance of utilizing multiple conformations obtained by alternative experimental methods, and analyzing results from both coarse-grained models and atomic simulations, for accurate assessment of motions accessible to proteins under native state conditions.
equilibrium dynamics; ensemble of conformations; inter-residue contact topology; crystal contacts; elastic network model; sugar-binding protein
Catalysis of ADP-ATP exchange by nucleotide exchange factors (NEFs) is central to the activity of Hsp70 molecular chaperones. Yet, the mechanism of interaction of this family of chaperones with NEFs is not well understood in the context of the sequence evolution and structural dynamics of Hsp70 ATPase domains. We studied the interactions of Hsp70 ATPase domains with four different NEFs on the basis of the evolutionary trace and co-evolution of the ATPase domain sequence, combined with elastic network modeling of the collective dynamics of the complexes. Our study reveals a subtle balance between the intrinsic (to the ATPase domain) and specific (to interactions with NEFs) mechanisms shared by the four complexes. Two classes of key residues are distinguished in the Hsp70 ATPase domain: (i) highly conserved residues, involved in nucleotide binding, which mediate, via a global hinge-bending, the ATPase domain opening irrespective of NEF binding, and (ii) not-conserved but co-evolved and highly mobile residues, engaged in specific interactions with NEFs (e.g., N57, R258, R262, E283, D285). The observed interplay between these respective intrinsic (pre-existing, structure-encoded) and specific (co-evolved, sequence-dependent) interactions provides us with insights into the allosteric dynamics and functional evolution of the modular Hsp70 ATPase domain.
The heat shock protein 70 (Hsp70) serves as a housekeeper in the cell, assisting in the correct folding, trafficking, and degradation of many proteins. The ATPase domain is the control unit of this molecular machine and its efficient functioning requires interactions with co-chaperones, including, in particular, the nucleotide exchange factors (NEFs). We examined the molecular motions of the ATPase domain in both NEF-bound and -unbound forms. We found that the NEF-binding surface enjoys large global movements prior to NEF binding, which presumably facilitates NEF recognition and binding. NEF binding stabilizes the ATPase domain in an open form and thereby facilitates the nucleotide exchange step of the chaperone cycle. A series of highly correlated amino acids were distinguished at the NEF-binding sites of the Hsp70 ATPase domain, which highlights the adaptability of the ATPase domain, both structurally and sequentially, to recognize NEFs. In contrast, the nucleotide-binding residues are tightly held near a global hinge center and are highly conserved. The contrasting properties of these two groups of residues point to an evolutionarily optimized balance between conserved/constrained and co-evolved/mobile amino acids, which enables the functional interactions of the modular Hps70 ATPase domains with NEFs.
The recent elucidation of the sodium/galactose symporter structure from the Vibrio parahaemolyticus bacterium, vSGLT, has revealed a similarity in the core architecture with transporters from different gene families, including the leucine transporter (LeuT). Even though several transporters sharing this core have been structurally determined over the past few years, vSGLT is the only one crystallized in the substrate-bound inward-facing conformation so far. In this study, we report the first insight into the dynamics and coordination of the galactose (Gal) and proposed Na+ ion in vSGLT using a series of molecular dynamics simulations with a total time of about 0.1 μs. Our study reveals new residues, not observed in the crystal structure, which closely interact with the Na+ ion or the substrate for extended times, and shows that in the crystallized conformation, a Na+ ion placed at the site equivalent to Na2 in LeuT can escape into the intracellular (IC) space in the absence of external forces. We have identified the highly conserved Asp189 as a likely binding residue on the pathway of Na+ into the IC cavity. The release of Gal, on the other hand, requires the rotation of the side chain of the inner hydrophobic gate, Tyr263, without a significant change in vSGLT backbone conformation. Our simulations further show that the crystal structure represents but one accessible binding pose of Gal and Na+ among an ensemble of microstates, and that the Gal undergoes versatile alternate interactions at the binding pocket.
The ATPase cycle of the heat shock protein 70 (HSP70) is largely dependent on the ability of its nucleotide binding domain (NBD), also called ATPase domain, to undergo structural changes between its open and closed conformations. We present here a combined study of the Hsp70 NBD sequence, structure and dynamic features to identify the residues that play a crucial role in mediating the allosteric signaling properties of the ATPase domain. Specifically, we identify the residues involved in the shortest-path communications of the domain modeled as a network of nodes (residues) and links (equilibrium interactions). By comparing the calculations on both closed and open conformation of Hsp70 NBD, we identified a subset of central residues located at the interface between the two lobes of the NBD near the nucleotide binding site, which form a putative communication pathway invariant to structural changes. Two pairs of residues forming contacts at the interface in the closed conformation of the NBD are observed to no longer interact in the open conformation, suggesting that these specific interactions may play a switch role in establishing the transition of the NBD between the two functional forms. Sequence co-evolution analysis and collective dynamics analysis with elastic network model further confirm the key roles of these residues in Hsp70 NBD dynamics and functions.
Biomolecular systems possess unique, structure-encoded dynamic properties that underlie their biological functions. Recent studies indicate that these dynamic properties are determined to a large extent by the topology of native contacts. In recent years, elastic network models used in conjunction with normal mode analyses have proven to be useful for elucidating the collective dynamics intrinsically accessible under native state conditions, including in particular the global modes of motions that are robustly defined by the overall architecture. With increasing availability of structural data for well-studied proteins in different forms (liganded, complexed, or free), there is increasing evidence in support of the correspondence between functional changes in structures observed in experiments and the global motions predicted by these coarse-grained analyses. These observed correlations suggest that computational methods may be advantageously employed for assessing functional changes in structure and allosteric mechanisms intrinsically favored by the native fold.
elastic network models; normal modes; principal component analysis; collective motions; allosteric changes in conformation; closed/open conformations
Comparison of elastic network model predictions with experimental data has provided important insights on the dominant role of the network of inter-residue contacts in defining the global dynamics of proteins. Most of these studies have focused on interpreting the mean-square fluctuations of residues, or deriving the most collective, or softest, modes of motions that are known to be insensitive to structural and energetic details. However, with increasing structural data, we are in a position to perform a more critical assessment of the structure-dynamics relations in proteins, and gain a deeper understanding of the major determinants of not only the mean-square fluctuations and lowest frequency modes, but the covariance or the cross-correlations between residue fluctuations and the shapes of higher modes. A systematic study of a large set of NMR-determined proteins is analyzed using a novel method based on entropy maximization to demonstrate that the next level of refinement in the elastic network model description of proteins ought to take into consideration properties such as contact order (or sequential separation between contacting residues) and the secondary structure types of the interacting residues, whereas the types of amino acids do not play a critical role. Most importantly, an optimal description of observed cross-correlations requires the inclusion of destabilizing, as opposed to exclusively stabilizing, interactions, stipulating the functional significance of local frustration in imparting native-like dynamics. This study provides us with a deeper understanding of the structural basis of experimentally observed behavior, and opens the way to the development of more accurate models for exploring protein dynamics.
As more protein structures are solved, we are able to perform a more critical assessment of the relationship between protein structure and dynamics, and to gain a deeper understanding of the major determinants of structural dynamics. Here we perform a systematic study on a set of proteins structurally determined by NMR spectroscopy. The dynamics are analyzed using elastic network models and a novel method based on entropy maximization to demonstrate that properties such as contact order and secondary structure do play a role in defining the experimentally observed covariance data. Most importantly, an optimal description of observed cross-correlations requires the inclusion of destabilizing, as well as stabilizing, interactions, stipulating the functional significance of local frustration in imparting native-like dynamics.
N-Acetyl-L-Glutamate Kinase (NAGK) is the structural paradigm for examining the catalytic mechanisms and dynamics of amino acid kinase family members. Given that the slow conformational dynamics of the NAGK (at the microseconds time scale or slower) may be rate-limiting, it is of importance to assess the mechanisms of the most cooperative modes of motion intrinsically accessible to this enzyme. Here, we present the results from normal mode analysis using an elastic network model representation, which shows that the conformational mechanisms for substrate binding by NAGK strongly correlate with the intrinsic dynamics of the enzyme in the unbound form. We further analyzed the potential mechanisms of allosteric signalling within NAGK using a Markov model for network communication. Comparative analysis of the dynamics of family members strongly suggests that the low-frequency modes of motion and the associated intramolecular couplings that establish signal transduction are highly conserved among family members, in support of the paradigm sequence→structure→dynamics→function.
During the last 20 years both the experimental and computational communities have provided strong evidence that proteins cannot be regarded as static entities, but as intrinsically flexible molecules that exploit their fluctuation dynamics for catalytic and ligand-binding events, as well as for allosteric regulation. This intrinsic dynamics is encoded in the protein structure and, therefore, those proteins with similar folding should share dynamic features essential to their biological function. In this work, we have applied an Elastic Network Model to predict the large-amplitude dynamics of different enzymes belonging to the same protein family (Amino Acid Kinase family). Subsequent comparison of the dynamics of these proteins reveals that this protein family follows the same dynamic pattern. The present results are strongly supported by experimental data and provide new insights into the performance of biological function by these enzymes. The investigation presented here provides us with a useful framework to identify dynamic fingerprints among proteins with structural similarities.
A 3-dimensional pharmacophore model was generated utilizing a set of known inhibitors of c-Myc-Max heterodimer formation. The model successfully identified a set of structurally diverse compounds with potential inhibitory activity against c-Myc. Nine compounds were tested in vitro, and four displayed affinities in the µM range and growth inhibitory activity against c-Myc-overexpressing cells. These studies demonstrate the applicability of pharmacophore modeling to the identification of novel and potentially more puissant inhibitors of the c-Myc oncoprotein.