The maturation of HIV‐1 protease precursor studied by discrete molecular dynamics

The equilibrium properties of a HIV‐1‐protease precursor are studied by means of an efficient molecular dynamics scheme, which allows for the simulation of the folding of the protein monomers and their dimerization into an active form and compare them with those of the mature protein. The results of the model provide, with atomic detail, an overall account of several experimental findings, including the NMR conformation of the mature dimer, the calorimetric properties of the system, the effects of the precursor tail on the dimerization constant, the secondary chemical shifts of the monomer, and the paramagnetic relaxation enhancement data associated with the conformations of the precursor. It is found that although the mature protein can dimerize in a unique, single way, the precursor populates several dimeric conformations in which monomers are always native‐like, but their binding can be non‐native. Proteins 2014; 82:633–639. © 2013 Wiley Periodicals, Inc.     

Intrinsic Disorder Mediates Cooperative Signal Transduction in STIM1

Intrinsically disordered domains have been reported to play important roles in signal transduction networks by introducing cooperativity into protein–protein interactions. Unlike intrinsically disordered domains that become ordered upon binding, the EF-SAM domain in the stromal interaction molecule (STIM) 1 is distinct in that it is ordered in the monomeric state and partially unfolded in its oligomeric state, with the population of the two states depending on the local Ca^2 + concentration. The oligomerization of STIM1, which triggers extracellular Ca^2 + influx, exhibits cooperativity with respect to the local endoplasmic reticulum Ca^2 + concentration. Although the physiological importance of the oligomerization reaction is well established, the mechanism of the observed cooperativity is not known. Here, we examine the response of the STIM1 EF-SAM domain to changes in Ca^2 + concentration using mathematical modeling based on in vitro experiments. We find that the EF-SAM domain partially unfolds and dimerizes cooperatively with respect to Ca^2 + concentration, with Hill coefficients and half-maximal activation concentrations very close to the values observed in vivo for STIM1 redistribution and extracellular Ca^2 + influx. Our mathematical model of the dimerization reaction agrees quantitatively with our analytical ultracentrifugation-based measurements and previously published free energies of unfolding. A simple interpretation of these results is that Ca^2 + loss effectively acts as a denaturant, enabling cooperative dimerization and robust signal transduction. We present a structural model of the Ca^2 +-unbound EF-SAM domain that is consistent with a wide range of evidence, including resistance to proteolytic cleavage of the putative dimerization portion.     

Insight into the Structure of the “Unstructured” Tau Protein

1. Download : Download high-res image (287KB)
2. Download : Download full-size image

     

Allosteric modulation balances thermodynamic stability and restores function of ΔF508 CFTR

Most cystic fibrosis is caused by a deletion of a single residue (F508) in CFTR (cystic fibrosis transmembrane conductance regulator) that disrupts the folding and biosynthetic maturation of the ion channel protein. Progress towards understanding the underlying mechanisms and overcoming the defect remains incomplete. Here, we show that the thermal instability of human ΔF508 CFTR channel activity evident in both cell-attached membrane patches and planar phospholipid bilayers is not observed in corresponding mutant CFTRs of several non-mammalian species. These more stable orthologs are distinguished from their mammalian counterparts by the substitution of proline residues at several key dynamic locations in first N-terminal nucleotide-binding domain (NBD1), including the structurally diverse region, the γ-phosphate switch loop, and the regulatory insertion. Molecular dynamics analyses revealed that addition of the prolines could reduce flexibility at these locations and increase the temperatures of unfolding transitions of ΔF508 NBD1 to that of the wild type. Introduction of these prolines experimentally into full-length human ΔF508 CFTR together with the already recognized I539T suppressor mutation, also in the structurally diverse region, restored channel function and thermodynamic stability as well as its trafficking to and lifetime at the cell surface. Thus, while cellular manipulations that circumvent its culling by quality control systems leave ΔF508 CFTR dysfunctional at physiological temperature, restoration of the delicate balance between the dynamic protein's inherent stability and channel activity returns a near-normal state.     

A Structural Model of the Pore-Forming Region of the Skeletal Muscle Ryanodine Receptor (RyR1)

Ryanodine receptors (RyRs) are ion channels that regulate muscle contraction by releasing calcium ions from intracellular stores into the cytoplasm. Mutations in skeletal muscle RyR (RyR1) give rise to congenital diseases such as central core disease. The absence of high-resolution structures of RyR1 has limited our understanding of channel function and disease mechanisms at the molecular level. Here, we report a structural model of the pore-forming region of RyR1. Molecular dynamics simulations show high ion binding to putative pore residues D4899, E4900, D4938, and D4945, which are experimentally known to be critical for channel conductance and selectivity. We also observe preferential localization of Ca²⁺ over K⁺ in the selectivity filter of RyR1. Simulations of RyR1-D4899Q mutant show a loss of preference to Ca²⁺ in the selectivity filter as seen experimentally. Electrophysiological experiments on a central core disease mutant, RyR1-G4898R, show constitutively open channels that conduct K⁺ but not Ca²⁺. Our simulations with G4898R likewise show a decrease in the preference of Ca²⁺ over K⁺ in the selectivity filter. Together, the computational and experimental results shed light on ion conductance and selectivity of RyR1 at an atomistic level.     

Topological determinants of protein domain swapping

Protein domain swapping has been repeatedly observed in a variety of proteins and is believed to result from destabilization due to mutations or changes in environment. Based on results from our studies and others, we propose that structures of the domain-swapped proteins are mainly determined by their native topologies. We performed molecular dynamics simulations of seven different proteins, known to undergo domain swapping experimentally, under mildly denaturing conditions and found in all cases that the domain-swapped structures can be recapitulated by using protein topology in a simple protein model. Our studies further indicated that, in many cases, domain swapping occurs at positions around which the protein tends to unfold prior to complete unfolding. This, in turn, enabled prediction of protein structural elements that are responsible for domain swapping. In particular, two distinct domain-swapped dimer conformations of the focal adhesion targeting domain of focal adhesion kinase were predicted computationally and were supported experimentally by data obtained from NMR analyses.     

Regulatory Insertion Removal Restores Maturation, Stability and Function of ΔF508 CFTR

The cystic fibrosis transmembrane conductance regulator (CFTR) epithelial anion channel is a large multidomain membrane protein that matures inefficiently during biosynthesis. Its assembly is further perturbed by the deletion of F508 from the first nucleotide-binding domain (NBD1) responsible for most cystic fibrosis. The mutant polypeptide is recognized by cellular quality control systems and is proteolyzed. CFTR NBD1 contains a 32-residue segment termed the regulatory insertion (RI) not present in other ATP-binding cassette transporters. We report here that RI deletion enabled F508 CFTR to mature and traffic to the cell surface where it mediated regulated anion efflux and exhibited robust single chloride channel activity. Long-term pulse-chase experiments showed that the mature ΔRI/ΔF508 had a T_1/2 of ∼ 14 h in cells, similar to the wild type. RI deletion restored ATP occlusion by NBD1 of ΔF508 CFTR and had a strong thermostabilizing influence on the channel with gating up to at least 40 °C. None of these effects of RI removal were achieved by deletion of only portions of RI. Discrete molecular dynamics simulations of NBD1 indicated that RI might indirectly influence the interaction of NBD1 with the rest of the protein by attenuating the coupling of the F508-containing loop with the F1-like ATP-binding core subdomain so that RI removal overcame the perturbations caused by F508 deletion. Restriction of RI to a particular conformational state may ameliorate the impact of the disease-causing mutation.     

Emergence of protein fold families through rational design

Diverse proteins with similar structures are grouped into families of homologs and analogs, if their sequence similarity is higher or lower, respectively, than 20%–30%. It was suggested that protein homologs and analogs originate from a common ancestor and diverge in their distinct evolutionary time scales, emerging as a consequence of the physical properties of the protein sequence space. Although a number of studies have determined key signatures of protein family organization, the sequence-structure factors that differentiate the two evolution-related protein families remain unknown. Here, we stipulate that subtle structural changes, which appear due to accumulating mutations in the homologous families, lead to distinct packing of the protein core and, thus, novel compositions of core residues. The latter process leads to the formation of distinct families of homologs. We propose that such differentiation results in the formation of analogous families. To test our postulate, we developed a molecular modeling and design toolkit, Medusa, to computationally design protein sequences that correspond to the same fold family. We find that analogous proteins emerge when a backbone structure deviates only 1–2 Å root-mean-square deviation from the original structure. For close homologs, core residues are highly conserved. However, when the overall sequence similarity drops to ~25%–30%, the composition of core residues starts to diverge, thereby forming novel families of protein homologs. This direct observation of the formation of protein homologs within a specific fold family supports our hypothesis. The conservation of amino acids in designed sequences recapitulates that of the naturally occurring sequences, thereby validating our computational design methodology.     

Clustering of identical oligomers in coding and noncoding DNA sequences.

We develop a quantitative method for analyzing repetitions of identical short oligomers in coding and noncoding DNA sequences. We analyze sequences presently available in the GenBank separately for primate, mammal, vertebrate, rodent, invertebrate and plant taxonomic partitions. We find that some oligomers "cluster" more than they would if randomly distributed, while other oligomers "repel" each other. To quantify this degree of clustering, we define clustering measures. We find that (i) clustering significantly differs in coding and noncoding DNA; (ii) in most cases, monomers, dimers and tetramers cluster in noncoding DNA but appear to repel each other in coding DNA. (iii) The degree of clustering for different sources (primates, invertebrates, and plants) is more conserved among these sources in the case of coding DNA than in the case of noncoding DNA. (iv) In contrast to other oligomers, we find that trimers always prefer to cluster. (v) Clustering of each particular oligomer is conserved within the same organism.