Hybrid QM/MM Free-Energy Evaluation of Drug-Resistant Mutational Effect on the Binding of an Inhibitor Indinavir to HIV-1 Protease.

A human immunodeficiency virus-1 (HIV-1) protease is a homodimeric aspartic protease essential for the replication of HIV. The HIV-1 protease is a target protein in drug discovery for antiretroviral therapy, and various inhibitor molecules of transition state analogues have been developed. However, serious drug-resistant mutants have emerged. For understanding the molecular mechanism of the drug resistance, an accurate examination of the impacts of the mutations on ligand binding and enzymatic activity is necessary. Here, we present a molecular simulation study on the ligand binding of indinavir, a potent transition state analogue inhibitor, to the wild-type protein and a V82T/I84V drug-resistant mutant of the HIV-1 protease. We employed a hybrid ab initio quantum mechanical/molecular mechanical (QM/MM) free-energy optimization technique which combines a highly accurate QM description of the ligand molecule and its interaction with statistically ample conformational sampling of the MM protein environment by long-time molecular dynamics simulations. Through the free-energy calculations of protonation states of catalytic groups at the binding pocket and of the ligand-binding affinity changes upon the mutations, we successfully reproduced the experimentally observed significant reduction of the binding affinity upon the drug-resistant mutations and elucidated the underlying molecular mechanism. The present study opens the way for understanding the molecular mechanism of drug resistance through the direct quantitative comparison of ligand binding and enzymatic reaction with the same accuracy.

A computational method to simulate global conformational changes of proteins induced by cosolvent.

A computational method to investigate the global conformational change of a protein is proposed by combining the linear response path following (LRPF) method and three-dimensional reference interaction site model (3D-RISM) theory, which is referred to as the LRPF/3D-RISM method. The proposed method makes it possible to efficiently simulate protein conformational changes caused by either solutions of varying concentrations or the presence of cosolvent species by taking advantage of the LRPF and 3D-RISM. The proposed method is applied to the urea-induced denaturation of ubiquitin. The LRPF/3D-RISM trajectories successfully simulate the early stage of the denaturation process within the simulation time of 300 ns, whereas no significant structural change is observed even in the 1 µs standard MD simulation. The obtained LRPF/3D-RISM trajectories reproduce the mechanism of the urea denaturation of ubiquitin reported in previous studies, and demonstrate the high efficiency of the method.

An Atomistic Model of a Precursor State of Light-Induced Channel Opening of Channelrhodopsin.

Channelrhodopsins (ChRs) are microbial light-gated ion channels with a retinal chromophore and are widely utilized in optogenetics to precisely control neuronal activity with light. Despite increasing understanding of their structures and photoactivation kinetics, the atomistic mechanism of light gating and ion conduction remains elusive. Here, we present an atomic structural model of a chimeric ChR in a precursor state of the channel opening determined by an accurate hybrid molecular simulation technique and a statistical theory of internal water distribution. The photoactivated structure features extensive tilt of the chromophore accompanied by redistribution of water molecules in its binding pocket, which is absent in previously known photoactivated structures of analogous photoreceptors, and widely agrees with structural and spectroscopic experimental evidence of ChRs. The atomistic model manifests a photoactivated ion-conduction pathway that is markedly different from a previously proposed one and successfully explains experimentally observed mutagenic effects on key channel properties.

QM/MM Geometry Optimization on Extensive Free-Energy Surfaces for Examination of Enzymatic Reactions and Design of Novel Functional Properties of Proteins.

Many remarkable molecular functions of proteins use their characteristic global and slow conformational dynamics through coupling of local chemical states in reaction centers with global conformational changes of proteins. To theoretically examine the functional processes of proteins in atomic detail, a methodology of quantum mechanical/molecular mechanical (QM/MM) free-energy geometry optimization is introduced. In the methodology, a geometry optimization of a local reaction center is performed with a quantum mechanical calculation on a free-energy surface constructed with conformational samples of the surrounding protein environment obtained by a molecular dynamics simulation with a molecular mechanics force field. Geometry optimizations on extensive free-energy surfaces by a QM/MM reweighting free-energy self-consistent field method designed to be variationally consistent and computationally efficient have enabled examinations of the multiscale molecular coupling of local chemical states with global protein conformational changes in functional processes and analysis and design of protein mutants with novel functional properties.

Crystal structure of the channelrhodopsin light-gated cation channel.

Channelrhodopsins (ChRs) are light-gated cation channels derived from algae that have shown experimental utility in optogenetics; for example, neurons expressing ChRs can be optically controlled with high temporal precision within systems as complex as freely moving mammals. Although ChRs have been broadly applied to neuroscience research, little is known about the molecular mechanisms by which these unusual and powerful proteins operate. Here we present the crystal structure of a ChR (a C1C2 chimaera between ChR1 and ChR2 from Chlamydomonas reinhardtii) at 2.3 Å resolution. The structure reveals the essential molecular architecture of ChRs, including the retinal-binding pocket and cation conduction pathway. This integration of structural and electrophysiological analyses provides insight into the molecular basis for the remarkable function of ChRs, and paves the way for the precise and principled design of ChR variants with novel properties.

X-ray Crystallographic Structure of Thermophilic Rhodopsin: IMPLICATIONS FOR HIGH THERMAL STABILITY AND OPTOGENETIC FUNCTION.

Thermophilic rhodopsin (TR) is a photoreceptor protein with an extremely high thermal stability and the first characterized light-driven electrogenic proton pump derived from the extreme thermophile Thermus thermophilus JL-18. In this study, we confirmed its high thermal stability compared with other microbial rhodopsins and also report the potential availability of TR for optogenetics as a light-induced neural silencer. The x-ray crystal structure of TR revealed that its overall structure is quite similar to that of xanthorhodopsin, including the presence of a putative binding site for a carotenoid antenna; but several distinct structural characteristics of TR, including a decreased surface charge and a larger number of hydrophobic residues and aromatic-aromatic interactions, were also clarified. Based on the crystal structure, the structural changes of TR upon thermal stimulation were investigated by molecular dynamics simulations. The simulations revealed the presence of a thermally induced structural substate in which an increase of hydrophobic interactions in the extracellular domain, the movement of extracellular domains, the formation of a hydrogen bond, and the tilting of transmembrane helices were observed. From the computational and mutational analysis, we propose that an extracellular LPGG motif between helices F and G plays an important role in the thermal stability, acting as a "thermal sensor." These findings will be valuable for understanding retinal proteins with regard to high protein stability and high optogenetic performance.

Key chemical factors of arginine finger catalysis of F1-ATPase clarified by an unnatural amino acid mutation.

A catalytically important arginine, called Arg finger, is employed in many enzymes to regulate their functions through enzymatic hydrolysis of nucleotide triphosphates. F1-ATPase (F1), a rotary motor protein, possesses Arg fingers which catalyze hydrolysis of adenosine triphosphate (ATP) for efficient chemomechanical energy conversion. In this study, we examined the Arg finger catalysis by single-molecule measurements for a mutant of F1 in which the Arg finger is substituted with an unnatural amino acid of a lysine analogue, 2,7-diaminoheptanoic acid (Lyk). The use of Lyk, of which the side chain is elongated by one CH2 unit so that its chain length to the terminal nitrogen of amine is set to be equal to that of arginine, allowed us to resolve key chemical factors in the Arg finger catalysis, i.e., chain length matching and chemical properties of the terminal groups. Rate measurements by single-molecule observations showed that the chain length matching of the side-chain length is not a sole requirement for the Arg finger to catalyze the ATP hydrolysis reaction step, indicating the crucial importance of chemical properties of the terminal guanidinium group in the Arg finger catalysis. On the other hand, the Lyk mutation prevented severe formation of an ADP inhibited state observed for a lysine mutant and even improved the avoidance of inhibition compared with the wild-type F1. The present study demonstrated that incorporation of unnatural amino acids can widely extend with its high "chemical" resolution biochemical approaches for elucidation of the molecular mechanism of protein functions and furnishing novel characteristics.

Molecular Mechanism of Wide Photoabsorption Spectral Shifts of Color Variants of Human Cellular Retinol Binding Protein II.

Color variants of human cellular retinol binding protein II (hCRBPII) created by protein engineering were recently shown to exhibit anomalously wide photoabsorption spectral shifts over ∼200 nm across the visible region. The remarkable phenomenon provides a unique opportunity to gain insight into the molecular basis of the color tuning of retinal binding proteins for understanding of color vision as well as for engineering of novel color variants of retinal binding photoreceptor proteins employed in optogenetics. Here, we report a theoretical investigation of the molecular mechanism underlying the anomalously wide spectral shifts of the color variants of hCRBPII. Computational modeling of the color variants with hybrid molecular simulations of free energy geometry optimization succeeded in reproducing the experimentally observed wide spectral shifts, and revealed that protein flexibility, through which the active site structure of the protein and bound water molecules is altered by remote mutations, plays a significant role in inducing the large spectral shifts.

Converting a light-driven proton pump into a light-gated proton channel.

There are two types of membrane-embedded ion transport machineries in nature. The ion pumps generate electrochemical potential by energy-coupled active ion transportation, while the ion channels produce action potential by stimulus-dependent passive ion transportation. About 80% of the amino acid residues of the light-driven proton pump archaerhodopsin-3 (AR3) and the light-gated cation channel channelrhodopsin (ChR) differ although they share the close similarity in architecture. Therefore, the question arises: How can these proteins function differently? The absorption maxima of ion pumps are red-shifted about 30-100 nm compared with ChRs, implying a structural difference in the retinal binding cavity. To modify the cavity, a blue-shifted AR3 named AR3-T was produced by replacing three residues located around the retinal (i.e., M128A, G132V, and A225T). AR3-T showed an inward H(+) flux across the membrane, raising the possibility that it works as an inward H(+) pump or an H(+) channel. Electrophysiological experiments showed that the reverse membrane potential was nearly zero, indicating light-gated ion channeling activity of AR3-T. Spectroscopic characterization of AR3-T revealed similar photochemical properties to some of ChRs, including an all-trans retinal configuration, a strong hydrogen bond between the protonated retinal Schiff base and its counterion, and a slow photocycle. From these results, we concluded that the functional determinant in the H(+) transporters is localized at the center of the membrane-spanning domain, but not in the cytoplasmic and extracellular domains.

A blue-shifted light-driven proton pump for neural silencing.

Ion-transporting rhodopsins are widely utilized as optogenetic tools both for light-induced neural activation and silencing. The most studied representative is Bacteriorhodopsin (BR), which absorbs green/red light (∼570 nm) and functions as a proton pump. Upon photoexcitation, BR induces a hyperpolarization across the membrane, which, if incorporated into a nerve cell, results in its neural silencing. In this study, we show that several residues around the retinal chromophore, which are completely conserved among BR homologs from the archaea, are involved in the spectral tuning in a BR homolog (HwBR) and that the combination mutation causes a large spectral blue shift (λmax = 498 nm) while preserving the robust pumping activity. Quantum mechanics/molecular mechanics calculations revealed that, compared with the wild type, the ß-ionone ring of the chromophore in the mutant is rotated ∼130° because of the lack of steric hindrance between the methyl groups of the retinal and the mutated residues, resulting in the breakage of the π conjugation system on the polyene chain of the retinal. By the same mutations, similar spectral blue shifts are also observed in another BR homolog, archearhodopsin-3 (also called Arch). The color variant of archearhodopsin-3 could be successfully expressed in the neural cells of Caenorhabditis elegans, and illumination with blue light (500 nm) led to the effective locomotory paralysis of the worms. Thus, we successfully produced a blue-shifted proton pump for neural silencing.

Adenosine triphosphate hydrolysis mechanism in kinesin studied by combined quantum-mechanical/molecular-mechanical metadynamics simulations.

Kinesin is a molecular motor that hydrolyzes adenosine triphosphate (ATP) and moves along microtubules against load. While motility and atomic structures have been well-characterized for various members of the kinesin family, not much is known about ATP hydrolysis inside the active site. Here, we study ATP hydrolysis mechanisms in the kinesin-5 protein Eg5 by using combined quantum mechanics/molecular mechanics metadynamics simulations. Approximately 200 atoms at the catalytic site are treated by a dispersion-corrected density functional and, in total, 13 metadynamics simulations are performed with their cumulative time reaching ~0.7 ns. Using the converged runs, we compute free energy surfaces and obtain a few hydrolysis pathways. The pathway with the lowest free energy barrier involves a two-water chain and is initiated by the Pγ-Oß dissociation concerted with approach of the lytic water to PγO3-. This immediately induces a proton transfer from the lytic water to another water, which then gives a proton to the conserved Glu270. Later, the proton is transferred back from Glu270 to HPO(4)2- via another hydrogen-bonded chain. We find that the reaction is favorable when the salt bridge between Glu270 in switch II and Arg234 in switch I is transiently broken, which facilitates the ability of Glu270 to accept a proton. When ATP is placed in the ADP-bound conformation of Eg5, the ATP-Mg moiety is surrounded by many water molecules and Thr107 blocks the water chain, which together make the hydrolysis reaction less favorable. The observed two-water chain mechanisms are rather similar to those suggested in two other motors, myosin and F1-ATPase, raising the possibility of a common mechanism.

Facile catch and release of fullerenes using a photoresponsive molecular tube.

A novel M2L2 molecular tube capable of binding fullerene C60 was synthesized from bispyridine ligands with embedded anthracene panels and Ag(I) hinges. Unlike previous molecular cages and capsules, this open-ended tubular host can accommodate a single molecule of various C60 derivatives with large substituents. The fullerene guest can then be released by using the ideal, noninvasive external stimulus, light.

The reaction mechanism of Claisen rearrangement obtained by transition state spectroscopy and single direct-dynamics trajectory.

Chemical bond breaking and formation during chemical reactions can be observed using "transition state spectroscopy". Comparing the measurement result of the transition state spectroscopy with the simulation result of single direct-dynamics trajectory, we have elucidated the reaction dynamics of Claisen rearrangement of allyl vinyl ether. Observed the reaction of the neat sample liquid, we have estimated the time constants of transformation from straight-chain structure to aromatic-like six-membered ring structure forming the C¹-C6 bond. The result clarifies that the reaction proceeds via three steps taking longer time than expected from the gas phase calculation. This finding provides new hypothesis and discussions, helping the development of the field of reaction mechanism analysis.

Liquid Structures and Diffusion Dynamics of Alkyl-Pyrene Liquids Studied by Molecular Dynamics Simulations.

Functional molecular liquids (FMLs) based on alkylated π-conjugated molecules have attracted attention as solvent-free and nonvolatile liquid materials with prominent optoelectronic features. Recently, novel FML compounds containing pyrene as the functional core were synthesized, and their rheological and photochemical properties were investigated. Although the molecules differ only in the number of alkyl chain substituents and their substitution positions, their viscosity coefficients are largely different beyond the Stokes-Einstein relation on the assumption of identical microscopic friction, indicating that local microscopic molecular interactions are crucial for the macroscopic rheological properties. Here, we report a theoretical study on the rheological properties of the alkyl-pyrene liquids by means of atomistic molecular dynamics (MD) simulations. We performed long-time MD simulations for tens of microseconds to obtain ample statistical samples of the alkyl-pyrene liquids and analyzed their liquid structures and diffusion dynamics based on spatiotemporal correlation functions. We found the formation of characteristic local liquid structures of π-π stacking of the pyrene moieties and locally anisotropic and anomalous diffusion dynamics, which remarkably vary depending on the alkyl substituent patterns. The present results provide an atomistic insight into the macroscopic rheological properties of alkyl-π FMLs and molecular design strategy for them.

Unique Vibrational Characteristics and Structures of the Photoexcited Retinal Chromophore in Ion-Pumping Rhodopsins.

Photoisomerization of an all-trans-retinal chromophore triggers ion transport in microbial ion-pumping rhodopsins. Understanding chromophore structures in the electronically excited (S1) state provides insights into the structural evolution on the potential energy surface of the photoexcited state. In this study, we examined the structure of the S1-state chromophore in Natronomonas pharaonis halorhodopsin (NpHR), a chloride ion-pumping rhodopsin, using time-resolved resonance Raman spectroscopy. The spectral patterns of the S1-state chromophore were completely different from those of the ground-state chromophore, resulting from unique vibrational characteristics and the structure of the S1 state. Mode assignments were based on a combination of deuteration shifts of the Raman bands and hybrid quantum mechanics-molecular mechanics calculations. The present observations suggest a weakened bond alternation in the π conjugation system. A strong hydrogen-out-of-plane bending band was observed in the Raman spectra of the S1-state chromophore in NpHR, indicating a twisted polyene structure. Similar frequency shifts for the C═N/C═C and C-C stretching modes of the S1-state chromophore in NpHR were observed in the Raman spectra of sodium ion-pumping and proton-pumping rhodopsins, suggesting that these unique features are common to the S1 states of ion-pumping rhodopsins.

Crucial role of protein flexibility in formation of a stable reaction transition state in an α-amylase catalysis.

Conformational flexibility of proteins provides enzymes with high catalytic activity. Although the conformational flexibility is known to be pivotal for the ligand binding and release, its role in the chemical reaction process of the reactive substrate remains unclear. We determined a transition state of an enzymatic reaction in a psychrophilic α-amylase by a hybrid molecular simulation that allows one to identify the optimal chemical state in an extensive conformational ensemble of protein. The molecular simulation uncovered that formation of the reaction transition state accompanies a large and slow movement of a loop adjacent to the catalytic site. Free energy calculations revealed that, although catalytic electrostatic potentials on the reactive moiety are formed by local and fast reorganization around the catalytic site, reorganization of the large and slow movement of the loop significantly contributes to reduction of the free energy barrier by stabilizing the local reorganization.

Molecular mechanism of ATP hydrolysis in F1-ATPase revealed by molecular simulations and single-molecule observations.

Enzymatic hydrolysis of nucleotide triphosphate (NTP) plays a pivotal role in protein functions. In spite of its biological significance, however, the chemistry of the hydrolysis catalysis remains obscure because of the complex nature of the reaction. Here we report a study of the molecular mechanism of hydrolysis of adenosine triphosphate (ATP) in F(1)-ATPase, an ATP-driven rotary motor protein. Molecular simulations predicted and single-molecule observation experiments verified that the rate-determining step (RDS) is proton transfer (PT) from the lytic water molecule, which is strongly activated by a metaphosphate generated by a preceding P(γ)-O(ß) bond dissociation (POD). Catalysis of the POD that triggers the chain activation of the PT is fulfilled by hydrogen bonds between Walker motif A and an arginine finger, which commonly exist in many NTPases. The reaction mechanism unveiled here indicates that the protein can regulate the enzymatic activity for the function in both the POD and PT steps despite the fact that the RDS is the PT step.

Cryo-EM structures of thermostabilized prestin provide mechanistic insights underlying outer hair cell electromotility.

Outer hair cell elecromotility, driven by prestin, is essential for mammalian cochlear amplification. Here, we report the cryo-EM structures of thermostabilized prestin (PresTS), complexed with chloride, sulfate, or salicylate at 3.52-3.63 Å resolutions. The central positively-charged cavity allows flexible binding of various anion species, which likely accounts for the known distinct modulations of nonlinear capacitance (NLC) by different anions. Comparisons of these PresTS structures with recent prestin structures suggest rigid-body movement between the core and gate domains, and provide mechanistic insights into prestin inhibition by salicylate. Mutations at the dimeric interface severely diminished NLC, suggesting that stabilization of the gate domain facilitates core domain movement, thereby contributing to the expression of NLC. These findings advance our understanding of the molecular mechanism underlying mammalian cochlear amplification.

A theoretical study on excited state double proton transfer reaction of a 7-azaindole dimer: an ab initio potential energy surface and its empirical valence bond model.

Double proton transfer (DPT) reaction of a 7-azaindole dimer in the first ππ* electronically excited state was studied theoretically. We investigated the reaction mechanism through constructing a full dimensional empirical valence bond potential energy function (PEF) based on potential energies evaluated by ab initio molecular orbital methods, and carrying out quantum dynamics calculations with the PEF. Potential energy surfaces of the DPT obtained at the multi-reference perturbation level of theory favors a concerted DPT mechanism, although a stepwise channel is suggested to open for an excited initial vibrational state. Reduced two dimensional quantum dynamics calculations for a reaction surface Hamiltonian of DPT coordinates were performed. Time constants of the reaction were evaluated to be on the order of picoseconds, which is consistent with experiments. On the other hand, the computed kinetic isotope effect deviates from experimental evidence, suggesting the importance of intermolecular stretching motion, which is not explicit in the present calculations for the quantum effect.

Ab Initio Evaluation of the Redox Potential of Cytochrome c.

Various biochemical activities of metabolism and biosynthesis are fulfilled by redox processes with explicit electron exchange, which furnish redox enzymes with high chemical reactivity. However, theoretical investigation of a redox process, which simultaneously involves a complex electronic change at a redox metal center and conformational reorganization of the surrounding protein environment coupled to the electronic change, requires computationally conflicting approaches, highly accurate quantum chemical calculations, and long-time molecular dynamics (MD) simulations, limiting the physicochemical understanding of biological redox processes. Here, we theoretically examined a redox process of cytochrome c by means of a hybrid molecular simulation technique, which enables one to consistently treat the redox center at the ab initio quantum chemistry level of theory and the protein reorganization with long-time MD simulations on the microsecond timescale. The calculations successfully evaluated a large absolute redox potential, 4.34 eV, with errors of only 0.03 to 0.34 eV to the experimental ones without any problem-specific empirical parameters. Through the long-time MD sampling, large and nonlinear reorganization of the protein environment was unveiled and the molecular determinants for the redox potential were identified. The present ab initio approach significantly expands the applicability of theoretical investigation to biological redox systems with more electronically complicated redox centers such as polynuclear transition metal complexes.