Search

We present a novel method for the calculation of the energy density of states D(E) for systems described by classical statistical mechanics. The method builds on an extension of a recently proposed strategy that allows the free-energy profile of a canonical system to be recovered within a preassigned accuracy [Proc. Natl. Acad. Sci. U.S.A. 99, 12562 (2002)]]. The method allows a good control over the error on the recovered system entropy. This fact is exploited to obtain D(E) more efficiently by combining measurements at different temperatures. The accuracy and efficiency of the method are tested for the two-dimensional Ising model (up to size 50 x 50) by comparison with both exact results and previous studies. This method is a general one and should be applicable to more realistic model systems.

Seven Brazilian Tay-Sachs disease cases were screened for the most frequent causative mutations. They all presented at least one copy of the IVS7+1g>c mutation. Three patients were homozygotes, three were compound heterozygotes, and in one case only the mother was tested and shown to carry the IVS7+1g>c mutation. In the second allele the compound heterozygotes presented: R178H (the DN allele), InsTATC1278 and an unidentified mutation. The IVS7+1g>c mutation has already been described in three Portuguese patients. In this study, all families were unaware of any Portuguese ancestry. Since Brazil was a Portuguese colony, the mutation most probably came from ancient common ancestry. The initial molecular analysis of Tay-Sachs disease patients in Brazil indicated a prevalence of the IVS7+1g>c mutation, possibly as a result of genetic drift.

A theoretical model for the folding of proteins containing disulfide bonds is introduced. The model exploits the knowledge of the native state to favor the progressive establishment of native interactions. At variance with traditional approaches based on native topology, not all native bonds are treated in the same way; in particular, a suitable energy term is introduced to account for the special strength of disulfide bonds, as well as their ability to undergo intramolecular reshuffling. The model thus possesses the minimal ingredients necessary to investigate the much debated issue of whether the refolding process occurs through partially structured intermediates with native or non-native disulfide bonds. This strategy is applied to a context of particular interest, the refolding process of hirudin, a thrombin-specific protease inhibitor, for which conflicting folding pathways have been proposed. We show that the only two parameters in the model (temperature and disulfide strength) can be tuned to reproduce well a set of experimental transitions between species with different number of formed disulfides. This model is then used to provide a characterization of the folding process and a detailed description of the species involved in the rate-limiting step of hirudin refolding., (Copyright 2003 Wiley-Liss, Inc.)

We characterize the equilibrium thermodynamics of a thick polymer confined in a spherical region of space. This is used to gain insight into the DNA packaging process. The experimental reference system for the present study is the recent characterization of the loading process of the genome inside the phi29 bacteriophage capsid. Our emphasis is on the modelling of double-stranded DNA as a flexible thick polymer (tube) instead of a beads-and-springs chain. By using finite-size scaling to extrapolate our results to genome lengths appropriate for phi29, we find that the thickness-induced force may account for up to half the one measured experimentally at high packing densities. An analogous agreement is found for the total work that has to be spent in the packaging process. Remarkably, such agreement can be obtained in the absence of any tunable parameters and is a mere consequence of the DNA thickness. Furthermore, we provide a quantitative estimate of how the persistence length of a polymer depends on its thickness. The expression accounts for the significant difference in the persistence lengths of single and double-stranded DNA (again with the sole input of their respective sections and natural nucleotide/base-pair spacing).

A variety of experimental and theoretical studies have established that the folding process of monomeric proteins is strongly influenced by the topology of the native state. In particular, folding times have been shown to correlate well with the contact order, a measure of contact locality. Our investigation focuses on identifying additional topologic properties that correlate with experimentally measurable quantities, such as folding rates and transition-state placement, for both two- and three-state folders. The validation against data from 40 experiments shows that a particular topological property that measures the interdependence of contacts, termed cliquishness or clustering coefficient, can account with statistically significant accuracy both for the transition state placement and especially for folding rates. The observed correlations can be further improved by optimally combining the distinct topological information captured by cliquishness and contact order., (Copyright 2003 Wiley-Liss, Inc.)

The relevance of various residue positions for the stability and the folding characteristics of the prion protein in its normal cellular form are investigated by using molecular dynamics simulations of models exploiting the topology of the native state. These models allow for reproducing the experimentally validated two-state behavior of the normal prion isoform. Highly significant correlations are found between the most topologically relevant sites in our analysis and the single point mutations known to be associated with the arousal of the genetic forms of prion disease. Insight into the conformational change is provided by comparing the folding process of cellular prion and doppel that share a similar native state topology: the folding pathways of the former can be grouped in two main classes according to which tertiary structure contacts are formed first enroute to the native state. For the latter a single class of pathways leads to the native state again through a two-state process. Our results are consistent and supportive of the recent experimental findings that doppel lacks the scrapie isoform and that such remarkably different behavior involves residues in the region containing the two beta-strands and the intervening helix.

We carry out a theoretical study of the vibrational and relaxation properties of naturally occurring proteins with the purpose of characterizing both the folding and equilibrium thermodynamics. By means of a suitable model, we provide a full characterization of the spectrum and eigenmodes of vibration at various temperatures by merely exploiting the knowledge of the protein native structure. It is shown that the rate at which perturbations decay at the folding transition correlates well with experimental folding rates. This validation is carried out on a list of about 30 two-state folders. Furthermore, the qualitative analysis of residues mean square displacements (shown to reproduce crystallographic data accurately) provides a reliable and statistically accurate method to identify crucial folding sites/contacts. This novel strategy is validated against clinical data for human immunodeficiency virus type 1 (HIV-1) protease. Finally, we compare the spectra and eigenmodes of vibration of natural proteins against randomly generated compact structures and regular random graphs. The comparison reveals a distinctive enhanced flexibility of natural structures accompanied by slow relaxation times at the folding temperature. The fact that these properties are connected intimately to the presence and assembly of secondary motifs hints at the special criteria adopted by evolution in the selection of viable folds.

An exactly solvable model based on the topology of a protein native state is applied to identify bottlenecks and key sites for the folding of human immunodeficiency virus type 1 (HIV-1) protease. The predicted sites are found to correlate well with clinical data on resistance to Food and Drug Administration-approved drugs. It has been observed that the effects of drug therapy are to induce multiple mutations on the protease. The sites where such mutations occur correlate well with those involved in folding bottlenecks identified through the deterministic procedure proposed in this study. The high statistical significance of the observed correlations suggests that the approach may be promisingly used in conjunction with traditional techniques to identify candidate locations for drug attacks.

A conceptual framework for understanding the protein folding problem has remained elusive in spite of many significant advances. We show that geometrical constraints imposed by chain connectivity, compactness, and the avoidance of steric clashes can be encompassed in a natural way using a three-body potential and lead to a selection in structure space, independent of chemical details. Strikingly, secondary motifs such as hairpins, sheets, and helices, which are the building blocks of protein folds, emerge as the chosen structures for segments of the protein backbone based just on elementary geometrical considerations., (Copyright 2002 Wiley-Liss, Inc.)

We introduce a simple theoretical approach for an equilibrium study of proteins with known native-state structures. We test our approach with results on well-studied globular proteins, chymotrypsin inhibitor (2ci2), barnase, and the alpha spectrin SH3 domain, and present evidence for a hierarchical onset of order on lowering the temperature with significant organization at the local level even at high temperatures. A further application to the folding process of HIV-1 protease shows that the model can be reliably used to identify key folding sites that are responsible for the development of drug resistance.

Drug resistance to HIV-1 protease involves the accumulation of multiple mutations in the protein. We investigate the role of these mutations by using molecular dynamics simulations that exploit the influence of the native-state topology in the folding process. Our calculations show that sites contributing to phenotypic resistance of FDA-approved drugs are among the most sensitive positions for the stability of partially folded states and should play a relevant role in the folding process. Furthermore, associations between amino acid sites mutating under drug treatment are shown to be statistically correlated. The striking correlation between clinical data and our calculations suggest a novel approach to the design of drugs tailored to bind regions crucial not only for protein function, but for folding as well., (Copyright 2001 Wiley-Liss, Inc.)

Menkes disease is a rare X-linked disorder related to a defect in the copper metabolism. According to the current literature, the most frequent neuroimaging findings are cortical atrophy, chronic subdural effusion or hygroma, and vascular abnormalities. White matter lesions may be present before other features of the disease and may evolve into atrophy. We hereby report a case of Menkes disease with typical history and progression, and an early phase imaging study with important white matter abnormalities, which could have lead to diagnostic difficulties.

The prediction of the three-dimensional structures of the native states of proteins from the sequences of their amino acids is one of the most important challenges in molecular biology. An essential task for solving this problem within coarse-grained models is the deduction of effective interaction potentials between the amino acids. Over the years, several techniques have been developed to extract potentials that are able to discriminate satisfactorily between the native and nonnative folds of a preassigned protein sequence. In general, when these potentials are used in actual dynamical folding simulations, they lead to a drift of the native structure outside the quasinative basin. In this article, we present and validate an approach to overcome this difficulty. By exploiting several numerical and analytical tools, we set up a rigorous iterative scheme to extract potentials satisfying a prerequisite of any viable potential: the stabilization of proteins within their native basin (less than 3-4 A RMSD). The scheme is flexible and is demonstrated to be applicable to a variety of parameterizations of the energy function, and it provides in each case the optimal potentials.

A simple and very efficient protein design strategy is proposed by developing some recently introduced theoretical tools which have been successfully applied to exactly solvable protein models. The design approach is implemented by using three amino acid classes and it is based on the minimization of an appropriate energy function. For a given native state the results of the design procedure are compared, through a statistical analysis, with the properties of an ensemble of sequences folding in the same conformation. If the success rate is computed on those sites designed with high confidence, it can be as high as 80%. The method is also able to identify key sites for the folding process: results for 2ci2 and barnase are in very good agreement with experimental results.

A novel scheme is introduced to capture the spatial correlations of consecutive amino acids in naturally occurring proteins. This knowledge-based strategy is able to carry out optimally automated subdivisions of protein fragments into classes of similarity. The goal is to provide the minimal set of protein oligomers (termed "oligons" for brevity) that is able to represent any other fragment. At variance with previous studies in which recurrent local motifs were classified, our concern is to provide simplified protein representations that have been optimised for use in automated folding and/or design attempts. In such contexts, it is paramount to limit the number of degrees of freedom per amino acid without incurring loss of accuracy of structural representations. The suggested method finds, by construction, the optimal compromise between these needs. Several possible oligon lengths are considered. It is shown that meaningful classifications cannot be done for lengths greater than six or smaller than four. Different contexts are considered for which oligons of length five or six are recommendable. With only a few dozen oligons of such length, virtually any protein can be reproduced within typical experimental uncertainties. Structural data for the oligons are made publicly available.

Optimal geometrical arrangements, such as the stacking of atoms, are of relevance in diverse disciplines. A classic problem is the determination of the optimal arrangement of spheres in three dimensions in order to achieve the highest packing fraction; only recently has it been proved that the answer for infinite systems is a face-centred-cubic lattice. This simply stated problem has had a profound impact in many areas, ranging from the crystallization and melting of atomic systems, to optimal packing of objects and the sub-division of space. Here we study an analogous problem--that of determining the optimal shapes of closely packed compact strings. This problem is a mathematical idealization of situations commonly encountered in biology, chemistry and physics, involving the optimal structure of folded polymeric chains. We find that, in cases where boundary effects are not dominant, helices with a particular pitch-radius ratio are selected. Interestingly, the same geometry is observed in helices in naturally occurring proteins.

A fascinating and open question challenging biochemistry, physics, and even geometry is the presence of highly regular motifs such as alpha helices in the folded state of biopolymers and proteins. Stimulating explanations ranging from chemical propensity to simple geometrical reasoning have been invoked to rationalize the existence of such secondary structures. We formulate a dynamical variational principle for selection in conformation space based on the requirement that the backbone of the native state of biologically viable polymers be rapidly accessible from the denatured state. The variational principle is shown to result in the emergence of helical order in compact structures.

A two amino acid (hydrophobic and polar) scheme is used to perform the design on target conformations corresponding to the native states of 20 single chain proteins. Strikingly, the percentage of successful identification of the nature of the residues benchmarked against naturally occurring proteins and their homologues is around 75%, independent of the complexity of the design procedure. Typically, the lowest success rate occurs for residues such as alanine that have a high secondary structure functionality. Using a simple lattice model, we argue that one possible shortcoming of the model studied may involve the coarse-graining of the 20 kinds of amino acids into just two effective types.

Matrices for electrokinetic separations, based on a unique class of mono- and disubstituted (on the amido nitrogen) acrylamides such as e.g., N-acryloylaminoethoxyethanol (AAEE) and acrylamido-N,N-diethoxyethanol, offer the following advantages: (i) strong resistance to alkaline hydrolysis (most zone separations occurring at basic pH values), (ii) high hydrophilicity and (iii) greater porosity, due to the higher molecular weight of the monomers. When compared with conventional poly(acrylamide), a poly(AAEE) matrix, when subjected to mild alkaline hydrolysis (0.1 N NaOH, 70 degrees C) appears to be 500 times more stable. Such stability is also confirmed under strong alkaline hydrolysis (1 N NaOH, 100 degrees C) as well as under mild and strong acidic hydrolysis. Mildly hydrolyzed poly(AAEE) matrices still perform extremely well in both conventional isoelectric focusing and immobilized pH gradients, techniques which are quite sensitive to traces of acrylate in the polymer coil. Conversely, mildly hydrolyzed poly(acrylamide) matrices, when used in isoelectric focusing, generate pH gradients between pH 4 and 5, having an inflection point (pH 4.6) equivalent to the pK value of acrylic acid. This novel class of monomers shows great promise for future applications in all electrokinetic methodologies.

1. The difference in planning and crafting strategy is the level of involvement in the process. It is the crafting of strategy that requires the involvement of state and local constituencies. 2. Creating vision enables constituencies to clarify and realize what they really want, independently of what presently seems possible at the state or local level. 3. By developing an ideal reality or vision first, a bridge can be built between the current and desired organization without loss of purpose. The purpose of the organization is defined by AAOHN's Mission Statement. 4. The bridge is built by using a decision making process similar to the nursing process.

Respiration rate, heart rate, and body temperature values were obtained from 14 cynomolgus monkeys (Macaca fascicularis) during neurosurgery under barbiturate anesthesia. Vital sign values markedly declined below baseline during the early stages of surgery, steadily increased as surgery progressed and neared completion, and finally returned to baseline by the end of the postsurgical recovery period. There was considerable variability among the 14 monkeys, but the ranking of each monkey relative to the others remained constant across the period of observation. The findings suggested that the cynomolgus monkey may be more sensitive to barbiturates than the rhesus monkey, and cynomolgus monkeys may exhibit considerable individual differences in their sensitivity to barbiturates.

Enzyme induction by isoniazid was studied by urinary D-glucaric acid estimation in slow and fast acetylators. Isoniazid administration increases significantively the D-glucaric acid elimination in the two classes of patients. In fast acetylators, the glucaricaciduria increases regularly up to 1.8 fold the physiological level in 30 days. In slow acetylators, after a progressive elevation during 20 days, the glucaricaciduria reaches quickly a 3 fold increase after 30 days. When a classic inducer such as phenobarbital is administrated in association with isoniazid, induction is stimulated in the two groups of patients. The estimation of plasmatic free isoniazid seems to indicate that the acetylation rate of isoniazid is not on the dependance of the induction process when the drug is administrated alone. In contrast, this rate increases when the association isoniazid-phenobarbital is administrated to slow acetylators. Possible consequences on the hepatic toxicity of isoniazid are discussed.