Search PubMed⌕ Search

Biomedical subjects

G Vriend

Publications and source records attributed to G Vriend.

At least 55 records · Page 3Linked to original sources

Extreme stabilization of a thermolysin-like protease by an engineered disulfide bond.

The thermal inactivation of broad specificity proteases such as thermolysin and subtilisin is initiated by partial unfolding processes that render the enzyme susceptible to autolysis. Previous studies have revealed that a surface-located region in the N-terminal domain of the thermolysin-like protease produced by Bacillus stearothermophilus is crucial for thermal stability. In this region a disulfide bridge between residues 8 and 60 was designed by molecular modelling, and the corresponding single and double cysteine mutants were constructed. The disulfide bridge was spontaneously formed in vivo and resulted in a drastic stabilization of the enzyme. This stabilization presents one of the very few examples of successful stabilization of a broad specificity protease by a designed disulfide bond. We propose that the success of the present stabilization strategy is the result of the localization and mutation of an area of the molecule involved in the partial unfolding processes that determine thermal stability.

Disulfides↗

CASP2 molecular docking predictions with the LIGIN software.

Seven docking predictions were made with the LIGIN program. In six cases the location of the binding pocket was identified correctly by systematically docking everywhere within the protein structure. In two cases the ligand was docked to within 1.8 A RMSD of the experimentally determined structure. LIGIN has not been optimized to deal with highly flexible ligands that dock at the surface of proteins. Consequently, in three cases the exposed part of the ligand was docked poorly, although the buried parts were docked well, and made similar atomic contacts with the protein as in the experimentally determined structure.

Amiloride↗

Prediction of protein conformational freedom from distance constraints.

A method is presented that generates random protein structures that fulfil a set of upper and lower interatomic distance limits. These limits depend on distances measured in experimental structures and the strength of the interatomic interaction. Structural differences between generated structures are similar to those obtained from experiment and from MD simulation. Although detailed aspects of dynamical mechanisms are not covered and the extent of variations are only estimated in a relative sense, applications to an IgG-binding domain, an SH3 binding domain, HPr, calmodulin, and lysozyme are presented which illustrate the use of the method as a fast and simple way to predict structural variability in proteins. The method may be used to support the design of mutants, when structural fluctuations for a large number of mutants are to be screened. The results suggest that motional freedom in proteins is ruled largely by a set of simple geometric constraints.

Magnetic Resonance Spectroscopy↗

ACGT and vicilin core sequences in a promoter domain required for seed-specific expression of a 2S storage protein gene are recognized by the opaque-2 regulatory protein.

The expression of Brazil nut storage albumin genes is highly regulated during seed development. Several sequences in the promoter of one of these genes show homologies with the target sites of the maize O2 bZIP regulatory protein. We therefore asked whether the O2 protein would recognize these promoter sequences. We show that the O2 protein binds to three different sequences (F1, F2 and F3). F1 and F3 are hybrid C/G and A/G boxes, respectively, that are homologous to the O2-binding site of a maize alpha-zein gene. F2 is a new O2-binding sequence related to the O2 target sites of the Coix alpha-coxin, the maize b-32 genes and the AP-1 pseudopalindrome. Molecular modelling showed that an Asn and a Ser in the 02 DNA binding domain make different base-specific contacts with each operator. 5' Promoter deletions of the be2S1 gene showed that the domain containing the O2 target sites F1 and F2 is required for detectable reporter gene expression in transgenic tobacco seeds. Moreover, the homologous coix O2 protein was shown to in situ transactivate the promoter region encompassing the three O2-binding sites F1, F2 and F3. Thus, these sites may be in vivo regulatory sequences mediating activation by bZIP regulatory proteins.

2S Albumins, Plant↗

Objectively judging the quality of a protein structure from a Ramachandran plot.

MOTIVATION: Statistical methods that compare observed and expected distributions of experimental observables provide powerful tools for the quality control of protein structures. The distribution of backbone dihedral angles ('Ramachandran plot') has often been used for such quality control, but without a firm statistical foundation. RESULTS: A new and-simple method is presented for judging the quality of a protein structure based on the distribution of backbone dihedral angles. Inputs to the method are 60 torsion angle distributions extracted from protein structures solved at high resolution; one for each combination of residue type and tri-state secondary structure. Output for a protein is a Ramachandran Z-score, expressing the quality of the Ramachandran plot relative to current state-of-the-art structures.

Algorithms↗

Model building of a thermolysin-like protease by mutagenesis.

The present study concerns the use of site-directed mutagenesis experiments to optimize a three-dimensional model of the neutral protease of Bacillus subtilis (NP-sub). An initial model of NP-sub was constructed using the crystal structures of the homologous neutral proteases of Bacillus thermoproteolyticus (thermolysin) and Bacillus cereus as templates. The largest portion of NP-sub could be modelled satisfactorily, using standard techniques, but several surface-located regions could only be modelled with a high degree of uncertainty. In order to make the model more reliable in these regions a 'model building by mutagenesis' approach was adopted. Mutations were designed such that their effect on thermal stability could indicate how their local environment should be modelled. This approach provided insight in the local structure of several regions in NP-sub that were hard to model on the basis of homology with the two known structures alone.

Amino Acid Sequence↗

Identification of class-determining residues in G protein-coupled receptors by sequence analysis.

G protein-coupled receptors (GPCRs) form a large superfamily of receptors that are characterised by a seven transmembrane helical motif. The functions they perform, such as binding ligands and G proteins, are related to the presence of certain amino acids in critical positions. We have developed a computational sequence pattern correlation technique for the recognition of such function-determining residues. The method searches for residues that are conserved in one class of proteins with a certain function but are different in other classes. The basic idea is that such residues are probably involved in this particular function. This technique was used to find residues that play a role in the binding of endogenous as well as exogenous ligands to various receptors. Many of the residues that were detected have been experimentally determined as important for ligand binding. More importantly, however, we also detected residues that are interesting targets for future mutation studies aimed at elucidating the sequence-function relationship in GPCRs. The information obtained may help improve three-dimensional GPCR models and can be useful for the study of receptor-ligand interactions.

Amino Acid Sequence↗

Molecular docking using surface complementarity.

A method is described to dock a ligand into a binding site in a protein on the basis of the complementarity of the intermolecular atomic contacts. Docking is performed by maximization of a complementarity function that is dependent on atomic contact surface area and the chemical properties of the contacting atoms. The generality and simplicity of the complementarity function ensure that a wide range of chemical structures can be handled. The ligand and the protein are treated as rigid bodies, but displacement of a small number of residues lining the ligand binding site can be taken into account. The method can assist in the design of improved ligands by indicating what changes in complementarity may occur as a result of the substitution of an atom in the ligand. The capabilities of the method are demonstrated by application to 14 protein-ligand complexes of known crystal structure.

Algorithms↗

Positioning hydrogen atoms by optimizing hydrogen-bond networks in protein structures.

A method is presented that positions polar hydrogen atoms in protein structures by optimizing the total hydrogen bond energy. For this goal, an empirical hydrogen bond force field was derived from small molecule crystal structures. Bifurcated hydrogen bonds are taken into account. The procedure also predicts ionization states of His, Asp, and Glu residues. During optimization, side-chain conformations of His, Gln, and Asn residues are allowed to change their last chi angle by 180 degrees to compensate for crystallographic misassignments. Crystal structure symmetry is taken into account where appropriate. The results can have significant implications for molecular dynamics simulations, protein engineering, and docking studies. The largest impact, however, is in protein structure verification: over 85% of protein structures tested can be improved by using our procedure.

Amino Acids↗

PRODRG, a program for generating molecular topologies and unique molecular descriptors from coordinates of small molecules.

A software package is described that operates on small molecules observed in the PDB collection of protein structures. Molecular topology files for many molecular modeling programs can be generated automatically. The three-dimensional coordinates of small molecules can be converted to molecular descriptor strings that encode them uniquely in order to enable small-molecule recognition, despite high variability in atom and molecule nomenclature. From this descriptor a plausible 3D structure can be regenerated using energy minimisation. Alternatively, an ensemble of structures can be generated using a distance-geometry-based algorithm.

Ligands↗

Improved coordinate reconstruction from stereo diagrams.

A program has been written that reconstructs three-dimensional coordinates for a protein structure given a stereo C alpha diagram. Initial three-dimensional coordinates are determined using an algorithm similar to the one used by Rossmann in the program STEREO. Thereafter the coordinates are refined such that the stereo image based on the reconstructed three-dimensional coordinates optimally fits the scanned stereo image while normal C alpha stereochemistry is enforced.

Computer Communication Networks↗

The PDBFINDER database: a summary of PDB, DSSP and HSSP information with added value.

MOTIVATION: The Protein Data Bank currently contains more than 4700 protein coordinate sets. It is often desirable to make a selection from these files based on a criterion like R-factor, experimental method, length of the amino acid sequence, or the number of homologous sequences in SWISSPROT. Doing this using the distributed form of the Protein Data Bank can be a tedious task, because (1) this requires reading one file for every single entry, and (2) not all of the information is present in a consistent computer readable way in all of the entries. RESULTS: The PDBFINDER database provides an easy to interpret file containing summary information about all Protein Data Bank files. Summary information from the DSSP (Definition of Secondary Structure of Proteins) and HSSP (Homology derived Secondary Structure of Proteins) databases is also included. Furthermore, where essential data were missing from the Protein Data Bank file, this information has been retrieved from the original literature. AVAILABILITY: The latest version of the PDBFINDER database can be downloaded by anonymous ftp from swift.embl-heidelberg.de, directory:/pdbfinder. CONTACT: E-mail address hooft@embl-heidelberg.de.

Amino Acid Sequence↗

Analysis of structural determinants of the stability of thermolysin-like proteases by molecular modelling and site-directed mutagenesis.

The thermolysin-like protease (TLP) produced by Bacillus stearothermophilus CU21 (TLP-ste) differs at 43 positions from the more thermally stable thermolysin (containing 316 residues in total). Of these differences, 26 were analysed by studying the effect of replacing residues in TLP-ste by the corresponding residues in thermolysin. Several stabilizing mutations were identified but, remarkably, considerable destabilizing mutational effects were also found. A Tyr-rich three residue insertion in TLP-ste (the only deletional/insertional difference between the two enzymes) appeared to make an important contribution to the stability of the enzyme. Mutations with large effects on stability were all localized in the beta-pleated N-terminal domain of TLP-ste, confirming observations that this domain has a lower intrinsic stability than the largely alpha-helical C-terminal domain. Rigidifying mutations such as Gly58-->Ala and Ala69-->Pro were among the most stabilizing ones. Apart from this observation, the analyses did not reveal general rules for stabilizing proteins. Instead, the results highlight the importance of context in evaluating the stability effects of mutations.

Amino Acid Sequence↗

The purification, characterization and analysis of primary and secondary-structure of prolyl oligopeptidase from human lymphocytes. Evidence that the enzyme belongs to the alpha/beta hydrolase fold family.

Prolyl oligopeptidase was isolated and purified to homogeneity from human lymphocytes, yielding a specific activity of 7780 mU/mg. The molecular mass using size-exclusion chromatography matches the 76 kDa obtained by SDS/PAGE. This provides evidence that prolyl oligopeptidase is a monomer. The isoelectric point is 4.8 as judged by isoelectric focusing in free solution. Di-isopropyl fluorophosphate and phenylmethylsulphonyl fluoride completely abolish the activity, classifying the enzyme as a serine proteinase. The inhibition by p-chloromercuribenzoic acid indicates the importance of a free sulfhydryl group near the active-site. alpha 1-Casein and ornithine decarboxylase, two proteins containing a PEST sequence, inhibit prolyl oligopeptidase, but were not hydrolyzed. This demonstrates that prolyl oligopeptidase is not participating in the metabolism of proteins according to a PEST-dependent pathway. alpha 1-Antitrypsin partially inhibits the enzyme but in contrast, aprotinin does not. Its inability to cleave corticotropin-releasing factor, ubiquitin, albumin and aprotinin, together with the hydrolysis of bradykinin between Pro7-Arg8 confirms the affinity of prolyl oligopeptidase for small peptides. Multiple sequence alignment does not reveal any similarity with proteases of known tertiary structure. Secondary-structure prediction displays striking similarity with dipeptidyl peptidase IV and acylaminoacyl peptidase. Two characteristic features of the members of the prolyl oligopeptidase family of serine proteases are high-lighted: the linear arrangement of the catalytic triad is nucleophile-acid-base and the proteolytic cleavage releasing the catalytically active C-terminal region of around 500 amino acids from the N-terminal sequence. Secondary structure prediction and comparison of the active-site of serine proteinases with known three-dimensional coordinates prove that Asp641 is the third member of the catalytic triad. The secondary structural organization of the protease domain of prolyl oligopeptidase is in accordance with the alpha/beta hydrolase fold.

Amino Acid Sequence↗

The essential dynamics of thermolysin: confirmation of the hinge-bending motion and comparison of simulations in vacuum and water.

Comparisons of the crystal structures of thermolysin and the thermolysin-like protease produced by B. cereus have recently led to the hypothesis that neutral proteases undergo a hinge-bending motion. We have investigated this hypothesis by analyzing molecular dynamics simulations of thermolysin in vacuum and water, using the essential dynamics method. This method is able to extract large concerted atomic motions of biological importance from a molecular dynamics trajectory. The analysis of the thermolysin trajectories indeed revealed a large rigid body hinge-bending motion of the N-terminal and C-terminal domains, similar to the motion hypothesized from the crystal structure comparisons. In addition, it appeared that the essential dynamics properties derived from the vacuum simulation were similar to those derived from the solvent simulation.

Computer Simulation↗

The use of position-specific rotamers in model building by homology.

In this study we concentrate on replacing side chains as a subtask of model building by homology. Two problems arise. How to determine potential low energy rotamers? And how to avoid the combinatorial explosion that results from the combination of many residues for which multiple good rotamers are predicted? We attempt to solve these problems by choosing position-specific rather than generalized rotamers and by sorting the residues that have to be modelled as a function of their freedom in rotamer space. The practical advantages of our method are the quality of the models for cases of high backbone similarity, the small amount of human intervention needed, and the fact that the method automatically estimates the reliability with which each residue has been modeled. Other methods described in this issue are probably more suitable if large backbone rearrangements or loop insertions and deletions need to be modeled.

Bacterial Proteins↗