Motivation

Cys450 side chain is involved in light state in bond formation with the cofactor. Cys450 can assume two conformational states, called here ON and OFF, and corresponding respectively, to the Sg being near or far from FMN.

The more populated conformer is the OFF, as showed by X-ray structures (OFF/ON=90/10 at low temperature and =70/30 at room temperature) but for the bond to form the Cys should switch the less favoured ON. The ON-OFF equilibrium is believed to be influenced by nearby residues. Moreover the steric hindrance of some nearby apolar residues is thought to have an effect also in the process of bond splitting. By looking at the structure and by examining previous mutation analysis we hypothesized that mutation favouring the OFF conformation of Cys450 in the dark state would likely result in favouring the splitting of the covalent adduct. On the other way round a strategy to increase the lifetime of covalent adduct could consist in mutating residue as to stabilize the ON conformation.

Since our simulation are capable of reproducing the experimental OFF/ON conformational equilibrium in the dark state (about 70/30) we decided to perform in-silico mutations, and to compute the OFF/ON ratios of mutants from MD simulations, as to have a guide to experimentally perform the mutations leading to lower OFF/ON ratios. Of course the OFF/ON conformational equilibrium is only indirectly connected to the energy barrier of the bond splitting, but some experimental evidences suggest a link between Cys conformational propensity and the splitting reaction.

We planned to perform QM calculations as to directly study the effect of specific mutations on the splitting reaction rate, but since those calcualtions require much more time we have at our disposal, we used for our mutations only chemical intuition and the results of the simulations of mutated systems.

Active site

The selected residues for mutations are ILE427 in purple on the next video, LEU453 in orange and VAL416 in yellow.

Another video using the surf function shows better the steric interactions between those residues. The residues have the same color as before, but we added the CYS450 in blue.

ILE427 to PHE

The isoleucine 427 is quite big. We tried in silico to change it to a cystein. As the cystein is a small residue that could lead some missfolding in vivo, we also planned to change it to valine in vivo. On the left, the ILE427 is colored in blue. On the right picture, the same residue is mutated to a PHE (I427F).

After 35 ns of simulation, we already see a clear difference in the ON/OFF ratio see the comparison with the wild type dark state result.

The calculated ratios are:

% of time the cystein is in the IN conformation = 57.2%.

% of time the cystein is in the OUT conformation = 42.8%.

These results are very encouraging because it proves that the cystein is "pushed" in the direction of the FMN by the steric interaction of the phenylalanine with the csytein.

Here is a plot of the evolution of dihedral angle during the simulation:

The plot of the dihedral angle validates the previous observations.

LEU453 to GLY

The light blue is the isoleucine and the dark blue is the new cystein, which is really much smaller, is in dark blue. This mutation gave good results in previous articles.

After 37 ns of simulation, we already see a small difference in the ON/OFF ratio see the comparison with the wild type dark state result.

The calculated ratios are:

% of time the cystein is in the IN conformation = 31.16%.

% of time the cystein is in the OUT conformation = 68.84%.

These results are less relevant because the cystein side chain have a little bit more place in the direction of the FMN than in the wild type. But, the ratios are not strongly influenced.

So, we could make an hypothesis: Maybe the steric pushing is more efficient to push the cystein side chain toward the C4A of the FMN than to let the cystein take the empty place. I might be interesting to do a simulation with both I427F and L453G to verify this.

Here is a plot of the evolution of dihedral angle during the simulation:

The plot of the dihedral angle is consistent with the previous observations.

VAL416 to ILE