scholarly journals Multi-State Design of Flexible Proteins Predicts Sequences Optimal for Conformational Change

2019 ◽  
Author(s):  
Marion Sauer ◽  
Alexander M. Sevy ◽  
James E. Crowe ◽  
Jens Meiler

AbstractComputational protein design of an ensemble of conformations for one protein – i.e., multi-state design – determines the side chain identity by optimizing the energetic contributions of that side chain in each of the backbone conformations. Sampling the resulting large sequence-structure search space limits the number of conformations and the size of proteins in multi-state design algorithms. Here, we demonstrated that the REstrained CONvergence (RECON) algorithm can simultaneously evaluate the sequence of large proteins that undergo substantial conformational changes, such as viral surface glycoproteins. Simultaneous optimization of side chain conformations across all conformations resulted in an increase of 30% to 40% in sequence conservation when compared to single-state designs. More importantly, the sampled sequence space of RECON designs resembled the evolutionary sequence space of functional proteins. This finding was especially true for sequence positions that require substantial changes in their local environment across an ensemble of conformations. To quantify this rewiring of contacts at a certain position in sequence and structure, we introduced a new metric designated ‘contact proximity deviation’ that enumerates contact map changes. This measure allows mapping of global conformational changes into local side chain proximity adjustments, a property not captured by traditional global similarity metrics such as RMSD or local similarity metrics such as changes in φ and ψ angles.Author SummaryMulti-state design can be used to engineer proteins that need to exist in multiple conformations or that bind to multiple partner molecules. In essence, multi-state design selects a compromise of protein sequences that allow for an ensemble of protein conformations, or states, associated with a particular biological function. In this paper, we used the REstrained CONvergence (RECON) algorithm with Rosetta to show that multi-state design of flexible proteins predicts sequences optimal for conformational change, mimicking mutation preferences sampled in evolution. Modeling optimal local side chain physicochemical environments within an ensemble selected significantly more native-like sequences than selections performed when all conformations states are designed independently. This outcome was particularly true for amino acids whose local side chain environment change between conformations. To quantify such contact map changes, we introduced a novel metric to show that sequence conservation is dependent on protein flexibility, i.e., changes in local side chain environments between stated limit the space of tolerated mutations. Additionally, such positions in sequence and structure are more likely to be energetically frustrated, at least in some states. Importantly, we showed that multi-state design over an ensemble of conformations (space) can explore evolutionary tolerated sequence space (time), thus enabling RECON to not only design proteins that require multiple states for function but also predict mutations that might be tolerated in native proteins but have not yet been explored by evolution. The latter aspect can be important to anticipate escape mutations, for example in pathogens or oncoproteins.

Molecules ◽  
2018 ◽  
Vol 23 (10) ◽  
pp. 2459
Author(s):  
Udeok Seo ◽  
Ku-Jin Kim ◽  
Beom Kang

In this paper, a novel method to compute side chain conformational variations for a protein molecule tunnel (or channel) is proposed. From the conformational variations, we compute the flexibly deformed shapes of the initial tunnel, and present a way to compute the maximum size of the ligand that can pass through the deformed tunnel. By using the two types of graphs corresponding to amino acids and their side chain rotamers, the suggested algorithm classifies amino acids and rotamers which possibly have collisions. Based on the divide and conquer technique, local side chain conformations are computed first, and then a global conformation is generated by combining them. With the exception of certain cases, experimental results show that the algorithm finds up to 327,680 valid side chain conformations from 128~1233 conformation candidates within three seconds.


1997 ◽  
Vol 52 (6) ◽  
pp. 749-756
Author(s):  
Zofia Plesnar ◽  
Stanisław Malanowski ◽  
Zenon Lotowski ◽  
Jacek W. Morzycki ◽  
Jadwiga Frelek ◽  
...  

Abstract The cryoscopic measurements show that title compounds are strongly associated in CHCl3 solutions. The association of the 20 R epimer is distinctly less pronounced than that of the 20 S epipmer. Self-association of the 20 S epimer leads to the formation of very large com­plexes. The 20 R epimer forms associates via water molecules. The dissimilarity may be ex­plained in terms of different accessibility of the lactam carbonyl groups in the two epimers for the association. It is proposed that the association process is controlled by the configura­tion at the carbon atom C(20) and conformation around the C(20)-C(22) bond. Populations of side chain conformations of both epimers were determined by means of proton nuclear magnetic resonance. It was found for the 20 R epimer that the t and the -g rotamers are almost equally populated, and the rotamer +g is excluded. For the 20 S epimer the +g rotamer predominates over the t one, and the -g rotamer is excluded. The NMR data are fully consistent with the results of the molecular modelling studies.


2002 ◽  
Vol 76 (9) ◽  
pp. 4456-4466 ◽  
Author(s):  
Jennifer A. Gruenke ◽  
R. Todd Armstrong ◽  
William W. Newcomb ◽  
Jay C. Brown ◽  
Judith M. White

ABSTRACT Influenza virus hemagglutinin undergoes a conformational change in which a loop-to-helix “spring-loaded” conformational change forms a coiled coil that positions the fusion peptide for interaction with the target bilayer. Previous work has shown that two proline mutations designed to disrupt this change disrupt fusion but did not determine the basis for the fusion defect. In this work, we made six additional mutants with single proline substitutions in the region that undergoes the spring-loaded conformational change and two additional mutants with double proline substitutions in this region. All double mutants were fusion inactive. We analyzed one double mutant, F63P/F70P, as an example. We observed that F63P/F70P undergoes key low-pH-induced conformational changes and binds tightly to target membranes. However, limited proteolysis and electron microscopy observations showed that the mutant forms a coiled coil that is only ∼50% the length of the wild type, suggesting that it is splayed in its N-terminal half. This work further supports the hypothesis that the spring-loaded conformational change is necessary for fusion. Our data also indicate that the spring-loaded conformational change has another role beyond presenting the fusion peptide to the target membrane.


2014 ◽  
Vol 143 (4) ◽  
pp. 449-464 ◽  
Author(s):  
Natascia Vedovato ◽  
David C. Gadsby

A single Na+/K+-ATPase pumps three Na+ outwards and two K+ inwards by alternately exposing ion-binding sites to opposite sides of the membrane in a conformational sequence coupled to pump autophosphorylation from ATP and auto-dephosphorylation. The larger flow of Na+ than K+ generates outward current across the cell membrane. Less well understood is the ability of Na+/K+ pumps to generate an inward current of protons. Originally noted in pumps deprived of external K+ and Na+ ions, as inward current at negative membrane potentials that becomes amplified when external pH is lowered, this proton current is generally viewed as an artifact of those unnatural conditions. We demonstrate here that this inward current also flows at physiological K+ and Na+ concentrations. We show that protons exploit ready reversibility of conformational changes associated with extracellular Na+ release from phosphorylated Na+/K+ pumps. Reversal of a subset of these transitions allows an extracellular proton to bind an acidic side chain and to be subsequently released to the cytoplasm. This back-step of phosphorylated Na+/K+ pumps that enables proton import is not required for completion of the 3 Na+/2 K+ transport cycle. However, the back-step occurs readily during Na+/K+ transport when external K+ ion binding and occlusion are delayed, and it occurs more frequently when lowered extracellular pH raises the probability of protonation of the externally accessible carboxylate side chain. The proton route passes through the Na+-selective binding site III and is distinct from the principal pathway traversed by the majority of transported Na+ and K+ ions that passes through binding site II. The inferred occurrence of Na+/K+ exchange and H+ import during the same conformational cycle of a single molecule identifies the Na+/K+ pump as a hybrid transporter. Whether Na+/K+ pump–mediated proton inflow may have any physiological or pathophysiological significance remains to be clarified.


Molecules ◽  
2018 ◽  
Vol 23 (9) ◽  
pp. 2217 ◽  
Author(s):  
Ryota Kirikoshi ◽  
Noriyoshi Manabe ◽  
Ohgi Takahashi

Spontaneous deamidation in the Asn-Gly-Arg (NGR) motif that yields an isoAsp-Gly-Arg (isoDGR) sequence has recently attracted considerable attention because of the possibility of application to dual tumor targeting. It is well known that Asn deamidation reactions in peptide chains occur via the five-membered ring succinimide intermediate. Recently, we computationally showed by the B3LYP density functional theory method, that inorganic phosphate and the Arg side chain can catalyze the NGR deamidation using a cyclic peptide, c[CH2CO–NGRC]–NH2. In this previous study, the tetrahedral intermediate of the succinimide formation was assumed to be readily protonated at the nitrogen originating from the Asn side chain by the solvent water before the release of an NH3 molecule. In the present study, we found a new mechanism for the decomposition of the tetrahedral intermediate that does not require the protonation by an external proton source. The computational method is the same as in the previous study. In the new mechanism, the release of an NH3 molecule occurs after a proton exchange between the peptide and the phosphate and conformational changes. The rate-determining step of the overall reaction course is the previously reported first step, i.e., the cyclization to form the tetrahedral intermediate.


1992 ◽  
Vol 285 (2) ◽  
pp. 419-425 ◽  
Author(s):  
U Christensen ◽  
L Mølgaard

The kinetics of a series of Glu-plasminogen ligand-binding processes were investigated at pH 7.8 and 25 degrees C (in 0.1 M-NaCl). The ligands include compounds analogous to C-terminal lysine residues and to normal lysine residues. Changes of the Glu-plasminogen protein fluorescence were measured in a stopped-flow instrument as a function of time after rapid mixing of Glu-plasminogen and ligand at various concentrations. Large positive fluorescence changes (approximately 10%) accompany the ligand-induced conformational changes of Glu-plasminogen resulting from binding at weak lysine-binding sites. Detailed studies of the concentration-dependencies of the equilibrium signals and the rate constants of the process induced by various ligands showed the conformational change to involve two sites in a concerted positive co-operative process with three steps: (i) binding of a ligand at a very weak lysine-binding site that preferentially, but not exclusively, binds C-terminal-type lysine ligands, (ii) the rate-determining actual-conformational-change step and (iii) binding of one more lysine ligand at a second weak lysine-binding site that then binds the ligand more tightly. Further, totally independent initial small negative fluorescence changes (approximately 2-4%) corresponding to binding at the strong lysine-binding site of kringle 1 [Sottrup-Jensen, Claeys, Zajdel, Petersen & Magnusson (1978) Prog. Chem. Fibrinolysis Thrombolysis 3, 191-209] were observed for the C-terminal-type ligands. The finding that the conformational change in Glu-plasminogen involves two weak lysine-binding sites indicates that the effect cannot be assigned to any single kringle and that the problem of whether kringle 4 or kringle 5 is responsible for the process resolves itself. Probably kringle 4 and 5 are both participating. The involvement of two lysine binding-sites further makes the high specificity of Glu-plasminogen effectors more conceivable.


2000 ◽  
Vol 47 (1) ◽  
pp. 65-78 ◽  
Author(s):  
J Mazerski ◽  
K Muchewicz

Imidazoacridinones (IAs) are a new group of highly active antitumor compounds. The intercalation of the IA molecule into DNA is the preliminary step in the mode of action of these compounds. There are no experimental data about the structure of an intercalation complex formed by imidazoacridinones. Therefore the design of new potentially better compounds of this group should employ the molecular modelling techniques. The results of molecular dynamics simulations performed for four IA analogues are presented. Each of the compounds was studied in two systems: i) in water, and ii) in the intercalation complex with dodecamer duplex d(GCGCGCGCGCGC)2. Significant differences in the conformation of the side chain in the two environments were observed for all studied IAs. These changes were induced by electrostatic as well as van der Waals interactions between the intercalator and DNA. Moreover, the results showed that the geometry of the intercalation complex depends on: i) the chemical constitution of the side chain, and ii) the substituent in position 8 of the ring system.


Sign in / Sign up

Export Citation Format

Share Document