Biomedical subjects

Andrei L Lomize

Publications and source records attributed to Andrei L Lomize.

6 recordsLinked to original sources

PubMed2006

OPM: orientations of proteins in membranes database.

SUMMARY: The Orientations of Proteins in Membranes (OPM) database provides a collection of transmembrane, monotopic and peripheral proteins from the Protein Data Bank whose spatial arrangements in the lipid bilayer have been calculated theoretically and compared with experimental data. The database allows analysis, sorting and searching of membrane proteins based on their structural classification, species, destination membrane, numbers of transmembrane segments and subunits, numbers of secondary structures and the calculated hydrophobic thickness or tilt angle with respect to the bilayer normal. All coordinate files with the calculated membrane boundaries are available for downloading. AVAILABILITY: http://opm.phar.umich.edu.

Computer Graphics↗

PubMed2006

Positioning of proteins in membranes: a computational approach.

A new computational approach has been developed to determine the spatial arrangement of proteins in membranes by minimizing their transfer energies from water to the lipid bilayer. The membrane hydrocarbon core was approximated as a planar slab of adjustable thickness with decadiene-like interior and interfacial polarity profiles derived from published EPR studies. Applicability and accuracy of the method was verified for a set of 24 transmembrane proteins whose orientations in membranes have been studied by spin-labeling, chemical modification, fluorescence, ATR FTIR, NMR, cryo-microscopy, and neutron diffraction. Subsequently, the optimal rotational and translational positions were calculated for 109 transmembrane, five integral monotopic and 27 peripheral protein complexes with known 3D structures. This method can reliably distinguish transmembrane and integral monotopic proteins from water-soluble proteins based on their transfer energies and membrane penetration depths. The accuracies of calculated hydrophobic thicknesses and tilt angles were approximately 1 A and 2 degrees, respectively, judging from their deviations in different crystal forms of the same proteins. The hydrophobic thicknesses of transmembrane proteins ranged from 21.1 to 43.8 A depending on the type of biological membrane, while their tilt angles with respect to the bilayer normal varied from zero in symmetric complexes to 26 degrees in asymmetric structures. Calculated hydrophobic boundaries of proteins are located approximately 5 A lower than lipid phosphates and correspond to the zero membrane depth parameter of spin-labeled residues. Coordinates of all studied proteins with their membrane boundaries can be found in the Orientations of Proteins in Membranes (OPM) database:http://opm.phar.umich.edu/.

Amino Acids↗

PubMed2005

Interactions of human melanocortin 4 receptor with nonpeptide and peptide agonists.

Specific interactions of human melanocortin-4 receptor (hMC4R) with its nonpeptide and peptide agonists were studied using alanine-scanning mutagenesis. The binding affinities and potencies of two synthetic, small-molecule agonists (THIQ, MB243) were strongly affected by substitutions in transmembrane alpha-helices (TM) 2, 3, 6, and 7 (residues Glu(100), Asp(122), Asp(126), Phe(261), His(264), Leu(265), and Leu(288)). In addition, a I129A mutation primarily affected the binding and potency of THIQ, while F262A, W258A, Y268A mutations impaired interactions with MB243. By contrast, binding affinity and potency of the linear peptide agonist NDP-MSH were substantially reduced only in D126A and H264A mutants. Three-dimensional models of receptor-ligand complexes with their agonists were generated by distance-geometry using the experimental, homology-based, and other structural constraints, including interhelical H-bonds and two disulfide bridges (Cys(40)-Cys(279), Cys(271)-Cys(277)) of hMC4R. In the models, all pharmacophore elements of small-molecule agonists are spatially overlapped with the corresponding key residues (His(6), d-Phe(7), Arg(8), and Trp(9)) of the linear peptide: their charged amine groups interact with acidic residues from TM2 and TM3, similar to His(6) and Arg(6) of NDP-MSH; their substituted piperidines mimic Trp(9) of the peptide and interact with TM5 and TM6, while the d-Phe aromatic rings of all three agonists contact with Leu(133), Trp(258), and Phe(261) residues.

Amino Acid Sequence↗

PubMed2004

Complex of an active mu-opioid receptor with a cyclic peptide agonist modeled from experimental constraints.

Site-directed mutagenesis and design of Zn(2+)-binding centers have been used to determine a set of specific tertiary interactions between the mu-opioid receptor, a rhodopsin-like G protein-coupled receptor (GPCR), and its cyclic peptide agonist ligand, Tyr(1)-c(S-Et-S)[d-Cys(2)-Phe(3)-d-Pen(4)]NH(2) (JOM6). The binding affinity of the tetrapeptide is strongly dependent on the nature of its first and third residues and on substitutions at positions 213, 216, 237, 300, 315, and 318 of the mu-opioid receptor. His(1) and His(3) analogues of the ligand were able to form metal-binding complexes with the V300C and G213C/T315C receptor mutants, respectively. Direct contact of the Phe(3) residue of JOM6 with Gly(213), Asp(216), Thr(315), and Trp(318) of the receptor was suggested by the binding affinities of His(3)-, Nle(3)-, Leu(3)-, Aci(3)-, Delta(E)Phe(3)-, and Delta(Z)Phe(3)-substituted peptides with the G213C/T315C, D216V, T315C, and W318L mutants. The improved binding affinity of the free carboxylate analogue of JOM6 for binding to the E229D mutant revealed an interaction between the C-terminal group of the peptide and Glu(229) of the receptor. The experimental constraints that were obtained were applied for distance geometry modeling of the mu-receptor in complex with the tetrapeptide agonist ligand, JOM6. The active conformation of the opioid receptor was calculated using the crystal structure of "inactive" rhodopsin and published engineered and intrinsic metal-binding sites and disulfide bonds that allow or facilitate activation of GPCRs. Interhelical H-bonds existing in the mu-receptor were applied as additional distance constraints. The calculated model of the receptor-ligand complex can serve as a prototype of the active state for all rhodopsin-like GPCRs. It displays a strongly shifted transmembrane helix 6 (TM6) and reorientation of the conserved Trp(293) residue in TM6 upon its interaction with the agonist. Importantly, the binding pockets of the active and inactive states are not identical, which implies distinct interaction modes of agonists and antagonists. In the active state, the binding pocket of the mu-receptor is complementary to the previously proposed receptor-bound conformation of JOM6.

Amino Acid Sequence↗

PubMed2004

Quantification of helix-helix binding affinities in micelles and lipid bilayers.

A theoretical approach for estimating association free energies of alpha-helices in nonpolar media has been developed. The parameters of energy functions have been derived from DeltaDeltaG values of mutants in water-soluble proteins and partitioning of organic solutes between water and nonpolar solvents. The proposed approach was verified successfully against three sets of published data: (1) dissociation constants of alpha-helical oligomers formed by 27 hydrophobic peptides; (2) stabilities of 22 bacteriorhodopsin mutants, and (3) protein-ligand binding affinities in aqueous solution. It has been found that coalescence of helices is driven exclusively by van der Waals interactions and H-bonds, whereas the principal destabilizing contributions are represented by side-chain conformational entropy and transfer energy of atoms from a detergent or lipid to the protein interior. Electrostatic interactions of alpha-helices were relatively weak but important for reproducing the experimental data. Immobilization free energy, which originates from restricting rotational and translational rigid-body movements of molecules during their association, was found to be less than 1 kcal/mole. The energetics of amino acid substitutions in bacteriorhodopsin was complicated by specific binding of lipid and water molecules to cavities created in certain mutants.

Amino Acid Substitution↗

PubMed2002

Interatomic potentials and solvation parameters from protein engineering data for buried residues.

Van der Waals (vdW) interaction energies between different atom types, energies of hydrogen bonds (H-bonds), and atomic solvation parameters (ASPs) have been derived from the published thermodynamic stabilities of 106 mutants with available crystal structures by use of an originally designed model for the calculation of free-energy differences. The set of mutants included substitutions of uncharged, inflexible, water-inaccessible residues in alpha-helices and beta-sheets of T4, human, and hen lysozymes and HI ribonuclease. The determined energies of vdW interactions and H-bonds were smaller than in molecular mechanics and followed the "like dissolves like" rule, as expected in condensed media but not in vacuum. The depths of modified Lennard-Jones potentials were -0.34, -0.12, and -0.06 kcal/mole for similar atom types (polar-polar, aromatic-aromatic, and aliphatic-aliphatic interactions, respectively) and -0.10, -0.08, -0.06, -0.02, and nearly 0 kcal/mole for different types (sulfur-polar, sulfur-aromatic, sulfur-aliphatic, aliphatic-aromatic, and carbon-polar, respectively), whereas the depths of H-bond potentials were -1.5 to -1.8 kcal/mole. The obtained solvation parameters, that is, transfer energies from water to the protein interior, were 19, 7, -1, -21, and -66 cal/moleA(2) for aliphatic carbon, aromatic carbon, sulfur, nitrogen, and oxygen, respectively, which is close to the cyclohexane scale for aliphatic and aromatic groups but intermediate between octanol and cyclohexane for others. An analysis of additional replacements at the water-protein interface indicates that vdW interactions between protein atoms are reduced when they occur across water.

Amino Acid Substitution↗