Thermodynamic stability of histone H3 is a necessary but not sufficient driving force for its evolutionary conservation

Determining the forces that conserve amino acid positions in proteins across species is a fundamental pursuit of molecular evolution. Evolutionary conservation is driven by either a protein's function or its thermodynamic stability. Highly conserved histone proteins offer a platform to evaluate these driving forces. While the conservation of histone H3 and H4 "tail" domains and surface residues are driven by functional importance, the driving force behind the conservation of buried histone residues has not been examined. Using a computational approach, we determined the thermodynamically preferred amino acids at each buried position in H3 and H4. In agreement with what is normally observed in proteins, we find a significant correlation between thermodynamic stability and evolutionary conservation in the buried residues in H4. In striking contrast, we find that thermodynamic stability of buried H3 residues does not correlate with evolutionary conservation. Given that these H3 residues are not post-translationally modified and only regulate H3-H3 and H3-H4 stabilizing interactions, our data imply an unknown function responsible for driving conservation of these buried H3 residues.     

Flanking bases influence the nature of DNA distortion by platinum 1,2-intrastrand (GG) cross-links

The differences in efficacy and molecular mechanisms of platinum anti-cancer drugs cisplatin (CP) and oxaliplatin (OX) are thought to be partially due to the differences in the DNA conformations of the CP and OX adducts that form on adjacent guanines on DNA, which in turn influence the binding of damage-recognition proteins that control downstream effects of the adducts. Here we report a comprehensive comparison of the structural distortion of DNA caused by CP and OX adducts in the TGGT sequence context using nuclear magnetic resonance (NMR) spectroscopy and molecular dynamics (MD) simulations. When compared to our previous studies in other sequence contexts, these structural studies help us understand the effect of the sequence context on the conformation of Pt-GG DNA adducts. We find that both the sequence context and the type of Pt-GG DNA adduct (CP vs. OX) play an important role in the conformation and the conformational dynamics of Pt-DNA adducts, possibly explaining their influence on the ability of many damage-recognition proteins to bind to Pt-DNA adducts.     

Automated minimization of steric clashes in protein structures

Molecular modeling of proteins including homology modeling, structure determination, and knowledge-based protein design requires tools to evaluate and refine three-dimensional protein structures. Steric clash is one of the artifacts prevalent in low-resolution structures and homology models. Steric clashes arise due to the unnatural overlap of any two nonbonding atoms in a protein structure. Usually, removal of severe steric clashes in some structures is challenging since many existing refinement programs do not accept structures with severe steric clashes. Here, we present a quantitative approach of identifying steric clashes in proteins by defining clashes based on the Van der Waals repulsion energy of the clashing atoms. We also define a metric for quantitative estimation of the severity of clashes in proteins by performing statistical analysis of clashes in high-resolution protein structures. We describe a rapid, automated, and robust protocol, Chiron, which efficiently resolves severe clashes in low-resolution structures and homology models with minimal perturbation in the protein backbone. Benchmark studies highlight the efficiency and robustness of Chiron compared with other widely used methods. We provide Chiron as an automated web server to evaluate and resolve clashes in protein structures that can be further used for more accurate protein design.     

Channel Gating Dependence on Pore Lining Helix Glycine Residues in Skeletal Muscle Ryanodine Receptor

Type 1 ryanodine receptors (RyR1s) release Ca²⁺ from the sarcoplasmic reticulum to initiate skeletal muscle contraction. The role of RyR1-G4934 and -G4941 in the pore-lining helix in channel gating and ion permeation was probed by replacing them with amino acid residues of increasing side chain volume. RyR1-G4934A, -G4941A, and -G4941V mutant channels exhibited a caffeine-induced Ca²⁺ release response in HEK293 cells and bound the RyR-specific ligand [³H]ryanodine. In single channel recordings, significant differences in the number of channel events and mean open and close times were observed between WT and RyR1-G4934A and -G4941A. RyR1-G4934A had reduced K⁺ conductance and ion selectivity compared with WT. Mutations further increasing the side chain volume at these positions (G4934V and G4941I) resulted in reduced caffeine-induced Ca²⁺ release in HEK293 cells, low [³H]ryanodine binding levels, and channels that were not regulated by Ca²⁺ and did not conduct Ca²⁺ in single channel measurements. Computational predictions of the thermodynamic impact of mutations on protein stability indicated that although the G4934A mutation was tolerated, the G4934V mutation decreased protein stability by introducing clashes with neighboring amino acid residues. In similar fashion, the G4941A mutation did not introduce clashes, whereas the G4941I mutation resulted in intersubunit clashes among the mutated isoleucines. Co-expression of RyR1-WT with RyR1-G4934V or -G4941I partially restored the WT phenotype, which suggested lessening of amino acid clashes in heterotetrameric channel complexes. The results indicate that both glycines are important for RyR1 channel function by providing flexibility and minimizing amino acid clashes.     

Interaction with the 5D3 monoclonal antibody is regulated by intramolecular rearrangements but not by covalent dimer formation of the human ABCG2 multidrug transporter

Human ABCG2 is a plasma membrane glycoprotein working as a homodimer or homo-oligomer. The protein plays an important role in the protection/detoxification of various tissues and may also be responsible for the multidrug-resistant phenotype of cancer cells. In our previous study we found that the 5D3 monoclonal antibody shows a function-dependent reactivity to an extracellular epitope of the ABCG2 transporter. In the current experiments we have further characterized the 5D3-ABCG2 interaction. The effect of chemical cross-linking and the modulation of extracellular S-S bridges on the transporter function and 5D3 reactivity of ABCG2 were investigated in depth. We found that several protein cross-linkers greatly increased 5D3 labeling in ABCG2 expressing HEK cells; however, there was no correlation between covalent dimer formation, the inhibition of transport activity, and the increase in 5D3 binding. Dithiothreitol treatment, which reduced the extracellular S-S bridge-forming cysteines of ABCG2, had no effect on transport function but caused a significant decrease in 5D3 binding. When analyzing ABCG2 mutants carrying Cys-to-Ala changes in the extracellular loop, we found that the mutant C603A (lacking the intermolecular S-S bond) showed comparable transport activity and 5D3 reactivity to the wild-type ABCG2. However, disruption of the intramolecular S-S bridge (in C592A, C608A, or C592A/C608A mutants) in this loop abolished 5D3 binding, whereas the function of the protein was preserved. Based on these results and ab initio folding simulations, we propose a model for the large extracellular loop of the ABCG2 protein.     

The coordinated evolution of yeast proteins is constrained by functional modularity

Functional modularity is a key attribute of cellular systems and has important roles in evolution. However, the extent to which functional modularity affects protein evolution is largely unknown. Here, we analyzed the evolution of both sequence and expression level of proteins in the yeast Saccharomyces cerevisiae and found that proteins within the same functional modules evolve at more similar rates than those between different modules. We also found stronger co-evolution of expression levels between proteins within functional modules than between them. These results suggest that a coordinated evolution of both the sequence and expression level of proteins is constrained by functional modularity.     

Ab initio folding of proteins with all-atom discrete molecular dynamics

Discrete molecular dynamics (DMD) is a rapid sampling method used in protein folding and aggregation studies. Until now, DMD was used to perform simulations of simplified protein models in conjunction with structure-based force fields. Here, we develop an all-atom protein model and a transferable force field featuring packing, solvation, and environment-dependent hydrogen bond interactions. We performed folding simulations of six small proteins (20-60 residues) with distinct native structures by the replica exchange method. In all cases, native or near-native states were reached in simulations. For three small proteins, multiple folding transitions are observed, and the computationally characterized thermodynamics are in qualitative agreement with experiments. The predictive power of all-atom DMD highlights the importance of environment-dependent hydrogen bond interactions in modeling protein folding. The developed approach can be used for accurate and rapid sampling of conformational spaces of proteins and protein-protein complexes and applied to protein engineering and design of protein-protein interactions.     

iFoldRNA: three-dimensional RNA structure prediction and folding

Three-dimensional RNA structure prediction and folding is of significant interest in the biological research community. Here, we present iFoldRNA, a novel web-based methodology for RNA structure prediction with near atomic resolution accuracy and analysis of RNA folding thermodynamics. iFoldRNA rapidly explores RNA conformations using discrete molecular dynamics simulations of input RNA sequences. Starting from simplified linear-chain conformations, RNA molecules (<50 nt) fold to native-like structures within half an hour of simulation, facilitating rapid RNA structure prediction. All-atom reconstruction of energetically stable conformations generates iFoldRNA predicted RNA structures. The predicted RNA structures are within 2-5 A root mean squre deviations (RMSDs) from corresponding experimentally derived structures. RNA folding parameters including specific heat, contact maps, simulation trajectories, gyration radii, RMSDs from native state, fraction of native-like contacts are accessible from iFoldRNA. We expect iFoldRNA will serve as a useful resource for RNA structure prediction and folding thermodynamic analyses. AVAILABILITY: http://iFoldRNA.dokhlab.org.     

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.