Molecular tailoring approach for geometry optimization of large molecules: energy evaluation and parallelization strategies

A linear-scaling scheme for estimating the electronic energy, gradients, and Hessian of a large molecule at ab initio level of theory based on fragment set cardinality is presented. With this proposition, a general, cardinality-guided molecular tailoring approach (CG-MTA) for ab initio geometry optimization of large molecules is implemented. The method employs energy gradients extracted from fragment wave functions, enabling computations otherwise impractical on PC hardware. Further, the method is readily amenable to large scale coarse-grain parallelization with minimal communication among nodes, resulting in a near-linear speedup. CG-MTA is applied for density-functional-theory-based geometry optimization of a variety of molecules including alpha-tocopherol, taxol, gamma-cyclodextrin, and two conformations of polyglycine. In the tests performed, energy and gradient estimates obtained from CG-MTA during optimization runs show an excellent agreement with those obtained from actual computation. Accuracy of the Hessian obtained employing CG-MTA provides good hope for the application of Hessian-based geometry optimization to large molecules.

Filamin A mutations cause periventricular heterotopia with Ehlers-Danlos syndrome

OBJECTIVE

To define the clinical, radiologic, and genetic features of periventricular heterotopia (PH) with Ehlers-Danlos syndrome (EDS).

METHODS

Exonic sequencing and single stranded conformational polymorphism (SSCP) analysis was performed on affected individuals. Linkage analysis using microsatellite markers on the X-chromosome was performed on a single pedigree. Western blotting evaluated for loss of filamin A (FLNA) protein and Southern blotting assessed for any potential chromosome rearrangement in this region.

RESULTS

The authors report two familial cases and nine additional sporadic cases of the EDS-variant form of PH, which is characterized by nodular brain heterotopia, joint hypermobility, and development of aortic dilatation in early adulthood. MRI typically demonstrated bilateral nodular PH, indistinguishable from PH due to FLNA mutations. Exonic sequencing or SSCP analyses of FLNA revealed a 2762 delG single base pair deletion in one affected female. Another affected female harbored a C116 single point mutation, resulting in an A39G change. A third affected female had a 4147 delG single base pair deletion. One pedigree with no detectable exonic mutation demonstrated positive linkage to the FLNA locus Xq28, an affected individual in this family also had no detectable FLNA protein, but no chromosomal rearrangement was detected.

CONCLUSION

These results suggest that the Ehlers-Danlos variant of periventricular heterotopia (PH), in part, represents an overlapping syndrome with X-linked dominant PH due to filamin A mutations.

Crystallographic structures of thrombin complexed with thrombin receptor peptides: existence of expected and novel binding modes

Many of the vital actions of thrombin on platelets and other cells appear to be mediated by the recently cloned seven-transmembrane-domain thrombin receptor. Thrombin activates this receptor by a novel proteolytic mechanism. The amino-terminal exodomain of the receptor contains the sequence LDPRSFLLRNPNDKYEPF. Structure-activity studies with mutant receptors and receptor peptides suggest that this sequence binds to thrombin at two sites: LDPR with the active center of thrombin and KYEPF with the fibrinogen recognition exosite of thrombin. Thrombin then cleaves the Arg41-Ser42 bond to unmask a new amino terminus, which functions as a tethered peptide ligand binding to as yet undefined sites within the body of the receptor to effect receptor activation. We have determined eight crystal structures of thrombin complexed with receptor-based peptides. Each of the two components of the bidentate docking model was captured in individual cocrystals. In one crystal type, the LDPR sequence docked in the active center of thrombin in a manner analogous to d-PheProArg chloromethyl ketone. In other crystals, the KYEPF sequence bound in the fibrinogen anion binding exosite of thrombin in a manner analogous to the DFEEI sequence of the carboxylate-terminal peptide of hirudin. Strikingly, however, generation of a single crystal that includes both components of the anticipated bidentate binding mode was not achieved, apparently because the peptides have a dominant solution S-like conformation that does not bind in a productive way at the active center. This peptide structure apparently favored a novel alternative mode of receptor peptide-thrombin interaction in which the receptor peptides formed an intermolecular bridge between neighboring thrombin molecules, resulting in an infinite peptide thrombin chain in crystals. In this structure, the KYEPF sequence docked in the expected manner at the exosite of one thrombin molecule, but the LDPR sequence docked in an unusual nonproductive mode with the active center of a neighboring molecule. Mutations that removed important determinants of the S-like receptor peptide structure underlying the bridging mode in the receptor itself did not significantly alter thrombin signaling. Additionally, a comparison of receptor density to the responsiveness of a cell did not support a role for receptor oligomerization in signaling. The physiological role for this unexpected intermolecular binding mode, if any, remains to be identified.(ABSTRACT TRUNCATED AT 400 WORDS)

Crystal structure of peanut lectin, a protein with an unusual quaternary structure

The x-ray crystal structure of the tetrameric T-antigen-binding lectin from peanut, M(r) 110,000, has been determined by using the multiple isomorphous replacement method and refined to an R value of 0.218 for 22,155 reflections within the 10- to 2.95-A resolution range. Each subunit has essentially the same characteristic tertiary fold that is found in other legume lectins. The structure, however, exhibits an unusual quaternary arrangement of subunits. Unlike other well-characterized tetrameric proteins with identical subunits, peanut lectin has neither 222 (D2) nor fourfold (C4) symmetry. A noncrystallographic twofold axis relates two halves of the molecule. The two monomers in each half are related by a local twofold axis. The mutual disposition of the axes is such that they do not lead to a closed point group. Furthermore, the structure of peanut lectin demonstrates that differences in subunit arrangement in legume lectins could be due to factors intrinsic to the protein molecule and, contrary to earlier suggestions, are not necessarily caused by interactions involving covalently linked sugar. The structure provides a useful framework for exploring the structural basis and the functional implications of the variability in the subunit arrangement in legume lectins despite all of them having nearly the same subunit structure, and also for investigating the general problem of "open" quaternary assembly in oligomeric proteins.

Self-assembled monolayers (SAMs) of alkoxycyanobiphenyl thiols on gold—A study of electron transfer reaction using cyclic voltammetry and electrochemical impedance spectroscopy

Self-assembled monolayers (SAMs) of liquid crystalline thiol-terminated alkoxycyanobiphenyl molecules with different alkyl chain lengths on Au surface have been studied for the first time using electrochemical techniques such as cyclic voltammetry (CV) and electrochemical impedance spectroscopy (EIS). The barrier property of the SAM-modified surfaces was evaluated using two different redox probes, namely potassium ferro/ferri cyanide and hexaammineruthenium(III) chloride. It was found that for short-length alkyl chain thiol (C5) the electron transfer reaction of hexaammineruthenium(III) chloride takes place through tunneling mechanism. In contrast, redox reaction of potassium ferro/ferri cyanide is almost completely blocked by the SAM-modified Au surface. From the impedance data, a surface coverage value of >99.9% was calculated for all the thiol molecules.

Enabling ab initio Hessian and frequency calculations of large molecules

A linear scaling method, termed as cardinality guided molecular tailoring approach, is applied for the estimation of the Hessian matrix and frequency calculations of spatially extended molecules. The method is put to test on a number of molecular systems largely employing the Hartree-Fock and density functional theory for a variety of basis sets. To demonstrate its ability for correlated methods, we have also performed a few test calculations at the Moller-Plesset second order perturbation theory. A comparison of central processing unit and memory requirements for medium-sized systems with those for the corresponding full ab initio computation reveals substantial gains with negligible loss of accuracy. The technique is further employed for a set of larger molecules, Hessian and frequency calculations of which are not possible on commonly available personal-computer-type hardware.

Vertical Alignment of Liquid Crystals Over a Functionalized Flexible Substrate

A simple and effective approach for vertical alignment of liquid crystals (LCs) over a functionalized transparent flexible substrate is described. Surface characterization of this commercially available plastic substrate through X-ray photoelectron spectroscopy (XPS) and attenuated total reflection infrared spectroscopy (ATR-IR) indicated that cellulose acetate is main component of the transparent substrate. This substrate was chemically functionalized with a suitable LC compound. A trimethoxysilane terminated new rod-shaped mesogen is synthesized and covalently attached to the pre-treated film through silane condensation reaction. LC functionalization of the polymer film is confirmed through contact angle (CA), atomic force microscopy (AFM), XPS and ATR-IR spectroscopy studies. Versatility of the LC modified flexible substrates for the alignment of bulk LC sample at substrate-LC interface was assessed for nematic (N) and smectic A (SmA) phases. Remarkably, LC functionalized cellulose acetate films were found to be highly efficient in assisting a perfect homeotropic alignment of LCs (for both, a room temperature N and a high temperature SmA phase) over the entire area of the LC sample under observation indicating their superior aligning ability in comparison to their unmodified and octadecyltrimethoxysilane (OTS) modified counterparts. The demonstrated method of surface modification of flexible polymer film is easy, surface modified substrates are stable for several months, retained their aligning ability intact and more importantly they are reusable with maximum delivery.

Many-body interaction analysis: Algorithm development and application to large molecular clusters

A completely automated algorithm for performing many-body interaction energy analysis of clusters (MBAC) [M. J. Elrodt and R. J. Saykally, Chem. Rev. 94, 1975 (1994); S. S. Xantheas, J. Chem. Phys. 104, 8821 (1996)] at restricted Hartree-Fock (RHF)/MA Plesset 2nd order perturbation theory (MP2)/density functional theory (DFT) level of theory is reported. Use of superior guess density matrices (DM's) for smaller fragments generated from DM of the parent system and elimination of energetically insignificant higher-body combinations, leads to a more efficient performance (speed-up up to 2) compared to the conventional procedure. MBAC approach has been tested out on several large-sized weakly bound molecular clusters such as (H(2)O)(n), n=8, 12, 16, 20 and hydrated clusters of amides and aldehydes. The MBAC results indicate that the amides interact more strongly with water than aldehydes in these clusters. It also reconfirms minimization of the basis set superposition error for large cluster on using superior quality basis set. In case of larger weakly bound clusters, the contributions higher than four body are found to be repulsive in nature and smaller in magnitude. The reason for this may be attributed to the increased random orientations of the interacting molecules separated from each other by large distances.

A novel form of pontocerebellar hypoplasia maps to chromosome 7q11-21

OBJECTIVE

To describe a novel form of pontocerebellar hypoplasia (PCH) and map its genetic locus.

BACKGROUND

PCH is a heterogeneous group of disorders that are characterized by abnormally small cerebellum and brainstem. Autosomal recessive inheritance has been implied in many cases, but no genetic loci have been mapped to date.

METHODS

The authors studied a consanguineous family from the Sultanate of Oman with three siblings with a novel form of PCH. The authors performed clinical studies and linkage analysis of this pedigree.

RESULTS

The clinical features of the affected children include developmental delay, progressive microcephaly with brachycephaly, seizures during the first year of life, hypotonia with hyperreflexia, short stature, and optic atrophy. Imaging studies showed a small pons and cerebellum, prominent sulci and lateral ventricles, and decreased cerebral white matter volume. A lack of dyskinesias distinguishes this pedigree from PCH type 2. Genetic studies of this family revealed evidence of significant linkage to chromosome 7q11-21 (maximum multipoint lod score 3.23).

CONCLUSIONS

This pedigree represents a novel form of autosomal recessive PCH, which the authors propose to call cerebellar atrophy with progressive microcephaly (CLAM). This disorder maps to chromosome 7q11-21, and this locus was named CLAM. This report represents the first identification of a genetic locus for PCH.

Facile and simultaneous synthesis of graphene quantum dots and reduced graphene oxide for bio-imaging and supercapacitor applications

A simple and facile method for the simultaneous preparation of graphene quantum dots (GQDs) having different emission colours, viz., yellow, green and blue, and reduced graphene oxide (RGO) utilized respectively for bio-imaging and supercapacitor applications is demonstrated.

Scanning Tunneling Microscopy, Fourier Transform Infrared Spectroscopy, and Electrochemical Characterization of 2-Naphthalenethiol Self-Assembled Monolayers on the Au Surface: A Study of Bridge-Mediated Electron Transfer in Ru(NH3)62+|Ru(NH3)63+ Redox Reactions

We have studied the structure, adsorption kinetics, and barrier properties of self-assembled monolayers of 2-naphthalenethiol on Au using electrochemical techniques, grazing-angle Fourier transform infrared (FTIR) spectroscopy, and scanning tunneling microscopy (STM). The results of cyclic voltammetric and impedance measurements using redox probes show that 2-naphthalenethiol on Au forms a stable and reproducible, but moderately blocking, monolayer. Annealing of the self-assembled monolayer (SAM)-modified surface at 72 +/- 2 degrees C remarkably improves the blocking property of the monolayer of 2-naphthalenethiol on Au. From the study of kinetics of SAM formation, we find that the self-assembly follows Langmuir adsorption isotherm. Our STM and FTIR results show that the molecules are adsorbed with the naphthalene ring tilted from the surface normal by forming a square root 3 x 3 R30 degrees overlayer structure. From our studies, we conclude that the electron-transfer reaction of ferro/ferricyanide in the freshly formed monolayer occurs predominantly through the pinholes and defects present in the monolayer. However, in the case of thermally annealed specimen, although the ferro/ferricyanide reaction is almost completely blocked, the electron-transfer reaction of hexaammineruthenium(III) chloride is not significantly inhibited. It is proposed that the electron-transfer reaction in the case of the ruthenium redox couple takes place by a tunneling mechanism through the high-electron-density aromatic naphthalene ring acting as a bridge between the monolayer-modified electrode and the ruthenium complex.

Contribution of regulatory T cells to alleviation of experimental allergic asthma after specific immunotherapy

BACKGROUND

Allergen-specific immunotherapy (SIT) has been used since 1911, yet its mechanism of action remains to be elucidated. There is evidence indicating that CD4(+)FOXP3(+) regulatory T cells (Treg cells) are induced during SIT in allergic patients. However, the contribution of these cells to SIT has not been evaluated in vivo.

OBJECTIVE

To evaluate the in vivo contribution of (i) CD4(+) CD25(+) T cells during SIT and of (ii) SIT-generated inducible FOXP3(+) Treg cells during allergen exposure to SIT-mediated suppression of asthmatic manifestations.

METHODS

We used a mouse model of SIT based on the classical OVA-driven experimental asthma. Treg cells were quantified by flow cytometry 24 and 96 h post SIT treatment. We depleted CD4(+) CD25(+) T cells prior to SIT, and CD4(+)FOXP3(+) T cells prior to allergen challenges to study their contribution to the suppression of allergic manifestations by SIT treatment.

RESULTS

Our data show that depletion of CD4(+)CD25(+) T cells at the time of SIT treatment reverses the suppression of airway hyperresponsiveness (AHR), but not of airway eosinophilia and specific IgE levels in serum. Interestingly, the number of CD4(+)CD25(+)FOXP3(+) T cells is transiently increased after SIT in the spleen and blood, suggesting the generation of inducible and presumably allergen-specific Treg cells during treatment. Depletion of CD4(+)FOXP3(+) Treg cells after SIT treatment partially reverses the SIT-induced suppression of airway eosinophilia, but not of AHR and serum levels of specific IgE.

CONCLUSION AND CLINICAL RELEVANCE

We conclude that SIT-mediated tolerance induction towards AHR requires CD4(+)CD25(+) T cells at the time of allergen injections. In addition, SIT generates CD4(+)CD25(+)FOXP3(+) T cells that contribute to the suppression of airway eosinophilia upon allergen challenges. Therefore, enhancing Treg cell number or their activity during and after SIT could be of clinical relevance to improve the therapeutic effects of SIT.

Detection of the neurotransmitter dopamine by a glassy carbon electrode modified with self-assembled perovskite LaFeO3 microspheres made up of nanospheres

Glassy carbon electrode modified with LaFeO₃ microspheres made up of nanospheres and having high electrocatalytic ability to detect the neurotransmitter dopamine.

Comparison of the structures of the cyclotheonamide A complexes of human alpha-thrombin and bovine beta-trypsin

Thrombin, a trypsin-like serine protease present in blood, plays a central role in the regulation of thrombosis and hemostasis. A cyclic pentapeptide, cyclotheonamide A (CtA), isolated from sponges of the genus Theonella, inhibits thrombin, trypsin, and certain other serine proteases. Enzyme inhibition data for CtA indicate that it is a moderate inhibitor of alpha-thrombin (K(i) = 1.0 nM), but substantially more potent toward trypsin (K(i) = 0.2 nM). The comparative study of the crystal structures of the CtA complexes of alpha-thrombin and beta-trypsin reported here focuses on structure-function relationships in general and the enhanced specificity of trypsin, in particular. The crystal structures of the CtA complexes of thrombin and trypsin were solved and refined at 1.7 and 2.0 A resolution, respectively. The structures show that CtA occupies the active site with the Pro-Arg motif positioned in the S2 and S1 binding sites. The alpha-keto group of CtA is involved in a tetrahedral intermediate hemiketal structure with Ser 195 OG of the catalytic triad and is positioned within bonding distance from, and orthogonal to, the re-face of the carbonyl of the arginine of CtA. As in other productive binding modes of serine proteases, the Ser 214-Gly 216 segment runs in a twisted antiparallel beta-strand manner with respect to the diaminopropionic acid (Dpr)-Arg segment of CtA. The Tyr 60A-Thr 60I insertion loop of thrombin makes a weak aromatic stacking interaction with the v-Tyr of CtA through Trp 60D. The Glu 39 Tyr and Leu 41 Phe substitutions in trypsin produce an enhanced aromatic interaction with D-Phe of CtA, which also leads to different orientations of the side chains of D-Phe and the v-Tyr. The comparison of the CtA complexes of thrombin and trypsin shows that the gross structural features of both in the active site region are the same, whereas the differences observed are mainly due to minor insertions and substitutions. In trypsin, the substitution of Ile 174-Arg 175 by Gly 174-Gln 175 makes the S3 aryl site more polar because the Arg 175 side chain is directed away from thrombin and into the solvent, whereas Gln 175 is not. Because the site is occupied by the Dpr group of CtA, the occupancy of the S3 site is better in trypsin than in thrombin. In trypsin, the D-Phe side chain of CtA fits between Tyr 39 and Phe 41 in a favorable manner, whereas in thrombin, these residues are Glu 39 and Leu 41. The higher degree of specificity for trypsin is most likely the result of these substitutions and the absence of the fairly rigid Tyr 60A-Thr 60I insertion loop of thrombin, which narrows access to the active site and forces less favorable orientations for the D-Phe and v-Tyr residues.

Simple and facile preparation of silver–polydopamine (Ag–PDA) core–shell nanoparticles for selective electrochemical detection of cysteine

A simple one-step method for the preparation of silver–polydopamine (Ag–PDA) core–shell nanoparticles is proposed and its application for non-enzymatic electrochemical detection of cysteine is demonstrated.