Elucidating the Membrane Binding Process of a Disordered Protein: Dynamic Interplay of Anionic Lipids and the Polybasic Region

Intrinsically disordered regions of proteins are responsible for many biological processes such as in the case of liver kinase B1 (LKB1)-a serine/threonine kinase relevant for cell proliferation and cell polarity. LKB1 becomes fully activated upon recruitment to the plasma membrane by binding of its disordered C-terminal polybasic motif consisting of eight lysines/arginines to phospholipids. Here, we present extensive molecular dynamics (MD) simulations of the polybasic motif interacting with a model membrane composed of 1-palmitoyl-2-oleoyl-sn-glycero-3-phosphocholine (POPC) and 1-palmitoyl-2-oleyl phosphatidic acid (PA) and cell culture experiments. Protein-membrane binding effects are due to the electrostatic interactions between the polybasic amino acids and PAs. For significant binding, the first three lysines turn out to be dispensable, which was also recapitulated in cell culture using transfected GFP-LKB1 variants. LKB1-membrane binding results in nonmonotonous changes in the structure of the protein as well as the membrane, in particular, accumulation of PAs and reduced thickness at the protein-membrane contact area. The protein-lipid binding turns out to be highly dynamic due to an interplay of PA-PA repulsion and protein-PA attraction. The thermodynamics of this interplay is captured by a statistical fluctuation model, which allows the estimation of both energies. Quantification of the significance of each polar amino acid in the polybasic provides detailed insights into the molecular mechanism of protein-membrane binding of LKB1. These results can likely be transferred to other proteins, which interact by intrinsically disordered polybasic regions with anionic membranes.

Conformational Selection of a Tryptophan Side Chain Drives the Generalized Increase in Activity of PET Hydrolases through a Ser/Ile Double Mutation

Poly(ethylene terephthalate) (PET) is the most common polyester plastic in the packaging industry and a major source of environmental pollution due to its single use. Several enzymes, termed PET hydrolases, have been found to hydrolyze this polymer at different temperatures, with the enzyme from Ideonella sakaiensis (IsPETase) having optimal catalytic activity at 30-35 °C. Crystal structures of IsPETase have revealed that the side chain of a conserved tryptophan residue within an active site loop (W185) shifts between three conformations to enable substrate binding and product release. This is facilitated by two residues unique to IsPETase, S214 and I218. When these residues are inserted into other PET hydrolases in place of the otherwise strictly conserved histidine and phenylalanine residues found at their respective positions, they enhance activity and decrease T _opt. Herein, we combine molecular dynamics and well-tempered metadynamics simulations to investigate dynamic changes of the S214/I218 and H214/F218 variants of IsPETase, as well as three other mesophilic and thermophilic PET hydrolases, at their respective temperature and pH optima. Our simulations show that the S214/I218 insertion both increases the flexibility of active site loop regions harboring key catalytic residues and the conserved tryptophan and expands the conformational plasticity of this tryptophan side chain, enabling the conformational transitions that allow for substrate binding and product release in IsPETase. The observed catalytic enhancement caused by this substitution in other PET hydrolases appears to be due to conformational selection, by capturing the conformational ensemble observed in IsPETase.

Increasing the O2 Resistance of the [FeFe]-Hydrogenase CbA5H through Enhanced Protein Flexibility

The high turnover rates of [FeFe]-hydrogenases under mild conditions and at low overpotentials provide a natural blueprint for the design of hydrogen catalysts. However, the unique active site (H-cluster) degrades upon contact with oxygen. The [FeFe]-hydrogenase fromClostridium beijerinckii (CbA5H) is characterized by the flexibility of its protein structure, which allows a conserved cysteine to coordinate to the active site under oxidative conditions. Thereby, intrinsic cofactor degradation induced by dioxygen is minimized. However, the protection from O₂ is only partial, and the activity of the enzyme decreases upon each exposure to O₂. By using site-directed mutagenesis in combination with electrochemistry, ATR-FTIR spectroscopy, and molecular dynamics simulations, we show that the kinetics of the conversion between the oxygen-protected inactive state (cysteine-bound) and the oxygen-sensitive active state can be accelerated by replacing a surface residue that is very distant from the active site. This sole exchange of methionine for a glutamate residue leads to an increased resistance of the hydrogenase to dioxygen. With our study, we aim to understand how local modifications of the protein structure can have a crucial impact on protein dynamics and how they can control the reactivity of inorganic active sites through outer sphere effects.

Short-Range Cooperative Slow-down of Water Solvation Dynamics Around SO4 2--Mg2+ Ion Pairs

The presence of ions affects the structure and dynamics of water on a multitude of length and time scales. In this context, pairs of Mg²⁺ and SO₄ ^2- ions in water constitute a prototypical system for which conflicting pictures of hydration geometries and dynamics have been reported. Key issues are the molecular pair and solvation shell geometries, the spatial range of electric interactions, and their impact on solvation dynamics. Here, we introduce asymmetric SO₄ ^2- stretching vibrations as new and most specific local probes of solvation dynamics that allow to access ion hydration dynamics at the dilute concentration (0.2 M) of a native electrolyte environment. Highly sensitive heterodyne 2D-IR spectroscopy in the fingerprint region of the SO₄ ^2- ions around 1100 cm^-1 reveals a specific slow-down of solvation dynamics for hydrated MgSO₄ and for Na₂SO₄ in the presence of Mg²⁺ ions, which manifests as a retardation of spectral diffusion compared to aqueous Na₂SO₄ solutions in the absence of Mg²⁺ ions. Extensive molecular dynamics and density functional theory QM/MM simulations provide a microscopic view of the observed ultrafast dephasing and hydration dynamics. They suggest a molecular picture where the slow-down of hydration dynamics arises from the structural peculiarities of solvent-shared SO₄ ^2--Mg²⁺ ion pairs.

Structure and Dynamic Properties of a Glycerol-Betaine Deep Eutectic Solvent: When Does a DES Become an Aqueous Solution?

Deep eutectic solvents (DESs) are an emerging class of green solvents with a wide spectrum of potential applications whose properties may be further tailored through the addition of water. Here, we study, through molecular dynamics, the influence of water on the properties of a betaine-glycerol-water (B:G:W) DES (1:2:ζ; ζ = 0 to 100), aiming at getting insight into the structural and dynamic crossover between a DES and an aqueous solution. The density, shear viscosity, and diffusion coefficients are found to exhibit a non-linear dependence of ζ, similar to that observed for the solvation layers' composition. Each Gly and Bet are replaced, respectively, by ∼3 and ∼5 water molecules, with the highest rates of depletion being found for Gly around Bet and Gly around Gly. Above ζ = 7 (70 mol %; 29.5 wt %), a major structural transformation occurs, with the complete disruption of the second Bet-Gly solvation layer and the formation of a new second layer at a shorter distance, accompanied by a sudden change in the rate of increase of the components' diffusion. Nonetheless, opposite to other DES, our results indicate a smooth crossover between a DES and an aqueous solution.

A Molecular Simulation Approach to Bond Reorganization in Epoxy Resins: From Curing to Deformation and Fracture

We model bond formation and dissociation processes in thermosetting polymer networks from molecular dynamics simulations. For this, a coarsened molecular mechanics model is derived from quantum calculations to provide effective interaction potentials that enable million-atoms scale simulations. The importance of bond (re)organization is demonstrated for (i) simulating epoxy resin formation-for which our approach leads to realistic network models which can now account for degrees of curing up to 98%. Moreover, (ii) we elucidate the competition of bond dissociation and bond reformation during plastic deformation and fracture. On this basis, we rationalize the molecular mechanisms that account for the irreversible nature of damaging epoxy polymers by mechanical load.

Application of Structure-based Methods to Analyze Resistance Mutations for Chemically Diverse Non-Nucleoside Reverse Transcriptase Inhibitors

BACKGROUND

Non-nucleoside reverse transcriptase inhibitors (NNRTIs) are used in combination with antiretroviral therapy to suppress viral loads in HIV patients. The chemical design of NNRTIs has changed in recent years in response to resistance-associated mutations (RAMs) and resistance. NNRTIs are chemically diverse compounds that bind an allosteric site of HIV RT. Resistance- associated mutations (RAMs) identified in HIV patients are associated with NNRTI resistance. RAMs confer amino acid changes that alter both structural and physiochemical properties of the allosteric site. Ultimately, these changes reduce NNRTI affinity. Previously, we used a combination of computational and experimental methods to analyze and validate RAMs for 3 diarylpyrimidine (DAPY) NNRTIs.

OBJECTIVE

The objective of this study is to apply these methods to other chemically diverse, non- DAPY NNRTIs.

MATERIALS AND METHODS

We selected MIV-150 (experimental microbicide) and doravirine for this study. A computational and molecular modeling strategy was used to evaluate the effects of RAMs. Calculated changes in drug affinity and stability (ΔS + ΔA) were used to determine overall resistance levels: susceptible, low, intermediate, and high. The ΔS + ΔA values for K101P suggest that this mutation confers intermediate/high-level resistance to MIV-150, but remains susceptible to doravirine. Based on the determined resistance levels, we analyzed the models and used Molecular Dynamics (MD) to compare the interactions of MIV-150/doravirine with RT wild-type (WT) and RT (K101P). From MD, we found that key interactions were lost with RT (K101P), but were retained with doravirine. To experimentally validate our findings, we conducted a fluorescence-based reverse transcription assay for MIV-150 with RT (WT) and RT (K101P). IC50 values determined in assays showed a 101-fold change in potency for MIV-150, but essentially no change for doravirine.

RESULTS

Our computational and experimental results are also consistent with antiviral data reported in the literature.

CONCLUSION

We believe that this approach is effective for analyzing mutations to determine resistance profiles for chemically diverse NNRTIs in development.

Simulated Protein Thermal Detection (SPTD) for Enzyme Thermostability Study and an Application Example for Pullulanase from Bacillus deramificans

BACKGROUND

The relationship between protein structure and its bioactivity is one of the fundamental problems for protein engineering and pharmaceutical design.

METHOD

A new method, called SPTD (Simulated Protein Thermal Detection), was proposed for studying and improving the thermal stability of enzymes. The method was based on the evidence observed by conducting the MD (Molecular Dynamics) simulation for all the atoms of an enzyme vibrating from the velocity at a room temperature (e.g., 25°C) to the desired working temperature (e.g., 65°C). According to the recorded MD trajectories and the coordinate deviations of the constituent residues under the two different temperatures, some new strategies have been found that are useful for both drug delivery and starch industry.

CONCLUSION

The SPTD technique presented in this paper may become a very useful tool for pharmaceutical design and protein engineering.

Atomic resolution monitoring of cation exchange in CdSe-PbSe heteronanocrystals during epitaxial solid-solid-vapor growth

Here, we show a novel solid-solid-vapor (SSV) growth mechanism whereby epitaxial growth of heterogeneous semiconductor nanowires takes place by evaporation-induced cation exchange. During heating of PbSe-CdSe nanodumbbells inside a transmission electron microscope (TEM), we observed that PbSe nanocrystals grew epitaxially at the expense of CdSe nanodomains driven by evaporation of Cd. Analysis of atomic-resolution TEM observations and detailed atomistic simulations reveals that the growth process is mediated by vacancies.