Conformational Analysis of Free and Bound Retinoic Acid.
ABSTRACT: The conformational profiles of unbound all-trans and 9-cis retinoic acid (RA) have been determined using classical and quantum mechanical calculations. Sixty-six all-trans-RA (ATRA) and forty-eight 9-cis-RA energy minimum conformers were identified via HF/6-31G* geometry optimizations in vacuo. Their relative conformational energies were estimated utilizing the M06, M06-2x and MP2 methods combined with the 6-311+G(d,p), aug-cc-pVDZ and aug-cc-pVTZ basis sets, as well as complete basis set MP2 extrapolations using the latter two basis sets. Single-point energy calculations performed with the M06-2x density functional were found to yield similar results to MP2/CBS for the low-energy retinoic acid conformations. Not unexpectedly, the conformational propensities of retinoic acid were governed by the orientation and arrangement of the torsion angles associated with the polyene tail. We also used previously reported QM/MM X-ray refinement results on four ATRA-protein crystal structures plus one newly refined 9-cis-RA complex (PDB ID 1XDK) in order to investigate the conformational preferences of bound retinoic acid. In the re-refined RA conformers the conjugated double bonds are nearly coplanar, which is consistent with the global minimum identified by the Omega/QM method rather than the corresponding crystallographically determined conformations given in the PDB. Consequently, a 91.3% average reduction of the local strain energy in the gas phase, as well as 92.1% in PCM solvent, was observed using the QM/MM refined structures versus the PDB deposited RA conformations. These results thus demonstrate that our QM/MM X-ray refinement approach can significantly enhance the quality of X-ray crystal structures refined by conventional refinement protocols, thereby providing reliable drug-target structural information for use in structure-based drug discovery applications.
Project description:Retinoic acid (RA) is a vitamin A derivative, which modifies the appearance of fine wrinkles and roughness of facial skin and treats acne and activates gene transcription by binding to heterodimers of the retinoic acid receptor (RAR) and the retinoic X receptor (RXR). There are series of protein bound RA complexes available in the protein databank (PDB), which provides a broad range of information about the different bioactive conformations of RA. To gain further insights into the observed bioactive RA conformations we applied quantum mechanic (QM)/molecular mechanic (MM) approaches to re-refine the available RA protein-ligand complexes. MP2 complete basis set (CBS) extrapolations single energy calculations are also carried out for both the experimental conformations and QM optimized geometries of RA in the gas as well as solution phase. The results demonstrate that the re-refined structures show better geometries for RA than seen in the originally deposited PDB structures through the use of QMs for the ligand in the X-ray refinement procedure. QM/MM re-refined conformations also reduced the computed strain energies found in the deposited crystal conformations for RA. Finally, the dependence of ligand strain on resolution is analyzed. It is shown that ligand strain is not converged in our calculations and is likely an artifact of the typical resolutions employed to study protein-ligand complexes.
Project description:The recognition and association of donepezil with acetylcholinesterase (AChE) has been extensively studied in the past several decades because of the former's use as a palliative treatment for mild Alzheimer disease. Herein we examine the conformational properties of donepezil and we re-examine the donepezil-AChE crystal structure using combined quantum mechanical/molecular mechanical (QM/MM) X-ray refinement tools. Donepezil's conformational energy surface was explored using the M06 suite of density functionals and with the MP2/complete basis set (CBS) method using the aug-cc-pVXZ (X = D and T) basis sets. The donepezil-AChE complex (PDB 1EVE) was also re-refined through a parallel QM/MM X-ray refinement approach based on an in-house ab initio code QUICK, which uses the message passing interface (MPI) in a distributed SCF algorithm to accelerate the calculation via parallelization. In the QM/MM re-refined donepezil structure, coordinate errors that previously existed in the PDB deposited geometry were improved leading to an improvement of the modeling of the interaction between donepezil and the aromatic side chains located in the AChE active site gorge. As a result of the re-refinement there was a 93% reduction in the donepezil conformational strain energy versus the original PDB structure. The results of the present effort offer further detailed structural and biochemical inhibitor-AChE information for the continued development of more effective and palliative treatments of Alzheimer disease.
Project description:Amino acid side chains adopt a discrete set of favorable conformations typically referred to as rotamers. The relative energies of rotamers partially determine which side chain conformations are more often observed in protein structures and accurate estimates of these energies are important for predicting protein structure and designing new proteins. Protein modelers typically calculate side chain rotamer energies by using molecular mechanics (MM) potentials or by converting rotamer probabilities from the protein database (PDB) into relative free energies. One limitation of the knowledge-based energies is that rotamer preferences observed in the PDB can reflect internal side chain energies as well as longer-range interactions with the rest of the protein. Here, we test an alternative approach for calculating rotamer energies. We use three different quantum mechanics (QM) methods (second order Møller-Plesset (MP2), density functional theory (DFT) energy calculation using the B3LYP functional, and Hartree-Fock) to calculate the energy of amino acid rotamers in a dipeptide model system, and then use these pre-calculated values in side chain placement simulations. Energies were calculated for over 36,000 different conformations of leucine, isoleucine, and valine dipeptides with backbone torsion angles from the helical and strand regions of the Ramachandran plot. In a subset of cases these energies differ significantly from those calculated with standard molecular mechanics potentials or those derived from PDB statistics. We find that in these cases the energies from the QM methods result in more accurate placement of amino acid side chains in structure prediction tests.
Project description:Development of the AMOEBA (Atomic Multipole Optimized Energetics for Biomolecular Simulation) force field for proteins is presented. The current version (AMOEBA-2013) utilizes permanent electrostatic multipole moments through the quadrupole at each atom, and explicitly treats polarization effects in various chemical and physical environments. The atomic multipole electrostatic parameters for each amino acid residue type are derived from high-level gas phase quantum mechanical calculations via a consistent and extensible protocol. Molecular polarizability is modeled via a Thole-style damped interactive induction model based upon distributed atomic polarizabilities. Inter- and intramolecular polarization is treated in a consistent fashion via the Thole model. The intramolecular polarization model ensures transferability of electrostatic parameters among different conformations, as demonstrated by the agreement between QM and AMOEBA electrostatic potentials, and dipole moments of dipeptides. The backbone and side chain torsional parameters were determined by comparing to gas-phase QM (RI-TRIM MP2/CBS) conformational energies of dipeptides and to statistical distributions from the Protein Data Bank. Molecular dynamics simulations are reported for short peptides in explicit water to examine their conformational properties in solution. Overall the calculated conformational free energies and J-coupling constants are consistent with PDB statistics and experimental NMR results, respectively. In addition, the experimental crystal structures of a number of proteins are well maintained during molecular dynamics (MD) simulation. While further calculations are necessary to fully validate the force field, initial results suggest the AMOEBA polarizable multipole force field is able to describe the structure and energetics of peptides and proteins, in both gas-phase and solution environments.
Project description:Macromolecular crystallographic refinement relies on sometimes dubious stereochemical restraints and rudimentary energy functionals to ensure the correct geometry of the model of the macromolecule and any covalently bound ligand(s). The ligand stereochemical restraint file (CIF) requires a priori understanding of the ligand geometry within the active site, and creation of the CIF is often an error-prone process owing to the great variety of potential ligand chemistry and structure. Stereochemical restraints have been replaced with more robust functionals through the integration of the linear-scaling, semiempirical quantum-mechanics (SE-QM) program DivCon with the PHENIX X-ray refinement engine. The PHENIX/DivCon package has been thoroughly validated on a population of 50 protein-ligand Protein Data Bank (PDB) structures with a range of resolutions and chemistry. The PDB structures used for the validation were originally refined utilizing various refinement packages and were published within the past five years. PHENIX/DivCon does not utilize CIF(s), link restraints and other parameters for refinement and hence it does not make as many a priori assumptions about the model. Across the entire population, the method results in reasonable ligand geometries and low ligand strains, even when the original refinement exhibited difficulties, indicating that PHENIX/DivCon is applicable to both single-structure and high-throughput crystallography.
Project description:The conformational preferences of the cationic nylon-3 ?NM [(3R,4)-diaminobutanoic acid, dAba] dipeptide in water were explored as the first step to understand the mode of action of polymers of ?NM against phylogenetically diverse and intrinsically drug-resistant pathogenic fungi. The CCSD(T), MP2, M06-2X, ?B97X-D, B2PLYP-D3BJ, and DSD-PBEP86-D3BJ levels of theory with various basis sets were assessed for relative energies of the 45 local minima of the cationic Ac-dAba-NHMe located at the SMD M06-2X/6-31+G(d) level of theory in water against the benchmark CCSD(T)/CBS-limit energies in water. The best performance was obtained at the double-hybrid DSD-PBEP86-D3BJ/def2-QZVP level of theory with RMSD = 0.12 kcal/mol in water. The M06-2X/def2-QZVP level of theory predicted reasonably the conformational preference with RMSD = 0.38 kcal/mol in water and may be an alternative level of theory with marginal deviations for the calculation of conformational energies of relatively longer cationic peptides in water. In particular, the H 14-helical structures appeared to be the most feasible conformations for the cationic Ac-dAba-NHMe populated at 48-64% by relative free energies in water. The hexamer built from the H 14-structure of the cationic Ac-dAba-NHMe adopted a left-handed 314-helix, which has a slightly narrower radius and a longer rise than the regular 314-helix of ?-peptides. Hence, the 314-helices of oligomers or polymers of the cationic dAba residues are expected to be the active conformation to exhibit the ability to bridge between charged lipid head groups that might cause a local depression or invagination of the membrane of fungi.
Project description:Presented is an extension of the CHARMM additive all-atom carbohydrate force field to enable the modeling of phosphate and sulfate linked to carbohydrates. The parameters are developed in a hierarchical fashion using model compounds containing the key atoms in the full carbohydrates. Target data for parameter optimization included full two-dimensional energy surfaces defined by the glycosidic dihedral angle pairs in the phosphate/sulfate model compound analogs of hexopyranose monosaccharide phosphates and sulfates, as determined by quantum mechanical (QM) MP2/cc-pVTZ single point energies on MP2/6-31+G(d) optimized structures. In order to achieve balanced, transferable dihedral parameters for the dihedral angles, surfaces for all possible anomeric and conformational states were included during the parametrization process. In addition, to model physiologically relevant systems both the mono- and di-anionic charged states were studied for the phosphates. This resulted in over 7000 MP2/cc-pVTZ//MP2/6-31G+(d) model compound conformational energies which, supplemented with QM geometries, were the main target data for the parametrization. Parameters were validated against crystals of relevant monosaccharide derivatives obtained from the Cambridge Structural Database (CSD) and larger systems, namely inositol-(tri/tetra/penta) phosphates non-covalently bound to the pleckstrin homology (PH) domain and oligomeric chondroitin sulfate in solution and in complex with cathepsin K protein.
Project description:We have found that refinement of protein NMR structures using Rosetta with experimental NMR restraints yields more accurate protein NMR structures than those that have been deposited in the PDB using standard refinement protocols. Using 40 pairs of NMR and X-ray crystal structures determined by the Northeast Structural Genomics Consortium, for proteins ranging in size from 5-22 kDa, restrained Rosetta refined structures fit better to the raw experimental data, are in better agreement with their X-ray counterparts, and have better phasing power compared to conventionally determined NMR structures. For 37 proteins for which NMR ensembles were available and which had similar structures in solution and in the crystal, all of the restrained Rosetta refined NMR structures were sufficiently accurate to be used for solving the corresponding X-ray crystal structures by molecular replacement. The protocol for restrained refinement of protein NMR structures was also compared with restrained CS-Rosetta calculations. For proteins smaller than 10 kDa, restrained CS-Rosetta, starting from extended conformations, provides slightly more accurate structures, while for proteins in the size range of 10-25 kDa the less CPU intensive restrained Rosetta refinement protocols provided equally or more accurate structures. The restrained Rosetta protocols described here can improve the accuracy of protein NMR structures and should find broad and general for studies of protein structure and function.
Project description:We evaluate the performance of the most widely used wavefunction, density functional theory, and semiempirical methods for the description of noncovalent interactions in a set of larger, mostly dispersion-stabilized noncovalent complexes (the L7 data set). The methods tested include MP2, MP3, SCS-MP2, SCS(MI)-MP2, MP2.5, MP2.X, MP2C, DFT-D, DFT-D3 (B3-LYP-D3, B-LYP-D3, TPSS-D3, PW6B95-D3, M06-2X-D3) and M06-2X, and semiempirical methods augmented with dispersion and hydrogen bonding corrections: SCC-DFTB-D, PM6-D, PM6-DH2 and PM6-D3H4. The test complexes are the octadecane dimer, the guanine trimer, the circumcoronene…adenine dimer, the coronene dimer, the guanine-cytosine dimer, the circumcoronene…guanine-cytosine dimer, and an amyloid fragment trimer containing phenylalanine residues. The best performing method is MP2.5 with relative root mean square deviation (rRMSD) of 4 %. It can thus be recommended as an alternative to the CCSD(T)/CBS (alternatively QCISD(T)/CBS) benchmark for molecular systems which exceed current computational capacity. The second best non-DFT method is MP2C with rRMSD of 8 %. A method with the most favorable "accuracy/cost" ratio belongs to the DFT family: BLYP-D3, with an rRMSD of 8 %. Semiempirical methods deliver less accurate results (the rRMSD exceeds 25 %). Nevertheless, their absolute errors are close to some much more expensive methods such as M06-2X, MP2 or SCS(MI)-MP2, and thus their price/performance ratio is excellent.
Project description:A re-refinement of 4xan, hen egg-white lysozyme (HEWL) with carboplatin crystallized in NaBr solution, has been made and is published here as an addendum to Tanley et al. [(2014), Acta Cryst. F70, 1135-1142]. This follows a previous re-refinement and PDB deposition (4yem) by Shabalin et al. [(2015), Acta Cryst. D71, 1965-1979]. The critical evaluation of the original PDB deposition (4xan), and the subsequent critical examination of the re-refined structure (4yem), has led to an improved model (PDB code 5hmj).