Energy-based de novo protein folding by conformational space annealing and an off-lattice united-residue force field: Application to the 10-55 fragment of staphylococcal protein A and to apo calbindin D9K

The conformational space annealing (CSA) method for global optimization has been applied to the 10-55 fragment of the B-domain of staphylococcal protein A (protein A) and to a 75-residue protein, apo calbindin D9K (PDB ID code 1CLB), by using the UNRES off-lattice united-residue force field. Although the potential was not calibrated with these two proteins, the native-like structures were found among the low-energy conformations, without the use of threading or secondary-structure predictions. This is because the CSA method can find many distinct families of low-energy conformations. Starting from random conformations, the CSA method found that there are two families of low-energy conformations for each of the two proteins, the native-like fold and its mirror image. The CSA method converged to the same low-energy folds in all cases studied, as opposed to other optimization methods. It appears that the CSA method with the UNRES force field, which is based on the thermodynamic hypothesis, can be used in prediction of protein structures in real time.

A quantitative method for evaluating the stabilities of nucleic acids

We report a general method for screening, in solution, the impact of deviations from canonical Watson-Crick composition on the thermodynamic stability of nucleic acid duplexes. We demonstrate how fluorescence resonance energy transfer (FRET) can be used to detect directly free energy differences between an initially formed “reference” duplex (usually a Watson-Crick duplex) and a related “test” duplex containing a lesion/alteration of interest (e.g., a mismatch, a modified, a deleted, or a bulged base, etc.). In one application, one titrates into a solution containing a fluorescently labeled, FRET-active, reference duplex, an unlabeled, single-stranded nucleic acid (test strand), which may or may not compete successfully to form a new duplex. When a new duplex forms by strand displacement, it will not exhibit FRET. The resultant titration curve (normalized fluorescence intensity vs. logarithm of test strand concentration) yields a value for the difference in stability (free energy) between the newly formed, test strand-containing duplex and the initial reference duplex. The use of competitive equilibria in this assay allows the measurement of equilibrium association constants that far exceed the magnitudes accessible by conventional titrimetric techniques. Additionally, because of the sensitivity of fluorescence, the method requires several orders of magnitude less material than most other solution methods. We discuss the advantages of this method for detecting and characterizing any modification that alters duplex stability, including, but not limited to, mutagenic lesions. We underscore the wide range of accessible free energy values that can be defined by this method, the applicability of the method in probing for a myriad of nucleic acid variations, such as single nucleotide polymorphisms, and the potential of the method for high throughput screening.

Randomised controlled trial of the READER method of critical appraisal in general practice

Objective: To evaluate the READER model for critical reading by comparing it with a free appraisal, and to explore what factors influence different components of the model.

Recent improvements in prediction of protein structure by global optimization of a potential energy function

Recent improvements of a hierarchical ab initio or de novo approach for predicting both α and β structures of proteins are described. The united-residue energy function used in this procedure includes multibody interactions from a cumulant expansion of the free energy of polypeptide chains, with their relative weights determined by Z-score optimization. The critical initial stage of the hierarchical procedure involves a search of conformational space by the conformational space annealing (CSA) method, followed by optimization of an all-atom model. The procedure was assessed in a recent blind test of protein structure prediction (CASP4). The resulting lowest-energy structures of the target proteins (ranging in size from 70 to 244 residues) agreed with the experimental structures in many respects. The entire experimental structure of a cyclic α-helical protein of 70 residues was predicted to within 4.3 Å α-carbon (Cα) rms deviation (rmsd) whereas, for other α-helical proteins, fragments of roughly 60 residues were predicted to within 6.0 Å Cα rmsd. Whereas β structures can now be predicted with the new procedure, the success rate for α/β- and β-proteins is lower than that for α-proteins at present. For the β portions of α/β structures, the Cα rmsd's are less than 6.0 Å for contiguous fragments of 30–40 residues; for one target, three fragments (of length 10, 23, and 28 residues, respectively) formed a compact part of the tertiary structure with a Cα rmsd less than 6.0 Å. Overall, these results constitute an important step toward the ab initio prediction of protein structure solely from the amino acid sequence.

On a psychophysical transformed-rule up and down method converging on a 75% level of correct responses

Tranformed-rule up and down psychophysical methods have gained great popularity, mainly because they combine criterion-free responses with an adaptive procedure allowing rapid determination of an average stimulus threshold at various criterion levels of correct responses. The statistical theory underlying the methods now in routine use is based on sets of consecutive responses with assumed constant probabilities of occurrence. The response rules requiring consecutive responses prevent the possibility of using the most desirable response criterion, that of 75% correct responses. The earliest transformed-rule up and down method, whose rules included nonconsecutive responses, did not contain this limitation but failed to become generally accepted, lacking a published theoretical foundation. Such a foundation is provided in this article and is validated empirically with the help of experiments on human subjects and a computer simulation. In addition to allowing the criterion of 75% correct responses, the method is more efficient than the methods excluding nonconsecutive responses in their rules.

Top-down free-energy minimization on protein potential energy landscapes

The hierarchical properties of potential energy landscapes have been used to gain insight into thermodynamic and kinetic properties of protein ensembles. It also may be possible to use them to direct computational searches for thermodynamically stable macroscopic states, i.e., computational protein folding. To this end, we have developed a top-down search procedure in which conformation space is recursively dissected according to the intrinsic hierarchical structure of a landscape's effective-energy barriers. This procedure generates an inverted tree similar to the disconnectivity graphs generated by local minima-clustering methods, but it fundamentally differs in the manner in which the portion of the tree that is to be computationally explored is selected. A key ingredient is a branch-selection algorithm that takes advantage of statistically predictive properties of the landscape to guide searches down the tree branches that are most likely to lead to the physically relevant macroscopic states. Using the computational folding of a β-hairpin-forming peptide as an example, we show that such predictive properties indeed exist and can be used for structure prediction by free-energy global minimization.

Dynamical hierarchy in transition states: Why and how does a system climb over the mountain?

How a reacting system climbs through a transition state during the course of a reaction has been an intriguing subject for decades. Here we present and quantify a technique to identify and characterize local invariances about the transition state of an N-particle Hamiltonian system, using Lie canonical perturbation theory combined with microcanonical molecular dynamics simulation. We show that at least three distinct energy regimes of dynamical behavior occur in the region of the transition state, distinguished by the extent of their local dynamical invariance and regularity. Isomerization of a six-atom Lennard–Jones cluster illustrates this: up to energies high enough to make the system manifestly chaotic, approximate invariants of motion associated with a reaction coordinate in phase space imply a many-body dividing hypersurface in phase space that is free of recrossings even in a sea of chaos. The method makes it possible to visualize the stable and unstable invariant manifolds leading to and from the transition state, i.e., the reaction path in phase space, and how this regularity turns to chaos with increasing total energy of the system. This, in turn, illuminates a new type of phase space bottleneck in the region of a transition state that emerges as the total energy and mode coupling increase, which keeps a reacting system increasingly trapped in that region.

Single photon emission computerized tomography imaging of amphetamine-induced dopamine release in drug-free schizophrenic subjects.

The dopamine hypothesis of schizophrenia proposes that hyperactivity of dopaminergic transmission is associated with this illness, but direct observation of abnormalities of dopamine function in schizophrenia has remained elusive. We used a newly developed single photon emission computerized tomography method to measure amphetamine-induced dopamine release in the striatum of fifteen patients with schizophrenia and fifteen healthy controls. Amphetamine-induced dopamine release was estimated by the amphetamine-induced reduction in dopamine D2 receptor availability, measured as the binding potential of the specific D2 receptor radiotracer [123I] (S)-(-)-3-iodo-2-hydroxy-6-methoxy-N-[(1-ethyl-2-pyrrolidinyl) methyl]benzamide ([123I]IBZM). The amphetamine-induced decrease in [123I]IBZM binding potential was significantly greater in the schizophrenic group (-19.5 +/- 4.1%) compared with the control group (-7.6 +/- 2.1%). In the schizophrenic group, elevated amphetamine effect on [123I]IBZM binding potential was associated with emergence or worsening of positive psychotic symptoms. This result suggests that psychotic symptoms elicited in this experimental setting in schizophrenic patients are associated with exaggerated stimulation of dopaminergic transmission. Such an observation would be compatible with an abnormal responsiveness of dopaminergic neurons in schizophrenia.

A gain-of-function of an amyotrophic lateral sclerosis-associated Cu,Zn-superoxide dismutase mutant: An enhancement of free radical formation due to a decrease in Km for hydrogen peroxide.

Cu,Zn-superoxide dismutase (SOD) is known to be a locus of mutation in familial amyotrophic lateral sclerosis (FALS). Transgenic mice that express a mutant Cu,Zn-SOD, Gly-93--> Ala (G93A), have been shown to develop amyotrophic lateral sclerosis (ALS) symptoms. We cloned the FALS mutant, G93A, and wild-type cDNA of human Cu,Zn-SOD, overexpressed them in Sf9 insect cells, purified the proteins, and studied their enzymic activities for catalyzing the dismutation of superoxide anions and the generation of free radicals with H2O2 as substrate. Our results showed that both enzymes contain one copper ion per subunit and have identical dismutation activity. However, the free radical-generating function of the G93A mutant, as measured by the spin trapping method, is enhanced relative to that of the wild-type enzyme, particularly at lower H2O2 concentrations. This is due to a small, but reproducible, decrease in the value of Km for H2O2 for the G93A mutant, while the kcat is identical for both enzymes. Thus, the ALS symptoms observed in G93A transgenic mice are not caused by the reduction of Cu,Zn-SOD activity with the mutant enzyme; rather, it is induced by a gain-of-function, an enhancement of the free radical-generating function. This is consistent with the x-ray crystallographic studies showing the active channel of the FALS mutant is slightly larger than that of the wild-type enzyme; thus, it is more accessible to H2O2. This gain-of-function, in part, may provide an explanation for the association between ALS and Cu,Zn-SOD mutants.

Thermodynamic parameters for the helix-coil transition of oligopeptides: molecular dynamics simulation with the peptide growth method.

The helix-coil transition equilibrium of polypeptides in aqueous solution was studied by molecular dynamics simulation. The peptide growth simulation method was introduced to generate dynamic models of polypeptide chains in a statistical (random) coil or an alpha-helical conformation. The key element of this method is to build up a polypeptide chain during the course of a molecular transformation simulation, successively adding whole amino acid residues to the chain in a predefined conformation state (e.g., alpha-helical or statistical coil). Thus, oligopeptides of the same length and composition, but having different conformations, can be incrementally grown from a common precursor, and their relative conformational free energies can be calculated as the difference between the free energies for growing the individual peptides. This affords a straightforward calculation of the Zimm-Bragg sigma and s parameters for helix initiation and helix growth. The calculated sigma and s parameters for the polyalanine alpha-helix are in reasonable agreement with the experimental measurements. The peptide growth simulation method is an effective way to study quantitatively the thermodynamics of local protein folding.

E-polynomial of the SL(3,C)-character variety of free groups.

We compute the E-polynomial of the character variety of representations of a rank r free group in SL(3,C). Expanding upon techniques of Logares, Muñoz and Newstead (Rev. Mat. Complut. 26:2 (2013), 635-703), we stratify the space of representations and compute the E-polynomial of each geometrically described stratum using fibrations. Consequently, we also determine the E-polynomial of its smooth, singular, and abelian loci and the corresponding Euler characteristic in each case. Along the way, we give a new proof of results of Cavazos and Lawton (Int. J. Math. 25:6 (2014), 1450058).

Tethered Monte Carlo: computing the effective potential without critical slowing down

We present Tethered Monte Carlo, a simple, general purpose method of computing the effective potential of the order parameter (Helmholtz free energy). This formalism is based on a new statistical ensemble, closely related to the micromagnetic one, but with an extended configuration space (through Creutz-like demons). Canonical averages for arbitrary values of the external magnetic field are computed without additional simulations. The method is put to work in the two-dimensional Ising model, where the existence of exact results enables us to perform high precision checks. A rather peculiar feature of our implementation, which employs a local Metropolis algorithm, is the total absence, within errors, of critical slowing down for magnetic observables. Indeed, high accuracy results are presented for lattices as large as L = 1024.

Tethered Monte Carlo: computing the effective potential without critical slowing down

We present Tethered Monte Carlo, a simple, general purpose method of computing the effective potential of the order parameter (Helmholtz free energy). This formalism is based on a new statistical ensemble, closely related to the micromagnetic one, but with an extended configuration space (through Creutz-like demons). Canonical averages for arbitrary values of the external magnetic field are computed without additional simulations. The method is put to work in the two-dimensional Ising model, where the existence of exact results enables us to perform high precision checks. A rather peculiar feature of our implementation, which employs a local Metropolis algorithm, is the total absence, within errors, of critical slowing down for magnetic observables. Indeed, high accuracy results are presented for lattices as large as L = 1024.

Free living astigmatid mites (Astigmatina): new taxa, rearing and use for mesostigmatid (Mesostigmata) predatory mite production

The cohort Astigmatina is divided in two major groups: Psoroptidia, composed mainly by feather and fur mites, and Non-psoroptidia, a dominant component of the acarofauna in ephemeral habitats. In these environments Astigmatina usually are saprophages or feed on fungi or bacteria. Astigmatina protonymphs undergo a complete reorganization of the body structure leading to the production of heteromorphic deutonymphs, generally specialized for dispersion through phoresy using arthropods and vertebrates as phoronts. Although most Astigmatina occur in natural environments, some species live in anthropic environments, such as food deposits, where some of them became pests; some Astigmatina infest subterraneous plant organs. Despite their economic and ecological importance, studies on the diversity and taxonomy of Astigmatina in Brazil have been rare over the last decades. The general objective of this thesis was to collaborate to the knowledge of the diversity and to evaluate the potential practical uses of these mites in Brazil. For this, new genera and species were described, method for rearing dust mites was studied and the efficiency of Astigmatina as prey for edaphic predators was evaluated. A new species of Thyreophagus (Astigmatina: Acaridae) was described based on specimens collected in Brazil, the association of three other species of this genus with stored food was reviewed and a key to all species of this genus was prepared. The genus Neotropacarus (Astigmatina: Acaridae), commonly found on plant leaves, was reviewed with the redescription of two species and description of new species collected in Brazil and from the Philippines. Two new genera and seven new species of Acaridae associated with the bee family Apidae was described and a key to Acaridae genera in subfamily Horstiinae was prepared. Several species of Astigmatina were evaluated as prey for predatory mites Stratiolaelaps scimitus (Womersley) (Mesostigmata: Laelapidae) and Protogamasellopsis zaheri Abo-Shnaf, Castilho and Moraes (Mesostigmata: Rhodacaridae), which oviposited on all evaluated astigmatids, with Tyrophagus putrescentiae (Schrank) and Aleuroglyphus ovatus (Tropeau) (Acaridae) being the most suitable prey. Seven foods and two development period, 30 and 60 days, after the introduction of 400 females of two important dust mite species, Blomia tropicalis van Bronswijk, de Cock e Oshima and Dermatophagoides pteronyssinus (Trouessart) were evaluate. With the most suitable foods, the population growth were higher than 20.2 and 15.3 for B. tropicalis and D. pteronyssinus, respectively.

Force-free twisted magnetospheres of neutron stars

Context. The X-ray spectra observed in the persistent emission of magnetars are evidence for the existence of a magnetosphere. The high-energy part of the spectra is explained by resonant cyclotron upscattering of soft thermal photons in a twisted magnetosphere, which has motivated an increasing number of efforts to improve and generalize existing magnetosphere models. Aims. We want to build more general configurations of twisted, force-free magnetospheres as a first step to understanding the role played by the magnetic field geometry in the observed spectra. Methods. First we reviewed and extended previous analytical works to assess the viability and limitations of semi-analytical approaches. Second, we built a numerical code able to relax an initial configuration of a nonrotating magnetosphere to a force-free geometry, provided any arbitrary form of the magnetic field at the star surface. The numerical code is based on a finite-difference time-domain, divergence-free, and conservative scheme, based of the magneto-frictional method used in other scenarios. Results. We obtain new numerical configurations of twisted magnetospheres, with distributions of twist and currents that differ from previous analytical solutions. The range of global twist of the new family of solutions is similar to the existing semi-analytical models (up to some radians), but the achieved geometry may be quite different. Conclusions. The geometry of twisted, force-free magnetospheres shows a wider variety of possibilities than previously considered. This has implications for the observed spectra and opens the possibility of implementing alternative models in simulations of radiative transfer aiming at providing spectra to be compared with observations.