N. Gautham

A new peptide docking strategy using a mean field technique with mutually orthogonal Latin square sampling

The theoretical prediction of the association of a flexible ligand with a protein receptor requires efficient sampling of the conformational space of the ligand. Several docking methodologies are currently available. We propose a new docking technique that performs well at low computational cost. The method uses mutually orthogonal Latin squares to efficiently sample the docking space. A variant of the mean field technique is used to analyze this sample to arrive at the optimum. The method has been previously applied to explore the conformational space of peptides and identify structures with low values for the potential energy. Here we extend this method to simultaneously identify both the low energy conformation as well as a ‘high-scoring’ docking mode. Application of the method to 56 protein–peptide complexes, in which the length of the peptide ligand ranges from three to seven residues, and comparisons with Autodock 3.05, showed that the method works well.

Hydrophobic clusters in protein structures

Globular proteins fold such that the hydrophobic groups are packed inside forming hydrophobic clusters, and the hydrophilic groups are present on the surface. In this article we analyze clusters of hydrophobic groups of atoms in 781 protein structures selected from the PDB. Our analysis showed that every structure consists of two types of clusters: at least one large cluster that forms the hydrophobic core and probably dictates the protein fold; and numerous smaller clusters, which might be involved in the stabilization of the fold. We also analyzed the preference of the hydrophobic groups in each of the amino acids toward forming hydrophobic clusters. We find that hydrophobic groups from the hydrophilic amino acids also contribute toward cluster formation. Proteins 2008.

Protein-ligand docking using mutually orthogonal Latin squares (MOLSDOCK).

The theoretical prediction of the association of a flexible ligand with a protein receptor requires efficient sampling of the conformational space of the ligand. Several docking methodologies are currently available. We have proposed a docking technique that performs well at low computational cost. The method uses mutually orthogonal Latin squares to efficiently sample the docking space. A variant of the mean field technique is used to analyze this sample to arrive at the optimum. The method has been previously applied to search through both the conformational space of a peptide as well its docking space. Here we extend this method to simultaneously identify both the low energy conformation as well as a high scoring docking mode for the small organic ligand molecules. Application of the method to 45 protein-ligand complexes, in which the number of rotatable torsions varies from 2 to 19, and comparisons with AutoDock 4.0, showed that the method works well.

Exploring the conformational space of protein loops using a mean field technique with MOLS sampling

We have recently developed a computational technique that uses mutually orthogonal Latin square sampling to explore the conformational space of oligopeptides in an exhaustive manner. In this article, we report its use to analyze the conformational spaces of 120 protein loop sequences in proteins, culled from the PDB, having the length ranging from 5 to 10 residues. The force field used did not have any information regarding the sequences or structures that flanked the loop. The results of the analyses show that the native structure of the loop, as found in the PDB falls at one of the low energy points in the conformational landscape of the sequences. Thus, a large portion of the structural determinants of the loop may be considered intrinsic to the sequence, regardless of either adjacent sequences or structures, or the interactions that the atoms of the loop make with other residues in the protein or in neighboring proteins. Proteins 2007.

Structure of the tetradecanucleotide d(CCCCGGTACCGGGG)2 as an A‐DNA duplex

The crystal structure of the tetradecanucleotide sequence d(CCCCGGTACCGGGG)(2) has been determined at 2.5 Å resolution in the tetragonal space group P4(1). This sequence was designed with the expectation of a four-way junction. However, the sequence crystallized as an A-DNA duplex and represents more than one full turn of the A-helix. The crystallographic asymmetric unit consists of one tetradecanucleotide duplex. The structural parameters of the A-type DNA duplex structure and the crystal-packing arrangement are described. One Mn(2+) ion was identified with direct coordination to the N7 position of G(13) and a water molecule at the major-groove side of the C(2)·G(13) base pair.

Hexammineruthenium(III) ion interactions with Z-DNA

The hexamer duplex d(CGCGCA).d(TGCGCG) was crystallized with hexammineruthenium(III) ions in an orthorhombic space group; the crystals diffracted to 1.54 A resolution. Strong ion interactions with the adenine base induce a tautomeric shift from the amino to the imino form. Consequently, the A.T base pairing is disrupted. This structural study may be relevant to metal toxicity.

Molecular dynamics simulation and molecular modelling studies on the insecticidal domain from jack bean urease

Urease is an enzyme which catalyses hydrolysis of urea to ammonia and carbon dioxide. Plant ureases but not bacterial ureases display insecticidal activity against insects with cathepsin B- and D-based digestive systems. It was found that a 10-kDa internal region of a plant urease exhibits this insecticidal activity, and more precisely it was predicted that a β-hairpin motif within this region may be responsible for this activity by functioning as a membrane pore former. We carried out molecular dynamics study of the jack bean urease insecticidal region in the presence of explicit water and in water/hexane interface, in order to gain more insight into the structural changes and behaviour of the hairpin region in the membrane-like interfacial environment. The results indicate that the hairpin anchors well in the polar/non-polar interface. Subsequent modelling study of the insecticidal region clearly suggests that this region can form a β-barrel assembly and thus supports our previous hypothesis that the insecticidal activity of plant ureases occurs through pore-forming β-barrel assembly.

Expression, purification, crystallization and preliminary X-ray diffraction analysis of carbonyl reductase from Candida parapsilosis ATCC 7330.

The NAD(P)H-dependent carbonyl reductase from Candida parapsilosis ATCC 7330 catalyses the asymmetric reduction of ethyl 4-phenyl-2-oxobutanoate to ethyl (R)-4-phenyl-2-hydroxybutanoate, a precursor of angiotensin-converting enzyme inhibitors such as Cilazapril and Benazepril. The carbonyl reductase was expressed in Escherichia coli and purified by GST-affinity and size-exclusion chromatography. Crystals were obtained by the hanging-drop vapour-diffusion method and diffracted to 1.86 Å resolution. The asymmetric unit contained two molecules of carbonyl reductase, with a solvent content of 48%. The structure was solved by molecular replacement using cinnamyl alcohol dehydrogenase from Saccharomyces cerevisiae as a search model.

Structure of d(TGCGCA)2 at 293 K: comparison of the effects of sequence and temperature

The crystal structure of a hexameric DNA fragment with the sequence d(TGCGCA)(2) has been solved and refined at 293 K at a resolution of 1.64 A. The molecule adopts a left-handed Z-type helical conformation which is common in alternating pyrimidine-purine sequences. The presence of A.T base pairs at the two terminals does not perturb the structure to any great degree. However, several sequence-specific microstructural changes are noticeable. The structure of the identical sequence determined at 120 K involving somewhat different crystallization conditions has been reported previously [Harper et al. (1998), Acta Cryst. D54, 1273-1284]. A comparison of the present structure with that at low temperature and with that of d(CGCGCG)(2) shows that the effect of the change in sequence is greater than the combined effect of changes in temperature and environment.

Interactions of Mn2+ with a non-self-complementary Z-type DNA duplex.

Crystal structures of the hexanucleotide d(CACGCG)·d(CGCGTG) were determined in two crystal lattices when different concentrations of the counterion Mn2+ were used in crystallization. The availability of Mn2+ during the crystallization process appears to play an important role in inducing different crystal packings that lead to crystals belonging to the two space groups P2(1) and P6(5). Analysis of the molecular interactions of Mn2+ with the Z-form duplexes shows direct coordination to the purine residues G and A.