Abstract:How is the native structure encoded in the amino acid sequence? For the traditional backbone centric view, the dominant forces are hydrogen bonds (backbone) and φ-ψ propensity. The role of hydrophobicity is non-specific. For the side-chain centric view, the dominant force of protein folding is hydrophobicity. In order to understand the balance between backbone and side-chain forces, we have studied the contributions of three components of a β-hairpin peptide: turn, backbone hydrogen bonding and side-chain interactions, of a 16-residue fragment of protein G. The peptide folds rapidly and cooperatively to a conformation with a defined secondary structure and a packed hydrophobic cluster of aromatic side-chains. Our strategy is to observe the structural stability of the β-hairpin under systematic perturbations of the turn region, backbone hydrogen bonds and the hydrophobic core formed by the side-chains, respectively. In our molecular dynamics simulations, the peptides are solvated. with explicit water molecules, and an all-atom force field (CFF91) is used. Starting from the original peptide (G41EWTYDDATKTFTVTE56), we carried out the following MD simulations. (1) unfolding at 350 K; (2) forcing the distance between the Cα atoms of ASP47 and LYS50 to be 8 Å; (3) deleting two turn residues (Ala48 and Thr49) to form a β-sheet complex of two short peptides, GEWTYDD and KTFTVTE; (4) four hydrophobic residues (W43, Y45, F52 and T53) are replaced by a glycine residue step-by-step; and (5) most importantly, four amide hydrogen atoms (T44, D46, T53, and T55, which are crucial for backbone hydrogen bonding), are substituted by fluorine atoms. The fluorination not only makes it impossible to form attractive hydrogen bonding between the two β-hairpin strands, but also introduces a repulsive force between the two strands due to the negative charges on the fluorine and oxygen atoms. Throughout all simulations, we observe that backbone hydrogen bonds are very sensitive to the perturbations and are easily broken. In contrast, the hydrophobic core survives most perturbations. In the decisive test of fluorination, the fluorinated peptide remains folded under our simulation conditions (5 ns, 278 K). Hydrophobic interactions keep the peptide folded, even with a repulsive force between the β-strands. Thus, our results strongly support a side-chain centric view for protein folding.

Local Structures in Proteins from Microsecond Molecular Dynamics Simulations: 2. The Role of Symmetry in GTPase Binding and Dimer Formation

Local Structures in Proteins from Microsecond Molecular Dynamics Simulations: A Symmetry-Based Perspective

Analysis of 15N-1H NMR Relaxation in Proteins by a Combined Experimental and Molecular Dynamics Simulation Approach: Picosecond-Nanosecond Dynamics of the Rho GTPase Binding Domain of Plexin-B1 in the Dimeric State Indicates Allosteric Pathways.

Exploring The Structure And Dynamics Of All-Atom Models For The Plexin Transmembrane Receptor Bound To Gtpases And To Lipid Bilayer

Investigation on the Interaction Between Plexin Intracellular Plus Transmembrane Domains with GTPases and with the Lipid Bilayer Using All-Atom Molecular Dynamics Simulations

Transient Structural Properties of the Rho GDP‐Dissociation Inhibitor

Intrinsically disordered membrane anchors of Rheb, RhoA and DiRas3 small GTPases: Molecular dynamics, membrane organization, and interactions

A putative structural mechanism underlying the antithetic effect of homologous RND1 and RhoD GTPases in mammalian plexin regulation

A Direct Coupling between Global and Internal Motions in a Single Domain Protein? MD Investigation of Extreme Scenarios

PDB-NMA of a protein homodimer reproduces distinct experimental motility asymmetry

Dynamic Protein-Protein Complexes: How Alternative Interactions Create Ensembles and How Solution Nmr and Md Simulations Can Characterize Them

Structure of an inactive conformation of GTP-bound RhoA GTPase

Local symmetry determines the phases of linear chains: a simple model for the self-assembly of peptide

Ras Isoforms Conformational Clustering And Community Networks Studies: Simulating Ras With Accelerated Molecular Dynamics

The dynamics of spatio-temporal Rho GTPase signaling: formation of signaling patterns

Structure and Dynamics of the Rhomboid Protease GlpG in Liposomes Studied by Solid-State NMR.

Prediction, Refinement, and Persistency of Transmembrane Helix Dimers in Lipid Bilayers Using Implicit and Explicit Solvent/lipid Representations: Microsecond Molecular Dynamics Simulations of ErbB1/B2 and EphA1.

Mechanical Unfolding of Segment-Swapped Protein G Dimer: Results from Replica Exchange Molecular Dynamics Simulations

SEC‐SAXS/MC Ensemble Structural Studies of the Microtubule Binding Protein Cdt1 Show Monomeric, Folded‐Over Conformations

Characterizing Conformational Ensemble of H‐, K‐, and N‐Ras Using Accelerated Molecular Dynamics and Wide‐Angle X‐ray Solution Scattering (WAXS)

Molecular dynamics simulations of a β-hairpin fragment of protein G: balance between side-chain and backbone forces