General trends of dihedral conformational transitions in a globular protein

Dihedral conformational transitions are analyzed systematically in a model globular protein, cytochrome P450cam, to examine their structural and chemical dependences through combined conventional molecular dynamics (cMD), accelerated molecular dynamics (aMD) and adaptive biasing force (ABF) simulations. The aMD simulations are performed at two acceleration levels, using dihedral and dual boost, respectively. In comparison with cMD, aMD samples protein dihedral transitions approximately two times faster on average using dihedral boost, and ～3.5 times faster using dual boost. In the protein backbone, significantly higher dihedral transition rates are observed in the bend, coil, and turn flexible regions, followed by the β bridge and β sheet, and then the helices. Moreover, protein side chains of greater length exhibit higher transition rates on average in the aMD‐enhanced sampling. Side chains of the same length (particularly N_χ = 2) exhibit decreasing transition rates with residues when going from hydrophobic to polar, then charged and aromatic chemical types. The reduction of dihedral transition rates is found to be correlated with increasing energy barriers as identified through ABF free energy calculations. These general trends of dihedral conformational transitions provide important insights into the hierarchical dynamics and complex free energy landscapes of functional proteins. Proteins 2016; 84:501–514. © 2016 Wiley Periodicals, Inc.     

Unconstrained enhanced sampling for free energy calculations of biomolecules: a review

Abstract

Free energy calculations are central to understanding the structure, dynamics and function of biomolecules. Yet insufficient sampling of biomolecular configurations is often regarded as one of the main sources of error. Many enhanced sampling techniques have been developed to address this issue. Notably, enhanced sampling methods based on biasing collective variables (CVs), including the widely used umbrella sampling, adaptive biasing force and metadynamics, have been discussed in a recent excellent review (Abrams and Bussi, Entropy, 2014). Here, we aim to review enhanced sampling methods that do not require predefined system-dependent CVs for biomolecular simulations and as such do not suffer from the hidden energy barrier problem as encountered in the CV-biasing methods. These methods include, but are not limited to, replica exchange/parallel tempering, self-guided molecular/Langevin dynamics, essential energy space random walk and accelerated molecular dynamics. While it is overwhelming to describe all details of each method, we provide a summary of the methods along with the applications and offer our perspectives. We conclude with challenges and prospects of the unconstrained enhanced sampling methods for accurate biomolecular free energy calculations.     

Superoxide dismutase: fluctuations in the structure and solvation of the active site channel studied by molecular dynamics simulation

The molecular dynamics (MD) simulation of superoxide dismutase (SOD) in water is carried out for a total of 23 ps. The simulation system is a 26 A sphere centered at the active site of SOD, including 1602 atoms from SOD and 1761 water molecules. There is no gross deviation from the x-ray structure for the average MD structure. The structure and potential fluctuations around the active site are examined. The results provide new insight to the interactions between SOD and its substrate superoxide.     

Helix–coil transitions in a simple polypeptide model

A simplified model of a polypeptide chain is described. Each residue is represented by a single interaction center. The energy of the chain and the force acting on each residue are given as a function of the residue coordinates. Terms to approximate the effect of solvent and the stabilization energy of helix formation are included. The model is used to study equilibrium and dynamical aspects of the helix–coil transition. The equilibrium properties examined include helix–coil equilibrium constants and their dependence on chain position. Dynamical properties are examined by a stochastic simulation of the Brownian motion of the chain in its solvent surroundings. Correlations in the motions of the residues are found to have an important influence on the helix–coil transition rates.     

Human alveolar macrophage fibronectin: synthesis, secretion, and ultrastructural localization during gelatin-coated latex particle binding

Human pulmonary alveolar macrophages synthesized and secreted several characteristic high molecular weight proteins for at least 7 d in vitro. Immunoprecipitates of medium and cell lysates from metabolically labeled cultures with specific anti-human plasma fibronectin IgG contained one major labeled polypeptide of molecular weight 440,000 (unreduced) or 220,000 (reduced). An identical polypeptide in conditioned medium from radiolabeled macrophages bound specifically to gelatin-Sepharose, demonstrating that alveolar macrophages synthesized and secreted a molecule immunologically and functionally similar to fibronectin. Fibronectin was the major newly synthesized and secreted polypeptide of freshly harvested alveolar macrophages. Pulse-chase experiments revealed that newly synthesized fibronectin was rapidly secreted into medium, approximately 50 percent appearing by 1 h and 80 percent by 8 h. Immunoperoxidase staining using antifibronectin F(ab’)(2)-peroxidase conjugates revealed the majority of immunoreactive fibronectin to be intracellular, localized to endoplasmic reticulum and Golgi apparatus. No extracellular matrix fibronectin was visualized, and cell surface staining was rarely seen, usually appearing only at sites where cells were closely apposed and not at sites of macrophage-substrate attachment. Similar immunostaining of fibroblast cultures revealed cell surface-associated fibrillar fibronectin. Ultrastructural localization of fibronectin during binding and phagocytosis of gelatin-coated and plain latex particles revealed fibronectin only on gelatin-latex beads and at their cell binding sites. Neigher plain latex beads nor their cell membrane binding sites stained for fibronectin. These results demonstrate that fibronectin is a major product of human alveolar macrophages, is rapidly secreted, and is localized at cell membrane binding sites for gelatin-coated particles. In view of the known binding properties of fibronectin, it may serve as an endogenous opsonic factor promoting the binding of staphylococcus, denatured collagen, fibrin, or other macromolecules to macrophages in the lower respiratory tract.     

Molecular dynamics of ferrocytochrome c. Magnitude and anisotropy of atomic displacements

Two computer simulations of the atomic motion in tuna ferrocytochrome c have been carried out. The average structures and the structural correlations of the magnitudes of the atomic position fluctuations are in substantial agreement with recent X-ray diffraction results, particularly for the protein interior. The simulations show, however, that the atomic displacements are quite anisotropic. The degree of anisotropy and the preferred directions of atomic displacement exhibit correlations with structural features of the protein.