Insights into the regioselectivity and RNA-binding affinity of HIV-1 nucleocapsid protein from linear-scaling quantum methods

Theoretical characterization of the shikimate 5-dehydrogenase reaction from Mycobacterium tuberculosis by hybrid QC/MM simulations and quantum chemical descriptors

In this study, we have investigated the enzyme shikimate 5-dehydrogenase from the causative agent of tuberculosis, Mycobacterium tuberculosis. We have employed a mixture of computational techniques, including molecular dynamics, hybrid quantum chemical/molecular mechanical potentials, relaxed surface scans, quantum chemical descriptors and free-energy simulations, to elucidate the enzyme's reaction pathway. Overall, we find a two-step mechanism, with a single transition state, that proceeds by an energetically uphill hydride transfer, followed by an energetically downhill proton transfer. Our mechanism and calculated free energy barrier for the reaction, 64.9 kJ mol^- 1, are in good agreement with those predicted from experiment. An analysis of quantum chemical descriptors along the reaction pathway indicated a possibly important, yet currently unreported, role of the active site threonine residue, Thr65.

Benchmarking quantum mechanical methods for calculating reaction energies of reactions catalyzed by enzymes

To assess the accuracy of different quantum mechanical methods for biochemical modeling, the reaction energies of 20 small model reactions (chosen to represent chemical steps catalyzed by commonly studied enzymes) were calculated. The methods tested included several popular Density Functional Theory (DFT) functionals, second-order Møller Plesset perturbation theory (MP2) and its spin-component scaled variant (SCS-MP2), and coupled cluster singles and doubles and perturbative triples (CCSD(T)). Different basis sets were tested. CCSD(T)/aug-cc-pVTZ results for all 20 reactions were used to benchmark the other methods. It was found that MP2 and SCS-MP2 reaction energy calculation results are similar in quality to CCSD(T) (mean absolute error (MAE) of 1.2 and 1.3 kcal mol−1, respectively). MP2 calculations gave a large error in one case, and are more subject to basis set effects, so in general SCS-MP2 calculations are a good choice when CCSD(T) calculations are not feasible. Results with different DFT functionals were of reasonably good quality (MAEs of 2.5–5.1 kcal mol−1), whereas popular semi-empirical methods (AM1, PM3, SCC-DFTB) gave much larger errors (MAEs of 11.6–14.6 kcal mol−1). These results should be useful in guiding methodological choices and assessing the accuracy of QM/MM calculations on enzyme-catalyzed reactions.

Specific Reaction Parametrization of the AM1/d Hamiltonian for Phosphoryl Transfer Reactions: H, O, and P Atoms

A semiempirical AM1/d Hamiltonian is developed to model phosphoryl transfer reactions catalyzed by enzymes and ribozymes for use in linear-scaling calculations and combined quantum mechanical/molecular mechanical simulations. The model, designated AM1/d-PhoT, is parametrized for H, O, and P atoms to reproduce high-level density-functional results from a recently constructed database of quantum calculations for RNA catalysis ( http://theory.chem.umn.edu/Database/QCRNA ), including geometries and relative energies of minima, transition states and reactive intermediates, dipole moments, proton affinities, and other relevant properties. The model is tested in the gas phase and in solution using a QM/MM potential. The results indicate that the method provides significantly higher accuracy than MNDO/d, AM1, and PM3 methods and, for the transphosphorylation reactions, is in close agreement with the density-functional calculations at the B3LYP/6-311++G(3df,2p) level with a reduction in computational cost of 3-4 orders of magnitude. The model is expected to have considerable impact on the application of semiempirical QM/MM methods to transphosphorylation reactions in solution, enzymes, and ribozymes and to ultimately facilitate the design of improved next-generation multiscale quantum models.

Computational enzymology

Techniques for modelling enzyme-catalyzed reaction mechanisms are making increasingly important contributions to biochemistry. They can address fundamental questions in enzyme catalysis and have the potential to contribute to practical applications such as drug development.

Biomolecular simulation and modelling: status, progress and prospects

Molecular simulation is increasingly demonstrating its practical value in the investigation of biological systems. Computational modelling of biomolecular systems is an exciting and rapidly developing area, which is expanding significantly in scope. A range of simulation methods has been developed that can be applied to study a wide variety of problems in structural biology and at the interfaces between physics, chemistry and biology. Here, we give an overview of methods and some recent developments in atomistic biomolecular simulation. Some recent applications and theoretical developments are highlighted.

Self-consistent polarization neglect of diatomic differential overlap: application to water clusters

Semiempirical self-consistent field (SCF) methods based on the neglect of diatomic differential overlap (NDDO) formalism have the ability to treat the formation and breaking of chemical bonds but have been found to poorly describe hydrogen bonding and weak electrostatic complexes. In contrast, most empirical potentials are not able to describe bond breaking and formation but have the ability to add missing elements of hydrogen bonding by using classical electrostatic interactions. We present a new method which combines aspects of both NDDO-based SCF techniques and classical descriptions of polarization to describe the diffuse nature of the electronic wavefunction in a self-consistent manner. We develop the "self-consistent polarization neglect of diatomic differential overlap" (SCP-NDDO) theory with the additional description of molecular dispersion developed as a second-order perturbation theory expression. The current study seeks to model water-water interactions as a test case. To this end, we have parametrized the method to accurate ab initio complete basis set limit estimates of small water cluster binding energies of Xantheas and co-workers [J. Chem. Phys. 116, 1493 (2002); 120, 823 (2004)]. Overall agreement with the ab initio binding energies (n=2-6, and 8) is achieved with a rms error of 0.19 kcal/mol. We achieve noticeable improvements in the structure, vibrational frequencies, and energetic predictions of water clusters (n< or =21) relative to standard NDDO-based methods.

Structural determination of large molecules through the reassembly of optimized fragments

The accurate determination of the optimized structures of large molecules is, computationally quite expensive. As an alternate to the conventional approaches to structural optimization, we have explored the accuracy and speed-up obtained when variants of the fragmentation optimization and recombination method (FORM) are used. Specifically, the method was applied to eight prototypical molecules -n-decane, hexa-alanine, a long conjugate hydrocarbon molecule, a large polar conjugated molecule, large (5,5) armchair single-walled carbon nanotubes, a salen-aluminum complex and a multiply H-bonded system (two conformers of vancomycin aglycon with Di-N-acetyl-l-Lys-d-Ala-d-Ala - without optimizing the structure of the whole molecules. We find that FORM can predict the structure of these molecules with an acceptable accuracy, all at a computational cost that is 2-11 times less than conventional quantum mechanical methods at the Hartree-Fock (HF), density functional theory (DFT) and MP2 level of accuracy. FORM may therefore represent a viable approach for the fast structural predictions of large molecules.

Extension of the Self-Consistent-Charge Density-Functional Tight-Binding Method: Third-Order Expansion of the Density Functional Theory Total Energy and Introduction of a Modified Effective Coulomb Interaction

The standard self-consistent-charge density-functional-tight-binding (SCC-DFTB) method (Phys. Rev. B 1998, 58, 7260) is derived by a second-order expansion of the density functional theory total energy expression, followed by an approximation of the charge density fluctuations by charge monopoles and an effective damped Coulomb interaction between the atomic net charges. The central assumptions behind this effective charge-charge interaction are the inverse relation of atomic size and chemical hardness and the use of a fixed chemical hardness parameter independent of the atomic charge state. While these approximations seem to be unproblematic for many covalently bound systems, they are quantitatively insufficient for hydrogen-bonding interactions and (anionic) molecules with localized net charges. Here, we present an extension of the SCC-DFTB method to incorporate third-order terms in the charge density fluctuations, leading to chemical hardness parameters that are dependent on the atomic charge state and a modification of the Coulomb scaling to improve the electrostatic treatment within the second-order terms. These modifications lead to a significant improvement in the description of hydrogen-bonding interactions and proton affinities of biologically relevant molecules.

Improvement of semiempirical response properties with charge-dependent response density

The present work outlines a new method for treatment of charge-dependent polarizability in semiempirical quantum models for use in combined quantum-mechanical/molecular mechanical simulations of biological reactions. The method addresses a major shortcoming in the performance of conventional semiempirical models for these simulations that is tied to the use of a localized minimal atomic-orbital basis set. The present approach has the advantages that it uses a density basis that retains a set of linear-response equations, does not increase the atomic-orbital basis, and avoids the problem of artificial charge transfer and scaling of the polarizability seen in related models that allow atomic charges to fluctuate. The model introduces four new atom-based parameters and has been tested with the modified neglect of differential overlap d-orbital Hamiltonian against 1132 molecules and ions and shown to decrease the dipole moment and polarizability errors by factors of 2 and 10, respectively, with respect to density-functional results. The method performs impressively for a variety of charge states (from 2+ to 2-), and offers a potentially powerful extension in the design of next generation semiempirical quantum models for accurate simulations of highly charged biological reactions.

Multilevel and density functional electronic structure calculations of proton affinities and gas-phase basicities involved in biological phosphoryl transfer

Five multilevel model chemistries (CBS-QB3, G3B3, G3MP2B3, MCG3/3, and MC-QCISD/3) and seven hybrid density functional methods (PBE0, B1B95, B3LYP, MPW1KCIS, PBE1KCIS, and MPW1B95) have been applied to the calculation of gas-phase basicity and proton affinity values for a series of 17 molecules relevant to the study of biological phosphoryl transfer. In addition, W1 calculations were performed on a subset of molecules. The accuracy of the methods was assessed and the nature of systematic errors was explored, leading to the introduction of a set of effective bond enthalpy and entropy correction terms. The multicoefficient correlation methods (MCG3/3 and MC-QCISD), with inclusion of specific zero-point scale factors, slightly outperform the other multilevel methods tested (CBS-QB3, G3B3, and G3MP2B3), with significantly less computational cost, and in the case of MC-QCISD, slightly less severe scaling. Four density functional methods, PBE1KCIS, MPW1B95, PBE0, and B1B95 perform nearly as well as the multilevel methods. These results provide an important set of benchmarks relevant to biological phosphoryl transfer reactions.

Smooth solvation method for d-orbital semiempirical calculations of biological reactions. 1. Implementation

The present paper describes the extension of a recently developed smooth conductor-like screening model for solvation to a d-orbital semiempirical framework (MNDO/d-SCOSMO) with analytic gradients that can be used for geometry optimizations, transition state searches, and molecular dynamics simulations. The methodology is tested on the potential energy surfaces for separating ions and the dissociative phosphoryl transfer mechanism of methyl phosphate. The convergence behavior of the smooth COSMO method with respect to discretization level is examined and the numerical stability of the energy and gradient are compared to that from conventional COSMO calculations. The present method is further tested in applications to energy minimum and transition state geometry optimizations of neutral and charged metaphosphates, phosphates, and phosphoranes that are models for stationary points in transphosphorylation reaction pathways of enzymes and ribozymes. The results indicate that the smooth COSMO method greatly enhances the stability of quantum mechanical geometry optimization and transition state search calculations that would routinely fail with conventional solvation methods. The present MNDO/d-SCOSMO method has considerable computational advantages over hybrid quantum mechanical/molecular mechanical methods with explicit solvation, and represents a potentially useful tool in the arsenal of multi-scale quantum models used to study biochemical reactions.

RNA recognition mechanism of the minimal active domain of the human immunodeficiency virus type-2 nucleocapsid protein

NCp8 of HIV-2 contains two CCHC-type zinc fingers connected by a linker, and is involved in many critical steps of the virus life cycle. It was previously shown that the first zinc finger flanked by the linker is the minimal active domain for specific binding to viral RNA. In our previous study, we determined the three-dimensional structure of NCp8-f1, including the minimal active domain, and found that a hydrogen bond between Asn(11) N(delta)H and Arg(27) O stabilized the conformation of the linker in the vicinity of the zinc finger [Kodera et al. (1998) Biochemistry 37, 17704-17713]. In this study, RNA binding activities of NCp8-f1 and three types of its mutant peptides were analysed by native PAGE assay. The activity and three-dimensional structure of NCp8-f1/N11A, in which alanine is substituted for Asn(11) thereby affecting the conformation of the linker, was analyzed and compared with those of NCp8-f1. We demonstrated that the existence of Arg(4) and/or Lys(5) and Arg(26) and/or Arg(27) were necessary for binding RNA. Furthermore, the linker's flexible orientation, which is controlled by the hydrogen bond between Asn(11) N(delta)H and Arg(27) O, appears to be a structural basis for NCp8 existing as a multi-functional protein.

Comparison of the Specificity of Interaction of Cellular and Viral Zinc-Binding Domains with 2-Mercaptobenzamide Thioesters

The interactions of two 2-mercaptobenzamide thioester compounds with six diverse zinc-binding domains (ZBDs) have been analyzed by UV/visible spectroscopy, NMR spectroscopy, and nucleic acid binding assays. These thioester compounds serve as useful tools for probing the intrinsic chemical stability of ZBDs that exist within a variety of cellular and viral proteins. In our studies, the classical (Cys(2)His(2)) zinc finger ZBDs, the interleaved RING like ZBDs of protein kinase C delta (Cys(2)HisCys and HisCys(3)), and the carboxyl-terminal (Cys(2)HisCys) ZBD of Mouse Mammary Tumor Virus nucleocapsid protein (MMTV NCp10) were resistant to reaction with the thioester compounds. In contrast, the thioester compounds were able to efficiently eject zinc from the amino-terminal (Cys(2)HisCys) ZBD of MMTV NCp10, a Cys(2)HisCys ZBD from Friend of GATA-1 (FOG-1), and from both Cys(4) ZBDs of GATA-1. In all cases, zinc ejection led to a loss of protein structure. Interestingly, GATA-1 was resistant to reaction with the thioester compounds when bound to its target DNA sequence. The electronic and steric screening was calculated for select ZBDs to further explore their reactivity. Based on these results, it appears that both first and second zinc-coordination shell interactions within ZBDs, as well as nucleic acid binding, play important roles in determining the chemical stability and reactivity of ZBDs. These studies not only provide information regarding the relative reactivity of cysteine residues within structural ZBDs but also are crucial for the design of future therapeutic agents that selectively target ZBDs, such as those that occur in the HIV-1 nucleocapsid protein.

QCRNA 1.0: a database of quantum calculations for RNA catalysis

This work outlines a new on-line database of quantum calculations for RNA catalysis (QCRNA) available via the worldwide web at http://theory.chem.umn.edu/QCRNA. The database contains high-level density functional calculations for a large range of molecules, complexes and chemical mechanisms important to phosphoryl transfer reactions and RNA catalysis. Calculations are performed using a strict, consistent protocol such that a wealth of cross-comparisons can be made to elucidate meaningful trends in biological phosphate reactivity. Currently, around 2000 molecules have been collected in varying charge states in the gas phase and in solution. Solvation was treated with both the PCM and COSMO continuum solvation models. The data can be used to study important trends in reactivity of biological phosphates, or used as benchmark data for the design of new semiempirical quantum models for hybrid quantum mechanical/molecular mechanical simulations.

Modelling enzyme reaction mechanisms, specificity and catalysis

Modern modelling methods can now give uniquely detailed understanding of enzyme-catalyzed reactions, including the analysis of mechanisms and the identification of determinants of specificity and catalytic efficiency. A new field of computational enzymology has emerged that has the potential to contribute significantly to structure-based design and to develop predictive models of drug metabolism and, for example, of the effects of genetic polymorphisms. This review outlines important techniques in this area, including quantum-chemical model studies and combined quantum-mechanics and molecular-mechanics (QM/MM) methods. Some recent applications to enzymes of pharmacological interest are also covered, showing the types of problems that can be tackled and the insight they can give.

Protic Conversion of Nitrile into Azavinylidene Complexes of Rhenium, a Mechanistic Theoretical Study

The mechanism of the protonation of the rhenium nitrile chloro-complexes [ReCl(NCCH3)(PH3)4] (2), taken as models of the real systems [ReCl(NCR)(dppe)(2)] (dppe = Ph2PCH2CH2PPh2), leading to the azavinylidene products [ReCl(NC(H)CH3)(PH3)4]+ (3) was investigated by theoretical methods at the B3LYP level of theory. Electrostatic and molecular orbital arguments and thermodynamic, kinetic, and steric factors are analyzed and indicate that the chlorine atom is the most probable site of the initial proton attack, although the direct protonation of the nitrile carbon atom is also possible as a concurrent process. For the cis-isomer of 2, the initially formed chloro-protonated species cis-[Re(ClH)(NCCH3)(PH3)4]+ further converts to the azavinylidene cis-3 via either an acid-independent 1,4-proton shift or an acid-base catalyzed pathway involving a second protonation of the nitrile carbon atom to give cis-[Re(ClH)(NC(H)CH3)(PH3)4]2+ followed by elimination of the proton from the chlorine atom.

Single point mutations in the zinc finger motifs of the human immunodeficiency virus type 1 nucleocapsid alter RNA binding specificities of the gag protein and enhance packaging and infectivity

A specific interaction between the nucleocapsid (NC) domain of the Gag polyprotein and the RNA encapsidation signal (Psi) is required for preferential incorporation of the retroviral genomic RNA into the assembled virion. Using the yeast three-hybrid system, we developed a genetic screen to detect human immunodeficiency virus type 1 (HIV-1) Gag mutants with altered RNA binding specificities. Specifically, we randomly mutated full-length HIV-1 Gag or its NC portion and screened the mutants for an increase in affinity for the Harvey murine sarcoma virus encapsidation signal. These screens identified several NC zinc finger mutants with altered RNA binding specificities. Furthermore, additional zinc finger mutants that also demonstrated this phenotype were made by site-directed mutagenesis. The majority of these mutants were able to produce normal virion-like particles; however, when tested in a single-cycle infection assay, some of the mutants demonstrated higher transduction efficiencies than that of wild-type Gag. In particular, the N17K mutant showed a seven- to ninefold increase in transduction, which correlated with enhanced vector RNA packaging. This mutant also packaged larger amounts of foreign RNA. Our results emphasize the importance of the NC zinc fingers, and not other Gag sequences, in achieving specificity in the genome encapsidation process. In addition, the described mutations may contribute to our understanding of HIV diversity resulting from recombination events between copackaged viral genomes and foreign RNA.

The contribution of phosphate-phosphate repulsions to the free energy of DNA bending

DNA bending is important for the packaging of genetic material, regulation of gene expression and interaction of nucleic acids with proteins. Consequently, it is of considerable interest to quantify the energetic factors that must be overcome to induce bending of DNA, such as base stacking and phosphate–phosphate repulsions. In the present work, the electrostatic contribution of phosphate–phosphate repulsions to the free energy of bending DNA is examined for 71 bp linear and bent-form model structures. The bent DNA model was based on the crystallographic structure of a full turn of DNA in a nucleosome core particle. A Green's function approach based on a linear-scaling smooth conductor-like screening model was applied to ascertain the contribution of individual phosphate–phosphate repulsions and overall electrostatic stabilization in aqueous solution. The effect of charge neutralization by site-bound ions was considered using Monte Carlo simulation to characterize the distribution of ion occupations and contribution of phosphate repulsions to the free energy of bending as a function of counterion load. The calculations predict that the phosphate–phosphate repulsions account for ∼30% of the total free energy required to bend DNA from canonical linear B-form into the conformation found in the nucleosome core particle.

Elimination of retroviral infectivity by N-ethylmaleimide with preservation of functional envelope glycoproteins

The zinc finger motifs in retroviral nucleocapsid (NC) proteins are essential for viral replication. Disruption of these Cys-X2-Cys-X4-His-X4-Cys zinc-binding structures eliminates infectivity. To determine if N-ethylmaleimide (NEM) can inactivate human immunodeficiency virus type 1 (HIV-1) or simian immunodeficiency virus (SIV) preparations by alkylating cysteines of NC zinc fingers, we treated infectious virus with NEM and evaluated inactivation of infectivity in cell-based assays. Inactivation was rapid and proportional to the NEM concentration. NEM treatment of HIV-1 or SIV resulted in extensive covalent modification of NC and other internal virion proteins. In contrast, viral envelope glycoproteins, in which the cysteines are disulfide bonded, remained intact and functional, as assayed by high-performance liquid chromatography, fusion-from-without analyses, and dendritic cell capture. Quantitative PCR assays for reverse transcription intermediates showed that NEM and 2,2'-dipyridyl disulfide (aldrithiol-2), a reagent which inactivates retroviruses through oxidation of cysteines in internal virion proteins such as NC, blocked HIV-1 reverse transcription prior to the formation of minus-strand strong-stop products. However, the reverse transcriptase from NEM-treated virions remained active in exogenous template assays, consistent with a role for NC in reverse transcription. Since disruption of NC zinc finger structures by NEM blocks early postentry steps in the retroviral infection cycle, virus preparations with modified NC proteins may be useful as vaccine immunogens and probes of the role of NC in viral replication.

Quantum descriptors for biological macromolecules from linear-scaling electronic structure methods

The characterization of electrostatic and chemical properties at the surface of biological macromolecules is of interest in elucidating the fundamental biological structure-function relationships as well as in problems of rational drug design. This paper presents a set of macromolecular quantum descriptors for the characterization of biological macromolecules in solution that can be obtained with modest computational cost from linear-scaling semi-empirical quantum/solvation methods. The descriptors discussed include: solvent-polarized electrostatic surface potential maps, equilibrated atomic charges, Fukui reactivity indices, approximate local hardness maps, and relative proton potentials. These properties are applied to study the conformational dependence of the electrostatic surface potential of the solvated phosphate-binding protein mutant (T141D), the regioselectivity of the zinc finger domains of HIV-1 nucleocapsid (NC) protein, and the order of pKa values of acidic residues in turkey ovomucoid third domain (OMTKY3) and of the zinc-binding residues in the carboxyl terminal zinc finger of NC. In all cases, insight beyond that obtainable from purely classical models is gained and can be used to rationalize the experimental observations. The macromolecular quantum descriptors presented here greatly extend the arsenal of tools for macromolecular characterization and offer promise in applications to modern structure-based drug design.