Inherent flexibility and protein function: The open/closed conformational transition in the N-terminal domain of calmodulin

The key to understand a protein's function often lies in its conformational dynamics. We develop a coarse-grained variational model to investigate the interplay between structural transitions, conformational flexibility, and function of the N-terminal calmodulin domain (nCaM). In this model, two energy basins corresponding to the "closed" apo conformation and "open" holo conformation of nCaM are coupled by a uniform interpolation parameter. The resulting detailed transition route from our model is largely consistent with the recently proposed EFbeta-scaffold mechanism in EF-hand family proteins. We find that the N-terminal parts of the calcium binding loops shows higher flexibility than the C-terminal parts which form this EFbeta-scaffold structure. The structural transition of binding loops I and II are compared in detail. Our model predicts that binding loop II, with higher flexibility and earlier structural change than binding loop I, dominates the open/closed conformational transition in nCaM.     

Using steered molecular dynamics to study the interaction between ADP and the nucleotide-binding domain of yeast Hsp70 protein Ssa1

Genetics experiments have identified six mutations located in the subdomain IA (A17V, R23H, G32D, G32S, R34K, V372I) of Ssa1 that influence propagation of the yeast [PSI⁺] prion. However, the underlining molecular mechanisms of these mutations are still unclear. The six mutation sites are present in the IA subdomain of the nucleotide-binding domain (NBD). The ATPase subdomain IA is a critical mediator of inter-domain allostery in Hsp70 molecular chaperones, so the mutation and changes in this subdomain may influence the function of the substrate-binding domain. In addition, ADP release is a rate-limiting step of the ATPase cycle and dysregulation of the ATPase cycle influences the propagation of the yeast [PSI⁺] prion. In this work, steered molecular dynamics (SMD) simulations were performed to explore the interaction between ADP and NBD. Results suggest that during the SMD simulations, hydrophobic interactions are predominant and variations in the binding state of ADP within the mutants is a potential reason for in vivo effects on yeast [PSI⁺] prion propagation. Additionally, we identify the primary residues in the ATPase domain that directly constitute the main hydrophobic interaction network and directly influence the ADP interaction state with the NBD of Ssa1. Furthermore, this in silico analysis reaffirms the importance of previously experimentally-determined residues in the Hsp70 ATPase domain involved in ADP binding and also identifies new residues potentially involved in this process.     

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

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

Inhibition of Hsp90 with synthetic macrolactones: synthesis and structural and biological evaluation of ring and conformational analogs of radicicol

A series of benzo-macrolactones of varying ring size and conformation has been prepared by chemical synthesis and evaluated by structural and biological techniques. Thus, 12- to 16-membered lactones were obtained by concise routes, involving ring-closing metathesis as a key step. In enzyme assays, the 13-, 15-, and 16-membered analogs are good inhibitors, suggesting that they can adopt the required conformation to fit in the ATP-binding site. This was confirmed by cocrystallization of 13-, 14-, and 15-membered lactones with the N-terminal domain of yeast Hsp90, showing that they bind similarly to the "natural" 14-membered radicicol. The most active compounds in the ATPase assays also showed the greatest growth-inhibitory potency in HCT116 human colon cancer cells and the established molecular signature of Hsp90 inhibition, i.e., depletion of client proteins with upregulation of Hsp70.     

Domain motion of individual F1-ATPase β-subunits during unbiased molecular dynamics simulations

F(1)-ATPase is the catalytic domain of F(1)F(o)-ATP synthase and consists of a hexameric arrangement of three noncatalytic α and three catalytic β subunits. We have used unbiased molecular dynamics simulations with a total simulation time of 900 ns to investigate the dynamic relaxation properties of isolated β-subunits as a step toward explaining the function of the integral F(1) unit. To this end, we simulated the open (β(E)) and the closed (β(TP)) conformations under unbiased conditions for up to 120 ns each using several samples. The simulations confirm that nucleotide-free β(E) retains its open configuration over the course of the simulations. The same is true when the neighboring α subunits are included. The nucleotide-depleted as well as the nucleotide-bound isolated β(TP) subunits show a significant trend toward the open conformation during our simulations, with one trajectory per case opening completely. Hence, our simulations suggest that the equilibrium conformation of a nucleotide-free β-subunit is the open conformation and that the transition from the closed to the open conformation can occur on a time scale of a few tens of nanoseconds.     

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

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

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

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

Quantum mechanical modeling of the structures, energetics and spectral properties of Iα and Iβ cellulose

Periodic planewave and molecular cluster density functional theory (DFT) calculations were performed on Iα and Iβ cellulose in four different conformations each. The results are consistent with the previous interpretation of experimental X-ray and neutron diffraction data that both Iα and Iβ cellulose are dominantly found in the tg conformation of the hydroxymethyl group with a H-bonding conformation termed “Network A”. Structural and energetic results of the periodic DFT calculations with dispersion corrections (DFT-D2) are consistent with observation suggesting that this methodology is accurate to within a few percent for modeling cellulose. The structural and energetic results were confirmed by comparison of calculated vibrational frequencies against observed infrared and Raman frequencies of Iα and Iβ cellulose. Structures extracted from the periodic DFT-D2 energy minimizations were used to calculate the ¹³C nuclear magnetic resonance chemical shifts (δ¹³C), and the tg/Network A conformations of both Iα and Iβ cellulose produced excellent correlations with observed δ¹³C values.     

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

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

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

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

Chaperone-assisted protein folding in the cell cytoplasm

Folding of polypeptides in the cell typically requires the assistance of a set of proteins termed molecular chaperones. Chaperones are an essential group of proteins necessary for cell viability under both normal and stress conditions. There are several chaperone systems which carry out a multitude of functions all aimed towards insuring the proper folding of target proteins. Chaperones can assist in the efficient folding of newly-translated proteins as these proteins are being synthesized on the ribosome and can maintain pre-existing proteins in a stable conformation. Chaperones can also promote the disaggregation of preformed protein aggregates. Many of the identified chaperones are also heat shock proteins. The general mechanism by which chaperones carry out their function usually involves multiple rounds of regulated binding and release of an unstable conformer of target polypeptides. The four main chaperone systems in the Escherichia coli cytoplasm are as follows. (1) Ribosome-associated trigger factor that assists in the folding of newly-synthesized nascent chains. (2) The Hsp 70 system consisting of DnaK (Hsp 70), its cofactor DnaJ (Hsp 40), and the nucleotide exchange factor GrpE. This system recognizes polypeptide chains in an extended conformation. (3) The Hsp 60 system, consisting of GroEL (Hsp 60) and its cofactor GroES (Hsp 10), which assists in the folding of compact folding intermediates that expose hydrophobic surfaces. (4) The Clp ATPases which are typically members of the Hsp 100 family of heat shock proteins. These ATPases can unfold proteins and disaggregate preformed protein aggregates to target them for degradation. Several advances have recently been made in characterizing the structure and function of all of these chaperone systems. These advances have provided us with a better understanding of the protein folding process in the cell.     

Hinged molecular capsules: synthesis and conformational control via temperature, pH, or solvent composition

Three new covalently linked molecular capsules were synthesized from their resorcinarene cavitand precursors in good yields. The capsules undergo reversible conformational switching between the closed "vase" form and the open "kite" form upon temperature or pH variation. The kite conformation obtained via either method in CDCl(3) switches to vase conformation upon addition of polar solvents such as acetone-d(6) or THF-d(8).     

Structural basis for matrix metalloproteinase 1-catalyzed collagenolysis

The proteolysis of collagen triple-helical structure (collagenolysis) is a poorly understood yet critical physiological process. Presently, matrix metalloproteinase 1 (MMP-1) and collagen triple-helical peptide models have been utilized to characterize the events and calculate the energetics of collagenolysis via NMR spectroscopic analysis of 12 enzyme-substrate complexes. The triple-helix is bound initially by the MMP-1 hemopexin-like (HPX) domain via a four amino acid stretch (analogous to type I collagen residues 782-785). The triple-helix is then presented to the MMP-1 catalytic (CAT) domain in a distinct orientation. The HPX and CAT domains are rotated with respect to one another compared with the X-ray "closed" conformation of MMP-1. Back-rotation of the CAT and HPX domains to the X-ray closed conformation releases one chain out of the triple-helix, and this chain is properly positioned in the CAT domain active site for subsequent hydrolysis. The aforementioned steps provide a detailed, experimentally derived, and energetically favorable collagenolytic mechanism, as well as significant insight into the roles of distinct domains in extracellular protease function.     

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

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

Solution conformation of a nitrobenzoxadiazole derivative of the polyene antibiotic nystatin: a FRET study

Nystatin is a polyene antibiotic frequently applied in the treatment of topical fungal infections. In this work, a 7-nitrobenz-2-oxa-1,3-diazole (NBD) hexanoyl amide derivative of nystatin was synthesized and its detailed photophysical characterization is presented. The average conformation of the labelled antibiotic in tetrahydrofuran, ethanol and methanol was determined by intramolecular (tetraene to NBD) fluorescence resonance energy transfer measurements. At variance with the literature [Can. J. Chem. 63 (1985) 77-85], it was concluded that there is no need to invoke a solvent-dependent conformational equilibrium between extended and closed conformers of the antibiotic, because the mean tetraene-to-NBD separating distance was found to remain constant (approximately 18 A) in all the solvents studied. In addition, the large solvent dependence of the fluorescence anisotropy observed for the non-derivatized nystatin, was rationalized on the basis of the prolate ellipsoidal geometry of the molecule. It was concluded that the rod shaped and amphipathic antibiotic remains monomeric in different solvents within the concentration range studied (2-20 microM).     

Probing Akt-inhibitor interaction by chemical cross-linking and mass spectrometry

The serine/threonine kinase Akt is a critical enzyme that regulates cell survival. As high Akt activity has been shown to contribute to the pathogenesis of various human malignancies, inhibition of Akt activation is a promising therapeutic strategy for cancers. We have previously demonstrated that changes in Akt interdomain arrangements from a closed to open conformation occur upon Akt-membrane interaction, which in turn allows Akt phosphorylation/activation. In the present study, we demonstrate a novel strategy to discern mechanisms for Akt inhibition based on Akt conformational changes using chemical cross-linking and ¹⁸O labeling mass spectrometry. By quantitative comparison of two interdomain cross-linked peptides, which represent the proximity of the domains involved, we found that the binding of Akt to an inhibitor (PI analog) caused the open interdomain conformation where the PH and regulatory domains moved away from the kinase domain, even before interacting with membranes, subsequently preventing translocation of Akt to the plasma membrane. In contrast, the interdomain conformation remained unchanged after incubating with another type of inhibitor (peptide TCL1). Subsequent interaction with unilamellar vesicles suggested that TCL1 impaired particularly the opening of the PH domain for exposing T308 for phosphorylation at the plasma membrane. This novel approach based on the conformation-based molecular interaction mechanism should be potentially useful for drug discovery efforts for specific Akt inhibitors or anti-tumor agents.     

NMR Analysis of Apo Glutamine‐Binding Protein Exposes Challenges in the Study of Interdomain Dynamics

Glutamine‐binding protein (GlnBP) displays an apo, “open” and a holo, “closed” crystal form, mutually related by a rigid‐body reorientation of its domains. A fundamental question about such large‐scale conformational transitions, whether the closed state exists in the absence of ligand, is controversial in the case of GlnBP. NMR observations have indicated no evidence of the closed form, whereas experimentally validated computations have suggested a remarkable ca. 40 % population. Herein, a paramagnetic NMR strategy designed to detect the putative apo‐closed species shows that a major population of the latter is highly improbable. Further, NMR residual dipolar couplings collected under three anisotropic conditions do not reveal differential domain alignment and establish that the average solution conformation is satisfied by the apo‐open crystal structure. Our results indicate that the computational prediction of large‐scale interdomain motions is not trivial and may lead to erroneous conclusions without proper experimental validation.     

Aggregation Pathways of the Amyloid β(1-42) Peptide Depend on Its Colloidal Stability and Ordered β-Sheet Stacking

Amyloid β (Aβ) fibrils are present as a major component in senile plaques, the hallmark of Alzheimer's disease (AD). Diffuse plaques (nonfibrous, loosely packed Aβ aggregates) containing amorphous Aβ aggregates are also formed in brain. This work examines the influence of Cu(2+) complexation by Aβ on the aggregation process in the context of charge and structural variations. Changes in the surface charges of Aβ molecules due to Cu(2+) binding, measured with a ζ-potential measurement device, were correlated with the aggregate morphologies examined by atomic force microscopy. As a result of the charge variation, the "colloid-like" stability of the aggregation intermediates, which is essential to the fibrillation process, is affected. Consequently, Cu(2+) enhances the amorphous aggregate formation. By monitoring variations in the secondary structures with circular dichroism spectroscopy, a direct transformation from the unstructured conformation to the β-sheet structure was observed for all types of aggregates observed (oligomers, fibrils, and/or amorphous aggregates). Compared to the Aβ aggregation pathway in the absence of Cu(2+) and taking other factors affecting Aβ aggregation (i.e., pH and temperature) into account, our investigation indicates that formations of amorphous and fibrous aggregates diverge from the same β-sheet-containing partially folded intermediate. This study suggests that the hydrophilic domain of Aβ also plays a role in the Aβ aggregation process. A kinetic model was proposed to account for the effects of the Cu(2+) binding on these two aggregation pathways in terms of charge and structural variations.