Watts, Charles R.

Matrix metaloproteinase-2 (MMP-2) is an extracellular Zn2+ protease specific to type I and IV collagens. Its expression is associated with several inflammatory, degenerative, and malignant diseases. Conformational properties, domain movements, and interactions between MMP-2 and its associated metal ions were characterized using a 1.0 µs molecular dynamics simulation. Dihedral principle component analysis revealed ten families of conformations with the greatest degree of variability occurring in the link region connecting the catalytic and hemopexin domains. Dynamic cross-correlation analysis indicated domain movements corresponding to the opening and closing of the hemopexin domain in relation to the fibronectin and catalytic domains facilitated by the link region. Interaction energies were calculated using the molecular mechanics Poisson Boltzman surface area-interaction entropy (MMPBSA-IE) analysis method and revealed strong binding energies for the catalytic Zn2+ ion 1, Ca2+ ion 1, and Ca2+ ion 3 with significant conformational stability at the binding sites of Zn2+ ion 1 and Ca2+ ion 1. Ca2+ ion 2 diffuses freely away from its crystallographically defined binding site. Zn2+ ion 2 plays a minor role in conformational stability of the catalytic domain while Ca2+ ion 3 is strongly attracted to the highly electronegative sidechains of the Asp residues around the central β-sheet core of the hemopexin domain; however, the interacting residue sidechain carboxyl groups are outside of Ca2+ ion 3′s coordination sphere.

Chlorotoxin (CTX) is a 36–amino acid peptide with eight Cys residues that forms four disulfide bonds. It has high affinity for the glioma-specific chloride channel and matrix metalloprotease-2. Structural and binding properties of CTX analogs with various Cys residue substitutions with l-α-aminobutyric acid (Abu) have been previously reported. Using 4.2 µs molecular dynamics, we compared the conformational and essential space sampling of CTX and analogs with selective substitution of the Cys residues and associated disulfide bonds with either Abu or Ser. The native and substituted peptides maintained a high degree of α-helix propensity from residues 8 through 21, with the exception of substitution of the Cys5–Cys28 residues with Ser and the Cys16–Cys33 residues with Abu. In agreement with previous circular dichroism spectropolarimetry results, the C-terminal β-sheet content varied less from residues 25 through 29 and 32 through 36 and was well conserved in most analogs. The Cys16–Cys33 and Cys20–Cys35 disulfide-bonded residues appear to be required to maintain the αβ motif of CTX. Selective substitution with the hydrophilic Ser, may mitigate the destabilizing effect of Cys16–Cys33 substitution through the formation of an inter residue H-bond from Ser16:OγH to Ser33:OγH bridged by a water molecule. All peptides shared considerable sampled conformational space, which explains the retained receptor binding of the non-native analogs.

Abstract The conformational space and structural ensembles of amyloid beta (Aβ) peptides and their oligomers in solution are inherently disordered and proven to be challenging to study. Optimum force field selection for molecular dynamics (MD) simulations and the biophysical relevance of results are still unknown. We compared the conformational space of the Aβ(1‐40) dimers by 300 ns replica exchange MD simulations at physiological temperature (310 K) using: the AMBER‐ff99sb‐ILDN, AMBER‐ff99sb*‐ILDN, AMBER‐ff99sb‐NMR, and CHARMM22* force fields. Statistical comparisons of simulation results to experimental data and previously published simulations utilizing the CHARMM22* and CHARMM36 force fields were performed. All force fields yield sampled ensembles of conformations with collision cross sectional areas for the dimer that are statistically significantly larger than experimental results. All force fields, with the exception of AMBER‐ff99sb‐ILDN (8.8 ± 6.4%) and CHARMM36 (2.7 ± 4.2%), tend to overestimate the α‐helical content compared to experimental CD (5.3 ± 5.2%). Using the AMBER‐ff99sb‐NMR force field resulted in the greatest degree of variance (41.3 ± 12.9%). Except for the AMBER‐ff99sb‐NMR force field, the others tended to under estimate the expected amount of β‐sheet and over estimate the amount of turn/bend/random coil conformations. All force fields, with the exception AMBER‐ff99sb‐NMR, reproduce a theoretically expected β‐sheet‐turn‐β‐sheet conformational motif, however, only the CHARMM22* and CHARMM36 force fields yield results compatible with collapse of the central and C‐terminal hydrophobic cores from residues 17‐21 and 30‐36. Although analyses of essential subspace sampling showed only minor variations between force fields, secondary structures of lowest energy conformers are different.

ABSTRACTReplica exchange molecular dynamics simulations (300 ns) were used to study the dimerization of amyloid β(1‐40) (Aβ(1‐40)) polypeptide. Configurational entropy calculations revealed that at physiological temperature (310 K, 37°C) dynamic dimers are formed by randomly docked monomers. Free energy of binding of the two chains to each other was −93.56 ± 6.341 kJ mol−1. Prevalence of random coil conformations was found for both chains with the exceptions of increased β‐sheet content from residues 16‐21 and 29‐32 of chain A and residues 15‐21 and 30‐33 of chain B with β‐turn/β‐bend conformations in both chains from residues 1‐16, 21‐29 of chain A, 1‐16, and 21‐29 of chain B. There is a mixed β‐turn/β‐sheet region from residues 33‐38 of both chains. Analysis of intra‐ and interchain residue distances shows that, although the individual chains are highly flexible, the dimer system stays in a loosely packed antiparallel β‐sheet configuration with contacts between residues 17‐21 of chain A with residues 17‐21 and 31‐36 of chain B as well as residues 31‐36 of chain A with residues 17‐21 and 31‐36 of chain B. Based on dihedral principal component analysis, the antiparallel β‐sheet‐loop‐β‐sheet conformational motif is favored for many low energy sampled conformations. Our results show that Aβ(1‐40) can form dynamic dimers in aqueous solution that have significant conformational flexibility and are stabilized by collapse of the central and C‐terminal hydrophobic cores with the expected β‐sheet‐loop‐β‐sheet conformational motif. Proteins 2017; 85:1024–1045. © 2017 Wiley Periodicals, Inc.