BOOMSLANG: a program for combinatorial structure generation

Probing the residual structure in avian prion hexarepeats by CD, NMR and MD techniques

Many proteins perform essential biological functions by means of regions that lacking specific organized structure exist as an ensemble of interconverting transient conformers. The characterization of such regions, including the description of their structural propensities, number of conformations and relative populations can provide useful insights. Prion diseases result from the conversion of a normal glycoprotein into a misfolded pathogenic isoform. The structures of mammal and chicken prion proteins show a similar fold with a globular domain and a flexible N-terminal portion that contains different repeated regions: octarepeats (PHGGGWGQ) in mammals and hexarepeats (PHNPGY) in chickens. The higher number of prolines in the hexarepeat region suggests that this region may retain a significant amount of residual secondary structure. Here, we report the CD, NMR and MD characterization of a peptide (2-HexaPY) composed of two hexarepeats. We combine experimental NMR data and MD to investigate at atomic level its ensemble-averaged structural properties, demonstrating how each residue of both repeats has a different quantified PPII propensity that shows a periodicity along the sequence. This feature explains the absence of cooperativity to stabilize a PPII conformation. Nonetheless, such residual structure can play a role in nucleating local structural transitions as well as modulating intra-molecular or inter-molecular interactions.

Modeling the interactions of phthalocyanines in water: from the Cu(II)-tetrasulphonate to the metal-free phthalocyanine

A quantum and statistical study on the effects of the ions Cu(2+) and SO(3)(-) in the solvent structure around the metal-free phthalocyanine (H(2)Pc) is presented. We developed an ab initio interaction potential for the system CuPc-H(2)O based on quantum chemical calculations and studied its transferability to the H(2)Pc-H(2)O and [CuPc(SO(3))(4)](4-)-H(2)O interactions. The use of the molecular dynamics technique allows the determination of energetic and structural properties of CuPc, H(2)Pc, and [CuPc(SO(3))(4)](4-) in water and the understanding of the keys for the different behaviors of the three phthalocyanine (Pc) derivatives in water. The inclusion of the Cu(2+) cation in the Pc structure reinforces the appearance of two axial water molecules and second-shell water molecules in the solvent structure, whereas the presence of SO(3)(-) anions implies a well defined hydration shell of about eight water molecules around them making the macrocycle soluble in water. Debye-Waller factors for axial water molecules have been obtained in order to examine the potential sensitivity of the extended x-ray absorption fine structure technique to detect the axial water molecules.

Generative models for chemical structures

We apply recently developed techniques for pattern recognition to construct a generative model for chemical structure. This approach can be viewed as ligand-based de novo design. We construct a statistical model describing the structural variations present in a set of molecules which may be sampled to generate new structurally similar examples. We prevent the possibility of generating chemically invalid molecules, according to our implicit hydrogen model, by projecting samples onto the nearest chemically valid molecule. By populating the input set with molecules that are active against a target, we show how new molecules may be generated that will likely also be active against the target.

Characterization of the tunneling conductance across DNA bases

Characterization of the electrical properties of the DNA bases (adenine, cytosine, guanine, and thymine), in addition to building the basic knowledge on these fundamental constituents of a DNA, is a crucial step in developing a DNA sequencing technology. We present a first-principles study of the current-voltage characteristics of nucleotidelike molecules of the DNA bases, placed in a 1.5 nm gap formed between gold nanoelectrodes. The quantum transport calculations in the tunneling regime are shown to vary strongly with the electrode-molecule geometry and the choice of the density-functional theory exchange-correlation functionals. Analysis of the results in the zero-bias limit indicates that distinguishable current-voltage characteristics of different DNA bases are dominated by the geometrical conformations of the bases and nanoelectrodes.

Modern computational chemistry and drug discovery: structure generating programs

During 1996 and 1997, the first reports were disclosed of active enzyme inhibitors based entirely on novel structures created by de novo methods. De novo methods have also been used to modify and significantly improve the binding affinity of an HIV protease inhibitor. Work continues in the improvement of methods for the de novo design of compounds which fit and chemically complement a binding site. De novo algorithms that generate only synthetically feasible structures have also been reported. In addition, methods are being developed for the automatic computer generation of virtual molecular libraries which can be searched to identify molecules to match a pharmacophore or fit into a binding site.