Electrostatic interactions contribute to the control of intramolecular thiol-disulfide isomerization in a protein.

The roles of structural factors and of electrostatic interactions with the environment on the outcome of thiol-disulfide exchange reactions were investigated in a mutated immunoglobulin domain (I27*) under mechanical stress. An extensive ensemble of molecular dynamics trajectories was generated by means of QM/MM simulations for a total sampling of 5.7 µs. A significant number of thiol-disulfide exchanges were observed, and the Cys32 thiolate preferred to attack Cys55 over Cys24, in agreement with previous experimental and computational studies. The structural features as well as electronic structures of the thiol-disulfide system along the reaction were analyzed, as were the electrostatic interactions with the environment. The previous findings of better accessibility of Cys55 were confirmed. Additionally, the reaction was found to be directed by the electrostatic interactions of the involved sulfur atoms with the molecular environment. The relationships of atomic charges, which stem from the electrostatic interactions, lead to the kinetic preference of the attack on Cys55. Further, QM/MM metadynamics simulations of thiol-disulfide exchange in a small model system with varied artificial external electric potentials revealed changes in reaction kinetics of the same magnitude as in I27*. Therefore, the electrostatic interactions are confirmed to play a role in the regioselectivity of the thiol-disulfide exchange reactions in the protein.

On the mechanism of spontaneous thiol-disulfide exchange in proteins.

The thiol-disulfide exchange reaction in model systems and small peptides was investigated by means of a combined QM/MM metadynamics scheme. The free energy landscapes of these systems were generated, providing the structures of reactants and products with atomic detail, as well as the heights of free energy barriers (or, activation energies) opposing the spontaneous exchange. A QM/MM scheme with purely classical water turned out to be an efficient and accurate compromise solution. The calculations yielded the expected symmetric trisulfide transition state at S-S distances of 2.7 Å, interestingly, with a slight deviation from linearity at an S-S-S angle of 165°. The structure of the transition state as well as the free energy barrier were very similar for the intramolecular thiol-disulfide reactions in model peptides. While CXC disulfide bonds were found sterically unfavorable, CXXC were favored over longer-range disulfide bonds along the peptide backbone, in line with the high abundance of CXXC motifs in redox proteins.

Orthogonal 19 F-Labeling for Solid-State NMR Spectroscopy Reveals the Conformation and Orientation of Short Peptaibols in Membranes.

Peptaibols are promising drug candidates in view of their interference with cellular membranes. Knowledge of their lipid interactions and membrane-bound structure is needed to understand their activity and should be, in principle, accessible by solid-state NMR spectroscopy. However, their unusual amino acid composition and noncanonical conformations make it very challenging to find suitable labels for NMR spectroscopy. Particularly in the case of short sequences, new strategies are required to maximize the structural information that can be obtained from each label. Herein, l-3-(trifluoromethyl)bicyclopent[1.1.1]-1-ylglycine, (R)- and (S)-trifluoromethylalanine, and 15 N-backbone labels, each probing a different direction in the molecule, have been combined to elucidate the conformation and membrane alignment of harzianin HK-VI. For the short sequence of 11 amino acids, 12 orientational constraints have been obtained by using 19 F and 15 N NMR spectroscopy. This strategy revealed a ß-bend ribbon structure, which becomes realigned in the membrane from a surface-parallel state towards a membrane-spanning state, with increasing positive spontaneous curvature of the lipids.

Structural Behavior of the Peptaibol Harzianin HK VI in a DMPC Bilayer: Insights from MD Simulations.

Microsecond molecular dynamics simulations of harzianin HK VI (HZ) interacting with a dimyristoylphosphatidylcholine bilayer were performed at the condition of low peptide-to-lipid ratio. Two orientations of HZ molecule in the bilayer were found and characterized. In the orientation perpendicular to the bilayer surface, HZ induces a local thinning of the bilayer. When inserted into the bilayer parallel to its surface, HZ is located nearly completely within the hydrophobic region of the bilayer. A combination of solid-state NMR and circular dichroism experiments found the latter orientation to be dominant. An extended sampling simulation provided qualitative results and showed the same orientation to be a global minimum of free energy. The secondary structure of HZ was characterized, and it was found to be located in the 310-helical family. The specific challenges of computer simulation of nonpolar peptides are discussed briefly.