Conformational variability of benzamidinium-based inhibitors

Accelerating reliable multiscale quantum refinement of protein-drug systems enabled by machine learning

Biomacromolecule structures are essential for drug development and biocatalysis. Quantum refinement (QR) methods, which employ reliable quantum mechanics (QM) methods in crystallographic refinement, showed promise in improving the structural quality or even correcting the structure of biomacromolecules. However, vast computational costs and complex quantum mechanics/molecular mechanics (QM/MM) setups limit QR applications. Here we incorporate robust machine learning potentials (MLPs) in multiscale ONIOM(QM:MM) schemes to describe the core parts (e.g., drugs/inhibitors), replacing the expensive QM method. Additionally, two levels of MLPs are combined for the first time to overcome MLP limitations. Our unique MLPs+ONIOM-based QR methods achieve QM-level accuracy with significantly higher efficiency. Furthermore, our refinements provide computational evidence for the existence of bonded and nonbonded forms of the Food and Drug Administration (FDA)-approved drug nirmatrelvir in one SARS-CoV-2 main protease structure. This study highlights that powerful MLPs accelerate QRs for reliable protein-drug complexes, promote broader QR applications and provide more atomistic insights into drug development.

Molecular Gas-Phase Conformational Ensembles

Accurately determining the global minima of a molecular structure is important in diverse scientific fields, including drug design, materials science, and chemical synthesis. Conformational search engines serve as valuable tools for exploring the extensive conformational space of molecules and for identifying energetically favorable conformations. In this study, we present a comparison of Auto3D, CREST, Balloon, and ETKDG (from RDKit), which are freely available conformational search engines, to evaluate their effectiveness in locating global minima. These engines employ distinct methodologies, including machine learning (ML) potential-based, semiempirical, and force field-based approaches. To validate these methods, we propose the use of collisional cross-section (CCS) values obtained from ion mobility-mass spectrometry studies. We hypothesize that experimental gas-phase CCS values can provide experimental evidence that we likely have the global minimum for a given molecule. To facilitate this effort, we used our gas-phase conformation library (GPCL) which currently consists of the full ensembles of 20 small molecules and can be used by the community to validate any conformational search engine. Further members of the GPCL can be readily created for any molecule of interest using our standard workflow used to compute CCS values, expanding the ability of the GPCL in validation exercises. These innovative validation techniques enhance our understanding of the conformational landscape and provide valuable insights into the performance of conformational generation engines. Our findings shed light on the strengths and limitations of each search engine, enabling informed decisions for their utilization in various scientific fields, where accurate molecular structure determination is crucial for understanding biological activity and designing targeted interventions. By facilitating the identification of reliable conformations, this study significantly contributes to enhancing the efficiency and accuracy of molecular structure determination, with particular focus on metabolite structure elucidation. The findings of this research also provide valuable insights for developing effective workflows for predicting the structures of unknown compounds with high precision.

Improving Small-Molecule Force Field Parameters in Ligand Binding Studies

Small molecules are major players of many chemical processes in diverse fields, from material science to biology. They are made by a combination of carbon and heteroatoms typically organized in system-specific structures of different complexity. This peculiarity hampers the application of standard force field parameters and their in silico study by means of atomistic simulations. Here, we combine quantum-mechanics and atomistic free-energy calculations to achieve an improved parametrization of the ligand torsion angles with respect to the state-of-the-art force fields in the paradigmatic molecular binding system benzamidine/trypsin. Funnel-Metadynamics calculations with the new parameters greatly reproduced the high-resolution crystallographic ligand binding mode and allowed a more accurate description of the binding mechanism, when the ligand might assume specific conformations to cross energy barriers. Our study impacts on future drug design investigations considering that the vast majority of marketed drugs are small-molecules.

Combining crystallography with quantum mechanics

In standard crystallographic refinement of biomacromolecules, the crystallographic raw data are supplemented by empirical restraints that ensure that the structure makes chemical sense. These restraints are typically accurate for amino acids and nucleic acids, but less so for cofactors, substrates, inhibitors, ligands and metal sites. In quantum refinement, this potential is replaced by more accurate quantum mechanical (QM) calculations. Several implementations have been presented, differing in the level of QM and whether it is used for the entire structure or only for a site of particular interest. It has been shown that the method can improve and correct errors in crystal structures and that it can be used to determine protonation and tautomeric states of various ligands and to decide what is really seen in the structure by refining different interpretations and using standard crystallographic and QM quality measures to decide which fits the structure best.

Intriguing role of water in protein-ligand binding studied by neutron crystallography on trypsin complexes

Hydrogen bonds are key interactions determining protein-ligand binding affinity and therefore fundamental to any biological process. Unfortunately, explicit structural information about hydrogen positions and thus H-bonds in protein-ligand complexes is extremely rare and similarly the important role of water during binding remains poorly understood. Here, we report on neutron structures of trypsin determined at very high resolutions ≤1.5 Å in uncomplexed and inhibited state complemented by X-ray and thermodynamic data and computer simulations. Our structures show the precise geometry of H-bonds between protein and the inhibitors N-amidinopiperidine and benzamidine along with the dynamics of the residual solvation pattern. Prior to binding, the ligand-free binding pocket is occupied by water molecules characterized by a paucity of H-bonds and high mobility resulting in an imperfect hydration of the critical residue Asp189. This phenomenon likely constitutes a key factor fueling ligand binding via water displacement and helps improving our current view on water influencing protein-ligand recognition.

Quantum Crystallography: Current Developments and Future Perspectives

Crystallography and quantum mechanics have always been tightly connected because reliable quantum mechanical models are needed to determine crystal structures. Due to this natural synergy, nowadays accurate distributions of electrons in space can be obtained from diffraction and scattering experiments. In the original definition of quantum crystallography (QCr) given by Massa, Karle and Huang, direct extraction of wavefunctions or density matrices from measured intensities of reflections or, conversely, ad hoc quantum mechanical calculations to enhance the accuracy of the crystallographic refinement are implicated. Nevertheless, many other active and emerging research areas involving quantum mechanics and scattering experiments are not covered by the original definition although they enable to observe and explain quantum phenomena as accurately and successfully as the original strategies. Therefore, we give an overview over current research that is related to a broader notion of QCr, and discuss options how QCr can evolve to become a complete and independent domain of natural sciences. The goal of this paper is to initiate discussions around QCr, but not to find a final definition of the field.

Capturing the photo-signaling state of a photoreceptor in a steady-state fashion by binding a transition metal complex

Binding a small molecule to proteins causes conformational changes, but often to a limited extent. Here, we demonstrate that the interaction of a CO-releasing molecule (CORM3) with a photoreceptor photoactive yellow protein (PYP) drives large structural changes in the latter. The interaction of CORM3 and a mutant of PYP, Met100Ala, not only trigger the isomerization of its chromophore, p-coumaric acid, from its anionic trans configuration to a protonated cis configuration, but also increases the content of β-sheet at the cost of α-helix and random coil in the secondary structure of the protein. The CORM3 derived Met100Ala is found to highly resemble the signaling state, which is one of the key photo-intermediates of this photoactive protein, in both protein local conformation and chromophore configuration. The organometallic reagents hold promise as protein engineering tools. This work highlights a novel approach to structurally accessing short lived intermediates of proteins in a steady-state fashion.

Quantum crystallography

Approximate wavefunctions can be improved by constraining them to reproduce observations derived from diffraction and scattering experiments. Conversely, charge density models, incorporating electron-density distributions, atomic positions and atomic motion, can be improved by supplementing diffraction experiments with quantum chemically calculated, tailor-made electron densities (form factors). In both cases quantum chemistry and diffraction/scattering experiments are combined into a single, integrated tool. The development of quantum crystallographic research is reviewed. Some results obtained by quantum crystallography illustrate the potential and limitations of this field.

Recent Progress in Treating Protein-Ligand Interactions with Quantum-Mechanical Methods

We review the first successes and failures of a “new wave” of quantum chemistry-based approaches to the treatment of protein/ligand interactions. These approaches share the use of “enhanced”, dispersion (D), and/or hydrogen-bond (H) corrected density functional theory (DFT) or semi-empirical quantum mechanical (SQM) methods, in combination with ensemble weighting techniques of some form to capture entropic effects. Benchmark and model system calculations in comparison to high-level theoretical as well as experimental references have shown that both DFT-D (dispersion-corrected density functional theory) and SQM-DH (dispersion and hydrogen bond-corrected semi-empirical quantum mechanical) perform much more accurately than older DFT and SQM approaches and also standard docking methods. In addition, DFT-D might soon become and SQM-DH already is fast enough to compute a large number of binding modes of comparably large protein/ligand complexes, thus allowing for a more accurate assessment of entropic effects.

Fragment quantum mechanical calculation of proteins and its applications

Conspectus The desire to study molecular systems that are much larger than what the current state-of-the-art ab initio or density functional theory methods could handle has naturally led to the development of novel approximate methods, including semiempirical approaches, reduced-scaling methods, and fragmentation methods. The major computational limitation of ab initio methods is the scaling problem, because the cost of ab initio calculation scales nth power or worse with system size. In the past decade, the fragmentation approach based on chemical locality has opened a new door for developing linear-scaling quantum mechanical (QM) methods for large systems and for applications to large molecular systems such as biomolecules. The fragmentation approach is highly attractive from a computational standpoint. First, the ab initio calculation of individual fragments can be conducted almost independently, which makes it suitable for massively parallel computations. Second, the electron properties, such as density and energy, are typically combined in a linear fashion to reproduce those for the entire molecular system, which makes the overall computation scale linearly with the size of the system. In this Account, two fragmentation methods and their applications to macromolecules are described. They are the electrostatically embedded generalized molecular fractionation with conjugate caps (EE-GMFCC) method and the automated fragmentation quantum mechanics/molecular mechanics (AF-QM/MM) approach. The EE-GMFCC method is developed from the MFCC approach, which was initially used to obtain accurate protein-ligand QM interaction energies. The main idea of the MFCC approach is that a pair of conjugate caps (concaps) is inserted at the location where the subsystem is divided by cutting the chemical bond. In addition, the pair of concaps is fused to form molecular species such that the overcounted effect from added concaps can be properly removed. By introducing the electrostatic embedding field in each fragment calculation and two-body interaction energy correction on top of the MFCC approach, the EE-GMFCC method is capable of accurately reproducing the QM molecular properties (such as the dipole moment, electron density, and electrostatic potential), the total energy, and the electrostatic solvation energy from full system calculations for proteins. On the other hand, the AF-QM/MM method was used for the efficient QM calculation of protein nuclear magnetic resonance (NMR) parameters, including the chemical shift, chemical shift anisotropy tensor, and spin-spin coupling constant. In the AF-QM/MM approach, each amino acid and all the residues in its vicinity are automatically assigned as the QM region through a distance cutoff for each residue-centric QM/MM calculation. Local chemical properties of the central residue can be obtained from individual QM/MM calculations. The AF-QM/MM approach precisely reproduces the NMR chemical shifts of proteins in the gas phase from full system QM calculations. Furthermore, via the incorporation of implicit and explicit solvent models, the protein NMR chemical shifts calculated by the AF-QM/MM method are in excellent agreement with experimental values. The applications of the AF-QM/MM method may also be extended to more general biological systems such as DNA/RNA and protein-ligand complexes.

Using quantum mechanical approaches to study biological systems

Conspectus Quantum mechanics (QM) has revolutionized our understanding of the structure and reactivity of small molecular systems. Given the tremendous impact of QM in this research area, it is attractive to believe that this could also be brought into the biological realm where systems of a few thousand atoms and beyond are routine. Applying QM methods to biological problems brings an improved representation to these systems by the direct inclusion of inherently QM effects such as polarization and charge transfer. Because of the improved representation, novel insights can be gleaned from the application of QM tools to biomacromolecules in aqueous solution. To achieve this goal, the computational bottlenecks of QM methods had to be addressed. In semiempirical theory, matrix diagonalization is rate limiting, while in density functional theory or Hartree-Fock theory electron repulsion integral computation is rate-limiting. In this Account, we primarily focus on semiempirical models where the divide and conquer (D&C) approach linearizes the matrix diagonalization step with respect to the system size. Through the D&C approach, a number of applications to biological problems became tractable. Herein, we provide examples of QM studies on biological systems that focus on protein solvation as viewed by QM, QM enabled structure-based drug design, and NMR and X-ray biological structure refinement using QM derived restraints. Through the examples chosen, we show the power of QM to provide novel insights into biological systems, while also impacting practical applications such as structure refinement. While these methods can be more expensive than classical approaches, they make up for this deficiency by the more realistic modeling of the electronic nature of biological systems and in their ability to be broadly applied. Of the tools and applications discussed in this Account, X-ray structure refinement using QM models is now generally available to the community in the refinement package Phenix. While the power of this approach is manifest, challenges still remain. In particular, QM models are generally applied to static structures, so ways in which to include sampling is an ongoing challenge. Car-Parrinello or Born-Oppenheimer molecular dynamics approaches address the short time scale sampling issue, but how to effectively use QM to study phenomenon covering longer time scales will be the focus of future research. Finally, how to accurately and efficiently include electron correlation effects to facilitate the modeling of, for example, dispersive interactions, is also a major hurdle that a broad range of groups are addressing The use of QM models in biology is in its infancy, leading to the expectation that the most significant use of these tools to address biological problems will be seen in the coming years. It is hoped that while this Account summarizes where we have been, it will also help set the stage for future research directions at the interface of quantum mechanics and biology.

Supramolecular association in proton-transfer adducts containing benzamidinium cations. III. Three molecular salts of 3-methoxy-, 4-methoxy- and 3,4,5-trimethoxybenzoates with benzamidine

Three molecular salts, benzamidinium 3-methoxybenzoate, C7H9N2(+)·C8H7O3(-), (I), benzamidinium 4-methoxybenzoate, C7H9N2(+)·C8H7O3(-), (II), and benzamidinium 3,4,5-trimethoxybenzoate monohydrate, C7H9N2(+)·C10H11O5(-)·H2O, (III), were formed from the proton-transfer reactions of 3-methoxy, 4-methoxy- and 3,4,5-trimethoxybenzoic acids with benzamidine (benzenecarboximidamide, benzam). Monoclinic salts (I) and (II) have a 1:1 ratio of cation to anion. In monoclinic salt (III), two cation-anion pairs and two water molecules constitute the asymmetric unit. In all three molecular salts, the amidinium fragments and the carboxylate groups are completely delocalized, and the delocalization favours the aggregation of the molecular components into nonplanar dimers with an R(2)(2)(8) graph-set motif by N(+)-H···O(-) (±) charge-assisted hydrogen bonding (CAHB). Of the three molecular salts, (I) and (II) show similar conformations of the anionic components and exhibit bidimensional isostructurality, which consists of alternating R(2)(2)(8) and R(6)(4)(16) rings resulting in a corrugated sheet propagated parallel to the crystallographic ab plane. In molecular salt (III), the R(2)(2)(8) synthon is retained but the supramolecular structure is different, due to the presence of three bulky methoxy substituents and a water molecule. The structures reported here further demonstrate the robustness of R(2)(2)(8) hydrogen-bonded synthons having the benzamidinium cation as a building block, whereas N(+)-H···O(-) hydrogen bonds external to the salt bridge contribute to the overall structure organization.

Automating crystallographic structure solution and refinement of protein-ligand complexes

High-throughput drug-discovery and mechanistic studies often require the determination of multiple related crystal structures that only differ in the bound ligands, point mutations in the protein sequence and minor conformational changes. If performed manually, solution and refinement requires extensive repetition of the same tasks for each structure. To accelerate this process and minimize manual effort, a pipeline encompassing all stages of ligand building and refinement, starting from integrated and scaled diffraction intensities, has been implemented in Phenix. The resulting system is able to successfully solve and refine large collections of structures in parallel without extensive user intervention prior to the final stages of model completion and validation.

Conformational analysis and parallel QM/MM X-ray refinement of protein bound anti-Alzheimer drug donepezil

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.

Supramolecular association in proton-transfer adducts containing benzamidinium cations. II. Concomitant polymorphs of the molecular salt of 2,6-dimethoxybenzoic acid with benzamidine

Two concomitant polymorphs of the molecular salt formed by 2,6-dimethoxybenzoic acid, C(9)H(10)O(4) (Dmb), with benzamidine, C(7)H(8)N(2) (benzenecarboximidamide, Benzam) from water solution have been identified. Benzamidinidium 2,6-dimethoxybenzoate, C(7)H(9)N(2)(+)·C(9)H(9)O(4)(-) (BenzamH(+)·Dmb(-)), was obtained through protonation at the imino N atom of Benzam as a result of proton transfer from the acidic hydroxy group of Dmb. In the monoclinic polymorph, (I) (space group P2(1)/n), the asymmetric unit consists of two Dmb(-) anions and two monoprotonated BenzamH(+) cations. In the orthorhombic polymorph, (II) (space group P2(1)2(1)2(1)), one Dmb(-) anion and one BenzamH(+) cation constitute the asymmetric unit. In both polymorphic salts, the amidinium fragments and carboxylate groups are completely delocalized. This delocalization favours the aggregation of the molecular components of these acid-base complexes into nonplanar dimers with an R(2)(2)(8) graph-set motif via N(+)-H···O(-) charge-assisted hydrogen bonding. Both the monoclinic and orthorhombic forms exhibit one-dimensional isostructurality, as the crystal structures feature identical hydrogen-bonding motifs consisting of dimers and catemers.

The central role of Gln63 for the hydrogen bonding network and UV-visible spectrum of the AppA BLUF domain

In blue-light sensing using flavin (BLUF) domains, the side-chain orientation of key residues close to the flavin chromophore is still under debate. We report quantum refinements of the wild-type AppA BLUF protein from Rhodobacter sphaeroides starting from two published X-ray structures (1YRX and 2IYG) with different arrangements of the residues around the chromophore. Quantum refinement uses the same experimental X-ray raw data as conventional refinement, but includes data from quantum mechanics/molecular mechanics (QM/MM) calculations as restraints, which is expected to be more reliable than the normally employed MM data. In addition to quantum refinement, pure QM/MM geometry optimizations are performed for the 1YRX and 2IYG structures and for five models derived therefrom. Vertical excitation energies are computed at the QM(DFT/MRCI)/MM level to assess the resulting structures. The experimental absorption maximum of the dark state of wild-type AppA is well reproduced for structures that contain the Gln63 residue in 1YRX-type orientation. The computed excitation energies are red-shifted for structures with a flipped Gln63 residue in 2IYG-type orientation. The calculated 1YRX- and 2IYG-type hydrogen-bonding networks are discussed in detail, particularly with regard to the orientation of the chromophore and the Gln63, Trp104, and Met106 residues.

PDB ligand conformational energies calculated quantum-mechanically

We present here a greatly updated version of an earlier study on the conformational energies of protein-ligand complexes in the Protein Data Bank (PDB) [Nicklaus et al. Bioorg. Med. Chem. 1995, 3, 411-428], with the goal of improving on all possible aspects such as number and selection of ligand instances, energy calculations performed, and additional analyses conducted. Starting from about 357,000 ligand instances deposited in the 2008 version of the Ligand Expo database of the experimental 3D coordinates of all small-molecule instances in the PDB, we created a "high-quality" subset of ligand instances by various filtering steps including application of crystallographic quality criteria and structural unambiguousness. Submission of 640 Gaussian 03 jobs yielded a set of about 415 successfully concluded runs. We used a stepwise optimization of internal degrees of freedom at the DFT level of theory with the B3LYP/6-31G(d) basis set and a single-point energy calculation at B3LYP/6-311++G(3df,2p) after each round of (partial) optimization to separate energy changes due to bond length stretches vs bond angle changes vs torsion changes. Even for the most "conservative" choice of all the possible conformational energies-the energy difference between the conformation in which all internal degrees of freedom except torsions have been optimized and the fully optimized conformer-significant energy values were found. The range of 0 to ~25 kcal/mol was populated quite evenly and independently of the crystallographic resolution. A smaller number of "outliers" of yet higher energies were seen only at resolutions above 1.3 Å. The energies showed some correlation with molecular size and flexibility but not with crystallographic quality metrics such as the Cruickshank diffraction-component precision index (DPI) and R(free)-R, or with the ligand instance-specific metrics such as occupancy-weighted B-factor (OWAB), real-space R factor (RSR), and real-space correlation coefficient (RSCC). We repeated these calculations with the solvent model IEFPCM, which yielded energy differences that were generally somewhat lower than the corresponding vacuum results but did not produce a qualitatively different picture. Torsional sampling around the crystal conformation at the molecular mechanics level using the MMFF94s force field typically led to an increase in energy.

Conformational Analysis of Free and Bound Retinoic Acid

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.

QM/MM refinement and analysis of protein bound retinoic acid

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.

Complete reconstruction of an enzyme-inhibitor binding process by molecular dynamics simulations

The understanding of protein-ligand binding is of critical importance for biomedical research, yet the process itself has been very difficult to study because of its intrinsically dynamic character. Here, we have been able to quantitatively reconstruct the complete binding process of the enzyme-inhibitor complex trypsin-benzamidine by performing 495 molecular dynamics simulations of free ligand binding of 100 ns each, 187 of which produced binding events with an rmsd less than 2 Å compared to the crystal structure. The binding paths obtained are able to capture the kinetic pathway of the inhibitor diffusing from solvent (S0) to the bound (S4) state passing through two metastable intermediate states S2 and S3. Rather than directly entering the binding pocket the inhibitor appears to roll on the surface of the protein in its transition between S3 and the final binding pocket, whereas the transition between S2 and the bound pose requires rediffusion to S3. An estimation of the standard free energy of binding gives ΔG° = -5.2 ± 0.4 kcal/mol (cf. the experimental value -6.2 kcal/mol), and a two-states kinetic model k(on) = (1.5 ± 0.2) × 10(8) M(-1) s(-1) and k(off) = (9.5 ± 3.3) × 10(4) s(-1) for unbound to bound transitions. The ability to reconstruct by simple diffusion the binding pathway of an enzyme-inhibitor binding process demonstrates the predictive power of unconventional high-throughput molecular simulations. Moreover, the methodology is directly applicable to other molecular systems and thus of general interest in biomedical and pharmaceutical research.

Challenges in the determination of the binding modes of non-standard ligands in X-ray crystal complexes

Despite its central role in structure based drug design the determination of the binding mode (position, orientation and conformation in addition to protonation and tautomeric states) of small heteromolecular ligands in protein:ligand complexes based on medium resolution X-ray diffraction data is highly challenging. In this perspective we demonstrate how a combination of molecular dynamics simulations and free energy (FE) calculations can be used to correct and identify thermodynamically stable binding modes of ligands in X-ray crystal complexes. The consequences of inappropriate ligand structure, force field and the absence of electrostatics during X-ray refinement are highlighted. The implications of such uncertainties and errors for the validation of virtual screening and fragment-based drug design based on high throughput X-ray crystallography are discussed with possible solutions and guidelines.

Supramolecular association in proton-transfer adducts containing benzamidinium cations. I. Four molecular salts with uracil derivatives

Four organic salts, namely benzamidinidium orotate (2,6-dioxo-1,2,3,6-tetrahydropyrimidine-4-carboxylate) hemihydrate, C(7)H(9)N(2)(+).C(5)H(3)N(2)O(4)(-).0.5H(2)O (BenzamH(+).Or(-)), (I), benzamidinium isoorotate (2,4-dioxo-1,2,3,4-tetrahydropyrimidine-5-carboxylate) trihydrate, C(7)H(9)N(2)(+).C(5)H(3)N(2)O(4)(-).3H(2)O (BenzamH(+).Isor(-)), (II), benzamidinium diliturate (5-nitro-2,6-dioxo-1,2,3,6-tetrahydropyrimidin-4-olate) dihydrate, C(7)H(9)N(2)(+).C(4)H(2)N(3)O(5)(-).2H(2)O (BenzamH(+).Dil(-)), (III), and benzamidinium 5-nitrouracilate (5-nitro-2,4-dioxo-1,2,3,4-tetrahydropyrimidin-1-ide), C(7)H(9)N(2)(+).C(4)H(2)N(3)O(4)(-) (BenzamH(+).Nit(-)), (IV), have been synthesized by a reaction between benzamidine (benzenecarboximidamide or Benzam) and the appropriate carboxylic acid. Proton transfer occurs to the benzamidine imino N atom. In all four acid-base adducts, the asymmetric unit consists of one tautomeric aminooxo anion (Or(-), Isor(-), Dil(-) and Nit(-)) and one monoprotonated benzamidinium cation (BenzamH(+)), plus one-half (which lies across a twofold axis), three and two solvent water molecules in (I), (II) and (III), respectively. Due to the presence of protonated benzamidine, these acid-base complexes form supramolecular synthons characterized by N(+)-H...O(-) and N(+)-H...N(-) (+/-)-charge-assisted hydrogen bonds (CAHB).