Formation Mechanism of Ion Channel in Channelrhodopsin-2: Molecular Dynamics Simulation and Steering Molecular Dynamics Simulations

Approaching the catalytic mechanism of protein lysine methyltransferases by biochemical and simulation techniques

Protein lysine methyltransferases (PKMTs) transfer up to three methyl groups to the side chains of lysine residues in proteins and fulfill important regulatory functions by controlling protein stability, localization and protein/protein interactions. The methylation reactions are highly regulated, and aberrant methylation of proteins is associated with several types of diseases including neurologic disorders, cardiovascular diseases, and various types of cancer. This review describes novel insights into the catalytic machinery of various PKMTs achieved by the combined application of biochemical experiments and simulation approaches during the last years, focusing on clinically relevant and well-studied enzymes of this group like DOT1L, SMYD1-3, SET7/9, G9a/GLP, SETD2, SUV420H2, NSD1/2, different MLLs and EZH2. Biochemical experiments have unraveled many mechanistic features of PKMTs concerning their substrate and product specificity, processivity and the effects of somatic mutations observed in PKMTs in cancer cells. Structural data additionally provided information about the substrate recognition, enzyme-substrate complex formation, and allowed for simulations of the substrate peptide interaction and mechanism of PKMTs with atomistic resolution by molecular dynamics and hybrid quantum mechanics/molecular mechanics methods. These simulation technologies uncovered important mechanistic details of the PKMT reaction mechanism including the processes responsible for the deprotonation of the target lysine residue, essential conformational changes of the PKMT upon substrate binding, but also rationalized regulatory principles like PKMT autoinhibition. Further developments are discussed that could bring us closer to a mechanistic understanding of catalysis of this important class of enzymes in the near future. The results described here illustrate the power of the investigation of enzyme mechanisms by the combined application of biochemical experiments and simulation technologies.

Optogenetics: Illuminating the Future of Hearing Restoration and Understanding Auditory Perception

Hearing loss is a prevalent sensory impairment significantly affecting communication and quality of life. Traditional approaches for hearing restoration, such as cochlear implants, have limitations in frequency resolution and spatial selectivity. Optogenetics, an emerging field utilizing light-sensitive proteins, offers a promising avenue for addressing these limitations and revolutionizing hearing rehabilitation. This review explores the methods of introducing Channelrhodopsin- 2 (ChR2), a key light-sensitive protein, into cochlear cells to enable optogenetic stimulation. Viral- mediated gene delivery is a widely employed technique in optogenetics. Selecting a suitable viral vector, such as adeno-associated viruses (AAV), is crucial in efficient gene delivery to cochlear cells. The ChR2 gene is inserted into the viral vector through molecular cloning techniques, and the resulting viral vector is introduced into cochlear cells via direct injection or round window membrane delivery. This allows for the expression of ChR2 and subsequent light sensitivity in targeted cells. Alternatively, direct cell transfection offers a non-viral approach for ChR2 delivery. The ChR2 gene is cloned into a plasmid vector, which is then combined with transfection agents like liposomes or nanoparticles. This mixture is applied to cochlear cells, facilitating the entry of the plasmid DNA into the target cells and enabling ChR2 expression. Optogenetic stimulation using ChR2 allows for precise and selective activation of specific neurons in response to light, potentially overcoming the limitations of current auditory prostheses. Moreover, optogenetics has broader implications in understanding the neural circuits involved in auditory processing and behavior. The combination of optogenetics and gene delivery techniques provides a promising avenue for improving hearing restoration strategies, offering the potential for enhanced frequency resolution, spatial selectivity, and improved auditory perception.

Simulation of the ligand-leaving process of the human heat shock protein

The human heat shock protein plays a critical role in various diseases and is an important target for pharmacological modulation. Simulation of conformational changes and free energy profiles of the human heat shock protein derived by the ligand-leaving process is a challenging issue. In this work, steered molecular dynamics simulation was adopted to simulate the ligand-leaving process. Two composite systems of heat shock protein NHSP90 and small molecules 6FJ and 6G7 are selected as research objects. The free energy during the leaving of ligand small molecules is calculated using conventional molecular dynamics simulation, steered molecular dynamics simulation (SMD), and the umbrella sampling method. We found that the a slower pulling velocity (0.001 nm ns-1) will result in 2.19 kcal mol-1, and the umbrella sampling method gives a value of 3.26 kcal mol-1 for the free energy difference for the two systems, which reasonably agrees with experimental results. A faster-pulling velocity (0.01 nm ns-1) leads to a large overestimation of free energy. At the same time, the conformational analysis indicated that the faster pulling velocity may lead to the conformational change of NHSP90, which was proved to be false by the slower pulling velocity and the umbrella sampling method.

The Mechanism of the Channel Opening in Channelrhodopsin-2: A Molecular Dynamics Simulation

Channelrhodopsin-2 (ChR2) has been one of the most important objects in the study of optogenetics. The retinal chromophore molecule absorbs photons and undergoes an isomerization reaction, which triggers the photocycle, resulting in a series of conformational changes. In this study, a series of intermediate structures (including D470, P500, P390-early, P390-late, and P520 states) of ChR2 in the photocycle were modeled, and molecular dynamics (MD) simulations were performed to elucidate the mechanism of ion channel opening of ChR2. The maximum absorption wavelength of these intermediates calculated by time-dependent density function theory (TD-DFT) is in general agreement with the experimental values, the distribution of water density gradually increases in the process of photocycle, and the radius of the ion channel is larger than 6 Å. All these results indicate that our structural models of the intermediates are reasonable. The evolution of protonation state of E90 during the photocycle is explained. E90 will deprotonate when the P390-early transforms into P390-late, in which the two conformations of P390-early and P390-late obtained from the simulations are consistent with the experimental descriptions. To validate the conductive P520 state, the potential mean force (PMF) of Na+ ions passing through the P520 intermediate was calculated by using steered molecular dynamics (SMD) simulation combined with umbrella sampling. The result shows that the Na+ ions passing through the channel with a very low energy barrier, especially in the central gate, is almost barrierless. This indicates that the channel is open in the P520 state.

Mechanistic basis of the increased methylation activity of the SETD2 protein lysine methyltransferase towards a designed super-substrate peptide

Protein lysine methyltransferases have important regulatory functions in cells, but mechanisms determining their activity and specificity are incompletely understood. Naturally, SETD2 introduces H3K36me3, but previously an artificial super-substrate (ssK36) was identified, which is methylated >100-fold faster. The ssK36-SETD2 complex structure cannot fully explain this effect. We applied molecular dynamics (MD) simulations and biochemical experiments to unravel the mechanistic basis of the increased methylation of ssK36, considering peptide conformations in solution, association of peptide and enzyme, and formation of transition-state (TS) like conformations of the enzyme-peptide complex. We observed in MD and FRET experiments that ssK36 adopts a hairpin conformation in solution with V35 and K36 placed in the loop. The hairpin conformation has easier access into the active site of SETD2 and it unfolds during the association process. Peptide methylation experiments revealed that introducing a stable hairpin conformation in the H3K36 peptide increased its methylation by SETD2. In MD simulations of enzyme-peptide complexes, the ssK36 peptide approached TS-like structures more frequently than H3K36 and distinct, substrate-specific TS-like structures were observed. Hairpin association, hairpin unfolding during association, and substrate-specific catalytically competent conformations may also be relevant for other PKMTs and hairpins could represent a promising starting point for SETD2 inhibitor development.

Time-resolved photoacoustics of channelrhodopsins: early energetics and light-driven volume changes

In biological photoreceptors, the energy stored in early transient species is a key feature to drive the photocycle or a chain of reactions. Time-resolved photoacoustics (PA) can explore the energy landscape of transient species formed within few ns after photoexcitation, as well as volumetric changes (ΔV) of these intermediates with respect to the parental state. In this work, PA identified these important parameters for several channelrhodopsins, namely CaChR1 from Chlamydomonas augustae and CrChR2 from Chlamydomonas reinhardtii and various variants. PA has access to the sub-ns formation of the early photoproduct P1 and to its relaxation, provided that this latter process occurs within a few μs. We found that ΔVP1 for CaChR1 is ca. 12 mL/mol, while it is much smaller for CrChR2 (4.7 mL/mol) and for H. salinarum bacteriorhodopsin (HsBR, ΔVK = 2.8 mL/mol). PA experiments on variants strongly indicate that part of this large ΔVP1 value for CaChR1 is caused by the protonation dynamics of the Schiff base counterion complex involving E169 and D299. PA data further show that the energy level of P1 is higher in CrChR2 (ca. 96 kJ/mol) than in CaChr1 (ca. 46 kJ/mol), comparable to the energy level of the K state of HsBR (60 kJ/mol). Instrumental to gain these molecular values from the raw PA data was the estimation of the quantum yield (Φ) for P1 formation via transient spectroscopy; for both channelrhodopsins, ΦP2 was evaluated as ca. 0.4.

Graphical Abstract

Modeling the syn-cycle in the light activated opening of the channelrhodopsin-2 ion channel

The ion channel of channelrhodopsin-2 (ChR2) is activated by absorbing light. The light stimulates retinal to isomerize to start the photocycle. There are two pathways for photocycles, which are caused by isomerization of the retinal from all-trans, 15-anti to 13-cis, 15-anti in the dark-adapted state (anti-cycle) and from 13-cis, 15-syn to all-trans, 15-syn in the light-adapted state (syn-cycle). In this work, the structure of the syn-cycle intermediate and mechanism of channel opening were studied by molecular dynamics (MD) and steered molecular dynamics (SMD) simulations. Due to the lack of crystal structure of intermediates in the syn-cycle of ChR2, the intermediate models were constructed from the homologous intermediates in the anti-cycle. The isomerization of retinal was shown to cause the central gate (CG) hydrogen bond network to rearrange, cutting the link between TM2 and TM7. TM2 is moved by the intrahelical hydrogen bond of E90 and K93, and induced the intracellular gate (ICG) to expand. The ion penetration pathway between TM1, TM2, TM3 and TM7 in the P500* state was observed by MD simulations. However, this channel is not fully opened compared with the homologous P500 state in the anti-cycle. In addition, the protons on Schiff bases were found to be unable to form hydrogen bonds with the counter residues (E123 and D253) in the P500* state, preventing an evolution of the P500* state to a P390-like state in the syn-cycle.

Modelling the syn-cycle is a series of operations on the ChR2 crystal structure (PDB ID: 6EID). By replacement and isomerization, we obtained P500* and P480 intermediates. A feasible explanation that no P390* was observed in experiment was inferred.

Molecular Dynamics Simulations in Designing DARPins as Phosphorylation-Specific Protein Binders of ERK2

Extracellular signal-regulated kinases 1 and 2 (ERK1/2) play key roles in promoting cell survival and proliferation through the phosphorylation of various substrates. Remarkable antitumour activity is found in many inhibitors that act upstream of the ERK pathway. However, drug-resistant tumour cells invariably emerge after their use due to the reactivation of ERK1/2 signalling. ERK1/2 inhibitors have shown clinical efficacy as a therapeutic strategy for the treatment of tumours with mitogen-activated protein kinase (MAPK) upstream target mutations. These inhibitors may be used as a possible strategy to overcome acquired resistance to MAPK inhibitors. Here, we report a class of repeat proteins-designed ankyrin repeat protein (DARPin) macromolecules targeting ERK2 as inhibitors. The structural basis of ERK2-DARPin interactions based on molecular dynamics (MD) simulations was studied. The information was then used to predict stabilizing mutations employing a web-based algorithm, MAESTRO. To evaluate whether these design strategies were successfully deployed, we performed all-atom, explicit-solvent molecular dynamics (MD) simulations. Two mutations, Ala → Asp and Ser → Leu, were found to perform better than the original sequence (DARPin E40) based on the associated energy and key residues involved in protein-protein interaction. MD simulations and analysis of the data obtained on these mutations supported our predictions.

The effect on ion channel of different protonation states of E90 in channelrhodopsin-2: a molecular dynamics simulation

Channelrhodopsin-2 (ChR2) is a cationic channel protein that has been extensively studied in optogenetics. The ion channel is opened via a series of proton transfers and H-bond changes during the photocycle but the detailed mechanism is still unknown. Molecular dynamics (MD) simulations with enhanced sampling were performed on the dark-adapted state (i.e., D470) and two photocycle intermediates (P₁⁵⁰⁰ and P₂³⁹⁰) to study the proton transfer path of the Schiff base and the subsequent conformational changes. The results suggest there are two possible proton transfer pathways from the Schiff base to proton acceptors (i.e., E123 or D253), depending on the protonation of E90. If E90 is protonated in the P₁⁵⁰⁰ state, the proton on the Schiff base will transfer to E123. The polyene chain of 13-cis retinal tilts and opens the channel that detours the blocking central gate (CG) and forms a narrow channel through the transmembrane helices (TM) 2, 3, 6 and 7. In contrast, if E90 deprotonates after retinal isomerization, the primary proton acceptor is D253, and an almost-open channel through TM1, 2, 3 and 7 is generated. The channel diameter is very close to the experimental value. The potential mean force (PMF) suggests that the free energy is extremely low for ions passing through this channel.

With E90 protonated, the proton acceptor of RSBH⁺ is E123 with a narrow channel along TM3; while with E90 deprotonated, proton transfer from RSBH⁺ to D253 generates an approximately open channel along TM2.

Molecular Dynamics Simulation of Transmembrane Transport of Chloride Ions in Mutants of Channelrhodopsin

Channelrhodopsins (ChRs) are light-gated transmembrane cation channels which are widely used for optogenetic technology. Replacing glutamate located at the central gate of the ion channel with positively charged amino acid residues will reverse ion selectivity and allow anion conduction. The structures and properties of the ion channel, the transport of chloride, and potential of mean force (PMF) of the chimera protein (C1C2) and its mutants, EK-TC, ER-TC and iChloC, were investigated by molecular dynamics simulation. The results show that the five-fold mutation in E122Q-E129R-E140S-D195N-T198C (iChloC) increases the flexibility of the transmembrane channel protein better than the double mutations in EK-TC and ER-TC, and results in an expanded ion channel pore size and decreased steric resistance. The iChloC mutant was also found to have a higher affinity for chloride ions and, based on surface electrostatic potential analysis, provides a favorable electrostatic environment for anion conduction. The PMF free energy curves revealed that high affinity Cl^- binding sites are generated near the central gate of the three mutant proteins. The energy barriers for the EK-TC and ER-TC were found to be much higher than that of iChloC. The results suggest that the transmembrane ion channel of iChloC protein is better at facilitating the capture and transport of chloride ions.