Coarse-grained free energy functions for studying protein conformational changes: a double-well network model

In this work, a double-well network model (DWNM) is presented for generating a coarse-grained free energy function that can be used to study the transition between reference conformational states of a protein molecule. Compared to earlier work that uses a single, multidimensional double-well potential to connect two conformational states, the DWNM uses a set of interconnected double-well potentials for this purpose. The DWNM free energy function has multiple intermediate states and saddle points, and is hence a “rough” free energy landscape. In this implementation of the DWNM, the free energy function is reduced to an elastic-network model representation near the two reference states. The effects of free energy function roughness on the reaction pathways of protein conformational change is demonstrated by applying the DWNM to the conformational changes of two protein systems: the coil-to-helix transition of the DB-loop in G-actin and the open-to-closed transition of adenylate kinase. In both systems, the rough free energy function of the DWNM leads to the identification of distinct minimum free energy paths connecting two conformational states. These results indicate that while the elastic-network model captures the low-frequency vibrational motions of a protein, the roughness in the free energy function introduced by the DWNM can be used to characterize the transition mechanism between protein conformations.     

On the roles of substrate binding and hinge unfolding in conformational changes of adenylate kinase

We characterized the conformational change of adenylate kinase (AK) between open and closed forms by conducting five all-atom molecular-dynamics simulations, each of 100 ns duration. Different initial structures and substrate binding configurations were used to probe the pathways of AK conformational change in explicit solvent, and no bias potential was applied. A complete closed-to-open and a partial open-to-closed transition were observed, demonstrating the direct impact of substrate-mediated interactions on shifting protein conformation. The sampled configurations suggest two possible pathways for connecting the open and closed structures of AK, affirming the prediction made based on available x-ray structures and earlier works of coarse-grained modeling. The trajectories of the all-atom molecular-dynamics simulations revealed the complexity of protein dynamics and the coupling between different domains during conformational change. Calculations of solvent density and density fluctuations surrounding AK did not show prominent variation during the transition between closed and open forms. Finally, we characterized the effects of local unfolding of an important hinge near Pro¹⁷⁷ on the closed-to-open transition of AK and identified a novel mechanism by which hinge unfolding modulates protein conformational change. The local unfolding of Pro¹⁷⁷ hinge induces alternative tertiary contacts that stabilize the closed structure and prevent the opening transition.     

Atomistic and coarse-grained analysis of double spectrin repeat units: the molecular origins of flexibility

Spectrin is an ubiquitous protein in metazoan cells, and its flexibility is one of the keys to maintaining cellular structure and organization. Both alpha-spectrin and beta-spectrin polypeptides consist primarily of triple coiled-coil modular repeat units, and two important factors that determine spectrin flexibility are the bending flexibility between two consecutive repeat units and the conformational flexibility of individual repeat units. Atomistic molecular dynamics (MD) simulations are used here to study double spectrin repeat units (DSRUs) from the human erythrocyte beta-spectrin (HEbeta89) and the chicken brain alpha-spectrin (CBalpha1617). From the results of MD simulations, a highly conserved Trp residue in the A-helix of most repeat units that has been suggested to be important in conferring stability to the coiled-coil structures is found not to have a significant effect on the conformational flexibility of individual repeat units. Characterization of the bending flexibility for two consecutive repeats of spectrin via atomistic simulations and coarse-grained (CG) modeling indicate that the bending flexibility is governed by the interactions between the AB-loop of the first repeat unit, the BC-loop of the second repeat unit and the linker region. Specifically, interactions between residues in these regions can lead to a strong directionality in the bending behavior of two repeat units. The biological implications of these finding are discussed.     

Coarse-grained modeling of the actin filament derived from atomistic-scale simulations

A coarse-grained (CG) procedure that incorporates the information obtained from all-atom molecular dynamics (MD) simulations is presented and applied to actin filaments (F-actin). This procedure matches the averaged values and fluctuations of the effective internal coordinates that are used to define a CG model to the values extracted from atomistic MD simulations. The fluctuations of effective internal coordinates in a CG model are computed via normal-mode analysis (NMA), and the computed fluctuations are matched with the atomistic MD results in a self-consistent manner. Each actin monomer (G-actin) is coarse-grained into four sites, and each site corresponds to one of the subdomains of G-actin. The potential energy of a CG G-actin contains three bonds, two angles, and one dihedral angle; effective harmonic bonds are used to describe the intermonomer interactions in a CG F-actin. The persistence length of a CG F-actin was found to be sensitive to the cut-off distance of assigning intermonomer bonds. Effective harmonic bonds for a monomer with its third nearest neighboring monomers are found to be necessary to reproduce the values of persistence length obtained from all-atom MD simulations. Compared to the elastic network model, incorporating the information of internal coordinate fluctuations enhances the accuracy and robustness for a CG model to describe the shapes of low-frequency vibrational modes. Combining the fluctuation-matching CG procedure and NMA, the achievable time- and length scales of modeling actin filaments can be greatly enhanced. In particular, a method is described to compute the force-extension curve using the CG model developed in this work and NMA. It was found that F-actin is easily buckled under compressive deformation, and a writhing mode is developed as a result. In addition to the bending and twisting modes, this novel writhing mode of F-actin could also play important roles in the interactions of F-actin with actin-binding proteins and in the force-generation process via polymerization.     

Dynamic ejection behaviour of water molecules passing through a nano-aperture nozzle

This study used molecular dynamics (MD) simulation to investigate the passage of water molecules through a composite graphene/Au nano-nozzle. Our focus was on the degree to which system temperature, extrusion speed, and nozzle diameter affect jet dynamics and the associated transient phenomena. Our findings show that high pressure and spatial confinement cause the nanojet from a small nozzle diameter (1.0?nm) to bend and twist, whereas the jets from a nozzle with a diameter of 1.5?nm present columns of greater stability. At 100?K, the H₂O nanojet froze at the outlet of the nozzle in the form of condensed icicles. At 500?K, the H₂O nanojet formed a loose spray and gaseous clusters. High extrusion speed of 55.824?m/s produced recirculating flow downstream from the nanojet with the appearance of an erupting volcano, which further prompted the jet column to thicken. Lower extrusion speeds produced jets with flow velocity insufficient to overcome the capillary force at the outlet of the nozzle, which subsequently manifests as unstable fluctuations in the flow rate.

• HIGHLIGHTS

• Water molecules through a composite graphene/Au nano-nozzle forming a nanojet is investigated.

• High pressure and spatial confinement cause the nanojet from a small nozzle diameter (≤1.0?nm) to bend and twist.

• High extrusion speed (≧55.824?m/s) produced recirculating flow downstream from the nanojet.

• Figure abstract: Schematic of the H₂O nano-jet through a nano-nozzle of graphene/Au

     

Synthesis, cytotoxicity and human telomerase inhibition activities of a series of 1,2-heteroannelated anthraquinones and anthra[1,2-d]imidazole-6,11-dione homologues

A series of 1,2-heteroannelated anthraquinones and anthra[1,2-d]imidazole-6,11-dione tetracyclic analogues with different side chain were prepared using an various synthetic route via acylation, cyclization, condensation, and intramolecular heterocyclization. Tetracyclic system containing alkyl and aryl, aromatic and heterocyclic, linear and cyclic, polar and apolar, and basic and acids residues were incorporated. They were evaluated for their effects on telomerase activity, hTERT expression, cell proliferations, and in vitro cytotoxicity against NCI’s 60 cell line human tumor screen. Compounds 4, 11, 12, 14, 15, 16, 17, 19, 20, 23, 25, and 26 were selected by the NCI for one dose screening program and further studies on 4, 23 and 25 where the curves cross these lines represent the interpolated values to cause 50% growth inhibition (GI₅₀), total growth inhibition (TGI) and 50% cell killing (LC₅₀), respectively. Further studies did not reveal any compound that showed potent and significant on telomerase inhibitory activity and hTERT repressing ability. Comparative testing of these compounds in the NCI’s screen revealed varying levels of potency and differential cytotoxicity, apparently related to the unsaturation levels in and substitution patterns on the core ring system. It appeared that addition of a fourth planar aromatic system to a tricyclic chromophore might enhances potent cytotoxic agents, at a level equivalent to a second side chain in one of the tricyclic series. Although the exact mechanism of how this pharmacophore contributes to its activity is still unclear, however, the group in the extended arm of the tetracyclic system might contribute to proper binding to the residues within the grove of G-quadruplex structure.     

Reticulon 3 interacts with NS4B of the hepatitis C virus and negatively regulates viral replication by disrupting NS4B self‐interaction

The non‐structural protein 4B (NS4B) of the hepatitis C virus (HCV) is an endoplasmic reticulum (ER) membrane protein comprising two consecutive amphipathic α‐helical domains (AH1 and AH2). Its self‐oligomerization via the AH2 domain is required for the formation of the membranous web that is necessary for viral replication. Previously, we reported that the host‐encoded ER‐associated reticulon 3 (RTN3) protein is involved in the formation of the replication‐associated membranes of (+)RNA enteroviruses during viral replication. In this study, we demonstrated that the second transmembrane region of RTN3 competed for, and bound to, the AH2 domain of NS4B, thus abolishing NS4B self‐interaction and leading to the downregulation of viral replication. This interaction was mediated by two crucial residues, lysine 52 and tyrosine 63, of AH2, and was regulated by the AH1 domain. The silencing of RTN3 in Huh7 and AVA5 cells harbouring an HCV replicon enhanced the replication of HCV, which was counteracted by the overexpression of recombinant RTN3. The synthesis of viral RNA was also increased in siRNA‐transfected human primary hepatocytes infected with HCV derived from cell culture. Our results demonstrated that RTN3 acted as a restriction factor to limit the replication of HCV.     

Protein-protein interaction site predictions with three-dimensional probability distributions of interacting atoms on protein surfaces

Protein-protein interactions are key to many biological processes. Computational methodologies devised to predict protein-protein interaction (PPI) sites on protein surfaces are important tools in providing insights into the biological functions of proteins and in developing therapeutics targeting the protein-protein interaction sites. One of the general features of PPI sites is that the core regions from the two interacting protein surfaces are complementary to each other, similar to the interior of proteins in packing density and in the physicochemical nature of the amino acid composition. In this work, we simulated the physicochemical complementarities by constructing three-dimensional probability density maps of non-covalent interacting atoms on the protein surfaces. The interacting probabilities were derived from the interior of known structures. Machine learning algorithms were applied to learn the characteristic patterns of the probability density maps specific to the PPI sites. The trained predictors for PPI sites were cross-validated with the training cases (consisting of 432 proteins) and were tested on an independent dataset (consisting of 142 proteins). The residue-based Matthews correlation coefficient for the independent test set was 0.423; the accuracy, precision, sensitivity, specificity were 0.753, 0.519, 0.677, and 0.779 respectively. The benchmark results indicate that the optimized machine learning models are among the best predictors in identifying PPI sites on protein surfaces. In particular, the PPI site prediction accuracy increases with increasing size of the PPI site and with increasing hydrophobicity in amino acid composition of the PPI interface; the core interface regions are more likely to be recognized with high prediction confidence. The results indicate that the physicochemical complementarity patterns on protein surfaces are important determinants in PPIs, and a substantial portion of the PPI sites can be predicted correctly with the physicochemical complementarity features based on the non-covalent interaction data derived from protein interiors.     

A structural and mechanistic study of the oxidation of methionine residues in hPTH(1-34) via experiments and simulations

The relationship between the conformational properties of 1-34 human parathyroid hormone [hPTH(1-34)] and the oxidation of its methionine residues, Met8 and Met18, by hydrogen peroxide is analyzed as a function of pH by measuring the rates of oxidation and by performing MD simulations with an explicit representation of water molecules. Between pH 4 and pH 8, both Met8 and Met18 have nearly pH independent rates of oxidation, and Met18 is oxidized at a rate that is 90-100% of that of freeMet and 10-20% faster than that of Met8. We also found that average 2SWCNs calculated from MD simulations correlate well to the rates of oxidation of Met8 and Met18. The use of 2SWCNs is based on the mechanism that we proposed, the water-mediated mechanism, in which water molecules stabilize the transition state via specific interactions, but the transfer of protons (acid-catalyzed mechanism) does not play a role [Chu, J. W., and Trout, B. L. (2004) J. Am. Chem. Soc. 126 (3), 900-908]. Only at very low pH values, pH 1 for the oxidation of freeMet, does the acid-catalyzed oxidation mechanism become important. For the oxidation of Met8 and Met18 in hPTH(1-34), the acid-catalyzed mechanism becomes significant at a higher pH value, pH 2, probably due to the proximity of nearby acidic residues to Met8 (Glu4) and Met18 (Glu22). In this study, we have demonstrated that the chemistry of oxidation and the structure of polypeptides can be correlated via a detailed understanding of the reaction mechanism, appropriate sampling of configurational space, and a suitable choice of a structural property, water coordination number.