scholarly journals Forced unfolding of protein-inspired single-chain random heteropolymers

Author(s):  
Zexiang Han ◽  
Shayna Hilburg ◽  
Alfredo Alexander-Katz

Synthetic random heteropolymers (RHPs) with high chemical heterogeneity can self-assemble into single-chain nanoparticles that exhibit features reminiscent of natural proteins, such as topological polymorphism. Using all-atom molecular dynamics simulations, this work investigates the structure and single-chain mechanical unfolding of a library of four-component RHPs in water, studying the effects of sequence, composition, configuration, and molecular weight. Results show that compactified RHPs can have highly dynamic unfolding behaviors which are dominated by complex side-chain interactions and prove markedly different from their homopolymer counterparts. For a given sequence and conformation, an RHP’s backbone topology can strongly impact its unfolding response, hinting at the importance of topological design in the nanoscale mechanics of heteropolymers. In addition, we identify enthalpically-driven reconfiguration upon unfolding, observing a solvent-shielding protection mechanism similar to protein stabilization by PEGylation. This work provides the first computational evidence for the force-induced unfolding of protein-inspired multicomponent heteropolymers.

2020 ◽  
Author(s):  
Julian Keupp ◽  
Johannes P. Dürholt ◽  
Rochus Schmid

The prototypical pillared layer MOFs, formed by a square lattice of paddle-<br>wheel units and connected by dinitrogen pillars, can undergo a breathing phase<br>transition by a “wine-rack” type motion of the square lattice. We studied this not<br>yet fully understood behavior using an accurate first principles parameterized force<br>field (MOF-FF) for larger nanocrystallites on the example of Zn 2 (bdc) 2 (dabco) [bdc:<br>benzenedicarboxylate, dabco: (1,4-diazabicyclo[2.2.2]octane)] and found clear indi-<br>cations for an interface between a closed and an open pore phase traveling through<br>the system during the phase transformation [Adv. Theory Simul. 2019, 2, 11]. In<br>conventional simulations in small supercells this mechanism is prevented by periodic<br>boundary conditions (PBC), enforcing a synchronous transformation of the entire<br>crystal. Here, we extend this investigation to pillared layer MOFs with flexible<br>side-chains, attached to the linker. Such functionalized (fu-)MOFs are experimen-<br>tally known to have different properties with the side-chains acting as fixed guest<br>molecules. First, in order to extend the parameterization for such flexible groups,<br>1a new parametrization strategy for MOF-FF had to be developed, using a multi-<br>structure force based fit method. The resulting parametrization for a library of<br>fu-MOFs is then validated with respect to a set of reference systems and shows very<br>good accuracy. In the second step, a series of fu-MOFs with increasing side-chain<br>length is studied with respect to the influence of the side-chains on the breathing<br>behavior. For small supercells in PBC a systematic trend of the closed pore volume<br>with the chain length is observed. However, for a nanocrystallite model a distinct<br>interface between a closed and an open pore phase is visible only for the short chain<br>length, whereas for longer chains the interface broadens and a nearly concerted trans-<br>formation is observed. Only by molecular dynamics simulations using accurate force<br>fields such complex phenomena can be studied on a molecular level.


2020 ◽  
Author(s):  
Kenneth Lucas ◽  
George Barnes

We present the results of direct dynamics simulations and DFT calculations aimed at elucidating the effect of \textit{O}-sulfonation on the collision induced dissociation for serine. Towards this end, direct dynamics simulations of both serine and sulfoserine were performed at multiple collision energies and theoretical mass spectra obtained. Comparisons to experimental results are favorable for both systems. Peaks related to the sulfo group are identified and the reaction dynamics explored. In particular, three significant peaks (m\z 106, 88, and 81) seen in the theoretical mass spectrum directly related to the sulfo group are analyzed as well as major peaks shared by both systems. Our analysis shows that the m\z 106 peaks result from intramolecular rearrangements, intermolecular proton transfer among complexes composed of initial fragmentation products, and at high energy side-chain fragmentation. The \mz 88 peak was found to contain multiple constitutional isomers, including a previously unconsidered, low energy structure. It was also seen that the RM1 semi empirical method was not able to obtain all of the major peaks seen in experiment for sulfoserine. In contrast, PM6 did obtain all major experimental peaks.


2020 ◽  
Vol 16 (4) ◽  
pp. 451-459 ◽  
Author(s):  
Fortunatus C. Ezebuo ◽  
Ikemefuna C. Uzochukwu

Background: Sulfotransferase family comprises key enzymes involved in drug metabolism. Oxamniquine is a pro-drug converted into its active form by schistosomal sulfotransferase. The conformational dynamics of side-chain amino acid residues at the binding site of schistosomal sulfotransferase towards activation of oxamniquine has not received attention. Objective: The study investigated the conformational dynamics of binding site residues in free and oxamniquine bound schistosomal sulfotransferase systems and their contribution to the mechanism of oxamniquine activation by schistosomal sulfotransferase using molecular dynamics simulations and binding energy calculations. Methods: Schistosomal sulfotransferase was obtained from Protein Data Bank and both the free and oxamniquine bound forms were subjected to molecular dynamics simulations using GROMACS-4.5.5 after modeling it’s missing amino acid residues with SWISS-MODEL. Amino acid residues at its binding site for oxamniquine was determined and used for Principal Component Analysis and calculations of side-chain dihedrals. In addition, binding energy of the oxamniquine bound system was calculated using g_MMPBSA. Results: The results showed that binding site amino acid residues in free and oxamniquine bound sulfotransferase sampled different conformational space involving several rotameric states. Importantly, Phe45, Ile145 and Leu241 generated newly induced conformations, whereas Phe41 exhibited shift in equilibrium of its conformational distribution. In addition, the result showed binding energy of -130.091 ± 8.800 KJ/mol and Phe45 contributed -9.8576 KJ/mol. Conclusion: The results showed that schistosomal sulfotransferase binds oxamniquine by relying on hybrid mechanism of induced fit and conformational selection models. The findings offer new insight into sulfotransferase engineering and design of new drugs that target sulfotransferase.


2021 ◽  
Vol 11 (1) ◽  
Author(s):  
Sara Bitam ◽  
Ahmad Elbahnsi ◽  
Geordie Creste ◽  
Iwona Pranke ◽  
Benoit Chevalier ◽  
...  

AbstractC407 is a compound that corrects the Cystic Fibrosis Transmembrane Conductance Regulator (CFTR) protein carrying the p.Phe508del (F508del) mutation. We investigated the corrector effect of c407 and its derivatives on F508del-CFTR protein. Molecular docking and dynamics simulations combined with site-directed mutagenesis suggested that c407 stabilizes the F508del-Nucleotide Binding Domain 1 (NBD1) during the co-translational folding process by occupying the position of the p.Phe1068 side chain located at the fourth intracellular loop (ICL4). After CFTR domains assembly, c407 occupies the position of the missing p.Phe508 side chain. C407 alone or in combination with the F508del-CFTR corrector VX-809, increased CFTR activity in cell lines but not in primary respiratory cells carrying the F508del mutation. A structure-based approach resulted in the synthesis of an extended c407 analog G1, designed to improve the interaction with ICL4. G1 significantly increased CFTR activity and response to VX-809 in primary nasal cells of F508del homozygous patients. Our data demonstrate that in-silico optimized c407 derivative G1 acts by a mechanism different from the reference VX-809 corrector and provide insights into its possible molecular mode of action. These results pave the way for novel strategies aiming to optimize the flawed ICL4–NBD1 interface.


2015 ◽  
Vol 1 (7) ◽  
pp. e1500263 ◽  
Author(s):  
Akihiko Nakamura ◽  
Takuya Ishida ◽  
Katsuhiro Kusaka ◽  
Taro Yamada ◽  
Shinya Fushinobu ◽  
...  

Hydrolysis of carbohydrates is a major bioreaction in nature, catalyzed by glycoside hydrolases (GHs). We used neutron diffraction and high-resolution x-ray diffraction analyses to investigate the hydrogen bond network in inverting cellulase PcCel45A, which is an endoglucanase belonging to subfamily C of GH family 45, isolated from the basidiomycete Phanerochaete chrysosporium. Examination of the enzyme and enzyme-ligand structures indicates a key role of multiple tautomerizations of asparagine residues and peptide bonds, which are finally connected to the other catalytic residue via typical side-chain hydrogen bonds, in forming the “Newton’s cradle”–like proton relay pathway of the catalytic cycle. Amide–imidic acid tautomerization of asparagine has not been taken into account in recent molecular dynamics simulations of not only cellulases but also general enzyme catalysis, and it may be necessary to reconsider our interpretation of many enzymatic reactions.


2000 ◽  
Vol 47 (1) ◽  
pp. 47-57 ◽  
Author(s):  
J Mazerski ◽  
I Antonini ◽  
S Martelli

Pyrimidoacridinetriones (PATs) are a new group of highly active antitumor compounds. It seems reasonable to assume that, like for some other acridine derivatives, intercalation into DNA is a necessary, however not a sufficient condition for antitumor activity of these compounds. Rational design of new compounds of this chemotype requires knowledge about the structure of the intercalation complex, as well as about interactions responsible for its stability. Computer simulation techniques such as molecular dynamics (MD) may provide valuable information about these problems. The results of MD simulations performed for three rationally selected PATs are presented in this paper. The compounds differ in the number and position of side chains. Each of the compounds was simulated in two systems: i) in water, and ii) in the intercalation complex with the dodecamer duplex d(GCGCGCGCGCGC)2. The orientation of the side chain in relation to the ring system is determined by the position of its attachment. Orientation of the ring system inside the intercalation cavity depends on the number and position of side chain(s). The conformations of the side chain(s) of all PATs studied in the intercalation complex were found to be very similar to those observed in water.


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