Fully quantum mechanical energy optimization for protein-ligand structure

We present a quantum mechanical approach to study protein-ligand binding structure with application to a Adipocyte lipid-binding protein complexed with Propanoic Acid. The present approach employs a recently develop molecular fractionation with a conjugate caps (MFCC) method to compute protein-ligand interaction energy and performs energy optimization using the quasi-Newton method. The MFCC method enables us to compute fully quantum mechanical ab initio protein-ligand interaction energy and its gradients that are used in energy minimization. This quantum optimization approach is applied to study the Adipocyte lipid-binding protein complexed with Propanoic Acid system, a complex system consisting of a 2057-atom protein and a 10-atom ligand. The MFCC calculation is carried out at the Hartree-Fock level with a 3-21G basis set. The quantum optimized structure of this complex is in good agreement with the experimental crystal structure. The quantum energy calculation is implemented in a parallel program that dramatically speeds up the MFCC calculation for the protein-ligand system. Similarly good agreement between MFCC optimized structure and the experimental structure is also obtained for the streptavidin-biotin complex. Due to heavy computational cost, the quantum energy minimization is carried out in a six-dimensional space that corresponds to the rigid-body protein-ligand interaction.     

Quantum study of HIV-1 protease-bridge water interaction

We present a fully quantum mechanical calculation for binding interaction between HIV-1 protease (PR) and the water molecule W301 which bridges the flaps of the protease with the inhibitors of PR. The quantum calculation is made possible by applying a recently developed molecular fractionation with conjugate caps (MFCC) method which divides a protein molecule into capped amino acid-based fragments and their conjugate caps. These individual fragments are properly treated to preserve the chemical property of bonds that are cut. Ab initio methods at HF, B3LYP, and MP2 levels with a fixed basis set 6-31+G* have been employed in the present calculation. The MFCC calculation produces a quantum mechanical interaction "map" representing interactions between individual residues of PR and W301. This enables a detailed quantitative analysis on binding of W301 to specific residues of PR at quantum mechanical level.     

Fractionation of peptide with disulfide bond for quantum mechanical calculation of interaction energy with molecules

We present a computational study of a recently developed molecular fractionation with conjugated caps (MFCC) method for application to peptide/protein that has disulfide bonds. Specifically, we employ the MFCC approach to generate peptide fragments in which a disulfide bond is cut and a pair of conjugated caps are inserted. The method is tested on two peptides interacting with a water molecule. The first is a dipeptide consisting of two cysteines (Cys-Cys) connected by a disulfide bond and the second is a seven amino acid peptide consisting of Gly-Cys-Gly-Gly-Gly-Cys-Gly with a disulfide cross link. One-dimensional peptide-water potential curves are computed using the MFCC method at various ab initio levels for a number of interaction geometries. The calculated interaction energies are found to be in excellent agreement with the results obtained from the corresponding full system ab initio calculations for both peptide/water systems. The current study provides further numerical support for the accuracy of the MFCC method in full quantum mechanical calculation of protein/peptide that contains disulfide bonds.     

Theoretical method for full ab initio calculation of DNA/RNA-ligand interaction energy

In this paper, we further develop the molecular fractionation with conjugate caps (MFCC) scheme for quantum mechanical computation of DNA-ligand interaction energy. We study three oligonuclear acid interaction systems: dinucleotide dCG/water, trinucleotide dCGT/water, and a Watson-Crick paired DNA segment, dCGT/dGCA. Using the basic MFCC approach, the nucleotide chains are cut at each phosphate group and a pair of conjugate caps (concaps) are inserted. Five cap molecules have been tested among which the dimethyl phosphate anion is proposed to be the standard concap for application. For each system, one-dimensional interaction potential curves are computed using the MFCC method and the calculated interaction energies are found to be in excellent agreement with corresponding results obtained from the full system ab initio calculations. The current study extends the application of the MFCC method to ab initio calculations for DNA- or RNA-ligand interaction energies.     

A new quantum method for electrostatic solvation energy of protein

A new method that incorporates the conductorlike polarizable continuum model (CPCM) with the recently developed molecular fractionation with conjugate caps (MFCC) approach is developed for ab initio calculation of electrostatic solvation energy of protein. The application of the MFCC method makes it practical to apply CPCM to calculate electrostatic solvation energy of protein or other macromolecules in solution. In this MFCC-CPCM method, calculation of protein solvation is divided into calculations of individual solvation energies of fragments (residues) embedded in a common cavity defined with respect to the entire protein. Besides computational efficiency, the current approach also provides additional information about contribution to protein solvation from specific fragments. Numerical studies are carried out to calculate solvation energies for a variety of peptides including alpha helices and beta sheets. Excellent agreement between the MFCC-CPCM result and those from the standard full system CPCM calculation is obtained. Finally, the MFCC-CPCM calculation is applied to several real proteins and the results are compared to classical molecular mechanics Poisson-Boltzmann (MM/PB) and quantum Divid-and-Conque Poisson-Boltzmann (D&C-PB) calculations. Large wave function distortion energy (solute polarization energy) is obtained from the quantum calculation which is missing in the classical calculation. The present study demonstrates that the MFCC-CPCM method is readily applicable to studying solvation of proteins.     

The generalized molecular fractionation with conjugate caps/molecular mechanics method for direct calculation of protein energy

A generalized molecular fractionation with conjugate caps/molecular mechanics (GMFCC/MM) scheme is developed for efficient linear-scaling quantum mechanical calculation of protein energy. In this GMFCC/MM scheme, the interaction energy between neighboring residues as well as between non-neighboring residues that are spatially in close contact are computed by quantum mechanics while the rest of the interaction energy is computed by molecular mechanics. Numerical studies are carried out to calculate torsional energies of six polypeptides using the GMFCC/MM approach and the energies are shown to be in general good agreement with the full system quantum calculation. Among those we tested is a polypeptide containing 396 atoms whose energies are computed at the MP26-31G* level. Our study shows that using GMFCC/MM, it is possible to perform high level ab initio calculation such as MP2 for applications such as structural optimization of protein complex and molecular dynamics simulation.     

Cation-pi interactions: an energy decomposition analysis and its implication in delta-opioid receptor-ligand binding

The nature and strength of the cation-pi interaction in protein-ligand binding are modeled by considering a series of nonbonded complexes involving N-substituted piperidines and substituted monocylic aromatics that mimic the delta-opioid receptor-ligand binding. High-level ab initio quantum mechanical calculations confirm the importance of such cation-pi interactions, whose intermolecular interaction energy ranges from -6 to -12 kcal/mol. A better understanding of the electrostatics, polarization, and other intermolecular interactions is obtained by appropriately decomposing the total interaction energy into their individual components. The energy decomposition analysis is also useful for parametrizing existing molecular mechanics force fields that could then account for energetic contributions arising out of cation-pi interactions in biomolecules. The present results further provide a framework for interpreting experimental results from point mutation reported for the delta-opioid receptor.     

Fully ab initio protein‐ligand interaction energies with dispersion corrected density functional theory

Dispersion corrected density functional theory (DFT‐D3) is used for fully ab initio protein‐ligand (PL) interaction energy calculation via molecular fractionation with conjugated caps (MFCC) and applied to PL complexes from the PDB comprising 3680, 1798, and 1060 atoms. Molecular fragments with n amino acids instead of one in the original MFCC approach are considered, thereby allowing for estimating the three‐body and higher many‐body terms. n > 1 is recommended both in terms of accuracy and efficiency of MFCC. For neutral protein side‐chains, the computed PL interaction energy is visibly independent of the fragment length n. The MFCC fractionation error is determined by comparison to a full‐system calculation for the 1060 atoms containing PL complex. For charged amino acid side‐chains, the variation of the MFCC result with n is increased. For these systems, using a continuum solvation model with a dielectricity constant typical for protein environments (? = 4) reduces both the variation with n and improves the stability of the DFT calculations considerably. The PL interaction energies for two typical complexes obtained ab initio for the first time are found to be rather large (?30 and ?54 kcal/mol). © 2012 Wiley Periodicals, Inc.     

An efficient approach for ab initio energy calculation of biopolymers

We present a new method for efficient total-energy calculation of biopolymers using the density-matrix (DM) scheme based on the molecular fractionation with conjugate caps (MFCC) approach. In this MFCC-DM method, a biopolymer such as a protein is partitioned into properly capped fragments whose density matrices are calculated by conventional ab initio methods which are then assembled to construct the full system density matrix. The assembled full density matrix is then employed to calculate the total energy and dipole moment of the protein using Hartree-Fock or density-functional theory methods. Using this MFCC-DM method, the self-consistent-field procedure for solving the full Hamiltonian problem is avoided and an efficient approach for ab initio energy calculation of biopolymers is achieved. Two implementations of the approach are presented in this paper. Systematic numerical studies are carried out on a series of extended polyglycines CH3CO-(GLY)n-NHCH3(n = 3-25) and excellent results are obtained.     

Reaction path potential for complex systems derived from combined ab initio quantum mechanical and molecular mechanical calculations

Combined ab initio quantum mechanical and molecular mechanical calculations have been widely used for modeling chemical reactions in complex systems such as enzymes, with most applications being based on the determination of a minimum energy path connecting the reactant through the transition state to the product in the enzyme environment. However, statistical mechanics sampling and reaction dynamics calculations with a combined ab initio quantum mechanical (QM) and molecular mechanical (MM) potential are still not feasible because of the computational costs associated mainly with the ab initio quantum mechanical calculations for the QM subsystem. To address this issue, a reaction path potential energy surface is developed here for statistical mechanics and dynamics simulation of chemical reactions in enzymes and other complex systems. The reaction path potential follows the ideas from the reaction path Hamiltonian of Miller, Handy and Adams for gas phase chemical reactions but is designed specifically for large systems that are described with combined ab initio quantum mechanical and molecular mechanical methods. The reaction path potential is an analytical energy expression of the combined quantum mechanical and molecular mechanical potential energy along the minimum energy path. An expansion around the minimum energy path is made in both the nuclear and the electronic degrees of freedom for the QM subsystem internal energy, while the energy of the subsystem described with MM remains unchanged from that in the combined quantum mechanical and molecular mechanical expression and the electrostatic interaction between the QM and MM subsystems is described as the interaction of the MM charges with the QM charges. The QM charges are polarizable in response to the changes in both the MM and the QM degrees of freedom through a new response kernel developed in the present work. The input data for constructing the reaction path potential are energies, vibrational frequencies, and electron density response properties of the QM subsystem along the minimum energy path, all of which can be obtained from the combined quantum mechanical and molecular mechanical calculations. Once constructed, it costs much less for its evaluation. Thus, the reaction path potential provides a potential energy surface for rigorous statistical mechanics and reaction dynamics calculations of complex systems. As an example, the method is applied to the statistical mechanical calculations for the potential of mean force of the chemical reaction in triosephosphate isomerase.     

Theoretical study on isotope and temperature effect in hydronium ion using ab initio path integral simulation

We have applied the ab initio path integral molecular dynamics simulation to study hydronium ion and its isotopes, which are the simplest systems for hydrated proton and deuteron. In this simulation, all the rotational and vibrational degrees of freedom are treated fully quantum mechanically, while the potential energies of the respective atomic configurations are calculated "on the fly" using ab initio quantum chemical approach. With the careful treatment of the ab initio electronic structure calculation by relevant choices in electron correlation level and basis set, this scheme is theoretically quite rigorous except for Born-Oppenheimer approximation. This accurate calculation allows a close insight into the structural shifts for the isotopes of hydronium ion by taking account of both quantum mechanical and thermal effects. In fact, the calculation is shown to be successful to quantitatively extract the geometrical isotope effect with respect to the Walden inversion. It is also shown that this leads to the isotope effect on the electronic structure as well as the thermochemical properties.     

A new method for direct calculation of total energy of protein

A new scheme is developed for efficient quantum mechanical calculation of total energy of protein based on a recently developed MFCC (molecular fractionation with conjugate caps) approach. In this scheme, the linear-scaling MFCC method is first applied to calculate total electron density of protein. The computed electron density is then employed for direct numerical integration in density functional theory (DFT) to yield total energy of protein, with the kinetic energy obtained by a proposed ansatz. Numerical studies are carried out to calculate torsional energies of two polypeptides using this approach and the energies are shown to be in good agreement with the corresponding full system DFT calculation.     

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.     

VISCANA: visualized cluster analysis of protein-ligand interaction based on the ab initio fragment molecular orbital method for virtual ligand screening

We have developed a visualized cluster analysis of protein-ligand interaction (VISCANA) that analyzes the pattern of the interaction of the receptor and ligand on the basis of quantum theory for virtual ligand screening. Kitaura et al. (Chem. Phys. Lett. 1999, 312, 319-324.) have proposed an ab initio fragment molecular orbital (FMO) method by which large molecules such as proteins can be easily treated with chemical accuracy. In the FMO method, a total energy of the molecule is evaluated by summation of fragment energies and interfragment interaction energies (IFIEs). In this paper, we have proposed a cluster analysis using the dissimilarity that is defined as the squared Euclidean distance between IFIEs of two ligands. Although the result of an ordered table by clustering is still a massive collection of numbers, we combine a clustering method with a graphical representation of the IFIEs by representing each data point with colors that quantitatively and qualitatively reflect the IFIEs. We applied VISCANA to a docking study of pharmacophores of the human estrogen receptor alpha ligand-binding domain (57 amino acid residues). By using VISCANA, we could classify even structurally different ligands into functionally similar clusters according to the interaction pattern of a ligand and amino acid residues of the receptor protein. In addition, VISCANA could estimate the correct docking conformation by analyzing patterns of the receptor-ligand interactions of some conformations through the docking calculation.     

Formal Estimation of Errors in Computed Absolute Interaction Energies of Protein-ligand Complexes

A largely unsolved problem in computational biochemistry is the accurate prediction of binding affinities of small ligands to protein receptors. We present a detailed analysis of the systematic and random errors present in computational methods through the use of error probability density functions, specifically for computed interaction energies between chemical fragments comprising a protein-ligand complex. An HIV-II protease crystal structure with a bound ligand (indinavir) was chosen as a model protein-ligand complex. The complex was decomposed into twenty-one (21) interacting fragment pairs, which were studied using a number of computational methods. The chemically accurate complete basis set coupled cluster theory (CCSD(T)/CBS) interaction energies were used as reference values to generate our error estimates. In our analysis we observed significant systematic and random errors in most methods, which was surprising especially for parameterized classical and semiempirical quantum mechanical calculations. After propagating these fragment-based error estimates over the entire protein-ligand complex, our total error estimates for many methods are large compared to the experimentally determined free energy of binding. Thus, we conclude that statistical error analysis is a necessary addition to any scoring function attempting to produce reliable binding affinity predictions.     

New method for direct linear-scaling calculation of electron density of proteins

A new scheme for direct linear-scaling quantum mechanical calculation of electron density of protein systems is developed. The new scheme gives much improved accuracy of electron density for proteins than the original MFCC (molecular fractionation with conjugate caps) approach in efficient linear-scaling calculation for protein systems. In this new approach, the error associated with each cut in the MFCC approach is estimated by computing the two neighboring amino acids in both cut and uncut calculations and is corrected. Numerical tests are performed on six oligopeptide taken from PDB (protein data bank), and the results show that the new scheme is efficient and accurate.     

An ab initio molecular orbital study of intramolecular H-bonding: 1,3 propanediol

An ab initio molecular orbital calculation has been carried out for three different conformations of 1,3 propanediol, one of which permits intramolecular H-bond studied by ab initio quantum mechanical methods. The ΔE for H-bonding formation is compated to be 0.9 kcal/mole and the charge redistributions and molecular orbital energy changes are compared to those found in intermolecular H-bonds.     

An extension of the grid empowered molecular simulator to quantum reactive scattering

Within the activities of the D37 COST Action, we have further developed the quantum dynamics framework of the grid empowered molecular simulator (GEMS) implemented on the segment of the European grid available to the COMPCHEM (computational chemistry) virtual organization. GEMS does now include in a full ab initio approach, the evaluation of the detailed quantum (both time dependent and time independent) dynamics of small systems starting from the calculation of the electronic structure properties as well as the direct calculation of thermalized properties. Illustrative, full dimensional applications of the extended simulator to the H + H(2) , N + N(2) , and O + O(2) systems are presented.     

Role of hydrogen-bond network in energy storage of bacteriorhodopsin's light-driven proton pump revealed by ab initio normal-mode analysis

Vibrational modes of the hydrogen-bond network in the binding site of bacteriorhodopsin (bR), a protein in halobacteria functioning as a light-driven proton pump, were investigated by an ab initio quantum mechanical/molecular mechanical (QM/MM) method. Normal-mode analysis calculations for O-D and N-D stretching modes of internal water molecules and the Schiff base of the retinal chromophore in the early intermediate state, K, reproduced well experimentally observed vibrational spectra. Supported by agreement with observed spectra, the QM/MM calculation suggests that weakened hydrogen bonds upon photoisomerization of the chromophore are an important means of energy storage in bR.