Understanding Cysteine Reactivity in Protein Environments with Electric Fields

Area of Science:

• Biochemistry
• Computational Biology
• Protein Chemistry

Background:

• Cysteine residues' reactivity is determined by their protonation state (thiol vs. thiolate).
• Predicting cysteine pKa values in proteins is difficult due to local environmental effects.
• Traditional methods like constant pH molecular dynamics (MD) can struggle with accurate pKa prediction.

Purpose of the Study:

• To develop a new computational approach for predicting cysteine reactivity.
• To accurately determine the protonation state and pKa of cysteine residues in various protein locations.
• To rationalize local pH and pKa effects within proteins.

Main Methods:

• Utilized electric field calculations focused on the thiol/thiolate group of cysteine.
• Validated the approach across different protein environments (active site, surface, interior).
• Applied the method to Cys-25 in papaya protease omega, a known challenge for other techniques.

Main Results:

• The electric field calculation method successfully predicted pKa shifts consistent with experimental data.
• Accurate predictions were achieved for cysteine residues in diverse protein microenvironments.
• Decomposition of electric fields provided insights into local pH and pKa modulation.

Conclusions:

• The novel electric field calculation approach offers a robust method for predicting cysteine pKa and reactivity in proteins.
• This method overcomes limitations of traditional techniques, particularly in heterogeneous protein environments.
• The approach provides a parameter-free way to understand local effects influencing cysteine ionization states.