Fragment molecular orbital study of the electronic excitations in the photosynthetic reaction center of Blastochloris viridis

All electron calculations were performed on the photosynthetic reaction center of Blastochloris viridis, using the fragment molecular orbital (FMO) method. The protein complex of 20,581 atoms and 77,754 electrons was divided into 1398 fragments, and the two‐body expansion of FMO/6‐31G* was applied to calculate the ground state. The excited electronic states of the embedded electron transfer system were separately calculated by the configuration interaction singles approach with the multilayer FMO method. Despite the structural symmetry of the system, asymmetric excitation energies were observed, especially on the bacteriopheophytin molecules. The asymmetry was attributed to electrostatic interaction with the surrounding proteins, in which the cytoplasmic side plays a major role. © 2009 Wiley Periodicals, Inc. J Comput Chem 2010     

Roles of K151 and D180 in L‐2‐haloacid dehalogenase from Pseudomonas sp. YL: Analysis by molecular dynamics and ab initio fragment molecular orbital calculations

L ‐2‐haloacid dehalogenase (L ‐DEX) catalyzes the hydrolytic dehalogenation of L ‐2‐haloalkanoic acids to produce the corresponding D ‐2‐hydroxyalkanoic acids. This enzyme is expected to be applicable to the bioremediation of environments contaminated with halogenated organic compounds. We analyzed the reaction mechanism of L ‐DEX from Pseudomonas sp. YL (L ‐DEX YL) by using molecular modeling. The complexes of wild‐type L ‐DEX YL and its K151A and D180A mutants with its typical substrate, L ‐2‐chloropropionate, were constructed by docking simulation. Subsequently, molecular dynamics (MD) and ab initio fragment molecular orbital (FMO) calculations of the complexes were performed. The ab initio FMO method was applied at the MP2/6‐31G level to estimate interfragment interaction energies. K151 and D180, which are experimentally shown to be important for enzyme activity, interact particularly strongly with L ‐2‐chloropropionate, catalytic water, nucleophile (D10), and with each other. Our calculations suggest that K151 stabilizes substrate orientation and balances the charge around the active site, while D180 stabilizes the rotation of the nucleophile D10, fixes catalytic water around D10, and prevents K151 from approaching D10. Further, D180 may activate catalytic water on its own or with K151, S175, and N177. These roles are consistent with the previous results. Thus, MD and ab initio FMO calculations are powerful tools for the elucidation of the mechanism of enzymatic reaction at the molecular level and can be applied to other catalytically important residues. The results obtained here will play an important role in elucidating the reaction mechanism and rational design of L ‐DEX YL with improved enzymatic activity or substrate specificity. © 2009 Wiley Periodicals, Inc. J Comput Chem, 2009     

All electron quantum chemical calculation of the entire enzyme system confirms a collective catalytic device in the chorismate mutase reaction

To elucidate the catalytic power of enzymes, we analyzed the reaction profile of Claisen rearrangement of Bacillus subtilis chorismate mutase (BsCM) by all electron quantum chemical calculations using the fragment molecular orbital (FMO) method. To the best of our knowledge, this is the first report of ab initio-based quantum chemical calculations of the entire enzyme system, where we provide a detailed analysis of the catalytic factors that accomplish transition-state stabilization (TSS). FMO calculations deliver an ab initio-level estimate of the intermolecular interaction between the substrate and the amino acid residues of the enzyme. To clarify the catalytic role of Arg90, we calculated the reaction profile of the wild-type BsCM as well as Lys90 and Cit90 mutant BsCMs. Structural refinement and the reaction path determination were performed at the ab initio QM/MM level, and FMO calculations were applied to the QM/MM refined structures. Comparison between three types of reactions established two collective catalytic factors in the BsCM reaction: (1) the hydrogen bonds connecting the Glu78-Arg90-substrate cooperatively control the stability of TS relative to the ES complex and (2) the positive charge on Arg90 polarizes the substrate in the TS region to gain more electrostatic stabilization.     

Assessment and acceleration of binding energy calculations for protein–ligand complexes by the fragment molecular orbital method

In the field of drug discovery, it is important to accurately predict the binding affinities between target proteins and drug applicant molecules. Many of the computational methods available for evaluating binding affinities have adopted molecular mechanics‐based force fields, although they cannot fully describe protein–ligand interactions. A noteworthy computational method in development involves large‐scale electronic structure calculations. Fragment molecular orbital (FMO) method, which is one of such large‐scale calculation techniques, is applied in this study for calculating the binding energies between proteins and ligands. By testing the effects of specific FMO calculation conditions (including fragmentation size, basis sets, electron correlation, exchange‐correlation functionals, and solvation effects) on the binding energies of the FK506‐binding protein and 10 ligand complex molecule, we have found that the standard FMO calculation condition, FMO2‐MP2/6‐31G(d), is suitable for evaluating the protein–ligand interactions. The correlation coefficient between the binding energies calculated with this FMO calculation condition and experimental values is determined to be R = 0.77. Based on these results, we also propose a practical scheme for predicting binding affinities by combining the FMO method with the quantitative structure–activity relationship (QSAR) model. The results of this combined method can be directly compared with experimental binding affinities. The FMO and QSAR combined scheme shows a higher correlation with experimental data (R = 0.91). Furthermore, we propose an acceleration scheme for the binding energy calculations using a multilayer FMO method focusing on the protein–ligand interaction distance. Our acceleration scheme, which uses FMO2‐HF/STO‐3G:MP2/6‐31G(d) at R_int = 7.0 Å, reduces computational costs, while maintaining accuracy in the evaluation of binding energy. © 2015 Wiley Periodicals, Inc.     

Ab initio quantum chemical calculation of electron density,electrostatic potential,and electric field of biomolecule based on fragment molecular orbital method

Efficient quantum chemical calculations of electrostatic properties, namely, the electron density (EDN), electrostatic potential (ESP), and electric field (EFL), were performed using the fragment molecular orbital (FMO) method. The numerical errors associated with the FMO scheme were examined at the HF, MP2, and RI‐MP2 levels of theory using 4 small peptides. As a result, the FMO errors in the EDN, ESP, and EFL were significantly smaller than the magnitude of the electron correlation effects, which indicated that the FMO method provides sufficiently accurate values of electrostatic properties. In addition, an attempt to reduce the computational effort was proposed by combining the FMO scheme and a point charge approximation. The error due to this approximation was examined using 2 proteins, prion protein and human immunodeficiency virus type 1 protease. As illustrative examples, the ESP values at the molecular surface of these proteins were calculated at the MP2 level of theory.     

Structural analysis of prion proteins by means of drift cell and traveling wave ion mobility mass spectrometry

The prion protein (PrP) is implicitly involved in the pathogenesis of transmissible spongiform encephalopathies (TSEs). The conversion of normal cellular PrP (PrP^C), a protein that is predominantly α-helical, to a β-sheet-rich isoform (PrP^Sc), which has a propensity to aggregate, is the key molecular event in prion diseases. During its short life span, PrP can experience two different pH environments; a mildly acidic environment, whilst cycling within the cell, and a neutral pH when it is glycosyl phosphatidylinositol (GPI)-anchored to the cell membrane. Ion mobility (IM) combined with mass spectrometry has been employed to differentiate between two conformational isoforms of recombinant Syrian hamster prion protein (SHaPrP). The recombinant proteins studied were α-helical SHaPrP(90-231) and β-sheet-rich SHaPrP(90-231) at pH 5.5 and pH 7.0. The recombinant proteins have the same nominal mass-to-charge ratio (m/z) but differ in their secondary and tertiary structures. A comparison of traveling-wave (T-Wave) ion mobility and drift cell ion mobility (DCIM) mass spectrometry estimated and absolute cross-sections showed an excellent agreement between the two techniques. The use of T-Wave ion mobility as a shape-selective separation technique enabled differentiation between the estimated cross-sections and arrival time distributions (ATDs) of α-helical SHaPrP(90-231) and β-sheet-rich SHaPrP(90-231) at pH 5.5. No differences in cross-section or ATD profiles were observed between the protein isoforms at pH 7.0. The findings have potential implications for a new ante-mortem screening assay, in bodily fluids, for prion misfolding diseases such as TSEs.     

Geometry optimizations of open-shell systems with the fragment molecular orbital method

The ability to perform geometry optimizations on large molecular systems is desirable for both closed- and open-shell species. In this work, the restricted open-shell Hartree-Fock (ROHF) gradients for the fragment molecular orbital (FMO) method are presented. The accuracy of the gradients is tested, and the ability of the method to reproduce adiabatic excitation energies is also investigated. Timing comparisons between the FMO method and full ab initio calculations are also performed, demonstrating the efficiency of the FMO method in modeling large open-shell systems.     

Energy gradients in combined fragment molecular orbital and polarizable continuum model (FMO/PCM) calculation

The analytic energy gradients for the combined fragment molecular orbital and polarizable continuum model (FMO/PCM) method are derived and implemented. Applications of FMO/PCM geometry optimization to polyalanine show that the structures obtained with the FMO/PCM method are very close to those obtained with the corresponding full ab initio PCM methods. FMO/PCM (RHF/6‐31G* level) is used to optimize the solution structure of the 304‐atom Trp‐cage miniprotein and the result is in agreement with NMR experiments. The key factors determining the relative stability of the α‐helix, β‐turn and the extended form in solution are elucidated for polyalanine. © 2009 Wiley Periodicals, Inc. J Comput Chem, 2010     

A configuration analysis for fragment interaction

We propose a modified version of configuration analysis (CA) for the fragment interaction in conjunction with Kitaura’s fragment molecular orbital (FMO) scheme. The proposal is abbreviated as CAFI. The MO sets of fragments are merged and then orthonormalized by the use of a weighted Löwdin orthonormalization. The energy calculation is performed with the concurrent electron relaxation functional (CERF). The relaxation energy is obtained in an orbital-wise fashion and is distinguished as the charge-transfer and the polarization. The utility of CAFI is demonstrated through test calculations on hydrogen-bonding systems.     

Extending the power of quantum chemistry to large systems with the fragment molecular orbital method

Following the brief review of the modern fragment-based methods and other approaches to perform quantum-mechanical calculations of large systems, the theoretical development of the fragment molecular orbital method (FMO) is covered in detail, with the emphasis on the physical properties, which can be computed with FMO. The FMO-based polarizable continuum model (PCM) for treating the solvent effects in large systems and the pair interaction energy decomposition analysis (PIEDA) are described in some detail, and a range of applications of FMO to biological studies is introduced. The factors determining the relative stability of polypeptide conformers (alpha-helix, beta-turn, and extended form) are elucidated using FMO/PCM and PIEDA, and the interactions in the Trp-cage miniprotein construct (PDB: 1L2Y) are analyzed using PIEDA.     

Application of fragment molecular orbital scheme to silicon-containing systems

The fragment molecular orbital (FMO) scheme has been successfully used for a variety of large-scale molecules such as proteins and nucleic acids so far. We have applied the FMO calculations to the silicon-containing systems like polysilanes. The error caused by the fragmentation was examined by the Hartree–Fock method and the second-order Møller-Plesset (MP2) perturbation method for the ground state energy. The dynamic polarizability as a linear response property was also evaluated with and without the fragmentation. A series of numerical comparisons showed that the FMO scheme is applicable to silicon-based molecules with reasonable accuracy. This implied a potential availability of FMO calculations for the issues relevant to nanoscience and nanotechnology.     

Environmental Effects on a Prion's Helix II Domain: Copper(II) and Membrane Interactions with PrP180–193 and Its Analogues

An abnormal interaction between copper and the prion protein is believed to play a pivotal role in the pathogenesis of prion diseases. Copper binding has been mainly attributed to the N‐terminal domain of the prion protein, but this hypothesis has recently been challenged in some papers which suggest that the C‐terminal domain might also compete for metal anchoring. In particular, the segment corresponding to the helix II region of the prion protein, namely PrP180–193, has been shown both to bind copper and to exhibit a copper‐enhanced cytotoxicity, as well as to interact with artificial membranes. The present work is aimed at extending these results by choosing the most representative model of this domain and by determining its copper affinity. With this aim, the different role played by the electrostatic properties of the C‐ and N‐termini of PrP180–193 (VNITIKQHTVTTTT) in determining its conformational behaviour, copper coordination and ability to perturb model membranes was investigated. Owing to the low solubility of PrP180–193, its copper affinity was evaluated by using the shorter PrPAc184–188NH₂ (IKQHT) analogue as a model. ESI‐MS, ESR, UV/Vis, and CD measurements were carried out on the copper(II )/PrPAc184–188NH₂ and copper(II )/PrP180–193NH₂ systems, and showed that PrPAc184–188NH₂ is a reliable model for the metal interaction with the helix II domain. The affinity of copper(II ) for the helix II fragment is higher than that for the octarepeat and PrP106–126 peptides. Finally, the different ability of PrP180–193 analogues to perturb the DPPC model membrane was assessed by DSC measurements. The possible biological consequences of these findings are also discussed briefly.     

Ab initio fragment molecular orbital (FMO) method applied to analysis of the ligand-protein interaction in a pheromone-binding protein

Full quantum computation of the electronic state of proteins has recently become possible by the advent of the ab initio fragment molecular orbital (FMO) method. We applied this method to the analysis of the interaction between the Bombyx mori pheromone-binding protein and its ligand, bombykol. The protein–ligand interaction of this molecular complex was minutely analyzed by the FMO method, and the analysis revealed several important interactions between the ligand and amino acid residues.     

Deciphering the molecular details for the binding of the prion protein to main ganglioside GM1 of neuronal membranes

The prion protein (PrP) resides in lipid rafts in?vivo, and lipids modulate misfolding of the protein to infectious isoforms. Here we demonstrate that binding of recombinant PrP to model raft membranes requires the presence of ganglioside GM1. A combination of liquid- and solid-state NMR revealed the binding sites of PrP to the saccharide head group of GM1. The binding epitope for GM1 was mapped to the folded C-terminal domain of PrP, and docking simulations identified key residues in the C-terminal region of helix C and the loop between strand S2 and helix B. Crucially, this region of PrP is linked to prion resistance in?vivo, and structural changes caused by lipid binding in this region may explain the requirement for lipids in the generation of infectious prions in?vitro.     

Binding interaction analysis of the active site and its inhibitors for neuraminidase (N1 subtype) of human influenza virus by the integration of molecular docking, FMO calculation and 3D-QSAR CoMFA modeling

Recently, the worldwide spread of A/H5N1 avian influenza with high virulence has highlighted the potential threat of human influenza pandemic. Tamiflu and Relenza are currently the only two anti-influenza drugs targeting the neuraminidase (NA) enzyme of human influenza virus. Reports of the emergence of drug resistance further make the development of new potent anti-influenza inhibitors a priority. The X-ray crystallographic study of A/H5N1 avian influenza NA subtypes (Russell, R. J. Nature 2006, 443, 45-49) has demonstrated that there exist two genetically distinct groups, group-1 (N1, N4, N5 and N8) and group-2 (N2, N3, N6, N7 and N9), whose conformations are substantially different. The detailed comparison of their active sites has established, heretofore, the most accurate and solid molecular basis of structure and mechanism for the development of new anti-influenza drugs. In the present study, a three-dimensional structure of N1 subtype of human influenza type A virus (N1hA) has been generated by homology modeling using the X-ray crystallographic structure of N1 subtype of avian influenza virus (N1aA) as the template. Binding interaction analysis between the active site and its inhibitors has been performed by combining ab initio fragment molecular orbital (FMO) calculations and three-dimensional quantitative structure-activity relationship with comparative molecular field analysis (3D-QSAR CoMFA) modeling. Integrated with docking-based 3D-QSAR CoMFA modeling, molecular surface property (electrostatic and steric) mapping and FMO pair interaction analysis, a set of new receptor-ligand binding models and bioaffinity predictive models for rational design and virtual screening of more potent inhibitors of N1hA are established. In addition, the flexibility of the loop-150 of N1hA and N1aA has been examined by a series of molecular dynamics simulations.     

Post-translational modifications in prion proteins

Prions are a novel class of infectious pathogens that cause a group of fatal prion diseases in which the benign cellular form of the prion protein (PrP(C)) is transformed into the disease-related scrapie variant (PrP(SC)). The two PrP isoforms differ in their structure and resistance to degradation. The molecular mechanism by which the PrP(SC) is formed and causes infectivity or neurodegeneration is not known. In a compelling and emerging view, post-translational modifications (or the lack thereof) play roles in the transformation of PrP(C) to PrP(SC). Human PrP contains two consensus sites for N-linked glycosylation, at Asn181 and Asn197. From the functional standpoint, glycosylation can modify either the conformation of PrP(C), or the stability of PrP(SC) and, hence, the rate of PrP(SC) clearance. So far the NMR structures of only recombinant, non-glycosylated prions are known, while the structure of the glycosylated form is estimated by molecular modeling. A number of native amino acid mutations in PrP can be mapped near the glycosylation sites. Normal prion protein has been demonstrated to be a copper binding protein, and increasing evidence has shown correlation between the level of PrP expression and tolerance to oxidative stress. Moreover, histochemistry for nitrotyrosine is used for detection of neuronal labeling, a sign of a peroxynitrite-mediated neuronal degradation and a marker for nitrative stress in scrapie-infected mouse brains. It is an intriguing proposition that the post-translational modifications alone, or in combination with amino acid changes, play dominant roles in the pathogenic transformation of PrP(C) to PrP(SC).     

DNA and estrogen receptor interaction revealed by fragment molecular orbital calculations

Molecular orbital calculations of the complex between DNA-ERE (estrogen response element) and ER (estrogen receptor)-DBD (DNA-binding domain) were performed using the fragment molecular orbital (FMO) method, which enables large-scale MO (molecular orbital) calculations by reducing the computational cost and by significantly increasing efficiency for parallel computation. Such a large system, which contains 3354 atoms, is impractical via conventional MO methods due to the immense computational cost. Details of the interaction between DNA-ERE and ER-DBD were revealed in this study as follows by using the FMO calculations to analyze the interfragment interaction energies (IFIEs) and the electrostatic potentials (ESPs). An area with a high positive ESP is identified on the DNA-binding side of ER-DBD and is the main driving force behind access to the DNA. The position of the ER-DBD monomer can be fixed on a phosphate group of DNA-ERE by the strong electrostatic interactions, whereas the rotation cannot be fixed. In contrast, both the position and rotation of the ER-DBD dimer can be fixed and can therefore form the stable (ER-DBD)2...DNA-ERE complex. Dimerization of the ER-DBD monomers, each of which have a charge of +5 , is mainly due to large attractive interaction energies of the second Zn fragments. The base pairs in the consensus sequence of DNA-ERE interact only with the recognition helix located in the major groove due to the large shielding effect of the phosphate groups of DNA. The recognition helix has weaker interactions with the base pairs than the electrostatic interactions with the phosphate groups. Thus, the DNA-binding machinery of the ER-DBD dimer, which can secure the recognition helix in the major groove of DNA, is crucial for interactions between the recognition helix and base pairs.     

Environmental effects on a prion's helix II domain: copper(II) and membrane interactions with PrP180-193 and its analogues

An abnormal interaction between copper and the prion protein is believed to play a pivotal role in the pathogenesis of prion diseases. Copper binding has been mainly attributed to the N-terminal domain of the prion protein, but this hypothesis has recently been challenged in some papers which suggest that the C-terminal domain might also compete for metal anchoring. In particular, the segment corresponding to the helix II region of the prion protein, namely PrP180-193, has been shown both to bind copper and to exhibit a copper-enhanced cytotoxicity, as well as to interact with artificial membranes. The present work is aimed at extending these results by choosing the most representative model of this domain and by determining its copper affinity. With this aim, the different role played by the electrostatic properties of the C- and N-termini of PrP180-193 (VNITIKQHTVTTTT) in determining its conformational behaviour, copper coordination and ability to perturb model membranes was investigated. Owing to the low solubility of PrP180-193, its copper affinity was evaluated by using the shorter PrPAc184-188NH2 (IKQHT) analogue as a model. ESI-MS, ESR, UV/Vis, and CD measurements were carried out on the copper(II)/PrPAc184-188NH2 and copper(II)/PrP180-193NH2 systems, and showed that PrPAc184-188NH2 is a reliable model for the metal interaction with the helix II domain. The affinity of copper(II) for the helix II fragment is higher than that for the octarepeat and PrP106-126 peptides. Finally, the different ability of PrP180-193 analogues to perturb the DPPC model membrane was assessed by DSC measurements. The possible biological consequences of these findings are also discussed briefly.     

Identification of the pathological prion protein allotypes in scrapie-infected heterozygous bank voles (Clethrionomys glareolus) by high-performance liquid chromatography-mass spectrometry

Cerebral formation of the pathological isoform of the prion protein (PrP) is a crucial molecular event in prion diseases. The bank vole (Clethrionomys glareolus) is a rodent species highly susceptible to natural scrapie. The PrP gene of bank vole is polymorphic (Met/Ile) at codon 109. Here we show that homozygous 109Met/Met voles have incubation times shorter than heterozygous 109Met/Ile voles after experimental challenge with three different scrapie isolates. An HPLC-MS/MS method was optimized and applied to investigate whether in heterozygous animals both PrP allotypes are able to undergo pathological conversion. The results demonstrate that both allotypes of the prion protein participate to pathological deposition.