Molecular Replacement Techniques Research Articles

The Automatic Solution of Macromolecular Crystal Structures via Molecular Replacement Techniques: REMO22 and Its Pipeline.

A description of REMO22, a new molecular replacement program for proteins and nucleic acids, is provided. This program, as with REMO09, can use various types of prior information through appropriate conditional distribution functions. Its efficacy in model searching has been validated through several test cases involving proteins and nucleic acids. Although REMO22 can be configured with different protocols according to user directives, it has been developed primarily as an automated tool for determining the crystal structures of macromolecules. To evaluate REMO22's utility in the current crystallographic environment, its experimental results must be compared favorably with those of the most widely used Molecular Replacement (MR) programs. To accomplish this, we chose two leading tools in the field, PHASER and MOLREP. REMO22, along with MOLREP and PHASER, were included in pipelines that contain two additional steps: phase refinement (SYNERGY) and automated model building (CAB). To evaluate the effectiveness of REMO22, SYNERGY and CAB, we conducted experimental tests on numerous macromolecular structures. The results indicate that REMO22, along with its pipeline REMO22 + SYNERGY + CAB, presents a viable alternative to currently used phasing tools.

Cyclic Automated Model Building (CAB) Applied to Nucleic Acids

Obtaining high-quality models for nucleic acid structures by automated model building programs (AMB) is still a challenge. The main reasons are the rather low resolution of the diffraction data and the large number of rotatable bonds in the main chains. The application of the most popular and documented AMB programs (e.g., PHENIX.AUTOBUILD, NAUTILUS and ARP/wARP) may provide a good assessment of the state of the art. Quite recently, a cyclic automated model building (CAB) package was described; it is a new AMB approach that makes the use of BUCCANEER for protein model building cyclic without modifying its basic algorithms. The applications showed that CAB improves the efficiency of BUCCANEER. The success suggested an extension of CAB to nucleic acids—in particular, to check if cyclically including NAUTILUS in CAB may improve its effectiveness. To accomplish this task, CAB algorithms designed for protein model building were modified to adapt them to the nucleic acid crystallochemistry. CAB was tested using 29 nucleic acids (DNA and RNA fragments). The phase estimates obtained via molecular replacement (MR) techniques were automatically submitted to phase refinement and then used as input for CAB. The experimental results from CAB were compared with those obtained by NAUTILUS, ARP/wARP and PHENIX.AUTOBUILD.

How far are we from automatic crystal structure solution via molecular-replacement techniques?

Although the success of molecular-replacement techniques requires the solution of a six-dimensional problem, this is often subdivided into two three-dimensional problems. REMO09 is one of the programs which have adopted this approach. It has been revisited in the light of a new probabilistic approach which is able to directly derive conditional distribution functions without passing through a previous calculation of the joint probability distributions. The conditional distributions take into account various types of prior information: in the rotation step the prior information may concern a non-oriented model molecule alone or together with one or more located model molecules. The formulae thus obtained are used to derive figures of merit for recognizing the correct orientation in the rotation step and the correct location in the translation step. The phases obtained by this new version of REMO09 are used as a starting point for a pipeline which in its first step extends and refines the molecular-replacement phases, and in its second step creates the final electron-density map which is automatically interpreted by CAB, an automatic model-building program for proteins and DNA/RNA structures.

Crystal structure determination and analysis of 11S coconut allergen: Cocosin

Allergy is an abnormal immune response against an innocuous target. Food allergy is an adverse reaction caused by common foods most well-known being those involving peanuts. Apart from mono sensitized food allergy, cross-reactivity with other food allergens is also commonly observed. To understand the phenomenon of cross-reactivity related to immune response, three dimensional structures of the allergens and their antigenic epitopes has to be analysed in detail. The X-ray crystal structure of Cocosin, a common 11S food allergen from coconut, has been determined at 2.2Å resolution using molecular replacement technique. The monomer of 52kDa is composed of two β-jelly roll domains, one with acidic and the other with basic character. The structure shows hexameric association with two trimers facing each other. Though the overall structure of Cocosin is similar to other 11S allergens, the occurrence of experimentally determined epitopes of the peanut allergen Ara h 3 at flexible as well as variable regions could be the reason for the clinically reported result of cross-reactivity that the peanut allergic patients are not sensitized with coconut allergen.

QSAR modeling and in silico design of small-molecule inhibitors targeting the interaction between E3 ligase VHL and HIF-1α.

Protein-protein interactions (PPIs) have attracted much attention recently because of their preponderant role in most biological processes. The prevention of the interaction between E3 ligase VHL and HIF-1[Formula: see text] may improve tolerance to hypoxia and ameliorate the prognosis of many diseases. To obtain novel potent inhibitors of VHL/HIF-1[Formula: see text] interaction, a series of hydroxyproline-based inhibitors were investigated for structural optimization using a combination of QSAR modeling and molecular docking. Here, 2D- and 3D-QSAR models were developed by genetic function approximation (GFA) and comparative molecular field analysis (CoMFA) and comparative molecular similarity index analysis (CoMSIA) methods, respectively. The top-ranked models with strict validation revealed satisfactory statistical parameters (CoMFA with [Formula: see text], 0.637; [Formula: see text], 0.955; [Formula: see text], 0.944; CoMSIA with [Formula: see text], 0.649; [Formula: see text], 0.954; [Formula: see text], 0.911; GFA with [Formula: see text], 0.721; [Formula: see text], 0.801; [Formula: see text], 0.861). The selected five 2D-QSAR descriptors were in good accordance with the 3D-QSAR results, and contour maps gave the visualization of feature requirements for inhibitory activity. A new diverse molecular database was created by molecular fragment replacement and BREED techniques for subsequent virtual screening. Eventually, 31 novel hydroxyproline derivatives stood out as potential VHL/HIF-1[Formula: see text] inhibitors with favorable predictions by the CoMFA, CoMSIA and GFA models. The reliability of this protocol suggests that it could also be applied to the exploration of lead optimization of other PPI targets.

Recombinant expression, purification and preliminary biophysical and structural studies of C-terminal carbohydrate recognition domain from human galectin-4

Galectin-4 (Gal4), a tandem-repeat type galectin, is expressed in healthy epithelium of the gastrointestinal tract. Altered levels of Gal4 expression are associated with different types of cancer, suggesting its usage as a diagnostic marker as well as target for drug development. The functional data available for this class of proteins suggest that the wide spectrum of cellular activities reported for Gal4 relies on distinct glycan specificity and structural characteristics of its two carbohydrate recognition domains. In the present work, two independent constructs for recombinant expression of the C-terminal domain of human galectin-4 (hGal4-CRD2) were developed. His6-tagged and untagged recombinant proteins were overexpressed in Escherichia coli, and purified by affinity chromatography followed by gel filtration. Correct folding and activity of hGal4-CRD2 were assessed by circular dichroism and fluorescence spectroscopies, respectively. Diffraction quality crystals were obtained by vapor-diffusion sitting drop setup and the crystal structure of CRD2 was solved by molecular replacement techniques at 1.78 Å resolution. Our work describes the development of important experimental tools that will allow further studies in order to correlate structure and binding properties of hGal4-CRD2 and human galectin-4 functional activities.

Author Index

Here, we report the crystallographic study of a lectin from Canavalia maritima seeds (ConM) and its relaxant activity on vascular smooth muscle, to provide new insights into the understanding of structure/function relationships of this class of proteins. ConM was crystallized and its structure determined by standard molecular replacement techniques. The amino acid residues, previously suggested incorrectly by manual sequencing, have now been determined as I17, I53, S129, S134, G144, S164, P165, S187, V190, S169, T196, and S202. Analysis of the structure indicated a dimer in the asymmetric unit, two metal binding sites per monomer, and loops involved in the molecular oligomerization. These confer 98% similarity between ConM and other previously described lectins, derived from Canavalia ensiformis and Canavalia brasiliensis. Our functional data indicate that ConM exerts a concentration-dependent relaxant action on isolated aortic rings that probably occurs via an interaction with a specific lectin-binding site on the endothelium, resulting in a release of nitric oxide.

Refining a model electron-density map via the Phantom Derivative method.

The Phantom Derivative (PhD) method [Giacovazzo (2015), Acta Cryst. A71, 483-512] has recently been described for ab initio and non-ab initio phasing. It is based on the random generation of structures with the same unit cell and the same space group as the target structure (called ancil structures), which are used to create derivatives devoid of experimental diffraction amplitudes. In this paper, the non-ab initio variant of the method was checked using phase sets obtained by molecular-replacement techniques as a starting point for phase extension and refinement. It has been shown that application of PhD is able to extend and refine phases in a way that is competitive with other electron-density modification techniques.

The X‐Ray Crystallographic Structure of Endoxylanase Xyn10C from Xanthomonas axonopodis at 1.6 Angstrom Resolution

The recent biochemical characterization of the xylanases of glycoside hydrolase family 10 (GH 10) has identified a distinctive endolytic mode of action, hydrolyzing the β-1,4 xylan chain at a specific site directed by the position of an α-1,2-linked glucuronate moiety. Xylanase C (XynC), the GH10 xylanase from Xanthomonas axonopodis, has been cloned, overexpressed and crystallized. Previous attempts to crystallize the protein using commercially available screens were unsuccessful, but analysis of the crystallization condition of a homologous protein allowed us to construct a novel grid screen and obtain small (< 0.5 mm) crystals of the protein in 24 hours via hanging drop vapor diffusion. Crystals grew in 100 mM sodium phosphate citrate, 200 mM NaCl, 25% (w/v) PEG 8000, 5% (v/v) glycerol, pH 5.5. The crystals belonged to the orthorhombic space group C2221 and had unit cell parameters of 71 Å, 80 Å and 161 Å and α=β=γ=90°. Diffraction data were collected to 1.6 Å resolution and an initial phasing model was created using molecular replacement techniques. Multiple rounds of model building and refinement were performed using Coot and Phenix, respectively. The protein was found to possess the typical (β/α)8 fold observed in all members of GH10. A large loop on the aglycone side of the active site was found to adopt a unique conformation. This loop was in a position that would prevent it from binding substrates larger than three xylose units in size and would be expected to produce a xylobiose and a free xylose upon cleavage of such a trimer. This work will help to identify the structural determinants that grant the exceptional specificity of this enzyme and the role it plays in the biological depolymerization and processing of glucuronoxylan.

Routine phasing of coiled-coil protein crystal structures with AMPLE.

Coiled-coil protein folds are among the most abundant in nature. These folds consist of long wound α-helices and are architecturally simple, but paradoxically their crystallographic structures are notoriously difficult to solve with molecular-replacement techniques. The program AMPLE can solve crystal structures by molecular replacement using ab initio search models in the absence of an existent homologous protein structure. AMPLE has been benchmarked on a large and diverse test set of coiled-coil crystal structures and has been found to solve 80% of all cases. Successes included structures with chain lengths of up to 253 residues and resolutions down to 2.9 Å, considerably extending the limits on size and resolution that are typically tractable by ab initio methodologies. The structures of two macromolecular complexes, one including DNA, were also successfully solved using their coiled-coil components. It is demonstrated that both the ab initio modelling and the use of ensemble search models contribute to the success of AMPLE by comparison with phasing attempts using single structures or ideal polyalanine helices. These successes suggest that molecular replacement with AMPLE should be the method of choice for the crystallo-graphic elucidation of a coiled-coil structure. Furthermore, AMPLE may be able to exploit the presence of a coiled coil in a complex to provide a convenient route for phasing.

The structure of rice weevil pectin methylesterase.

Rice weevils (Sitophilus oryzae) use a pectin methylesterase (EC 3.1.1.11), along with other enzymes, to digest cell walls in cereal grains. The enzyme is a right-handed β-helix protein, but is circularly permuted relative to plant and bacterial pectin methylesterases, as shown by the crystal structure determination reported here. This is the first structure of an animal pectin methylesterase. Diffraction data were collected to 1.8 Å resolution some time ago for this crystal form, but structure solution required the use of molecular-replacement techniques that have been developed and similar structures that have been deposited in the last 15 years. Comparison of the structure of the rice weevil pectin methylesterase with that from Dickeya dandantii (formerly Erwinia chrysanthemi) indicates that the reaction mechanisms are the same for the insect, plant and bacterial pectin methylesterases. The similarity of the structure of the rice weevil enzyme to the Escherichia coli lipoprotein YbhC suggests that the evolutionary origin of the rice weevil enzyme was a bacterial lipoprotein, the gene for which was transferred to a primitive ancestor of modern weevils and other Curculionidae. Structural comparison of the rice weevil pectin methylesterase with plant and bacterial enzymes demonstrates that the rice weevil protein is circularly permuted relative to the plant and bacterial molecules.

Crystal structure of Rv2258c, a putative SAM-dependent methyltransferase

The Mycobacterium tuberculosis protein Rv2258c belongs to a large family of putative S-adenosylmethionine (SAM)-dependent methyltransferases in mycobacteria. As part of a laboratory-scale structural genomics project on conserved hypothetical proteins in pathogenic bacteria, we have determined the crystal structure of Rv2258c from M. tuberculosis H37Rv at 1.8 Å resolution. The crystals of apo Rv2258c belong to space group C2, with unit cell parameters a = 109.2 Å, b = 140.6 Å, c = 97.2 Å, and β = 98.50. Assuming that three monomers are present in the asymmetric unit, the Matthews parameter and the solvent fraction are 3.32 Å3 Da-1 and 63.0%, respectively. The crystal structure of Rv2258c has been solved by a combination of the molecular replacement and single wavelength anomalous diffraction technique, using platinum atoms as anomalous scattering centers and the coordinates of the RebM from Lechevalieria aerocolonigenes (PDB entry 3BUS) as the search model. In addition, two binary complex structures, Rv2258c-SAM and Rv2258c-S-adenosylhomocysteine (SAH), have been determined by molecular replacement using the model of apo Rv2258c.

Crystallization and preliminary X-ray crystallographic analysis of importin-α from Neurospora crassa.

Importin-α recognizes cargo proteins that contain classical nuclear localization sequences (NLS) and, in complex with importin-β, is able to translocate nuclear proteins through the nuclear pore complex. The filamentous fungus Neurospora crassa is a well studied organism that has been widely used as a model organism for fundamental aspects of eukaryotic biology, and is important for understanding the specific mechanisms of protein transport to the cell nucleus. In this work, the crystallization and preliminary X-ray diffraction analysis of importin-α from N. crassa (IMPα-Nc) complexed with a classical NLS peptide (SV40 NLS) are reported. IMPα-Nc-SV40 NLS crystals diffracted X-rays to 2.0 Å resolution and the structure was solved by molecular-replacement techniques, leading to a monomeric structure. The observation of the electron-density map indicated the presence of SV40 NLSs interacting at both the minor and major NLS-binding sites of the protein.

Applications of molecular replacement to G protein-coupled receptors.

G protein-coupled receptors (GPCRs) are a large class of integral membrane proteins involved in regulating virtually every aspect of human physiology. Despite their profound importance in human health and disease, structural information regarding GPCRs has been extremely limited until recently. With the advent of a variety of new biochemical and crystallographic techniques, the structural biology of GPCRs has advanced rapidly, offering key molecular insights into GPCR activation and signal transduction. To date, almost all GPCR structures have been solved using molecular-replacement techniques. Here, the unique aspects of molecular replacement as applied to individual GPCRs and to signaling complexes of these important proteins are discussed.

Rapid sampling of local minima in protein energy surface and effective reduction through a multi-objective filter.

BackgroundMany problems in protein modeling require obtaining a discrete representation of the protein conformational space as an ensemble of conformations. In ab-initio structure prediction, in particular, where the goal is to predict the native structure of a protein chain given its amino-acid sequence, the ensemble needs to satisfy energetic constraints. Given the thermodynamic hypothesis, an effective ensemble contains low-energy conformations which are similar to the native structure. The high-dimensionality of the conformational space and the ruggedness of the underlying energy surface currently make it very difficult to obtain such an ensemble. Recent studies have proposed that Basin Hopping is a promising probabilistic search framework to obtain a discrete representation of the protein energy surface in terms of local minima. Basin Hopping performs a series of structural perturbations followed by energy minimizations with the goal of hopping between nearby energy minima. This approach has been shown to be effective in obtaining conformations near the native structure for small systems. Recent work by us has extended this framework to larger systems through employment of the molecular fragment replacement technique, resulting in rapid sampling of large ensembles.MethodsThis paper investigates the algorithmic components in Basin Hopping to both understand and control their effect on the sampling of near-native minima. Realizing that such an ensemble is reduced before further refinement in full ab-initio protocols, we take an additional step and analyze the quality of the ensemble retained by ensemble reduction techniques. We propose a novel multi-objective technique based on the Pareto front to filter the ensemble of sampled local minima.Results and conclusionsWe show that controlling the magnitude of the perturbation allows directly controlling the distance between consecutively-sampled local minima and, in turn, steering the exploration towards conformations near the native structure. For the minimization step, we show that the addition of Metropolis Monte Carlo-based minimization is no more effective than a simple greedy search. Finally, we show that the size of the ensemble of sampled local minima can be effectively and efficiently reduced by a multi-objective filter to obtain a simpler representation of the probed energy surface.

Elucidating the ensemble of functionally-relevant transitions in protein systems with a robotics-inspired method

BackgroundMany proteins tune their biological function by transitioning between different functional states, effectively acting as dynamic molecular machines. Detailed structural characterization of transition trajectories is central to understanding the relationship between protein dynamics and function. Computational approaches that build on the Molecular Dynamics framework are in principle able to model transition trajectories at great detail but also at considerable computational cost. Methods that delay consideration of dynamics and focus instead on elucidating energetically-credible conformational paths connecting two functionally-relevant structures provide a complementary approach. Effective sampling-based path planning methods originating in robotics have been recently proposed to produce conformational paths. These methods largely model short peptides or address large proteins by simplifying conformational space.MethodsWe propose a robotics-inspired method that connects two given structures of a protein by sampling conformational paths. The method focuses on small- to medium-size proteins, efficiently modeling structural deformations through the use of the molecular fragment replacement technique. In particular, the method grows a tree in conformational space rooted at the start structure, steering the tree to a goal region defined around the goal structure. We investigate various bias schemes over a progress coordinate for balance between coverage of conformational space and progress towards the goal. A geometric projection layer promotes path diversity. A reactive temperature scheme allows sampling of rare paths that cross energy barriers.Results and conclusionsExperiments are conducted on small- to medium-size proteins of length up to 214 amino acids and with multiple known functionally-relevant states, some of which are more than 13Å apart of each-other. Analysis reveals that the method effectively obtains conformational paths connecting structural states that are significantly different. A detailed analysis on the depth and breadth of the tree suggests that a soft global bias over the progress coordinate enhances sampling and results in higher path diversity. The explicit geometric projection layer that biases the exploration away from over-sampled regions further increases coverage, often improving proximity to the goal by forcing the exploration to find new paths. The reactive temperature scheme is shown effective in increasing path diversity, particularly in difficult structural transitions with known high-energy barriers.

Crystallization and preliminary X-ray crystallographic analysis of the oxysterol-binding protein Osh3 from Saccharomyces cerevisiae.

Oxysterol-binding protein (OSBP) related proteins (ORPs) are conserved from yeast to humans and are implicated in regulation of sterol homeostasis and in signal transduction pathways. Osh3 of Saccharomyces cerevisiae is a pleckstrin-homology (PH) domain-containing ORP member that regulates phosphoinositide metabolism at endoplasmic reticulum-plasma membrane contact sites. The N-terminal PH domain of Osh3 was purified and crystallized as a lysozyme fusion and the resulting crystal diffracted to 2.3 Å resolution. The crystal belonged to the monoclinic space group C2, with unit-cell parameters a=98.03, b=91.31, c=84.13 Å, β=81.41°. With two molecules in the asymmetric unit, the Matthews coefficient was 3.13 Å3 Da(-1). Initial attempts to solve the structure by molecular-replacement techniques using T4 lysozyme as a search model were successful. The C-terminal OSBP-related domain (OBD) of Osh3 was crystallized by the vapour-diffusion method and the resulting crystal diffracted to 1.5 Å resolution. The crystal was orthorhombic, belonging to space group P2(1)2(1)2(1), with unit-cell parameters a=41.57, b=87.52, c=100.58 Å. With one molecule in the asymmetric unit, the Matthews coefficient was 2.01 Å3 Da(-1). Initial attempts to solve the structure by the single-wavelength anomalous dispersion technique using bromine were successful.

Purification, crystallization and preliminary X-ray crystallographic analysis of the ETS domain of human Ergp55 in complex with the cfos promoter DNA sequence.

The Ergp55 protein belongs to the Ets family of transciption factors. The Ets transcription factors are involved in various developmental processes and the regulation of cancer metabolism. They contain a highly similar DNA-binding domain known as the ETS domain and have diverse functions in oncogenesis and physiology. The Ets transcription factors differ in their DNA-binding preference at the ETS site and the mechanisms by which they target genes are not clearly understood. To understand its DNA-binding mechanism, the ETS domain of Ergp55 was expressed and purified. The ETS domain was crystallized in the native form and in complex forms with DNA sequences from the E74 and cfos promoters. An X-ray diffraction data set was collected from an ETS-cfos DNA complex crystal at a wavelength of 0.9725 Å on the BM14 synchrotron beamline at the ESRF, France. The ETS-cfos DNA complex crystal belonged to space group C222(1), with four molecules in the asymmetric unit. For structure analysis, initial phases for the ETS-cfos DNA complex were obtained by the molecular-replacement technique with Phaser in the CCP4 suite using the coordinates of Fli-1 protein (PDB entry 1fli) and cfos DNA (PDB entry 1bc7) as search models. Structure analysis of the ETS-cfos DNA complex may possibly explain the DNA-binding specificity and its mechanism of interaction with the ETS domain of Ergp55.

Crystal Structure of Cryptosporidium parvum Pyruvate Kinase

Pyruvate kinase plays a critical role in cellular metabolism of glucose by serving as a major regulator of glycolysis. This tetrameric enzyme is allosterically regulated by different effector molecules, mainly phosphosugars. In response to binding of effector molecules and substrates, significant structural changes have been identified in various pyruvate kinase structures. Pyruvate kinase of Cryptosporidium parvum is exceptional among known enzymes of protozoan origin in that it exhibits no allosteric property in the presence of commonly known effector molecules. The crystal structure of pyruvate kinase from C. parvum has been solved by molecular replacement techniques and refined to 2.5 Å resolution. In the active site a glycerol molecule is located near the γ-phosphate site of ATP, and the protein structure displays a partially closed active site. However, unlike other structures where the active site is closed, the α6' helix in C. parvum pyruvate kinase unwinds and assumes an extended conformation. In the crystal structure a sulfate ion is found at a site that is occupied by a phosphate of the effector molecule in many pyruvate kinase structures. A new feature of the C. parvum pyruvate kinase structure is the presence of a disulfide bond cross-linking the two monomers in the asymmetric unit. The disulfide bond is formed between cysteine residue 26 in the short N-helix of one monomer with cysteine residue 312 in a long helix (residues 303–320) of the second monomer at the interface of these monomers. Both cysteine residues are unique to C. parvum, and the disulfide bond remained intact in a reduced environment. However, the significance of this bond, if any, remains unknown at this time.

Crystallization and preliminary X-ray diffraction analysis of three myotoxic phospholipases A2 from Bothrops brazili venom.

Two myotoxic and noncatalytic Lys49-phospholipases A(2) (braziliantoxin-II and MT-II) and a myotoxic and catalytic phospholipase A(2) (braziliantoxin-III) from the venom of the Amazonian snake Bothrops brazili were crystallized. The crystals diffracted to resolutions in the range 2.56-2.05 Å and belonged to space groups P3(1)21 (braziliantoxin-II), P6(5)22 (braziliantoxin-III) and P2(1) (MT-II). The structures were solved by molecular-replacement techniques. Both of the Lys49-phospholipases A(2) (braziliantoxin-II and MT-II) contained a dimer in the asymmetric unit, while the Asp49-phospholipase A(2) braziliantoxin-III contained a monomer in its asymmetric unit. Analysis of the quaternary assemblies of the braziliantoxin-II and MT-II structures using the PISA program indicated that both models have a dimeric conformation in solution. The same analysis of the braziliantoxin-III structure indicated that this protein does not dimerize in solution and probably acts as a monomer in vivo, similar to other snake-venom Asp49-phospholipases A(2).