Assessing the acid-base and conformational properties of histidine residues in human prion protein (125-228) by means of pK(a) calculations and molecular dynamics simulations

A thorough study of the acid-base behavior of the four histidines and the other titratable residues of the structured domain of human prion protein (125-228) is presented. By using multi-tautomer electrostatic calculations, average titration curves have been built for all titratable residues, using the whole bundles of NMR structures determined at pH 4.5 and 7.0. According to our results, (1) only histidine residues are likely to be involved in the first steps of the pH-driven conformational transition of prion protein; (2) the pK(a)'s of His140 and His177 are approximately 7.0, whereas those of His155 and His187 are < 5.5. 10-ns long molecular dynamics simulations have been performed on five different models, corresponding to the most significant combinations of histidine protonation states. A critical comparison between the available NMR structures and our computational results (1) confirms that His155 and His187 are the residues whose protonation is involved in the conformational rearrangement of huPrP in mildly acidic condition, and (2) shows how their protonation leads to the destructuration of the C-terminal part of HB and to the loss of the last turn of HA that represent the crucial microscopic steps of the rearrangement.     

Conformational studies of human islet amyloid peptide using molecular dynamics and simulated annealing methods

Molecular dynamics simulations and simulated annealing in vacuum, model aqueous solution, and simulated membrane were used to analyze the conformational preferences of a segment spanning 20–29 residues of human islet amyloid polypeptide, [referred to as IAPP^H(20–29)]. Molecular dynamics simulations were conducted at 300 K on IAPP^H(20–29). The minimum energy conformers obtained in model aqueous solution and vacuum exhibited similar structures. Even in the absence of any constraints on peptide bonds, trans conformation was preferred consistently by all the peptide bonds. Analysis of the minimum energy conformers indicated that IAPP^H(20–29) showed a strong preference for turn structures in all the environments. These turn structures were stabilized by the formation of hydrogen bonds between the backbone amide and carbonyl groups. A good agreement was found between the results obtained from the molecular dynamics simulation and solid-state nmr experimental studies. © 1998 John Wiley & Sons, Inc. Biopoly 45: 9–20, 1998     

Computer simulations of cyclic enkephalin analogues

Molecular dynamics simulations and energy minimization studies of cyclic enkephalin analogues incorporating retro-inverso modifications have been carried out. The dynamic trajectories are analyzed in terms of the relative mobility of the 14-membered rings, conformational transitions among equilibrium states, and hydrogen-bonding patterns. The cyclization of the molecules reduces the motion of the ring structures substantially. Time-correlated conformational transitions resulting in the reorientation of peptide units are observed. Hydrogen bonds form principally C7 structures. Because of the incorporation of retro-inverso residues, C6 and C8 structures are also formed. Starting conformations for energy minimizations were obtained from the molecular dynamics simulations and from a systematic search of the conformational space available to the molecules. Several minimum energy backbone and side-chain conformations were found for each analogue. The effect of retro-inverso residues on hydrogen-bonding patterns and backbone conformations is discussed.     

Amyloid beta-protein monomer folding: free-energy surfaces reveal alloform-specific differences

Alloform-specific differences in structural dynamics between amyloid β-protein (Aβ) 40 and Aβ42 appear to underlie the pathogenesis of Alzheimer's disease. To elucidate these differences, we performed microsecond timescale replica-exchange molecular dynamics simulations to sample the conformational space of the Aβ monomer and constructed its free-energy surface. We find that neither peptide monomer is unstructured, but rather that each may be described as a unique statistical coil in which five relatively independent folding units exist, comprising residues 1-5, 10-13, 17-22, 28-37, and 39-42, which are connected by four turn structures. The free-energy surfaces of both peptides are characterized by two large basins, comprising conformers with either substantial α-helix or β-sheet content. Conformational transitions within and between these basins are rapid. The two additional hydrophobic residues at the Aβ42 C-terminus, Ile41 and Ala42, significantly increase contacts within the C-terminus, and between the C-terminus and the central hydrophobic cluster (Leu17-Ala21). As a result, the β-structure of Aβ42 is more stable than that of Aβ40, and the conformational equilibrium in Aβ42 shifts towards β-structure. These results suggest that drugs stabilizing α-helical Aβ conformers (or destabilizing the β-sheet state) would block formation of neurotoxic oligomers. The atomic-resolution conformer structures determined in our simulations may serve as useful targets for this purpose. The conformers also provide starting points for simulations of Aβ oligomerization—a process postulated to be the key pathogenetic event in Alzheimer's disease.     

The sequence TGAAKAVALVL from glyceraldehyde-3-phosphate dehydrogenase displays structural ambivalence and interconverts between alpha-helical and beta-hairpin conformations mediated by collapsed conformational states.

The peptide TGAAKAVALVL from glyceraldehyde-3-phosphate dehydrogenase adopts a helical conformation in the crystal structure and is a site for two hydrated helical segments, which are thought to be helical folding intermediates. Overlapping sequences of four to five residues from the peptide, sample both helical and strand conformations in known protein structures, which are dissimilar to glyceraldehyde-3-phosphate dehydrogenase suggesting that the peptide may have a structural ambivalence. Molecular dynamics simulations of the peptide sequence performed for a total simulation time of 1.2 micros, starting from the various initial conformations using GROMOS96 force field under NVT conditions, show that the peptide samples a large number of conformational forms with transitions from alpha-helix to beta-hairpin and vice versa. The peptide, therefore, displays a structural ambivalence. The mechanism from alpha-helix to beta-hairpin transition and vice versa reveals that the compact bends and turns conformational forms mediate such conformational transitions. These compact structures including helices and hairpins have similar hydrophobic radius of gyration (Rgh) values suggesting that similar hydrophobic interactions govern these conformational forms. The distribution of conformational energies is Gaussian with helix sampling lowest energy followed by the hairpins and coil. The lowest potential energy of the full helix may enable the peptide to take up helical conformation in the crystal structure of the glyceraldehyde-3-phosphate dehydrogenase, even though the peptide has a preference for hairpin too. The relevance of folding and unfolding events observed in our simulations to hydrophobic collapse model of protein folding are discussed.     

Caveolin-1 hydrophobic segment peptides insertion into membrane mimetic systems: role of proline residue

Caveolin-1 has a segment of hydrophobic amino acids comprising approximately residues 103-122 that are anchored to the membrane with cholesterol-rich domains. Previously, we reported that changing the Pro(110) residue to Ala (the P110A mutant) prevents not only the localization of the protein into lipid rafts but also the formation and functioning of caveolae. The conformational state of caveolin-1 can be shifted toward the transmembrane arrangement by this single amino acid mutation. To model the conformation, and extent of membrane insertion of this segment into membrane-mimetic environments, we have prepared a peptide corresponding to this hydrophobic segment of caveolin-1 having the sequence KKKKLSTIFGIPMALIWGIYFAILKKKKK-amide and the mutated version, KKKKLSTIFGIAMALIWGIYFAILKKKKK-amide. These peptides contain flanking Lys residues to facilitate purification and handling of the peptide. Circular dichroism measurements demonstrated that the mutated peptide has increased helical content compared with the wild type both in the presence and absence of lipid. The fluorescence emission from the Trp residues in the peptide showed significant blue shifts in the presence of liposomes, however the presence of cholesterol in hydrated vesicle bilayers decreases its helical content. Our overall findings support our studies with the intact protein in cells and suggest that the peptide of WT caveolin-1 hydrophobic segment has an intrinsic preference not to maintain its conformation as a rigid transmembrane helix. Substituting the Pro residue with an Ala allows the peptide to exist in a more hydrophobic environment likely as a consequence of a change in its conformation to a straight hydrophobic helix that traverses the membrane.     

Protein conformational dynamics of crambin in crystal, solution and in the trajectories of molecular dynamics simulations

Atomic displacement parameters — B factors of the eight crambin crystal structures obtained at 0.54–1.5 Å resolution and temperatures of 100–293 K have been analyzed. The comparable contributions to the B factor values are the intramolecular motions which are modeled by the harmonic vibration calculations and derived from the molecular dynamics simulation (MD) as well as rigid body changes in the position of a protein molecule as a whole. In solution for the average NMR structure of crambin the amplitudes of the backbone atomic fluctuations of the most residues of the segments with the regular backbone conformations are close to the amplitudes of the small scale harmonic vibrations. For the same residues the probability of the medium scale fluctuations fixed by the hydrogen exchange method is very low. The restricted conformational mobility of those segments is coupled with the depressed amplitudes of the fluctuation changes of the tertiary structure registered by the residue accessibility changes in an ensemble of NMR structures that forms the average NMR structure of crambin. The amplitudes of temperature fluctuations of backbone atoms and the tertiary structure raise in the segment with the irregular conformations, turn and loops. In the same segments the amplitudes of the calculated harmonic vibrations also increase, but to a lesser extent and especially in the interhelical loop with the most strong and complicated fluctuation changes of the backbone conformation. In solution for the NMR structure in this loop the conformational transitions occur between the conformational substates separated by the energy barriers, but they are not observed even in the long 100 ns trajectories from the MD simulation of crambin. These strong local fluctuation changes of the structure may play a key role in the protein functioning and modern performance improvements in the MD simulation techniques are oriented to increase the probability of protein appearance in the trajectories from the MD simulations.     

Elucidation of interactions regulating conformational stability and dynamics of SARS-CoV-2 S-protein

The ongoing COVID-19 pandemic caused by the new coronavirus, SARS-CoV-2, calls for urgent developments of vaccines and antiviral drugs. The spike protein of SARS-CoV-2 (S-protein), which consists of trimeric polypeptide chains with glycosylated residues on the surface, triggers the virus entry into a host cell. Extensive structural and functional studies on this protein have rapidly advanced our understanding of the S-protein structure at atomic resolutions, although most of these structural studies overlook the effect of glycans attached to the S-protein on the conformational stability and functional motions between the inactive down and active up forms. Here, we performed all-atom molecular dynamics simulations of both down and up forms of a fully glycosylated S-protein in solution as well as targeted molecular dynamics simulations between them to elucidate key interdomain interactions for stabilizing each form and inducing the large-scale conformational transitions. The residue-level interaction analysis of the simulation trajectories detects distinct amino acid residues and N-glycans as determinants on conformational stability of each form. During the conformational transitions between them, interdomain interactions mediated by glycosylated residues are switched to play key roles on the stabilization of another form. Electrostatic interactions, as well as hydrogen bonds between the three receptor binding domains, work as driving forces to initiate the conformational transitions toward the active form. This study sheds light on the mechanisms underlying conformational stability and functional motions of the S-protein, which are relevant for vaccine and antiviral drug developments.     

Synthesis of substrates and inhibitors of botulinum neurotoxin type A metalloprotease.

Botulinum neurotoxin (BoNT) metalloproteases and related proteases are the most selective proteases known. X-ray crystal structures suggest that the active sites of the native enzymes exist in catalytically incompetent forms that must be activated by substrate binding. In order to characterize the postulated substrate-induced conformational changes for enzyme activation, we synthesized a series of transition-state analog inhibitors in which the dipeptide cleavage site is replaced by tetrahedral intermediate analogs within the minimal substrate peptide sequence. In this paper, we report our efforts to design inhibitors of BoNT/A metalloprotease. We confirm that an effective substrate sequence for BoNT/A metalloprotease is a 17-mer peptide corresponding to residues 187-203 of SNAP-25. A more stable substrate, Nle202SNAP-25 [187-203] was synthesized in order to develop an assay for proteolytic activity of BoNT/A metalloprotease that can be used to monitor time-dependent inhibition. Alpha-thiol amide analogs of Gln-197 were incorporated via solid-phase peptide synthesis into both 17-mer minimal peptide substrate sequences. The synthesis, characterization and inhibition kinetics for the alpha-thiol amide analogs of holotoxin A substrate are described. These substrate-derived inhibitors were shown to be submicromolar inhibitors of BoNT/A catalytic activity.     

Human EGF-derived direct and reverse short linear motifs: conformational dynamics insight into the receptor-binding residues

Short linear motifs (SLiMs) have been recognized to perform diverse functions in a variety of regulatory proteins through the involvement in protein–protein interactions, signal transduction, cell cycle regulation, protein secretion, etc. However, detailed molecular mechanisms underlying their functions including roles of definite amino acid residues remain obscure. In our previous studies, we demonstrated that conformational dynamics of amino acid residues in oligopeptides derived from regulatory proteins such as alpha-fetoprotein (AFP), carcino-embryonic antigen (CEA), and pregnancy specific β1-glycoproteins (PSGs) contributes greatly to their biological activities. In the present work, we revealed the 22-member linear modules composed of direct and reverse AFP_14–20-like heptapeptide motifs linked by CxxGY/FxGx consensus motif within epidermal growth factor (EGF), growth factors of EGF family and numerous regulatory proteins containing EGF-like modules. We showed, first, the existence of similarity in amino acid signatures of both direct and reverse motifs in terms of their physicochemical properties. Second, molecular dynamics (MD) simulation study demonstrated that key receptor-binding residues in human EGF in the aligned positions of the direct and reverse motifs may have similar distribution of conformational probability densities and dynamic behavior despite their distinct physicochemical properties. Third, we found that the length of a polypeptide chain (from 7 to 53 residues) has no effect, while disulfide bridging and backbone direction significantly influence the conformational distribution and dynamics of the residues. Our data may contribute to the atomic level structure–function analysis and protein structure decoding; additionally, they may provide a basis for novel protein/peptide engineering and peptide-mimetic drug design.     

Enhanced sampling of peptide and protein conformations using replica exchange simulations with a peptide backbone biasing-potential

During replica exchange molecular dynamics (RexMD) simulations, several replicas of a system are simulated at different temperatures in parallel allowing for exchange between replicas at frequent intervals. This technique allows significantly improved sampling of conformational space and is increasingly being used for structure prediction of peptides and proteins. A drawback of the standard temperature RexMD is the rapid increase of the replica number with increasing system size to cover a desired temperature range. In an effort to limit the number of replicas, a new Hamiltonian-RexMD method has been developed that is specifically designed to enhance the sampling of peptide and protein conformations by applying various levels of a backbone biasing potential for each replica run. The biasing potential lowers the barrier for backbone dihedral transitions and promotes enhanced peptide backbone transitions along the replica coordinate. The application on several peptide cases including in all cases explicit solvent indicates significantly improved conformational sampling when compared with standard MD simulations. This was achieved with a very modest number of 5-7 replicas for each simulation system making it ideally suited for peptide and protein folding simulations as well as refinement of protein model structures in the presence of explicit solvent.     

Conformational choreography of a molecular switch region in myelin basic protein--molecular dynamics shows induced folding and secondary structure type conversion upon threonyl phosphorylation in both aqueous and membrane-associated environments

The 18.5 kDa isoform of myelin basic protein is essential to maintaining the close apposition of myelin membranes in central nervous system myelin, but its intrinsic disorder (conformational dependence on environment), a variety of post-translational modifications, and a diversity of protein ligands (e.g., actin and tubulin) all indicate it to be multifunctional. We have performed molecular dynamics simulations of a conserved central segment of 18.5 kDa myelin basic protein (residues Glu80-Gly103, murine sequence numbering) in aqueous and membrane-associated environments to ascertain the stability of constituent secondary structure elements (α-helix from Glu80-Val91 and extended poly-proline type II from Thr92-Gly103) and the effects of phosphorylation of residues Thr92 and Thr95, individually and together. In aqueous solution, all four forms of the peptide bent in the middle to form a hydrophobic cluster. The phosphorylated variants were stabilized further by electrostatic interactions and formation of β-structures, in agreement with previous spectroscopic data. In simulations performed with the peptide in association with a dimyristoylphosphatidylcholine bilayer, the amphipathic α-helical segment remained stable and membrane-associated, although the degree of penetration was less in the phosphorylated variants, and the tilt of the α-helix with respect to the plane of the membrane also changed significantly with the modifications. The extended segment adjacent to this α-helix represents a putative SH3-ligand and remained exposed to the cytoplasm (and thus accessible to binding partners). The results of these simulations demonstrate how this segment of the protein can act as a molecular switch: an amphipathic α-helical segment of the protein is membrane-associated and presents a subsequent proline-rich segment to the cytoplasm for interaction with other proteins. Phosphorylation of threonyl residues alters the degree of membrane penetration of the α-helix and the accessibility of the proline-rich ligand and can stabilize a β-bend. A bend in this region of 18.5 kDa myelin basic protein suggests that the N- and C-termini of the proteins can interact with different leaflets of the myelin membrane and explain how a single protein can bring them close together.     

Simulation of the activation of alpha-chymotrypsin: analysis of the pathway and role of the propeptide

Alpha-chymotrypsin undergoes a reversible conformational change from an inactive chymotrypsinogen-like structure at high pH to an active conformation at neutral pH. In order to gain insight into this process on a structural level, we applied molecular dynamics and targeted molecular dynamics simulations in aqueous environment on the activation and inactivation processes of three different types of chymotrypsin. These are the wild-type bovine chymotrypsin containing the propeptide and the bovine and rat chymotrypsin lacking the propeptide. From these simulations, the importance of the propeptide and of the sequence differences between the rat and bovine variants from the viewpoint of activation could be evaluated and compared with previous fluorescence stopped flow results. The obtained results show the unambiguous influence of the propeptide on the explored conformational space, whereas the sequence differences between bovine and rat chymotrypsin play a minor role. The main features of activation are present in both the wild type and the variant lacking the propeptide, despite the fact that different parts of the conformational space were explored. The comparison of all trajectories shows that particular amino acid residues, such as 17, 18, 19, 187, 217, 218, and 223, undergo large dihedral transitions during the activation process, suggesting a role as hinge residues during the conformational change.     

A thermodynamic approach to the conformational preferences of the 180–195 segment derived from the human prion protein α2‐helix

On consideration that intrinsic structural weakness could affect the segment spanning the α2‐helical residues 173–195 of the PrP, we have investigated the conformational stabilities of some synthetic Ala‐scanned analogs of the peptide derived from the 180–195 C‐terminal sequence, using a novel approach whose theoretical basis originates from protein thermodynamics. Even though a quantitative comparison among peptides could not be assessed to rank them according to the effect caused by single amino acid substitution, as a general trend, all peptides invariably showed an appreciable preference for an α‐type organization, consistently with the fact that the wild‐type sequence is organized as an α‐helix in the native protein. Moreover, the substitution of whatever single amino acid in the wild‐type sequence reduced the gap between the α‐ and the β‐propensity, invariably enhancing the latter, but in any case this gap was larger than that evaluated for the full‐length α2‐helix‐derived peptide. It appears that the low β‐conformation propensity of the 180–195 region depends on the simultaneous presence of all of the Ala‐scanned residues, indirectly confirming that the N‐terminal 173–179 segment could play a major role in determining the chameleon conformational behavior of the entire 173–195 region in the PrP. Copyright © 2008 European Peptide Society and John Wiley & Sons, Ltd.     

Molecular dynamics simulation studies of the structural response of an isolated Aβ1–42 monomer localized in the vicinity of the hydrophilic TiO2 surface

We have probed the effect of a model hydrophilic surface, rutile TiO₂, on the full-length amyloid beta (Aβ_1–42) monomer using molecular dynamics simulations. The rutile surface brings about sharp changes in the peptide’s intrinsic behavior in a distance-dependent manner. The intrinsic collapse of the peptide is disrupted, while the β-sheet propensity is sharply enhanced with increased proximity to the surface. The results may have implications for Aβ self-assembly and fibrillogenesis on hydrophilic surfaces and should be taken into consideration in the design of novel nanomaterials for perturbing amyloidogenic behavior.     

Enhanced sampling protocol to elucidate fusion peptide opening of SARS-CoV-2 spike protein

Large-scale conformational transitions in the spike protein S2 domain are required during host-cell infection of the severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) virus. Although conventional molecular dynamics simulations have been extensively used to study therapeutic targets of SARS-CoV-2, it is still challenging to gain molecular insight into the key conformational changes because of the size of the spike protein and the long timescale required to capture these transitions. In this work, we have developed an efficient simulation protocol that leverages many short simulations, a dynamic selection algorithm, and Markov state models to interrogate the structural changes of the S2 domain. We discovered that the conformational flexibility of the dynamic region upstream of the fusion peptide in S2 is coupled to the proteolytic cleavage state of the spike protein. These results suggest that opening of the fusion peptide likely occurs on a submicrosecond timescale after cleavage at the S2′ site. Building on the structural and dynamical information gained to date about S2 domain dynamics, we provide proof of principle that a small molecule bound to a seam neighboring the fusion peptide can slow the opening of the fusion peptide, leading to a new inhibition strategy for experiments to confirm. In aggregate, these results will aid the development of drug cocktails to inhibit infections caused by SARS-CoV-2 and other coronaviruses.