Serpin latency transition at atomic resolution.

Protease inhibition by serpins requires a large conformational transition from an active, metastable state to an inactive, stable state. Similar reactions can also occur in the absence of proteases, and these latency transitions take hours, making their time scales many orders of magnitude larger than are currently accessible using conventional molecular dynamics simulations. Using a variational path sampling algorithm, we simulated the entire serpin active-to-latent transition in all-atom detail with a physically realistic force field using a standard computing cluster. These simulations provide a unifying picture explaining existing experimental data for the latency transition of the serpin plasminogen activator inhibitor-1 (PAI-1). They predict a long-lived intermediate that resembles a previously proposed, partially loop-inserted, prelatent state; correctly predict the effects of PAI-1 mutations on the kinetics; and provide a potential means to identify ligands able to accelerate the latency transition. Interestingly, although all of the simulated PAI-1 variants readily access the prelatent intermediate, this conformation is not populated in the active-to-latent transition of another serpin, α1-antitrypsin, which does not readily go latent. Thus, these simulations also help elucidate why some inhibitory serpin families are more conformationally labile than others.

Unfolding thermodynamics of cysteine-rich proteins and molecular thermal-adaptation of marine ciliates.

Euplotes nobilii and Euplotes raikovi are phylogenetically closely allied species of marine ciliates, living in polar and temperate waters, respectively. Their evolutional relation and the sharply different temperatures of their natural environments make them ideal organisms to investigate thermal-adaptation. We perform a comparative study of the thermal unfolding of disulfide-rich protein pheromones produced by these ciliates. Recent circular dichroism (CD) measurements have shown that the two psychrophilic (E. nobilii) and mesophilic (E. raikovi) protein families are characterized by very different melting temperatures, despite their close structural homology. The enhanced thermal stability of the E. raikovi pheromones is realized notwithstanding the fact that these proteins form, as a rule, a smaller number of disulfide bonds. We perform Monte Carlo (MC) simulations in a structure-based coarse-grained (CG) model to show that the higher stability of the E. raikovi pheromones is due to the lower locality of the disulfide bonds, which yields a lower entropy increase in the unfolding process. Our study suggests that the higher stability of the mesophilic E. raikovi phermones is not mainly due to the presence of a strongly hydrophobic core, as it was proposed in the literature. In addition, we argue that the molecular adaptation of these ciliates may have occurred from cold to warm, and not from warm to cold. To provide a testable prediction, we identify a point-mutation of an E. nobilii pheromone that should lead to an unfolding temperature typical of that of E. raikovi pheromones.