RESUMEN

Coarse-grained molecular dynamics simulations enable the modeling of increasingly complex systems at millisecond timescales. The transferable coarse-grained force field Martini 3 has shown great promise in modeling a wide range of biochemical processes, yet folded proteins in Martini 3 are not stable without the application of external bias potentials, such as elastic networks or Go̅-like models. We herein develop an algorithm, called OLIVES, which identifies native contacts with hydrogen bond capabilities in coarse-grained proteins and use it to implement a novel Go̅-like model for Martini 3. We show that the protein structure instability originates in part from the lack of hydrogen bond energy in the coarse-grained force field representation. By using realistic hydrogen bond energies obtained from literature ab initio calculations, it is demonstrated that protein stability can be recovered by the reintroduction of a coarse-grained hydrogen bond network and that OLIVES removes the need for secondary structure restraints. OLIVES is validated against known protein complexes and at the same time addresses the open question of whether there is a need for protein quaternary structure bias in Martini 3 simulations. It is shown that OLIVES can reduce the number of bias terms, hereby speeding up Martini 3 simulations of proteins by up to ≈30% on a GPU architecture compared to the established Go̅MARTINI Go̅-like model.

RESUMEN

Halogenated bisphenol A (BPA) derivatives are produced during disinfection treatment of drinking water or are synthesized as flame retardants (TCBPA or TBBPA). BPA is considered as an endocrine disruptor especially on human follicle-stimulating hormone receptor (FSHR). Using a global experimental approach, we assessed the effect of halogenated BPA derivatives on FSHR activity and estimated the risk of halogenated BPA derivatives to the reproductive health of exposed populations. For the first time, we show that FSHR binds halogenated BPA derivatives, at 10 nM, a concentration lower than those requires to modulate the activity of nuclear receptors and/or steroidogenesis enzymes. Indeed, bioluminescence assays show that FSHR response is lowered up to 42.36 % in the presence of BPA, up to 32.79 % by chlorinated BPA derivatives and up to 27.04 % by brominated BPA derivatives, at non-cytotoxic concentrations and without modification of basal receptor activity. Moreover, molecular docking, molecular dynamics simulations, and site-directed mutagenesis experiments demonstrate that the halogenated BPA derivatives bind the FSHR transmembrane domain reducing the signal transduction efficiency which lowers the cellular cAMP production and in fine disrupts the physiological effect of FSH. The potential reproductive health risk of exposed individuals was estimated by comparing urinary concentrations (through a collection of human biomonitoring data) with the lowest effective concentrations derived from in vitro cell assays. Our results suggest a potentially high concern for the risk of inhibition of the FSHR pathway. This global approach based on FSHR activity could enable the rapid characterization of the toxicity of halogenated BPA derivatives (or other compounds) and assess the associated risk of exposure to these halogenated BPA derivatives.

RESUMEN

Conventional bulk protein structure determination methods are not suitable for understanding the distinct and diverse interactions of proteins with interfaces. Notably, interfacial activation is a feature common to many lipases involving movement of a helical "lid" region upon contact with a hydrophobic surface to expose the catalytic site. Here we use the surface specificity of vibrational sum frequency generation spectroscopy (VSFG) spectroscopy to directly probe the conformation of Thermomyces lanuginosus lipase (TLL) at hydrophobic interfaces. The TLL-catalyzed reaction at the air/water interface is monitored by VSFG spectroscopy, showing loss of ester carbonyl modes and appearance of carboxylate stretching modes of the fatty acid products. Furthermore, comparison of experimental and calculated VSFG spectra of the amide I band of TLL allows us to discern the subtle structural changes involved with lid-opening at a hydrophobic surface. Finally, we report a likely orientation of this lid-open state, which interacts with the surface through a loop region away from the lid and active site. This experimental framework for probing protein structure and function at interfaces addresses a significant problem in protein science that is not only impeding the design of better enzymes for biotechnology applications but also drug discovery targeting membrane associated proteins.

Asunto(s)

RESUMEN

The regulation of dopamine (DA) removal from the synaptic cleft is a crucial process in neurotransmission and is facilitated by the sodium- and chloride-coupled dopamine transporter DAT. Psychostimulant drugs, cocaine, and amphetamine, both block the uptake of DA, while amphetamine also triggers the release of DA. As a result, they prolong or even amplify neurotransmitter signaling. Atypical inhibitors of DAT lack cocaine-like rewarding effects and offer a promising strategy for the treatment of drug use disorders. Here, we present the 3.2 Å resolution cryo-electron microscopy structure of the Drosophila melanogaster dopamine transporter (dDAT) in complex with the atypical non-competitive inhibitor AC-4-248. The inhibitor partially binds at the central binding site, extending into the extracellular vestibule, and locks the transporter in an outward open conformation. Our findings propose mechanisms for the non-competitive inhibition of DAT and attenuation of cocaine potency by AC-4-248 and provide a basis for the rational design of more efficacious atypical inhibitors.

RESUMEN

Vortioxetine (VTX) is a recently approved antidepressant that targets a variety of serotonin receptors. Here, we investigate the drug's molecular mechanism of operation at the serotonin 5-HT3 receptor (5-HT3R), which features two properties: VTX acts differently on rodent and human 5-HT3R, and VTX appears to suppress any subsequent response to agonists. Using a combination of cryo-EM, electrophysiology, voltage-clamp fluorometry and molecular dynamics, we show that VTX stabilizes a resting inhibited state of the mouse 5-HT3R and an agonist-bound-like state of human 5-HT3R, in line with the functional profile of the drug. We report four human 5-HT3R structures and show that the human receptor transmembrane domain is intrinsically fragile. We also explain the lack of recovery after VTX administration via a membrane partition mechanism.

Asunto(s)

RESUMEN

Antibiotic resistance of bacteria is considered one of the most alarming developments in modern medicine. While varied pathways for bacteria acquiring antibiotic resistance have been identified, there still are open questions concerning the mechanisms underlying resistance. Here, we show that alpha phenol-soluble modulins (PSMαs), functional bacterial amyloids secreted by Staphylococcus aureus, catalyze hydrolysis of ß-lactams, a prominent class of antibiotic compounds. Specifically, we show that PSMα2 and, particularly, PSMα3 catalyze hydrolysis of the amide-like bond of the four membered ß-lactam ring of nitrocefin, an antibiotic ß-lactam surrogate. Examination of the catalytic activities of several PSMα3 variants allowed mapping of the active sites on the amyloid fibrils' surface, specifically underscoring the key roles of the cross-α fibril organization, and the combined electrostatic and nucleophilic functions of the lysine arrays. Molecular dynamics simulations further illuminate the structural features of ß-lactam association upon the fibril surface. Complementary experimental data underscore the generality of the functional amyloid-mediated catalytic phenomenon, demonstrating hydrolysis of clinically employed ß-lactams by PSMα3 fibrils, and illustrating antibiotic degradation in actual S. aureus biofilms and live bacteria environments. Overall, this study unveils functional amyloids as catalytic agents inducing degradation of ß-lactam antibiotics, underlying possible antibiotic resistance mechanisms associated with bacterial biofilms.

Asunto(s)

RESUMEN

The amyloid aggregation of α-synuclein (αS), related to Parkinson's disease, can be catalyzed by lipid membranes. Despite the importance of lipid surfaces, the 3D-structure and orientation of lipid-bound αS is still not known in detail. Here, we report interface-specific vibrational sum-frequency generation (VSFG) experiments that reveal how monomeric αS binds to an anionic lipid interface over a large range of αS-lipid ratios. To interpret the experimental data, we present a frame-selection method ("ViscaSelect") in which out-of-equilibrium molecular dynamics simulations are used to generate structural hypotheses that are compared to experimental amide-I spectra via excitonic spectral calculations. At low and physiological αS concentrations, we derive flat-lying helical structures as previously reported. However, at elevated and potentially disease-related concentrations, a transition to interface-protruding αS structures occurs. Such an upright conformation promotes lateral interactions between αS monomers and may explain how lipid membranes catalyze the formation of αS amyloids at elevated protein concentrations.

Asunto(s)

RESUMEN

Sucrose import from photosynthetic tissues into the phloem is mediated by transporters from the low-affinity sucrose transporter family (SUC/SUT family). Furthermore, sucrose redistribution to other tissues is driven by phloem sap movement, the product of high turgor pressure created by this import activity. Additionally, sink organs such as fruits, cereals and seeds that accumulate high concentrations of sugar also depend on this active transport of sucrose. Here we present the structure of the sucrose-proton symporter, Arabidopsis thaliana SUC1, in an outward open conformation at 2.7 Å resolution, together with molecular dynamics simulations and biochemical characterization. We identify the key acidic residue required for proton-driven sucrose uptake and describe how protonation and sucrose binding are strongly coupled. Sucrose binding is a two-step process, with initial recognition mediated by the glucosyl moiety binding directly to the key acidic residue in a stringent pH-dependent manner. Our results explain how low-affinity sucrose transport is achieved in plants, and pinpoint a range of SUC binders that help define selectivity. Our data demonstrate a new mode for proton-driven symport with links to cation-driven symport and provide a broad model for general low-affinity transport in highly enriched substrate environments.

Asunto(s)

RESUMEN

Noncanonical DNA structures, termed G-quadruplexes, are present in human genomic DNA and are important elements in many DNA metabolic processes. Multiple sites in the human genome have G-rich DNA stretches able to support formation of several consecutive G-quadruplexes. One of those sites is the telomeric overhang region that has multiple repeats of TTAGGG and is tightly associated with both cancer and aging. We investigated the folding of consecutive G-quadruplexes in both potassium- and sodium-containing solutions using single-molecule FRET spectroscopy, circular dichroism, thermal melting and molecular dynamics simulations. Our observations show coexistence of partially and fully folded DNA, the latter consisting of consecutive G-quadruplexes. Following the folding process over hours in sodium-containing buffers revealed fast G-quadruplex folding but slow establishment of thermodynamic equilibrium. We find that full consecutive G-quadruplex formation is inhibited by the many DNA structures randomly nucleating on the DNA, some of which are off-path conformations that need to unfold to allow full folding. Our study allows describing consecutive G-quadruplex formation in both nonequilibrium and equilibrium conditions by a unified picture, where, due to the many possible DNA conformations, full folding with consecutive G-quadruplexes as beads on a string is not necessarily achieved.

Asunto(s)

RESUMEN

Amyloid fibrils are protein aggregates formed by protein assembly through cross ß structures. Inhibition of amyloid fibril formation may contribute to therapy against amyloid-related disorders like Parkinson's, Alzheimer's, and type 2 diabetes. Here we report that several fluorinated sulfonamide compounds, previously shown to inhibit human carbonic anhydrase, also inhibit the fibrillation of different proteins. Using a range of spectroscopic, microscopic and chromatographic techniques, we found that the two fluorinated sulfonamide compounds completely inhibit insulin fibrillation over a period of 16 h and moderately suppress α-synuclein and Aß fibrillation. In addition, these compounds decreased cell toxicity of insulin incubated under fibrillation-inducing conditions. We ascribe these effects to their ability to maintain insulin in the native monomeric state. Molecular dynamic simulations suggest that these compounds inhibit insulin self-association by interacting with residues at the dimer interface. This highlights the general anti-aggregative properties of aromatic sulfonamides and suggests that sulfonamide compounds which inhibit carbonic anhydrase activity may have potential as therapeutic agents against amyloid-related disorders.

Asunto(s)

RESUMEN

Independent force field validation is an essential practice to keep track of developments and for performing meaningful Molecular Dynamics simulations. In this work, atomistic force fields for intrinsically disordered proteins (IDP) are tested by simulating the archetypical IDP α-synuclein in solution for 2.5 µs. Four combinations of protein and water force fields were tested: ff19SB/OPC, ff19SB/TIP4P-D, ff03CMAP/TIP4P-D, and a99SB-disp/TIP4P-disp, with four independent repeat simulations for each combination. We compare our simulations to the results of a 73 µs simulation using the a99SB-disp/TIP4P-disp combination, provided by D. E. Shaw Research. From the trajectories, we predict a range of experimental observations of α-synuclein and compare them to literature data. This includes protein radius of gyration and hydration, intramolecular distances, NMR chemical shifts, and 3 J-couplings. Both ff19SB/TIP4P-D and a99SB-disp/TIP4P-disp produce extended conformational ensembles of α-synuclein that agree well with experimental radius of gyration and intramolecular distances while a99SB-disp/TIP4P-disp reproduces a balanced α-synuclein secondary structure content. It was found that ff19SB/OPC and ff03CMAP/TIP4P-D produce overly compact conformational ensembles and show discrepancies in the secondary structure content compared to the experimental data.

Asunto(s)

RESUMEN

Phosphatidylinositol 4,5-bisphosphate (PI(4,5)P2) lipids have been shown to stabilize an active conformation of class A G-protein coupled receptors (GPCRs) through a conserved binding site, not present in class B GPCRs. For class B GPCRs, previous molecular dynamics (MD) simulation studies have shown PI(4,5)P2 interacting with the Glucagon receptor (GCGR), which constitutes an important target for diabetes and obesity therapeutics. In this work, we applied MD simulations supported by native mass spectrometry (nMS) to study lipid interactions with GCGR. We demonstrate how tail composition plays a role in modulating the binding of PI(4,5)P2 lipids to GCGR. Specifically, we find the PI(4,5)P2 lipids to have a higher affinity toward the inactive conformation of GCGR. Interestingly, we find that in contrast to class A GPCRs, PI(4,5)P2 appear to stabilize the inactive conformation of GCGR through a binding site conserved across class B GPCRs but absent in class A GPCRs. This suggests differences in the regulatory function of PI(4,5)P2 between class A and class B GPCRs.

Asunto(s)

RESUMEN

G-protein coupled receptors (GPCRs) are important pharmacological targets. Despite substantial progress, important questions still remain concerning the details of activation: how can a ligand act as an agonist in one receptor but as an antagonist in a homologous receptor, and how can agonists activate a receptor despite lacking polar functional groups able to interact with helix 5 as is the case for the related adrenergic receptors? Studying vortioxetine (VXT), an important multimodal antidepressant drug, may elucidate both questions. Herein, we present a thorough in silico analysis of VXT binding to 5-HT1A, 5-HT1B, and 5-HT7 receptors and compare it with available experimental data. We are able to rationalize the differential mode of action of VXT at different receptors, but also, in the case of the 5-HT1A receptor, we observe the initial steps of activation that inform about an activation mechanism that does not involve polar interaction with helix 5. The results extend our current understanding of agonist and antagonist action at aminergic GPCRs.

Asunto(s)

RESUMEN

Sugars are essential sources of energy and carbon and also function as key signalling molecules in plants. Sugar transport proteins (STP) are proton-coupled symporters responsible for uptake of glucose from the apoplast into plant cells. They are integral to organ development in symplastically isolated tissues such as seed, pollen and fruit. Additionally, STPs play a vital role in plant responses to stressors such as dehydration and prevalent fungal infections like rust and mildew. Here we present a structure of Arabidopsis thaliana STP10 in the inward-open conformation at 2.6 Å resolution and a structure of the outward-occluded conformation at improved 1.8 Å resolution, both with glucose and protons bound. The two structures describe key states in the STP transport cycle. Together with molecular dynamics simulations that establish protonation states and biochemical analysis, they pinpoint structural elements, conserved in all STPs, that clarify the basis of proton-to-glucose coupling. These results advance our understanding of monosaccharide uptake, which is essential for plant organ development, and set the stage for bioengineering strategies in crops.

Asunto(s)

RESUMEN

C99, a naturally occurring peptide, is a precursor of the amyloid ß-peptide (Aß) and plays an important role in the so-called amyloidogenic pathway of degradation of amyloid precursor protein. While the effect of C99's dimerization is not clearly determined, it has been hypothesized that the dimerization protects C99 from being cleaved further. Cholesterol (CHOL) is known to interact with C99 and its presence in high concentrations has been linked to an increase in the production of Aß; however, to what extent this is correlated, and how, has not yet been determined. In this study, we systematically examine the effect of increasing cholesterol concentration on the homodimerization propensity of C99, combining unbiased atomistic molecular dynamics simulations with biased simulations using a coarse grained resolution. Through the use of umbrella sampling, we show how the presence of high levels of CHOL destabilizes the interaction between two C99 monomers. The interaction pattern between the two C99s has shifted several residues, from the N-terminal end of the transmembrane region toward the corresponding C-terminal in the presence of CHOL. The umbrella sampling shows that the presence of high levels of CHOL led to a decrease of the disassociation energy by approximately 3 kJ/mol. In conclusion, this suggests that increasing CHOL destabilizes the interaction between the two C99 monomers, which may possibly cause an increase in the production of Aß42.

Asunto(s)

RESUMEN

Islet amyloid polypeptide (IAPP) is a proposed cause of the decreased beta-cell mass in patients with type-II diabetes. The molecular composition of the cell-membrane is important for regulating IAPP cytotoxicity and aggregation. Cholesterol is present at high concentrations in the pancreatic beta-cells, and in-vitro experiments have indicated that it affects the amyloid formation of IAPP either by direct interactions or by changing the properties of the membrane. In this study we apply atomistic, unbiased molecular dynamics simulations at a microsecond timescale to investigate the effect of cholesterol on membrane bound IAPP. Simulations were performed with various combinations of cholesterol, phosphatidylcholine (PC) and phosphatidylserine (PS) lipids. In all simulations, the helical structure of monomer IAPP was stabilized by the membrane. We found that cholesterol decreased the insertion depth of IAPP compared to pure phospholipid membranes, while PS lipids counteract the effect of cholesterol. The aggregation propensity has previously been proposed to correlate with the insertion depth of IAPP, which we found to decrease with the increased ordering of the lipids induced by cholesterol. Cholesterol is depleted in the vicinity of IAPP, and thus our results suggest that the effect of cholesterol is indirect.

RESUMEN

Nanodisc technology is increasingly being applied for structural and biophysical studies of membrane proteins. In this work, we present a general protocol for constructing molecular models of nanodiscs for molecular dynamics simulations. The protocol is written in python and based on geometric equations, making it fast and easy to modify, enabling automation and customization of nanodiscs in silico. The novelty being the ability to construct any membrane scaffold protein (MSP) variant fast and easy given only an input sequence. We validated and tested the protocol by simulating seven different nanodiscs of various sizes and with different membrane scaffold proteins, both circularized and noncircularized. The structural and biophysical properties were analyzed and shown to be in good agreement with previously reported experimental data and simulation studies.

Asunto(s)

RESUMEN

The human dopamine transporter (hDAT) is one in three members of the monoamine transporter family (MAT). hDAT is essential for regulating the dopamine concentration in the synaptic cleft through dopamine reuptake into the presynaptic neuron; thereby controlling hDAT dopamine signaling. Dysfunction of the transporter is linked to several psychiatric disorders. hDAT and the other MATs have been shown to form oligomers in the plasma membrane, but only limited data exists on which dimeric and higher order oligomeric states are accessible and energetically favorable. In this work, we present several probable dimer conformations using computational coarse-grained self-assembly simulations and assess the relative stability of the different dimer conformations using umbrella sampling replica exchange molecular dynamics. Overall, the dimer conformations primarily involve TM9 and/or TM11 and/or TM12 at the interface. Furthermore, we show that a palmitoyl group (palm) attached to hDAT on TM12 modifies the free energy of separation for interfaces involving TM12, suggesting that S-palmitoylation may change the relative abundance of dimers involving TM12 in a biological context. Finally, a comparison of the identified interfaces of hDAT and palmitoylated hDAT to the human serotonin transporter interfaces and the leucine transporter interface, suggests similar dimer conformations across these protein family.

Asunto(s)

RESUMEN

Kinetic properties and crystal structures of the Na+,K+-ATPase in complex with cardiotonic steroids (CTS) revealed significant differences between CTS subfamilies (Laursen et al.). Thus, we found beneficial effects of K+ on bufadienolide binding, which strongly contrasted with the well-known antagonism between K+ and cardenolides. In order to understand this peculiarity of bufalin interactions, we used docking and molecular dynamics simulations of the complexes involving Na+,K+-ATPase, bufadienolides (bufalin, cinobufagin), and ions (K+, Na+, Mg2+). The results revealed that bufadienolide binding is affected by (i) electrostatic attraction of the lactone ring by a cation and (ii) the ability of a cation to stabilize and "shape" the site constituted by transmembrane helices of the α-subunit (αM1-6). The latter effect was due to varying coordination patterns involving amino acid residues from helix bundles αM1-4 and αM5-10. Substituents on the steroid core of a bufadienolide add to and modify the cation effects. The above rationale is fully consistent with the ion effects on the kinetics of Na+,K+-ATPase/bufadienolide interactions.

Asunto(s)