Characterizing protein protonation microstates using Monte Carlo sampling

Proteins are polyelectrolytes with acidic or basic amino acids making up ≈25% of the residues. The protonation state of all Asp, Glu, Arg, Lys, His and other protonatable residues, cofactors and ligands define each protonation microstate. As all of these residues will not be fully ionized or neutral, proteins exist in a mixture of microstates. The microstate distribution changes with pH. As the protein environment modifies the proton affinity of each site the distribution may also change in different reaction intermediates or as ligands are bound. Particular protonation microstates may be required for function, while others exist simply because there are many states with similar energy. Here, the protonation microstates generated in Monte Carlo sampling in MCCE are characterized in HEW lysozyme as a function of pH and bacterial photosynthetic reaction centers (RCs) in different reaction intermediates. The lowest energy and highest probability microstates are compared. The ∆G, ∆H and ∆S between the four protonation states of Glu35 and Asp52 in lysozyme are shown to be calculated with reasonable precision. A weighted Pearson correlation analysis identifies coupling between residue protonation states in RCs and how they change when the quinone in the QB site is reduced.

Unravelling Constant pH Molecular Dynamics in Oligopeptides with Explicit Solvation Model

An accurate description of the protonation state of amino acids is essential to correctly simulate the conformational space and the mechanisms of action of proteins or other biochemical systems. The pH and the electrochemical environments are decisive factors to define the effective pKa of amino acids and, therefore, the protonation state. However, they are poorly considered in Molecular Dynamics (MD) simulations. To deal with this problem, constant pH Molecular Dynamics (cpHMD) methods have been developed in recent decades, demonstrating a great ability to consider the effective pKa of amino acids within complex structures. Nonetheless, there are very few studies that assess the effect of these approaches in the conformational sampling. In a previous work of our research group, we detected strengths and weaknesses of the discrete cpHMD method implemented in AMBER when simulating capped tripeptides in implicit solvent. Now, we progressed this assessment by including explicit solvation in these peptides. To analyze more in depth the scope of the reported limitations, we also carried out simulations of oligopeptides with distinct positions of the titratable amino acids. Our study showed that the explicit solvation model does not improve the previously noted weaknesses and, furthermore, the separation of the titratable amino acids in oligopeptides can minimize them, thus providing guidelines to improve the conformational sampling in the cpHMD simulations.

pH-dependent polymorphism of the structure of SARS-CoV-2 nsp7

The solution structure of SARS-CoV-2 nonstructural protein 7 (nsp7) at pH 7.0 has been determined by NMR spectroscopy. nsp7 is conserved in the coronavirinae subfamily and is an essential co-factor of the viral RNA-dependent RNA polymerase for active and processive replication. Similar to the previously deposited structures of SARS-CoV-1 nsp7 at acidic and basic conditions, SARS-CoV-2 nsp7 has a helical bundle folding at neutral pH. Remarkably, the α4 helix shows gradual dislocation from the core α2-α3 structure as pH increases from 6.5 to 7.5. The protonation state of residue H36 contributes to the change of nsp7’s intramolecular interactions, and thus, to the structural variation near-neutral pH. Spin-relaxation results revealed that all three loop regions in nsp7 possess dynamic properties associated with this structural variation.

Effects of Cardiolipin on the Conformational Dynamics of Membrane-Anchored Bcl-xL

The anti-apoptotic protein Bcl-xL regulates apoptosis by preventing the permeation of the mitochondrial outer membrane by pro-apoptotic pore-forming proteins, which release apoptotic factors into the cytosol that ultimately lead to cell death. Two different membrane-integrated Bcl-xL constructs have been identified: a membrane-anchored and a membrane-inserted conformation. Here, we use molecular dynamics simulations to study the effect of the mitochondrial specific lipid cardiolipin and the protein protonation state on the conformational dynamics of membrane-anchored Bcl-xL. The analysis reveals that the protonation state of the protein and cardiolipin content of the membrane modulate the orientation of the soluble head region (helices α1 through α7) and hence the exposure of its BH3-binding groove, which is required for its interaction with pro-apoptotic proteins.