Biomedical subjects

Mihaly Mezei

Publications and source records attributed to Mihaly Mezei.

10 recordsLinked to original sources

PubMed2006

TRAJELIX: a computational tool for the geometric characterization of protein helices during molecular dynamics simulations.

We have developed a computer program with the necessary mathematical formalism for the geometric characterization of distorted conformations of alpha-helices proteins, such as those that can potentially be sampled during typical molecular dynamics simulations. This formalism has been incorporated into TRAJELIX, a new module within the SIMULAID framework (http://inka.mssm.edu/~mezei/simulaid/) that is capable of monitoring distortions of alpha-helices in terms of their displacement, global and local tilting, rotation around their axes, compression/extension, winding/unwinding, and bending. Accurate evaluation of these global and local structural properties of the helix can help study possible intramolecular and intermolecular changes in the helix packing of alpha-helical membrane proteins, as shown here in an application to the interacting helical domains of rhodopsin dimers. Quantification of the dynamic structural behavior of alpha-helical membrane proteins is critical for our understanding of signal transduction, and may enable structure-based design of more specific and efficient drugs.

Algorithms↗

PubMed2005

Origin of the sequence-dependent polyproline II structure in unfolded peptides.

Recent studies have indicated that the unfolded states of polypeptides contain a substantial amount of polyproline type II (P(II)) structure. This energetically and structurally preorganized state may contribute to the reduction of the folding search, as well as to the recognition of intrinsically unstructured proteins and polyproline ligands. Using Monte Carlo simulations of natively unfolded peptides in the presence of explicit aqueous solvation, we observe that residue-specific P(II) conformational propensity is the result of the modulation of polypeptide backbone hydration by a proximal side-chain. Such a mechanism may be unique among those that contribute to the modulation of secondary structures in proteins. The calculated conformational propensities should prove useful for the development of a configurational P(II) scale necessary for the prediction and design of natural-like polypeptides.

Kinetics↗

PubMed2005

Effect of cholesterol on the properties of phospholipid membranes. 4. Interatomic voids.

The properties of the interatomic voids present in fully hydrated dimyristoylphosphatidylcholine (DMPC)-cholesterol mixed membranes of different compositions are analyzed in detail using a generalized variant of the Voronoi-Delaunay method on the basis of computer simulation results. The systems investigated are chosen from both sides of the DMPC-cholesterol miscibility gap; the pure DMPC bilayer has also been included in the analysis as a reference system. The results obtained show that the empty space is organized in a more compact way, forming larger voids in the presence than in the absence of cholesterol. The voids located in the region of the rigid cholesterol rings become, on average, less spherical, oriented more parallel with the membrane normal axis with increasing cholesterol concentration, whereas an opposite effect of cholesterol is observed in the middle of the membrane among the chain terminal methyl groups. In general, the preferential orientation of the voids is found to strongly correlate with that of the molecules in the hydrocarbon phase of the membranes. The membranes are found to contain rather large voids, the volume of which can be an order of magnitude larger than the largest spherical cavities present in the systems. These voids are elongated or branching channels rather than big empty holes. The voids located among the DMPC and cholesterol molecules are lying preferably parallel with the membrane normal axis. The existence of such empty channels can be of great importance in the cross-membrane permeation of small, uncharged penetrants, in particular, of polar molecules.

Cholesterol↗

PubMed2005

Interfacial water as a "hydration fingerprint" in the noncognate complex of BamHI.

The molecular code of specific DNA recognition by proteins as a paradigm in molecular biology remains an unsolved puzzle primarily because of the subtle interplay between direct protein-DNA interaction and the indirect contribution from water and ions. Transformation of the nonspecific, low affinity complex to a specific, high affinity complex is accompanied by the release of interfacial water molecules. To provide insight into the conversion from the loose to the tight form, we characterized the structure and energetics of water at the protein-DNA interface of the BamHI complex with a noncognate sequence and in the specific complex. The fully hydrated models were produced with Grand Canonical Monte Carlo simulations. Proximity analysis shows that water distributions exhibit sequence dependent variations in both complexes and, in particular, in the noncognate complex they discriminate between the correct and the star site. Variations in water distributions control the number of water molecules released from a given sequence upon transformation from the loose to the tight complex as well as the local entropy contribution to the binding free energy. We propose that interfacial waters can serve as a "hydration fingerprint" of a given DNA sequence.

Base Sequence↗

PubMed2004

A novel method reveals that solvent water favors polyproline II over beta-strand conformation in peptides and unfolded proteins: conditional hydrophobic accessible surface area (CHASA).

In aqueous solution, the ensemble of conformations sampled by peptides and unfolded proteins is largely determined by their interaction with water. It has been a long-standing goal to capture these solute-water energetics accurately and efficiently in calculations. Historically, accessible surface area (ASA) has been used to estimate these energies, but this method breaks down when applied to amphipathic peptides and proteins. Here we introduce a novel method in which hydrophobic ASA is determined after first positioning water oxygens in hydrogen-bonded orientations proximate to all accessible peptide/protein backbone N and O atoms. This conditional hydrophobic accessible surface area is termed CHASA. The CHASA method was validated by predicting the polyproline-II (P(II)) and beta-strand conformational preferences of non-proline residues in the coil library (i.e., non-alpha-helix, non-beta-strand, non-beta-turn library derived from X-ray elucidated structures). Further, the method successfully rationalizes the previously unexplained solvation energies in polyalanyl peptides and compares favorably with published experimentally determined P(II) residue propensities. We dedicate this paper to Frederic M. Richards.

Hydrogen Bonding↗

PubMed2004

Polyproline II helix is the preferred conformation for unfolded polyalanine in water.

Does aqueous solvent discriminate among peptide conformers? To address this question, we computed the solvation free energy of a blocked, 12-residue polyalanyl-peptide in explicit water and analyzed its solvent structure. The peptide was modeled in each of 4 conformers: alpha-helix, antiparallel beta-strand, parallel beta-strand, and polyproline II helix (P(II)). Monte Carlo simulations in the canonical ensemble were performed at 300 K using the CHARMM 22 forcefield with TIP3P water. The simulations indicate that the solvation free energy of P(II) is favored over that of other conformers for reasons that defy conventional explanation. Specifically, in these 4 conformers, an almost perfect correlation is found between a residue's solvent-accessible surface area and the volume of its first solvent shell, but neither quantity is correlated with the observed differences in solvation free energy. Instead, solvation free energy tracks with the interaction energy between the peptide and its first-shell water. An additional, previously unrecognized contribution involves the conformation-dependent perturbation of first-shell solvent organization. Unlike P(II), beta-strands induce formation of entropically disfavored peptide:water bridges that order vicinal water in a manner reminiscent of the hydrophobic effect. The use of explicit water allows us to capture and characterize these dynamic water bridges that form and dissolve during our simulations.

Computer Simulation↗

PubMed2004

Unfolded state of polyalanine is a segmented polyproline II helix.

Definition of the unfolded state of proteins is essential for understanding their stability and folding on biological timescales. Here, we find that under near physiological conditions the configurational ensemble of the unfolded state of the simplest protein structure, polyalanine alpha-helix, cannot be described by the commonly used Flory random coil model, in which configurational probabilities are derived from conformational preferences of individual residues. We utilize novel effectively ergodic sampling algorithms in the presence of explicit aqueous solvation, and observe water-mediated formation of polyproline II helical (P(II)) structure in the natively unfolded state of polyalanine, and its facilitation of alpha-helix formation in longer peptides. The segmented P(II) helical coil preorganizes the unfolded state ensemble for folding pathway entry by reducing the conformational space available to the diffusive search. Thus, as much as half of the folding search in polyalanine is accomplished by preorganization of the unfolded state.

Computer Simulation↗

PubMed2003

MC-PHS: a Monte Carlo implementation of the primary hydration shell for protein folding and design.

A primary hydration shell (PHS) approach is developed for Monte Carlo simulations of conformationally rich macromolecular systems in an environment that efficiently captures principal solvation effects. It has been previously demonstrated that molecular dynamics using PHS is an efficient method to study peptide structure and dynamics in aqueous solution. Here, we extend the PHS approach to Monte Carlo simulations, whereby a stable shell of water molecules is maintained with a flexible, nonspherical, half-harmonic potential, tuned to maintain a constant restraining energy, with the difference between the restraint and shell energies used to dynamically adjust the shell radius. Examination of the shell and system size dependence of the restraining potential reveals its robustness. Moreover, its suitability for biomolecular simulations is evaluated using small spheres of water, hydration properties of small biological molecules, and configurational sampling of beta-hairpin pentapeptide YPGDV. This method, termed MC-PHS, appears to provide efficient representation of dominant solvation effects and should prove useful in the study of protein folding and design.

Amino Acids↗

PubMed2003

A new method for mapping macromolecular topography.

A new method, using circular variance, is introduced for mapping macromolecular topography. Circular variance, generally used to measures angular spread, can be used to characterize of molecular structures based on a simple idea. It will be shown that the circular variance of vectors drawn from some origin to a set of points is well correlated with the degree to which that origin is inside/outside the chosen points. In addition, it has continuous derivatives that are also easy to compute. This concept will be shown to be useful for: (i) distinguishing between atoms near the surface of a macromolecule and those in either the deep interior or remote exterior; (ii) identifying invaginations (even shallow ones); and (iii) detecting linker regions that interconnect two domains.

Algorithms↗

PubMed2003

A novel fingerprint for the characterization of protein folds.

A novel fingerprint, defined without the use of distances, is introduced to characterize protein folds. It is of the form of binary matrices whose elements are defined by angles between the C=O direction, the backbone axis and the line connecting the alpha-carbons of the various residues. It is shown that matches in the fingerprint matrices correspond to low r.m.s.d.

Peptide Mapping↗