Molecular Dynamics of Channelrhodopsin at the Early Stages of Channel Opening.

Channelrhodopsin (ChR) is a light-gated cation channel that responds to blue light. Since ChR can be readily expressed in specific neurons to precisely control their activities by light, it has become a powerful tool in neuroscience. Although the recently solved crystal structure of a chimeric ChR, C1C2, provided the structural basis for ChR, our understanding of the molecular mechanism of ChR still remains limited. Here we performed electrophysiological analyses and all-atom molecular dynamics (MD) simulations, to investigate the importance of the intracellular and central constrictions of the ion conducting pore observed in the crystal structure of C1C2. Our electrophysiological analysis revealed that two glutamate residues, Glu122 and Glu129, in the intracellular and central constrictions, respectively, should be deprotonated in the photocycle. The simulation results suggested that the deprotonation of Glu129 in the central constriction leads to ion leakage in the ground state, and implied that the protonation of Glu129 is important for preventing ion leakage in the ground state. Moreover, we modeled the 13-cis retinal bound; i.e., activated C1C2, and performed MD simulations to investigate the conformational changes in the early stage of the photocycle. Our simulations suggested that retinal photoisomerization induces the conformational change toward channel opening, including the movements of TM6, TM7 and TM2. These insights into the dynamics of the ground states and the early photocycle stages enhance our understanding of the channel function of ChR.

Structural basis for dynamic mechanism of proton-coupled symport by the peptide transporter POT.

Proton-dependent oligopeptide transporters (POTs) are major facilitator superfamily (MFS) proteins that mediate the uptake of peptides and peptide-like molecules, using the inwardly directed H(+) gradient across the membrane. The human POT family transporter peptide transporter 1 is present in the brush border membrane of the small intestine and is involved in the uptake of nutrient peptides and drug molecules such as ß-lactam antibiotics. Although previous studies have provided insight into the overall structure of the POT family transporters, the question of how transport is coupled to both peptide and H(+) binding remains unanswered. Here we report the high-resolution crystal structures of a bacterial POT family transporter, including its complex with a dipeptide analog, alafosfalin. These structures revealed the key mechanistic and functional roles for a conserved glutamate residue (Glu310) in the peptide binding site. Integrated structural, biochemical, and computational analyses suggested a mechanism for H(+)-coupled peptide symport in which protonated Glu310 first binds the carboxyl group of the peptide substrate. The deprotonation of Glu310 in the inward open state triggers the release of the bound peptide toward the intracellular space and salt bridge formation between Glu310 and Arg43 to induce the state transition to the occluded conformation.

Molecular dynamics simulation of Autotaxin: roles of the nuclease-like domain and the glycan modification.

Autotaxin (ATX) is a secreted lysophospholipase D that produces lysophosphatidic acid, a lipid mediator that activates G protein-coupled receptors to evoke various cellular responses. The nuclease-like domain of ATX and the Asn524-linked glycan are reportedly critical for the catalytic activity. Recently, the crystal structures of ATX were determined, but the means by which the nuclease-like domain and the N-glycosylation participate in the catalytic activity still remain undetermined. To address this question, we conducted molecular dynamics simulations of ATX. The simulation trajectories starting from the full-length structure and from structures lacking the nuclease-like domain and/or the glycan were compared. The results suggested that an allosteric interaction pathway, formed by the catalytic domain, including the two insertion regions, the essential glycan modification, and the nuclease-like domain, may stabilize the proper location of the catalytic threonine residue. The results complement those from previous biochemical experiments.

Quantitative Prediction of Regioselectivity Toward Cytochrome P450/3A4 Using Machine Learning Approaches.

In the drug discovery process, it is important to know the properties of both drug candidates and their metabolites. Fast and precise prediction of metabolites is essential. However, it has been difficult to predict metabolites because of the complexity of the mechanism of cytochrome P450/3A4 (CYP 3A4), which is the main metabolite enzyme of drugs. In this study, we focus on the regioselectivity of CYP 3A4, i.e., the selectivity of metabolic sites. We have developed a model to predict the regioselectivity of drug candidates by using machine learning (ML) approaches.