Search

Multistate Bennett acceptance ratio (MBAR) works as a method to analyze molecular dynamics (MD) simulation data after the simulations have been finished. It is widely used to estimate free-energy changes between different states and averaged properties at the states of interest. MBAR allows us to treat a wide range of states from those at different temperature/pressure to those with different model parameters. Due to the broad applicability, the MBAR equations are rather difficult to apply for free-energy calculations using different types of MD simulations including enhanced conformational sampling methods and free-energy perturbation. In this review, we first summarize the basic theory of the MBAR equations and categorize the representative usages into the following four: (i) perturbation, (ii) scaling, (iii) accumulation, and (iv) full potential energy. For each, we explain how to prepare input data using MD simulation trajectories for solving the MBAR equations. MBAR is also useful to estimate reliable free-energy differences using MD trajectories based on a semi-empirical quantum mechanics/molecular mechanics (QM/MM) model and ab initio QM/MM energy calculations on the MD snapshots. We also explain how to use the MBAR software in the GENESIS package, which we call mbar_analysis, for the four representative cases. The proposed estimations of free-energy changes and thermodynamic averages are effective and useful for various biomolecular systems.

Fibroblast growth factor receptor 3 (FGFR3) is a member of receptor tyrosine kinases, which is involved in skeletal cell growth, differentiation, and migration. FGFR3 transduces biochemical signals from the extracellular ligand-binding domain to the intracellular kinase domain through the conformational changes of the transmembrane (TM) helix dimer. Here, we apply generalized replica exchange with solute tempering method to wild type (WT) and G380R mutant (G380R) of FGFR3. The dimer interface in G380R is different from WT and the simulation results are in good agreement with the solid-state nuclear magnetic resonance (NMR) spectroscopy. TM helices in G380R are extended more than WT, and thereby, G375 in G380R contacts near the N-termini of the TM helix dimer. Considering that both G380R and G375C show the constitutive activation, the formation of the N-terminal contacts of the TM helices can be generally important for the activation mechanism. © 2019 Wiley Periodicals, Inc.

Significance Molecular recognition via protein–ligand binding is essential for biomolecular functions. Current understanding of ligand recognition is mostly brought from a bound state structure resolved by X-ray crystallography. Using molecular dynamics simulation with enhanced sampling techniques, we here characterize multiple bound poses and transient encounter complexes of an inhibitor kinase on the statistically converged free-energy landscapes. Our simulations propose parallel binding pathways connected through encounter complexes, highlighting the importance of early binding processes. Determination of transient encounter complexes is still difficult experimentally, and the present finding helps to understand the process and widens the possibility of drug compound design., Modern drug discovery increasingly focuses on the drug-target binding kinetics which depend on drug (un)binding pathways. The conventional molecular dynamics simulation can observe only a few binding events even using the fastest supercomputer. Here, we develop 2D gREST/REUS simulation with enhanced flexibility of the ligand and the protein binding site. Simulation (43 μs in total) applied to an inhibitor binding to c-Src kinase covers 100 binding and unbinding events. On the statistically converged free-energy landscapes, we succeed in predicting the X-ray binding structure, including water positions. Furthermore, we characterize hidden semibound poses and transient encounter complexes on the free-energy landscapes. Regulatory residues distant from the catalytic core are responsible for the initial inhibitor uptake and regulation of subsequent bindings, which was unresolved by experiments. Stabilizing/blocking of either the semibound poses or the encounter complexes can be an effective strategy to optimize drug-target residence time.

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.

GENeralized-Ensemble SImulation System (GENESIS) is a software package for molecular dynamics (MD) simulation of biological systems. It is designed to extend limitations in system size and accessible time scale by adopting highly parallelized schemes and enhanced conformational sampling algorithms. In this new version, GENESIS 1.1, new functions and advanced algorithms have been added. The all-atom and coarse-grained potential energy functions used in AMBER and GROMACS packages now become available in addition to CHARMM energy functions. The performance of MD simulations has been greatly improved by further optimization, multiple time-step integration, and hybrid (CPU + GPU) computing. The string method and replica-exchange umbrella sampling with flexible collective variable choice are used for finding the minimum free-energy pathway and obtaining free-energy profiles for conformational changes of a macromolecule. These new features increase the usefulness and power of GENESIS for modeling and simulation in biological research. © 2017 Wiley Periodicals, Inc.

Rhodopsin is a G-protein coupled receptor functioning as a photoreceptor for vision through photoactivation of a covalently bound ligand of a retinal protonated Schiff base chromophore. Despite the availability of structural information on the inactivated and activated forms of the receptor, the transition processes initiated by the photoabsorption have not been well understood. Here we theoretically examined the photoactivation processes by means of molecular dynamics (MD) simulations and ab initio quantum mechanical/molecular mechanical (QM/MM) free energy geometry optimizations which enabled accurate geometry determination of the ligand molecule in ample statistical conformational samples of the protein. Structures of the intermediate states of the activation process, blue-shifted intermediate and Lumi, as well as the dark state first generated by MD simulations and then refined by the QM/MM free energy geometry optimizations were characterized by large displacement of the β-ionone ring of retinal along with change in the hydrogen bond of the protonated Schiff base. The ab initio calculations of vibrational and electronic spectroscopic properties of those states well reproduced the experimental observations and successfully identified the molecular origins underlying the spectroscopic features. The structural evolution in the formation of the intermediates provides a molecular insight into the efficient activation processes of the receptor.

Molecular recognition underpins all specific protein-ligand interactions and is essential for biomolecular functions. The prediction of canonical binding poses and distinguishing binders from nonbinders are much sought after goals. Here, we apply the generalized replica exchange with solute tempering method, gREST, combined with a flat-bottom potential to evaluate binder and nonbinder interactions with a T4 lysozyme Leu99Ala mutant. The buried hydrophobic cavity and possibility of coupled conformational changes in this protein make binding predictions difficult. The present gREST simulations, enabling enhanced flexibilities of the ligand and protein residues near the binding site, sample bindings in multiple poses, and correct portrayal of X-ray structures. The free-energy profiles of binders (benzene, ethylbenzene, and

In this study, a replica-exchange method was developed to overcome conformational sampling difficulties in computer simulations of spin glass or other systems with rugged free-energy landscapes. This method was then applied to the protein-folding problem in combination with molecular dynamics (MD) simulation. Owing to its simplicity and sampling efficiency, the replica-exchange method has been applied to many other biological problems and has been continuously improved. The method has often been combined with other sampling techniques, such as umbrella sampling, free-energy perturbation, metadynamics, and Gaussian accelerated MD (GaMD). In this chapter, we first summarize the original replica-exchange molecular dynamics (REMD) method and discuss how new algorithms related to the original method are implemented to add new features. Heterogeneous and flexible structures of an N-glycan in a solution are simulated as an example of applications by REMD, replica exchange with solute tempering, and GaMD. The sampling efficiency of these methods on the N-glycan system and the convergence of the free-energy changes are compared. REMD simulation protocols and trajectory analysis using the GENESIS software are provided to facilitate the practical use of advanced simulation methods.

Replica-exchange molecular dynamics (REMD) and their variants have been widely used in simulations of the biomolecular structure and dynamics. Replica exchange with solute tempering (REST) is one of the methods where temperature of a pre-defined solute molecule is exchanged between replicas, while solvent temperatures in all the replicas are kept constant. REST greatly reduces the number of replicas compared to the temperature REMD, while replicas at low temperatures are often trapped under their conditions, interfering with the conformational sampling. Here, we introduce a new scheme of REST, referred to as generalized REST (gREST), where the solute region is defined as a part of a molecule or a part of the potential energy terms, such as the dihedral-angle energy term or Lennard-Jones energy term. We applied this new method to folding simulations of a β-hairpin (16 residues) and a Trp-cage (20 residues) in explicit water. The protein dihedral-angle energy term is chosen as the solute region in the simulations. gREST reduces the number of replicas necessary for good random walks in the solute-temperature space and covers a wider conformational space compared to the conventional REST2. Considering the general applicability, gREST should become a promising tool for the simulations of protein folding, conformational dynamics, and an in silico drug design.

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.

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.

Molecular dynamics (MD) simulations and excitation energy calculations of C1C2 chimera channelrhodopsin, a light-gated ion channel protein utilized as a biotechnological tool for optogenetics, based on a protein structure determined recently by X-ray crystallography were performed to investigate its structural and spectral properties. The MD simulations showed stability of hydrogen-bonds responsible for the channel gating observed in the crystallographic structural model. Analysis of electrostatic contribution of the surrounding protein groups to the absorption energy proposes several site-specific mutations that shift absorption maxima significantly, and provides a clear and controlled guide for engineering design of color variant proteins utilized in optogenetics.

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.

Microbial opsins with a bound chromophore function as photosensitive ion transporters and have been employed in optogenetics for the optical control of neuronal activity. Molecular engineering has been utilized to create colour variants for the functional augmentation of optogenetics tools, but was limited by the complexity of the protein–chromophore interactions. Here we report the development of blue-shifted colour variants by rational design at atomic resolution, achieved through accurate hybrid molecular simulations, electrophysiology and X-ray crystallography. The molecular simulation models and the crystal structure reveal the precisely designed conformational changes of the chromophore induced by combinatory mutations that shrink its π-conjugated system which, together with electrostatic tuning, produce large blue shifts of the absorption spectra by maximally 100 nm, while maintaining photosensitive ion transport activities. The design principle we elaborate is applicable to other microbial opsins, and clarifies the underlying molecular mechanism of the blue-shifted action spectra of microbial opsins recently isolated from natural sources., Retinal-bound opsins are widely used tools for optical control of neuronal activity in vivo, so called optogenetics. Here, using molecular simulations, biochemistry, electrophysiology and X-ray crystallography, the authors present new molecular design principles for the generation of blue-shifted variants of microbial rhodopsins.

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.

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. Background: Light-driven proton pumps are utilized to control the neural activity. Results: We have succeeded to produce a blue-shifted proton pump. The rotation of the β-ionone ring contributes to the spectral shift. Conclusion: The designed color variant provides a tool that allows the control of neural activity by blue light. Significance: The knowledge will help to understand the color-tuning mechanism and can be utilized for optogenetics.

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.

The structure of the putative signaling state in the photocycle of photoactive yellow protein, pB, is investigated by molecular dynamics (MD), steered MD, and free energy calculations. The present work suggests that the structural change toward pB involves a translational displacement of the N-terminal region. The structural rearrangement around alpha3, which is induced by the water molecular migration to the interior part of the protein during the formation of pB', makes the joint region between beta1 and the N-terminal domain flexible and induces the structural displacement of beta1 and the N-terminal domain.

The mechanism of the proton transfer and the concomitant molecular structural and hydrogen bond rearrangements after the photoisomerization of the chromophore in the photocycle of photoactive yellow protein are theoretically investigated by using the QM/MM method and molecular dynamics calculations. The free energy surface along this proton-transfer process is determined. This work suggests the important role of the water molecular migration into the moiety of chromophore, which facilitates proton transfer by the hydrogen bond rearrangement and the hydration of the pB' state.