Probing the Flexibility of Large Conformational Changes in Protein Structures through Local Perturbations

Protein conformational changes and dynamic behavior are fundamental for such processes as catalysis, regulation, and substrate recognition. Although protein dynamics have been successfully explored in computer simulation, there is an intermediate-scale of motions that has proven difficult to simulate—the motion of individual segments or domains that move independently of the body the protein. Here, we introduce a molecular-dynamics perturbation method, the Rotamerically Induced Perturbation (RIP), which can generate large, coherent motions of structural elements in picoseconds by applying large torsional perturbations to individual sidechains. Despite the large-scale motions, secondary structure elements remain intact without the need for applying backbone positional restraints. Owing to its computational efficiency, RIP can be applied to every residue in a protein, producing a global map of deformability. This map is remarkably sparse, with the dominant sites of deformation generally found on the protein surface. The global map can be used to identify loops and helices that are less tightly bound to the protein and thus are likely sites of dynamic modulation that may have important functional consequences. Additionally, they identify individual residues that have the potential to drive large-scale coherent conformational change. Applying RIP to two well-studied proteins, Dihdydrofolate Reductase and Triosephosphate Isomerase, which possess functionally-relevant mobile loops that fluctuate on the microsecond/millisecond timescale, the RIP deformation map identifies and recapitulates the flexibility of these elements. In contrast, the RIP deformation map of α-lytic protease, a kinetically stable protein, results in a map with no significant deformations. In the N-terminal domain of HSP90, the RIP deformation map clearly identifies the ligand-binding lid as a highly flexible region capable of large conformational changes. In the Estrogen Receptor ligand-binding domain, the RIP deformation map is quite sparse except for one large conformational change involving Helix-12, which is the structural element that allosterically links ligand binding to receptor activation. RIP analysis has the potential to discover sites of functional conformational changes and the linchpin residues critical in determining these conformational states

Cooperative Transition between Open and Closed Conformations in Potassium Channels

Potassium (K+) ion channels switch between open and closed conformations. The nature of this important transition was revealed by comparing the X-ray crystal structures of the MthK channel from Methanobacterium thermoautotrophicum, obtained in its open conformation, and the KcsA channel from Streptomyces lividans, obtained in its closed conformation. We analyzed the dynamic characteristics and energetics of these homotetrameric structures in order to study the role of the intersubunit cooperativity in this transition. For this, elastic models and in silico alanine-scanning mutagenesis were used, respectively. Reassuringly, the calculations manifested motion from the open (closed) towards the closed (open) conformation. The calculations also revealed a network of dynamically and energetically coupled residues. Interestingly, the network suggests coupling between the selectivity filter and the gate, which are located at the two ends of the channel pore. Coupling between these two regions was not observed in calculations that were conducted with the monomer, which emphasizes the importance of the intersubunit interactions within the tetrameric structure for the cooperative gating behavior of the channel

Deuteration protects asparagine residues against racemization

Racemization in proteins and peptides at sites of L-asparaginyl and L-aspartyl residues contributes to their spontaneous degradation, especially in the biological aging process. Amino acid racemization involves deprotonation of the alpha carbon and replacement of the proton in the opposite stereoconfiguration; this reaction is much faster for aspartate/asparagine than for other amino acids because these residues form a succinimide ring in which resonance stabilizes the carbanion resulting from proton loss. To determine if the replacement of the hydrogen atom on the alpha carbon with a deuterium atom might decrease the rate of racemization and thus stabilize polypeptides, we synthesized a hexapeptide, VYPNGA, in which the three carbon-bound protons in the asparaginyl residue were replaced with deuterium atoms. Upon incubation of this peptide in pH 7.4 buffer at 37 °C, we found that the rate of deamidation via the succinimide intermediate was unchanged by the presence of the deuterium atoms. However, the accumulation of the D-aspartyl and D-isoaspartyl forms resulting from racemization and hydrolysis of the succinimide was decreased more than five-fold in the deuterated peptide over a 20 day incubation at physiological temperature and pH. Additionally, we found that the succinimide intermediate arising from the degradation of the deuterated asparaginyl peptide was slightly less likely to open to the isoaspartyl configuration than was the protonated succinimide. These findings suggest that the kinetic isotope effect resulting from the presence of deuteriums in asparagine residues can limit the accumulation of at least some of the degradation products that arise as peptides and proteins age

Protein motions, dynamic effects and thermal stability in dihydrofolate reductase from the hyperthermophile thermotoga maritima

Dihydrofolate reductase (DHFR) has long been used as a model system in studies of the relationship between enzyme structure and catalysis. DHFR from the hyperthermophilic bacterium Thermotoga maritima (TmDHFR) is substantially different to other chromosomal DHFRs. It is dimeric where most others are monomeric, it lacks the conformational behaviour of monomeric DHFRs, and the kinetics of the catalysed reaction are significantly different. Experimental and computational studies of TmDHFR and comparison to other DHFRs have yielded deep insights into the role of enzyme motions and dynamics in catalysis. Mutational studies and formation of hybrids between TmDHFR and a monomeric homologue have demonstrated that dimerisation is required for extreme thermostability, but also leads to an inability to adequately close the active site with detrimental effects for the speed of the catalysed reaction. However, in common with other DHFRs there is no involvement of large-scale enzyme motions in the chemical reaction itself and dynamic coupling to the reaction coordinate is efficiently minimised. Studies of DHFRs from hyperthermophilic organisms and comparisons to their mesophilic counterparts remain a rich source of information on the fundamental nature of enzyme catalysis