Biomolecular electrostatics with the linearized Poisson-Boltzmann equation

4openopenFOGOLARI F.; ZUCCATO P.; ESPOSITO G.; VIGLINO PFogolari, Federico; Zuccato, P.; Esposito, Gennaro; Viglino, Paol

Scoring predictive models using a reduced representation of proteins: model and energy definition

BACKGROUND: Reduced representations of proteins have been playing a keyrole in the study of protein folding. Many such models are available, with different representation detail. Although the usefulness of many such models for structural bioinformatics applications has been demonstrated in recent years, there are few intermediate resolution models endowed with an energy model capable, for instance, of detecting native or native-like structures among decoy sets. The aim of the present work is to provide a discrete empirical potential for a reduced protein model termed here PC2CA, because it employs a PseudoCovalent structure with only 2 Centers of interactions per Amino acid, suitable for protein model quality assessment. RESULTS: All protein structures in the set top500H have been converted in reduced form. The distribution of pseudobonds, pseudoangle, pseudodihedrals and distances between centers of interactions have been converted into potentials of mean force. A suitable reference distribution has been defined for non-bonded interactions which takes into account excluded volume effects and protein finite size. The correlation between adjacent main chain pseudodihedrals has been converted in an additional energetic term which is able to account for cooperative effects in secondary structure elements. Local energy surface exploration is performed in order to increase the robustness of the energy function. CONCLUSION: The model and the energy definition proposed have been tested on all the multiple decoys' sets in the Decoys'R'us database. The energetic model is able to recognize, for almost all sets, native-like structures (RMSD less than 2.0 Å). These results and those obtained in the blind CASP7 quality assessment experiment suggest that the model compares well with scoring potentials with finer granularity and could be useful for fast exploration of conformational space. Parameters are available at the url:

Molecular models for intrastrand DNA G-quadruplexes

<p>Abstract</p> <p>Background</p> <p>Independent surveys of human gene promoter regions have demonstrated an overrepresentation of G₃X_<it>n</it>1G3X_<it>n</it>2G₃X_<it>n</it>3G₃motifs which are known to be capable of forming intrastrand quadruple helix structures. In spite of the widely recognized importance of G-quadruplex structures in gene regulation and growing interest around this unusual DNA structure, there are at present only few such structures available in the Nucleic Acid Database. In the present work we generate by molecular modeling feasible G-quadruplex structures which may be useful for interpretation of experimental data.</p> <p>Results</p> <p>We have used all quadruplex DNA structures deposited in the Nucleic Acid Database in order to select a list of fragments entailing a strand of three adjacent G's paired with another strand of three adjacent G's separated by a loop of one to four residues. These fragments were further clustered and representative fragments were finally selected. Further fragments were generated by assemblying the two strands of each fragment with loops from different fragments whenever the anchor G's were superimposable. The fragments were used to assemble G quadruplex based on a superimposability criterion.</p> <p>Conclusion</p> <p>Molecular models have been generated for a large number of G₃X_<it>n</it>1G₃X_<it>n</it>2G3X_<it>n</it>3G₃sequences. For a given sequence not all topologies are possible with the available repertoire of fragments due to steric hindrance and low superimposability. Since all molecular models are generated by fragments coming from observed quadruplex structures, molecular models are in principle reliable and may be used for interpretation of experimental data. Some examples of applications are given.</p

Properties of Some Variants of Human β2-Microglobulin and Amyloidogenesis

Three variants of human beta(2)-microglobulin (beta(2)-m) were compared with wild-type protein. For two variants, namely the mutant R3Abeta(2)-m and the form devoid of the N-terminal tripeptide (DeltaN3beta(2)-m), a reduced unfolding free energy was measured compared with wild-type beta(2)-m, whereas an increased stability was observed for the mutant H31Ybeta(2)-m. The solution structure could be determined by (1)H NMR spectroscopy and restrained modeling only for R3Abeta(2)-m that showed the same conformation as the parent species, except for deviations at the interstrand loops. Analogous conclusions were reached for H31Ybeta(2)-m and DeltaN3beta(2)-m. Precipitation and unfolding were observed over time periods shorter than 4-6 weeks with all the variants and, sometimes, with wild-type protein. The rate of structured protein loss from solution as a result of precipitation and unfolding always showed pseudo-zeroth order kinetics. This and the failure to observe an unfolded species without precipitation suggest that a nucleated conformational conversion scheme should apply for beta(2)-m fibrillogenesis. The mechanism is consistent with the previous and present results on beta(2)-m amyloid transition, provided a nucleated oligomeric species be considered the stable intermediate of fibrillogenesis, the monomeric intermediate being the necessary transition step along the pathway from the native protein to the nucleated oligomer

Structural and Folding Dynamic Properties of the T70N Variant of Human Lysozyme

Definition of the transition mechanism from the native globular protein into fibrillar polymer was greatly improved by the biochemical and biophysical studies carried out on the two amyloidogenic variants of human lysozyme, I56T and D67H. Here we report thermodynamic and kinetic data on folding as well as structural features of a naturally occurring variant of human lysozyme, T70N, which is present in the British population at an allele frequency of 5% and, according to clinical and histopathological data, is not amyloidogenic. This variant is less stable than the wild-type protein by 3.7 kcal/mol, but more stable than the pathological, amyloidogenic variants. Unfolding kinetics in guanidine are six times faster than in the wild-type, but three and twenty times slower than in the amyloidogenic variants. Enzyme catalytic parameters, such as maximal velocity and affinity, are reduced in comparison to the wild-type. The solution structure, determined by 1H NMR and modeling calculations, exhibits a more compact arrangement at the interface between the beta-sheet domain and the subsequent loop on one side and part of the alpha domain on the other side, compared with the wild-type protein. This is the opposite of the conformational variation shown by the amyloidogenic variant D67H, but it accounts for the reduced stability and catalytic performance of T70N

Homology model building of the thyroid transcription factor 1 homeodomain

A possible structure for the homeodomain of rat thyroid transcription factor 1 is proposed on the basis of the homology with other homeodomains whose structures have been solved by X-ray crystallography. A structure very similar to the reference ones is feasible and may account for the observed DNA-binding specificity. Structural features of the model, which are likely to be shared by other homeodomains, are discussed and may help researchers in the field to analyse their own experimental data

Uptake and lifetime of fluoride ion in rats by 19F NMR

19F nuclear magnetic resonance (NMR) was utilized to obtain information on the uptake and half-life time of fluoride ion in rats. Changes in tissue fluoride level after acute loading were monitored over time in blood and tissue homogenates obtained from liver and brain. The rate of fluoride elimination from various tissues was roughly similar, following in all cases a first-order kinetic rate law. The F 12 concentration in brain was about 20% of that found in liver, indicating a reduced fluoride diffusion across the blood-brain barrier. In vivo F 12 spectra were obtained in rat brain in few minutes with a good signal-to-noise ratio; this confirms the possibility of extending the use of F 12 as a probe of biomolecules to in vivo applications