Computational method to reduce the search space for directed protein evolution

We introduce a computational method to optimize the in vitro evolution of proteins. Simulating evolution with a simple model that statistically describes the fitness landscape, we find that beneficial mutations tend to occur at amino acid positions that are tolerant to substitutions, in the limit of small libraries and low mutation rates. We transform this observation into a design strategy by applying mean-field theory to a structure-based computational model to calculate each residue's structural tolerance. Thermostabilizing and activity-increasing mutations accumulated during the experimental directed evolution of subtilisin E and T4 lysozyme are strongly directed to sites identified by using this computational approach. This method can be used to predict positions where mutations are likely to lead to improvement of specific protein properties.

A detergent-free method for purifying caveolae membrane from tissue culture cells.

Current methods for purifying caveolae from tissue culture cells take advantage of the Triton X-100 insolubility of this membrane domain. To circumvent the use of detergents, we have developed a method that depends upon the unique buoyant density of caveolae membrane. The caveolae fractions that we obtain are highly enriched in caveolin. As a consequence we are able to identify caveolae-associated proteins that had previously gone undetected. Moreover, resident caveolae proteins that are soluble in Triton X-100 are retained during the isolation.

Blood glutathione synthesis rates in healthy adults receiving a sulfur amino acid-free diet

The availability of cysteine is thought to be the rate limiting factor for synthesis of the tripeptide glutathione (GSH), based on studies in rodents. GSH status is compromised in various disease states and by certain medications leading to increased morbidity and poor survival. To determine the possible importance of dietary cyst(e)ine availability for whole blood glutathione synthesis in humans, we developed a convenient mass spectrometric method for measurement of the isotopic enrichment of intact GSH and then applied it in a controlled metabolic study. Seven healthy male subjects received during two separate 10-day periods an l-amino acid based diet supplying an adequate amino acid intake or a sulfur amino acid (SAA) (methionine and cysteine) free mixture (SAA-free). On day 10, l-[1-13C]cysteine was given as a primed, constant i.v. infusion (3μmol⋅kg−1⋅h−1) for 6 h, and incorporation of label into whole blood GSH determined by GC/MS selected ion monitoring. The fractional synthesis rate (mean ± SD; day-1) of whole blood GSH was 0.65 ± 0.13 for the adequate diet and 0.49 ± 0.13 for the SAA-free diet (P < 0.01). Whole blood GSH was 1,142 ± 243 and 1,216 ± 162 μM for the adequate and SAA-free periods (P > 0.05), and the absolute rate of GSH synthesis was 747 ± 216 and 579 ± 135 μmol⋅liter−1⋅day−1, respectively (P < 0.05). Thus, a restricted dietary supply of SAA slows the rate of whole blood GSH synthesis and diminishes turnover, with maintenance of the GSH concentration in healthy subjects.

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.