High-resolution crystallography and drug design

The relationship between Hirshfeld potential and cytotoxic activity: a study along a series of flavonoid and chromanone derivatives

The present study examines a series of six biologically-active flavonoid and chromanone derivatives by X-ray crystal structure analysis: (E)-3-benzylidene-2-phenylchroman-4-one, C₂₂H₁₆O₂, I, (E)-3-(4-methylbenzylidene)-2-phenylchroman-4-one, C₂₃H₁₈O₂, II, (E)-3-(3-methylbenzylidene)-2-phenylchroman-4-one, C₂₃H₁₈O₂, III, (E)-3-(4-methoxybenzylidene)-2-phenylchroman-4-one, C₂₃H₁₈O₃, IV, (E)-3-benzylidenechroman-4-one, C₁₆H₁₂O₂, V, and (E)-3-(4-methoxybenzylidene)chroman-4-one, C₁₇H₁₄O₃, VI. The cytotoxic activities of the presented crystal structures have been determined, together with their intermolecular interaction preferences and Hirshfeld surface characteristics. An inverse relationship was found between the contribution of C...C close contacts to the Hirshfeld surface and cytotoxic activity against the WM-115 cancer line. Dependence was also observed between the logP value and the percentage contribution of C...H contacts to the Hirshfeld surface.

Combined quantum mechanics/molecular mechanics (QM/MM) methods to understand the charge density distribution of estrogens in the active site of estrogen receptors

The ligand binding to protein and host–guest interactions are ubiquitous for molecular recognition. In drug design, the ligand binding to the active site of proteins is influenced by the charge density distribution and the electrostatic interactions of ligands and the nearby amino acids of the protein. The charge density analyses of ligand–protein complexes need accurate positions of hydrogen atoms and their valence electron distribution and the fine structure of proteins. Such information cannot be obtained from the conventional protein X-ray crystallography analysis in the resolution range of 1.5 to 3.5 Å. This can be realized from QM/MM based structure and charge density analysis of estrogens with the estrogen receptor. The charge density properties such as electron density, Laplacian of electron density and electrostatic properties of estrogens in the presence of active site amino acid residues have been determined and compared with the isolated estrogen molecules from theory and experimental. The present study reveals the chemical bonding nature of estrogen molecules and the strength of the intermolecular interactions in the active site of estrogen receptor, and also the importance of π⋯π interactions between the estrogens and Phe404 amino acid residue and protonation state of His524 amino acid residue have been identified using electrostatic potential maps. The difference in the electrostatic potential map of estrogens displays the hormone dependent actions of estrogen receptor. This method is very helpful to derive the charge density distribution of macromolecules to understand their biological recognition and interactions.

The ligand binding to protein and host–guest interactions are ubiquitous for molecular recognition.

Tailored Robust Hydrogel Composite Membranes for Continuous Protein Crystallization with Ultrahigh Morphology Selectivity

The tailored and robust hydrogel composite membranes (HCMs) with diverse ion adsorption and interfacial nucleation property are prepared and successfully used in the continuous lysozyme crystallization. Beyond the heterogeneous supporter, the HCMs functioning as an interface ion concentration controller and nucleation generator are demonstrated. By constructing accurately controlled nucleation and growth circumstances in the HCM-equipped membrane crystallizer, the target desired morphology (hexagon cube) and brand-new morphology (multiple flower shape) that differ from the ones created in the conventional crystallizer are continuously and repetitively generated with ultrahigh morphology selectivity. These tailored robust HCMs show great potential for improving current approaches to continuous protein crystallization with specific crystal targets from laboratorial research to actual engineering applications.

Conformational energy range of ligands in protein crystal structures: The difficult quest for accurate understanding

In this review, we address a fundamental question: What is the range of conformational energies seen in ligands in protein-ligand crystal structures? This value is important biophysically, for better understanding the protein-ligand binding process; and practically, for providing a parameter to be used in many computational drug design methods such as docking and pharmacophore searches. We synthesize a selection of previously reported conflicting results from computational studies of this issue and conclude that high ligand conformational energies really are present in some crystal structures. The main source of disagreement between different analyses appears to be due to divergent treatments of electrostatics and solvation. At the same time, however, for many ligands, a high conformational energy is in error, due to either crystal structure inaccuracies or incorrect determination of the reference state. Aside from simple chemistry mistakes, we argue that crystal structure error may mainly be because of the heuristic weighting of ligand stereochemical restraints relative to the fit of the structure to the electron density. This problem cannot be fixed with improvements to electron density fitting or with simple ligand geometry checks, though better metrics are needed for evaluating ligand and binding site chemistry in addition to geometry during structure refinement. The ultimate solution for accurately determining ligand conformational energies lies in ultrahigh-resolution crystal structures that can be refined without restraints.

Measuring the orientation of taurine in the active site of the non-heme Fe(II)/α-ketoglutarate-dependent taurine hydroxylase (TauD) using electron spin echo envelope modulation (ESEEM) spectroscopy

The position and orientation of taurine near the non-heme Fe(II) center of the α-ketoglutarate (α-KG)-dependent taurine hydroxylase (TauD) was measured using Electron Spin Echo Envelope Modulation (ESEEM) spectroscopy. TauD solutions containing Fe(II), α-KG, and natural abundance taurine or specifically deuterated taurine were prepared anaerobically and treated with nitric oxide (NO) to make an S = 3/2 {FeNO}(7) complex that is suitable for robust analysis with EPR spectroscopy. Using ratios of ESEEM spectra collected for TauD samples having natural abundance taurine or deuterated taurine, (1)H and (14)N modulations were filtered out of the spectra and interactions with specific deuterons on taurine could be studied separately. The Hamiltonian parameters used to calculate the amplitudes and line shapes of frequency spectra containing isolated deuterium ESEEM were obtained with global optimization algorithms. Additional statistical analysis was performed to validate the interpretation of the optimized parameters. The strongest (2)H hyperfine coupling was to a deuteron on the C1 position of taurine and was characterized by an effective dipolar distance of 3.90 ± 0.25 Å from the {FeNO}(7) paramagnetic center. The principal axes of this C1-(2)H hyperfine coupling and nuclear quadrupole interaction tensors were found to make angles of 26 ± 5 and 52 ± 17°, respectively, with the principal axis of the {FeNO}(7) zero-field splitting tensor. These results are discussed within the context of the orientation of substrate taurine prior to the initiation of hydrogen abstraction.

An overview of biological macromolecule crystallization

The elucidation of the three dimensional structure of biological macromolecules has provided an important contribution to our current understanding of many basic mechanisms involved in life processes. This enormous impact largely results from the ability of X-ray crystallography to provide accurate structural details at atomic resolution that are a prerequisite for a deeper insight on the way in which bio-macromolecules interact with each other to build up supramolecular nano-machines capable of performing specialized biological functions. With the advent of high-energy synchrotron sources and the development of sophisticated software to solve X-ray and neutron crystal structures of large molecules, the crystallization step has become even more the bottleneck of a successful structure determination. This review introduces the general aspects of protein crystallization, summarizes conventional and innovative crystallization methods and focuses on the new strategies utilized to improve the success rate of experiments and increase crystal diffraction quality.

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.

Complexomics study of two Helicobacter pylori strains of two pathological origins: potential targets for vaccine development and new insight in bacteria metabolism

Helicobacter pylori infection plays a causal role in the development of gastric mucosa-associated lymphoid tissue (MALT) lymphoma (LG-MALT) and duodenal ulcer (DU). Although many virulence factors have been associated with DU, many questions remain unanswered regarding the evolution of the infection toward this exceptional event, LG-MALT. The present study describes and compares the complexome of two H. pylori strains, strain J99 associated with DU and strain B38 associated with LG-MALT, using the two-dimensional blue native/SDS-PAGE method. It was possible to identify 90 different complexes (49 and 41 in the B38 and J99 strains, respectively); 12 of these complexes were common to both strains (seven and five in the membrane and cytoplasm, respectively), reflecting the variability of H. pylori strains. The 44 membrane complexes included numerous outer membrane proteins, such as the major adhesins BabA and SabA retrieved from a complex in the B38 strain, and also proteins from the hor family rarely studied. BabA and BabB adhesins were found to interact independently with HopM/N in the B38 and J99 strains, respectively. The 46 cytosolic complexes essentially comprised proteins involved in H. pylori physiology. Some orphan proteins were retrieved from heterooligomeric complexes, and a function could be proposed for a number of them via the identification of their partners, such as JHP0119, which may be involved in the flagellar function. Overall, this study gave new insights into the membrane and cytoplasm structure, and those which could help in the design of molecules for vaccine and/or antimicrobial agent development are highlighted.

Internal motion in protein crystal structures

The binding states of the substrates and the environment have significant influence on protein motion. We present the analysis of such motion derived from anisotropic atomic displacement parameters (ADPs) in a set of atomic resolution protein structures. Local structural motion caused by ligand binding as well as functional loops showing cooperative patterns of motion could be inferred. The results are in line with proposed protonation states, hydrogen bonding patterns and the location of distinctly flexible regions: we could locate the mobile active site loop in a virus integrase, distinguish the subdomains in RNAse A and hydroxynitrile lyase, and reconstruct the molecular architecture in a xylanase. We demonstrate that the ADP-based motion analysis provides information at high level of detail and that the structural changes needed for substrate attachment or release may be derived from single X-ray structures.

Environmental effects on charge densities of biologically active molecules: do molecule crystal environments indeed approximate protein surroundings?

In the present paper, we investigate whether crystal and enzyme environments influence the electron density (ED) of active compounds in a similar manner. This supposition is essential for high-resolution X-ray studies, which use the EDs obtained from crystals of the pure active compound as approximations for the ED of the active compound in its complex with the target enzyme. The EDs of such complexes determine the molecular recognition process between the targeted enzyme and active compound and are, hence, extremely useful tools for rational drug design. The approximation of such EDs by data obtained from crystals of the pure active compound is needed since high-resolution X-ray experiments of the target-ligand complexes are still extremely demanding. Quantum mechanical/molecular mechanical (QM/MM) and pure QM calculations are employed to determine the EDs of two inhibitors, the reversible trans-4-(aminomethyl)cyclohexane-1-carboxylic acid (AMCHA) and the irreversible E64c in four different environments (the enzyme-inhibitor complex, crystals of the pure compounds, a continuum solvation model, and the gas phase). Our investigation shows that the environment inside of the crystal of the pure active compound generally influences the ED of an active compound in a very similar way as the enzyme surrounding in the complex between the active compound and target enzyme. However, this does not hold any more if the geometrical arrangement of the inhibitor in the enzyme differs significantly from that in the crystal. While EDs computed for gas-phase environments deviate strongly from those in crystal and protein surroundings, polar solvent environments provide rather similar electron distributions. Thus, such continuum solvation models are very well suited to compute density databases which are to be employed for the determination of the ED of macromolecules.

X-ray crystallography of chemical compounds

AIMS

Accurate knowledge of molecular structure is a prerequisite for rational drug design. This review examines the role of X-ray crystallography in providing the required structural information and advances in the field of X-ray crystallography that enhance or expand its role.

MAIN METHODS

X-ray crystallography of new drugs candidates and intermediates can provide valuable information of new syntheses and parameters for quantitative structure activity relationships (QSAR).

KEY FINDINGS

Crystallographic studies play a vital role in many disciplines including materials science, chemistry, pharmacology, and molecular biology. X-ray crystallography is the most comprehensive technique available to determine molecular structure. A requirement for the high accuracy of crystallographic structures is that a 'good crystal' must be found, and this is often the rate-limiting step. In the past three decades developments in detectors, increases in computer power, and powerful graphics capabilities have contributed to a dramatic increase in the number of materials characterized by X-ray crystallography. More recently the advent of high-throughput crystallization techniques has enhanced our ability to produce that one good crystal required for crystallographic analysis.

SIGNIFICANCE

Continuing advances in all phases of a crystallographic study have expanded the ranges of samples which can be analyzes by X-ray crystallography to include larger molecules, smaller or weakly diffracting crystals, and twinned crystals.

Testing geometrical discrimination within an enzyme active site: constrained hydrogen bonding in the ketosteroid isomerase oxyanion hole

Enzymes are classically proposed to accelerate reactions by binding substrates within active-site environments that are structurally preorganized to optimize binding interactions with reaction transition states rather than ground states. This is a remarkably formidable task considering the limited 0.1-1 A scale of most substrate rearrangements. The flexibility of active-site functional groups along the coordinate of substrate rearrangement, the distance scale on which enzymes can distinguish structural rearrangement, and the energetic significance of discrimination on that scale remain open questions that are fundamental to a basic physical understanding of enzyme active sites and catalysis. We bring together 1.2-1.5 A resolution X-ray crystallography, (1)H and (19)F NMR spectroscopy, quantum mechanical calculations, and transition-state analogue binding measurements to test the distance scale on which noncovalent forces can constrain the structural relaxation or translation of side chains and ligands along a specific coordinate and the energetic consequences of such geometric constraints within the active site of bacterial ketosteroid isomerase (KSI). Our results strongly suggest that packing and binding interactions within the KSI active site can constrain local side-chain reorientation and prevent hydrogen bond shortening by 0.1 A or less. Further, this constraint has substantial energetic effects on ligand binding and stabilization of negative charge within the oxyanion hole. These results provide evidence that subtle geometric effects, indistinguishable in most X-ray crystallographic structures, can have significant energetic consequences and highlight the importance of using synergistic experimental approaches to dissect enzyme function.

Fast electron density methods in the life sciences--a routine application in the future?

The understanding of mutual recognition of biologically interacting systems on an atomic scale is of paramount importance in the life sciences. Electron density distributions that can be obtained from a high resolution X-ray diffraction experiment can provide--in addition to steric information--electronic properties of the species involved in these interactions. In recent years experimental ED methods have seen several favourable developments towards successful application in the life sciences. Experimental and methodological advances have made possible on the one hand high-speed X-ray diffraction experiments, and have allowed on the other hand the quantitative derivation of bonding, non-bonding and atomic electronic properties. This has made the investigation of a large number of molecules possible, and moreover, molecules with 200 or more atoms can be subject of experimental ED studies, as has been demonstrated by the example of vitamin B12. Supported by the experimentally verified transferability concept of submolecular electronic properties, a key issue in Bader's The Quantum Theory of Atoms in Molecules, activities have emerged to establish databases for the additive generation of electron densities of macromolecules from submolecular building blocks. It follows that the major aims of any experimental electron density work in the life sciences, namely the generation of electronic information for a series of molecules in a reasonable time and the study of biological macromolecules (proteins, polynucleotides), are within reach in the near future.

Expression screening of integral membrane proteins from Helicobacter pylori 26695

The efficiency of Helicobacter pylori as a mucosal pathogen is caused by unique soluble and integral membrane proteins, which allow its survival at acidic pH and successful colonization of the gastric environment. With about one-fourth of the H. pylori's proteome comprising integral membrane proteins, the need for solution of their three-dimensional (3D) structures becomes persistent as it can potentially drive the generation of more effective drugs. This study presents a medium-throughput approach for cloning and expression screening of integral membrane proteins from H. pylori (26695) using Escherichia coli as the expression host. One-hundred sixteen H. pylori targets were cloned into two different vector systems and heterologously expressed in E. coli. Eighty-four percent of these proteins displayed medium to high expression. No clear-cut correlation was found between expression levels and number of putative transmembrane spans, predicted functionality, and molecular mass. Nonetheless, expression of transporters and hypothetical proteins < or =40 kDa with two to four transmembrane spans displayed generally high expression levels. To statistically strengthen the quality of the data from the medium-throughput approach, a comparison with data derived from robotic-based methodologies was conducted. Optimization of expression and solubilization conditions for selected targets was also performed. Seventeen targets have been purified and subjected to crystallization so far. Eighteen percent of these targets (2/17) produced crystals under specific sets of crystallization conditions.