Implementation and application of basis set superposition error-correction schemes to the theoretical modeling of weak intermolecular interactions

This thesis deals with the so-called Basis Set Superposition Error (BSSE) from both a methodological and a practical point of view. The purpose of the present thesis is twofold: (a) to contribute step ahead in the correct characterization of weakly bound complexes and, (b) to shed light the understanding of the actual implications of the basis set extension effects in the ab intio calculations and contribute to the BSSE debate. The existing BSSE-correction procedures are deeply analyzed, compared, validated and, if necessary, improved. A new interpretation of the counterpoise (CP) method is used in order to define counterpoise-corrected descriptions of the molecular complexes. This novel point of view allows for a study of the BSSE-effects not only in the interaction energy but also on the potential energy surface and, in general, in any property derived from the molecular energy and its derivatives A program has been developed for the calculation of CP-corrected geometry optimizations and vibrational frequencies, also using several counterpoise schemes for the case of molecular clusters. The method has also been implemented in Gaussian98 revA10 package. The Chemical Hamiltonian Approach (CHA) methodology has been also implemented at the RHF and UHF levels of theory for an arbitrary number interacting systems using an algorithm based on block-diagonal matrices. Along with the methodological development, the effects of the BSSE on the properties of molecular complexes have been discussed in detail. The CP and CHA methodologies are used for the determination of BSSE-corrected molecular complexes properties related to the Potential Energy Surfaces and molecular wavefunction, respectively. First, the behaviour of both BSSE-correction schemes are systematically compared at different levels of theory and basis sets for a number of hydrogen-bonded complexes. The Complete Basis Set (CBS) limit of both uncorrected and CP-corrected molecular properties like stabilization energies and intermolecular distances has also been determined, showing the capital importance of the BSSE correction. Several controversial topics of the BSSE correction are addressed as well. The application of the counterpoise method is applied to internal rotational barriers. The importance of the nuclear relaxation term is also pointed out. The viability of the CP method for dealing with charged complexes and the BSSE effects on the double-well PES blue-shifted hydrogen bonds is also studied in detail. In the case of the molecular clusters the effect of high-order BSSE effects introduced with the hierarchical counterpoise scheme is also determined. The effect of the BSSE on the electron density-related properties is also addressed. The first-order electron density obtained with the CHA/F and CHA/DFT methodologies was used to assess, both graphically and numerically, the redistribution of the charge density upon BSSE-correction. Several tools like the Atoms in Molecules topologycal analysis, density difference maps, Quantum Molecular Similarity, and Chemical Energy Component Analysis were used to deeply analyze, for the first time, the BSSE effects on the electron density of several hydrogen bonded complexes of increasing size. The indirect effect of the BSSE on intermolecular perturbation theory results is also pointed out It is shown that for a BSSE-free SAPT study of hydrogen fluoride clusters, the use of a counterpoise-corrected PES is essential in order to determine the proper molecular geometry to perform the SAPT analysis.

Analytical potentials for triatomic molecules from spectroscopic data V. Application to HOX (X=F, Cl, Br, I)

Analytical potential energy functions are reported for HOX (X=F, Cl, Br, I). The surface for HOF predicts two metastable minima as well as the equilibrium configuration. These correspond to HFO (bent) and OHF (linear). Ab initio calculations performed for the HOF surface confirm these predictions. Comparisons are drawn between the two sets of results, and a vibrational analysis is undertaken for the hydrogen bonded OHF species. For HOCl, one further minimum is predicted, corresponding to HClO (bent), the parameters for which compare favourably with those reported from ab initio studies. In contrast, only the equilibrium configurations are predicted to be stable for HOBr and HOI.

Amino acid pairing preferences in parallel beta-sheets in proteins

Statistical approaches have been applied to examine amino acid pairing preferences within parallel beta-sheets. The main chain hydrogen bonding pattern in parallel beta-sheets means that, for each residue pair, only one of the residues is involved in main chain hydrogen bonding with the strand containing the partner residue. We call this the hydrogen bonded (HB) residue and the partner residue the non-hydrogen bonded (nHB) residue, and differentiate between the favorability of a pair and that of its reverse pair, e.g. Asn(HB)-Thr(nHB)versus Thr(HB)-Asn(nHB). Significantly (p < or = 0.000001) favoured pairings were rationalised using stereochemical arguments. For instance, Asn(HB)-Thr(nHB) and Arg(HB)-Thr(nHB) were favoured pairs, where the residues adopted favoured chi1 rotamer positions that allowed side-chain interactions to occur. In contrast, Thr(HB)-Asn(nHB) and Thr(HB)-Arg(nHB) were not significantly favoured, and could only form side-chain interactions if the residues involved adopted less favourable chi1 conformations. The favourability of hydrophobic pairs e.g. Ile(HB)-Ile(nHB), Val(HB)-Val(nHB) and Leu(HB)-Ile(nHB) was explained by the residues adopting their most preferred chi1 and chi2 conformations, which enabled them to form nested arrangements. Cysteine-cysteine pairs are significantly favoured, although these do not form intrasheet disulphide bridges. Interactions between positively and negatively charged residues were asymmetrically preferred: those with the negatively charged residue at the HB position were more favoured. This trend was accounted for by the presence of general electrostatic interactions, which, based on analysis of distances between charged atoms, were likely to be stronger when the negatively charged residue is the HB partner. The Arg(HB)-Asp(nHB) interaction was an exception to this trend and its favorability was rationalised by the formation of specific side-chain interactions. This research provides rules that could be applied to protein structure prediction, comparative modelling and protein engineering and design. The methods used to analyse the pairing preferences are automated and detailed results are available (http://www.rubic.rdg.ac.uk/betapairprefsparallel/).

Amino acid pairing preferences in parallel beta-sheets in proteins

Statistical approaches have been applied to examine amino acid pairing preferences within parallel beta-sheets. The main chain hydrogen bonding pattern in parallel beta-sheets means that, for each residue pair, only one of the residues is involved in main chain hydrogen bonding with the strand containing the partner residue. We call this the hydrogen bonded (HB) residue and the partner residue the non-hydrogen bonded (nHB) residue, and differentiate between the favourability of a pair and that of its reverse pair, e.g. Asn(HB)-Thr(nHB) versus Thr(HB)-Asn(nHB). Significantly (p <= 0.000001) favoured pairings were rationalised using stereochemical arguments. For instance, Asn(HB)-Thr(nHB) and Arg(HB)-Thr(nHB) were favoured pairs, where the residues adopted favoured chi(1) rotamer positions that allowed side-chain interactions to occur. In contrast, Thr(HB)-Asn(nHB) and Thr(HB)-Arg(nHB) were not significantly favoured, and could only form side-chain interactions if the residues involved adopted less favourable chi(1) conformations. The favourability of hydrophobic pairs e.g. Ile(HB)-Ile(nHB), Val(HB)-Val(nHB) and Leu(HB)-Ile(nHB) was explained by the residues adopting their most preferred chi(1) and chi(2) conformations, which enabled them to form nested arrangements. Cysteine-cysteine pairs are significantly favoured, although these do not form intrasheet disulphide bridges. Interactions between positively and negatively charged residues were asymmetrically preferred: those with the negatively charged residue at the HB position were more favoured. This trend was accounted for by the presence of general electrostatic interactions, which, based on analysis of distances between charged atoms, were likely to be stronger when the negatively charged residue is the HB partner. The Arg(HB)-Asp(nHB) interaction was an exception to this trend and its favourability was rationalised by the formation of specific side-chain interactions. This research provides rules that could be applied to protein structure prediction, comparative modelling and protein engineering and design. The methods used to analyse the pairing preferences are automated and detailed results are available (http:// www.rubic.rdg.ac.uk/betapairprefsparallel/). (c) 2005 Elsevier Ltd. All rights reserved.

Stabilization of a helical water chain in a metal-organic host of a trinuclear Schiff base complex

Three heterometallic trinuclear Schiff base complexes, [{GuL(1)(H2O)}(2)Ni(CN)(4)]center dot 4H(2)O (1), [{CuL2(H2O)}(2)Ni(CN)(4)] (2), and [{CuL3(H2O)}(2)Ni(CN)(4)] (3) (HL1 = 7-amino-4-methyl-5-azahept-3-en-2-one, HL2 = 7-methylamino-4-methyl-5-azahept-3-en-2-one, and HL3 = 7-dimethylamino-4-methyl-5-azahept-3-en-2-one), were synthesized. All three complexes were characterized by elemental analysis, IR and UV spectroscopies, and thermal analysis. Two of them (1 and 3) were also characterized by single crystal X-ray crystallography. Complex 1 forms a hydrogen-bonded one-dimensional metal-organic framework that stabilizes a helical water chain into its cavity, but when any of the amine hydrogen atoms of the Schiff base are replaced by methyl groups, as in L 2 and L 3, the water chain, vanishes, showing explicitly the importance of the host-guest H-bonding interactions for the stabilization of a water cluster.

Low-temperature determination of theophylline dimethyl sulfoxide solvate

The title solvate, C7H8N4O2 center dot C2H6OS, was obtained unintentionally from a cocrystal screen involving theophylline and isophthalic acid. One molecule each of theophylline and dimethyl sulfoxide is present in the asymmetric unit. The packing consists of molecular sheets lying parallel to the ( 040) series of lattice planes, in which each theophylline molecule is hydrogen bonded to one dimethyl sulfoxide molecule through an N-H center dot center dot center dot O [2.7658 (15) angstrom] hydrogen bond. This particular hydrogen-bond donor was found to be used in this type of interaction in a variety of other crystal structures of theophylline.

Stepwise self-assembly of a tripeptide from molecular dimers to supramolecular beta-sheets in crystals and amyloid-like fibrils in the solid state

The terminally protected tripeptide Boc-Ala(1)-Leu(2)-Ala(3)-OMe 1 forms antiparallel hydrogen-bonded dimers of two different conformers in the asymmetric unit and the individual dimers then self-associate to form supramolecular beta-sheet structures in crystals and amyloid-like fibrils in the solid state.

Supramolecular helix and beta-sheet through self-assembly of two isomeric tetrapeptides in crystals and formation of filaments and ribbons in the solid state

Single crystal X-ray diffraction studies show that the beta-turn structure of tetrapeptide I, Boc-Gly-Phe-Aib-Leu-OMe (Aib: alpha-amino isobutyric acid) self-assembles to a supramolecular helix through intermolecular hydrogen bonding along the crystallographic a axis. By contrast the beta-turn structure of an isomeric tetrapeptide II, Boc-Gly-Leu-Aib-Phe-OMe self-assembles to a supramolecular beta-sheet-like structure via a two-dimensional (a, b axis) intermolecular hydrogen bonding network and pi-pi interactions. FT-IR studies of the peptides revealed that both of them form intermolecularly hydrogen bonded supramolecular structures in the solid state. Field emission scanning electron micrographs (FE-SEM) of the dried fibrous materials of the peptides show different morphologies, non-twisted filaments in case of peptide I and non-twisted filaments and ribbon-like structures in case of peptide II.

Design of a new foldamer turn-linker-turn in acyclic hexapeptides and formation of channels through self-assembly

Single crystal X-ray diffraction studies and solvent dependent NMR titration reveal that the designed pepticles I and 11, Boc-Xx(1)-Aib(2)-Yy(3)-NH(CH2)(2)NH-Yy(3)-Aib(2)-Xx(1)-Boc, where Xx and Yy are lie and Leu in peptide I and Leu and Val in peptide 11, respectively, fold into a turn-linker-turn (T-L-T) conformation both in the solid state and in solution. In the crystalline state the T-L-T foldamers; of peptide I and II self-assemble to form a three-dimensional framework of channels. The insides of the channels are hydrophilic and found to contain solvent CHCl3 hydrogen bonded to exposed C=O of Aib located at the turn regions. (c) 2008 Elsevier B.V. All rights reserved.

A new molecular scaffold for the formation of supramolecular peptide double helices: the crystallographic insight

A series of water-soluble synthetic dipeptides (1-3) with an N-terminally located beta-alanine residue, beta-alanyl-L-valine (1), beta-alanyl-L-isoleucine (2), and beta-alanyl-L-phenylalanine (3, form hydrogen-bonded supramolecular double helices with a pitch length of 1 nm, whereas the C-terminally positioned beta-alanine containing dipeptide (4), L-phenylalanyl-beta-alanine, does not form a supramolecular double helical structure. beta-Ala-Xaa (Xaa = Val/Ile/Phe) can be regarded as a new motif for the formation of supramolecular double helical structures in the solid state.

First crystallographic signature of an acyclic peptide nanorod: molecular mechanism of nanorod formation by a self-assembled tetrapeptide

A terminally protected acyclic tetrapeptide Boc-Aib-Val-Aib-beta-Ala-OMe 1 (Aib: alpha-aminoisobutyric acid, beta-Ala: beta-Alanine) self-assembles into a continuous hydrogen-bonded supramolecular helix with an average diameter of 10Angstrom (1nm) starting from a double bend molecular conformation in crystals and further self-assembly of this supramolecular architecture leads to the formation of polydisperse nanorods of diameters 10-40 nm.

Solvothermal syntheses, crystal structures and properties of five new thloantimonates(III) containing the [Sb4S7](2-) anion

Five new thioantimonates have been synthesized in the presence of organic amines under solvothermal conditions and their structures determined by single-crystal X-ray diffraction. All of the compounds are layered and contain antimony-sulphide anions of stoichiometry [Sb4S7](2-), but the structure of the anion formed is dependent on the amine used in synthesis. (H3N(CH2)(4)NH3)[Sb4S7] (1) contains [Sb4S7](2-) double chains directed along [010]. Weak interchain Sb-S interactions between neighbouring chains cause the double chains to pack into layers in the ab plane. In the [001] direction, the layers of double chains alternate with doubly protonated diaminobutane molecules to which the chains are hydrogen bonded. Compounds of general formula (TH)(2)[Sb4S7] (T= CH3(CH2)(2)NH2 (2), (CH3)(2)CHNH2 (3), CH3(CH2)(3)NH2 (4) and CH3(CH2)(4)NH2 (5)) adopt a more complex structure in which [Sb3S8](7-) units are linked by Sb-3(3-) pyramids to form chains, which in turn are bridged by sulphur atoms to create sheets containing large heterorings. Pairs of such sheets form double layers of four atoms thickness that are stacked along [001]. Protonated amine molecules are located between anionic antimony-sulphide layers to which they are hydrogen bonded. Thermal analysis reveals that the decomposition temperature of materials containing [Sb4S7](2-) anions is dependent both on the structure of the anion, the lowest decomposition temperature being that of the low-dimensional phase (1) and on the identity of the amine, the decomposition temperature decreasing with an increasing number of carbon atoms and decreasing density. (c) 2005 Elsevier Inc. All rights reserved.

A short water-soluble self-assembling peptide forms amyloid-like fibrils

A water-soluble tripeptide Val-Ile-Ala (VIA) 1, bearing sequence identity with the C-terminal portion of the Alzheimer A beta-peptide (A beta(40-42)), self-assembles, in crystalline form, to produce an intermolecularly hydrogen bonded supramolecular beta-sheet structure which self-associates to form straight, unbranched nanofibrils exhibiting amyloid-like behavior; in contrast, the synthetic tripeptide Ala-Val-Ile (AVI) 2 self-assembles to produce a beta-sheet structure that forms branched nanofibrils which do not show any characteristic features of amyloid-like fibrils.

The role of terminal tyrosine residues in the formation of tripeptide nanotubes: a crystallographic insight

Terminally protected acyclic tripeptides containing tyrosine residues at both termini self-assemble into nanotubes in crystals through various non-covalent interactions including intermolecular hydrogen bonds. The nanotube has an average internal diameter of 5 angstrom (0.5 nm) and the tubular ensemble is developed through the hydrogen-bonded phenolic-OH side chains of tyrosine (Tyr) residues [Org. Lett. 2004, 6, 4463]. We have synthesized and studied several tripeptides 3-6 to probe the role of tyrosine residues in nanotube structure formation. These peptides either have only one Tyr residue at N- or C-termini or they have one or two terminally located phenylalanine (Phe) residues. These tripeptides failed to form any kind of nanotubular structure in the solid state. Single crystal X-ray diffraction studies of these peptides 3-6 clearly demonstrate that substitution of any one of the terminal Tyr residues in the Boc-Tyr-X-Tyr-OMe (X=VaI or Ile) sequence disrupts the formation of the nanotubular structure indicating that the presence of two terminally located Tyr residues is vital for nanotube formation. (c) 2006 Elsevier Ltd. All rights reserved.

Supramolecular beta-sheet and nanofibril formation by self-assembling tripeptides containing an N-terminally located gamma-aminobutyric acid residue

Three terminally protected tripeptides Boc-gamma-Abu-Val-Leu-OMe 1, Boc-gamma-Abu-Leu-Phe-OMe 2 and Boc-gamma-Abu-Val-Tyr-OMe 3 (gamma-Abu = gamma-aminobutyric acid) each containing an N-terminally positioned gamma-aminobutyric acid residue have been synthesized, purified and studied. FT-IR studies of all these peptides revealed that these peptides form intermolecularly hydrogen bonded supramolecular beta-sheet structures. Peptides 1, 2 and 3 adopt extended backbone beta-strand molecular structures in crystals. Crystal packing of all these peptides demonstrates that these beta-strand structures self-assemble to form intermolecularly H-bonded parallel beta-sheet structures. Peptide 3 uses a side chain tyrosyl -OH group as an additional hydrogen bonding functionality in addition to the backbone CONH groups to pack in crystals. Transmission electron microscopic studies of all peptides indicate that they self-assemble to form nanofibrillar structures of an average diameter of 65 nm. These peptide fibrils exhibit amyloid-like behavior as they bind to a physiological dye Congo red and show a characteristic green-gold birefringence under polarizing microscope.