Multiple copy sampling in protein loop modeling: Computational efficiency and sensitivity to dihedral angle perturbations

Multiple copy sampling and the bond scaling‐relaxation technique are combined to generate 3‐dimensional conformations of protein loop segments. The computational efficiency and sensitivity to initial loop copy dispersion are analyzed. The multicopy loop modeling method requires approximately 20‐50% of the computational time required by the single‐copy method for the various protein segments tested. An analytical formula is proposed to estimate the computational gain prior to carrying out a multicopy simulation. When 7‐residue loops within flexible proteins are modeled, each multicopy simulation can sample a set of loop conformations with initial dispersions up to ±15° for backbone and ±3° for side‐chain rotatable dihedral angles. The dispersions are larger for shorter and smaller for longer and/or surface loops. The degree of convergence of loop copies during a simulation can be used to complement commonly used target functions (such as potential energy) for distinguishing between native and misfolded conformations. Furthermore, this convergence also reflects the conformational flexibility of the modeled protein segment. Application to simultaneously building all 6 hypervariable loops of an antibody is discussed.

[1]  Jiří Novotný,et al.  Structure of antibody hypervariable loops reproduced by a conformational search algorithm , 1988, Nature.

[2]  Harold A. Scheraga,et al.  Standard‐geometry chains fitted to X‐ray derived structures: Validation of the rigid‐geometry approximation. I. Chain closure through a limited search of “loop” conformations , 1991 .

[3]  R. Elber,et al.  Computational studies of ligand diffusion in globins: I. Leghemoglobin , 1991, Proteins.

[4]  S. Nakanishi,et al.  Different binding epitopes on the NK1 receptor for substance P and a non-peptide antagonist , 1993, Nature.

[5]  M. Karplus,et al.  Enhanced sampling in molecular dynamics: use of the time-dependent Hartree approximation for a simulation of carbon monoxide diffusion through myoglobin , 1990 .

[6]  A. Lesk,et al.  Canonical structures for the hypervariable regions of immunoglobulins. , 1987, Journal of molecular biology.

[7]  John E. Straub,et al.  ENERGY EQUIPARTITIONING IN THE CLASSICAL TIME-DEPENDENT HARTREE APPROXIMATION , 1991 .

[8]  I. Kuntz,et al.  Protein docking and complementarity. , 1991, Journal of molecular biology.

[9]  R. Bruccoleri,et al.  Criteria that discriminate between native proteins and incorrectly folded models , 1988, Proteins.

[10]  C. DeLisi,et al.  Computing the structure of bound peptides. Application to antigen recognition by class I major histocompatibility complex receptors. , 1993, Journal of molecular biology.

[11]  Ron Elber,et al.  The thermal equilibrium aspects of the time dependent Hartree and the locally enhanced sampling approximations: Formal properties, a correction, and computational examples for rare gas clusters , 1993 .

[12]  R. Bruccoleri,et al.  On the attribution of binding energy in antigen-antibody complexes McPC 603, D1.3, and HyHEL-5. , 1989, Biochemistry.

[13]  M. Karplus,et al.  Multiple copy simultaneous search and construction of ligands in binding sites: application to inhibitors of HIV-1 aspartic proteinase. , 1993, Journal of medicinal chemistry.

[14]  C. Levinthal,et al.  Predicting antibody hypervariable loop conformation. I. Ensembles of random conformations for ringlike structures , 1987, Biopolymers.

[15]  M. Karplus,et al.  Functionality maps of binding sites: A multiple copy simultaneous search method , 1991, Proteins.

[16]  R. Elber,et al.  Modeling side chains in peptides and proteins: Application of the locally enhanced sampling and the simulated annealing methods to find minimum energy conformations , 1991 .

[17]  T. A. Jones,et al.  Using known substructures in protein model building and crystallography. , 1986, The EMBO journal.

[18]  J. Greer Comparative modeling methods: Application to the family of the mammalian serine proteases , 1990, Proteins.

[19]  M. Karplus,et al.  CHARMM: A program for macromolecular energy, minimization, and dynamics calculations , 1983 .

[20]  D J Kyle,et al.  A proposed model of bradykinin bound to the rat B2 receptor and its utility for drug design. , 1994, Journal of medicinal chemistry.

[21]  Robert Langridge,et al.  A distance geometry study of ring systems , 1983 .

[22]  S Vajda,et al.  Determining protein loop conformation using scaling‐relaxation techniques , 1993, Protein science : a publication of the Protein Society.

[23]  Rakefet Rosenfeld,et al.  Theoretical analysis of the multicopy sampling method in molecular modeling , 1993 .

[24]  Qiang Zheng,et al.  Loop closure via bond scaling and relaxation , 1993, J. Comput. Chem..

[25]  Gennady M Verkhivker,et al.  Microscopic modeling of ligand diffusion through the protein leghemoglobin : computer simulations and experiments , 1992 .

[26]  M Karplus,et al.  Modeling of globular proteins. A distance-based data search procedure for the construction of insertion/deletion regions and Pro----non-Pro mutations. , 1990, Journal of molecular biology.

[27]  Robert E. Bruccoleri,et al.  Chain closure with bond angle variations , 1985 .

[28]  N. Go,et al.  Ring Closure and Local Conformational Deformations of Chain Molecules , 1970 .