Search

The long-term goal of theoretically characterizing kinetics of long-range proton transfer in biomolecular pumps such as cytochrome c oxidase, which might exhibit interesting kinetic properties as required by their biological function are presented. The results suggest that coordinates based on the modified center of excess charge can be used as the reaction coordinate for computing a meaningful potential of mean force for long-range proton transfer.

ABCG2 is a multidrug transporter that protects tissues from xenobiotics, affects drug pharmacokinetics, and contributes to multidrug resistance of cancer cells. Here, we present tetracyclic fumitremorgin C analog Ko143 derivatives, evaluate their in vitro modulation of purified ABCG2, and report four high-resolution cryo-EM structures and computational analyses to elucidate their interactions with ABCG2. We found that Ko143 derivatives that are based on a ring-opened scaffold no longer inhibit ABCG2-mediated transport activity. In contrast, closed-ring, tetracyclic analogs were highly potent inhibitors. Strikingly, the least potent of these compounds, MZ82, bound deeper into the central ABCG2 cavity than the other inhibitors and it led to partial closure of the transmembrane domains and increased flexibility of the nucleotide-binding domains. Minor structural modifications can thus convert a potent inhibitor into a compound that induces conformational changes in ABCG2 similar to those observed during binding of a substrate. Molecular dynamics simulations and free energy binding calculations further supported the correlation between reduced potency and distinct binding pose of the compounds. We introduce the highly potent inhibitor AZ99 that may exhibit improved in vivo stability.

A detailed understanding of ligand-protein interaction is essential for developing rational drug-design strategies. In recent years, technological advances in cryo-electron microscopy (cryo-EM) brought a new era to the structural determination of biological macromolecules and assemblies at high resolution, marking cryo-EM as a promising tool for studying ligand-protein interactions. However, even in high-resolution cryo-EM results, the densities for the bound small-molecule ligands are often of lower quality due to their relatively dynamic and flexible nature, frustrating their accurate coordinate assignment. To address the challenge of ligand modeling in cryo-EM maps, here we report the development of GOLEM (Genetic Optimization of Ligands in Experimental Maps), an automated and robust ligand docking method that predicts a ligand's pose and conformation in cryo-EM maps. GOLEM employs a Lamarckian genetic algorithm to perform a hybrid global/local search for exploring the ligand's conformational, orientational, and positional space. As an important feature, GOLEM explicitly considers water molecules and places them at optimal positions and orientations. GOLEM takes into account both molecular energetics and the correlation with the cryo-EM maps in its scoring function to optimally place the ligand. We have validated GOLEM against multiple cryo-EM structures with a wide range of map resolutions and ligand types, returning ligand poses in excellent agreement with the densities. As a VMD plugin, GOLEM is free of charge and accessible to the community. With these features, GOLEM will provide a valuable tool for ligand modeling in cryo-EM efforts toward drug discovery.

Biological membranes play key roles in cellular compartmentalization, structure, and its signaling pathways. At varying temperatures, individual membrane lipids sample from different configurations, a process that frequently leads to higher-order phase behavior and phenomena. Here, we present a persistent homology (PH)-based method for quantifying the structural features of individual and bulk lipids, providing local and contextual information on lipid tail organization. Our method leverages the mathematical machinery of algebraic topology and machine learning to infer temperature-dependent structural information on lipids from static coordinates. To train our model, we generated multiple molecular dynamics trajectories of dipalmitoyl-phosphatidylcholine membranes at varying temperatures. A fingerprint was then constructed for each set of lipid coordinates by PH filtration, in which interaction spheres were grown around the lipid atoms while tracking their intersections. The sphere filtration formed a simplicial complex that captures enduring key topological features of the configuration landscape using homology, yielding persistence data . Following fingerprint extraction for physiologically relevant temperatures, the persistence data were used to train an attention-based neural network for assignment of effective temperature values to selected membrane regions. Our persistence homology-based method captures the local structural effects, via effective temperature, of lipids adjacent to other membrane constituents, e.g., sterols and proteins. This topological learning approach can predict lipid effective temperatures from static coordinates across multiple spatial resolutions. The tool, called MembTDA, can be accessed at https://github.com/hyunp2/Memb-TDA.

Cholesterol (CHL) plays an integral role in modulating the function and activity of various mammalian membrane proteins. Due to the slow dynamics of lipids, conventional computational studies of protein-CHL interactions rely on either long-time scale atomistic simulations or coarse-grained approximations to sample the process. A highly mobile membrane mimetic (HMMM) has been developed to enhance lipid diffusion and thus used to facilitate the investigation of lipid interactions with peripheral membrane proteins and, with customized in silico solvents to replace phospholipid tails, with integral membrane proteins. Here, we report an updated HMMM model that is able to include CHL, a nonphospholipid component of the membrane, henceforth called HMMM-CHL. To this end, we had to optimize the effect of the customized solvents on CHL behavior in the membrane. Furthermore, the new solvent is compatible with simulations using force-based switching protocols. In the HMMM-CHL, both improved CHL dynamics and accelerated lipid diffusion are integrated. To test the updated model, we have applied it to the characterization of protein-CHL interactions in two membrane protein systems, the human β 2 -adrenergic receptor (β 2 AR) and the mitochondrial voltage-dependent anion channel 1 (VDAC-1). Our HMMM-CHL simulations successfully identified CHL binding sites and captured detailed CHL interactions in excellent consistency with experimental data as well as other simulation results, indicating the utility of the improved model in applications where an enhanced sampling of protein-CHL interactions is desired.

Calculating the binding free energy of integral transmembrane (TM) proteins is crucial for understanding the mechanisms by which they recognize one another and reversibly associate. The glycophorin A (GpA) homodimer, composed of two α-helical segments, has long served as a model system for studying TM protein reversible association. The present work establishes a methodological framework for calculating the binding affinity of the GpA homodimer in the heterogeneous environment of a membrane. Our investigation carefully considered a variety of protocols, including the appropriate choice of the force field, rigorous standardization reflecting the experimental conditions, sampling algorithm, anisotropic environment, and collective variables, to accurately describe GpA dimerization via molecular dynamics-based approaches. Specifically, two strategies were explored: (i) an unrestrained potential mean force (PMF) calculation, which merely enhances sampling along the separation of the two binding partners without any restraint, and (ii) a so-called "geometrical route", whereby the α-helices are progressively separated with imposed restraints on their orientational, positional, and conformational degrees of freedom to accelerate convergence. Our simulations reveal that the simplified, unrestrained PMF approach is inadequate for the description of GpA dimerization. Instead, the geometrical route, tailored specifically to GpA in a membrane environment, yields excellent agreement with experimental data within a reasonable computational time. A dimerization free energy of -10.7 kcal/mol is obtained, in fairly good agreement with available experimental data. The geometrical route further helps elucidate how environmental forces drive association before helical interactions stabilize it. Our simulations also brought to light a distinct, long-lived spatial arrangement that potentially serves as an intermediate state during dimer formation. The methodological advances in the generalized geometrical route provide a powerful tool for accurate and efficient binding-affinity calculations of intricate TM protein complexes in inhomogeneous environments.

The increasing resistance of bacteria to commercially available antibiotics threatens patient safety in healthcare settings. Perturbation of ion homeostasis has emerged as a potential therapeutic strategy to fight against antibacterial resistance and other channelopathies. This study reports the development of 8-aminoquinoline (QN) derivatives and their transmembrane Zn 2+ transport activities. Our findings showed that a potent QN-based Zn 2+ transporter exhibits promising antibacterial properties against Gram-positive bacteria with reduced hemolytic activity and cytotoxicity to mammalian cells. Furthermore, this combination showed excellent in vivo efficacy against Staphylococcus aureus . Interestingly, this combination prevented bacterial resistance and restored susceptibility of gentamicin and methicillin-resistant S. aureus to commercially available β-lactam and other antibiotics that had lost their activity against the drug-resistant bacterial strain. Our findings suggest that the transmembrane transport of Zn 2+ by QN derivatives could be a promising strategy to combat bacterial infections and restore the activity of other antibiotics.

Modeling and simulation of small molecules such as drugs and biological cofactors have been both a major focus of computational chemistry for decades and a growing need among computational biophysicists who seek to investigate the interaction of different types of ligands with biomolecules. Of particular interest in this regard are quantum mechanical (QM) calculations that are used to more accurately describe such small molecules, which can be of heterogeneous structures and chemistry, either in purely QM calculations or in hybrid QM/molecular mechanics (MM) simulations. QM programs are also used to develop MM force field parameters for small molecules to be used along with established force fields for biomolecules in classical simulations. With this growing need in mind, here we report a set of software tools developed and closely integrated within the broadly used molecular visualization/analysis program, VMD, that allow the user to construct, modify, and parametrize small molecules and prepare them for QM, hybrid QM/MM, or classical simulations. The tools also provide interactive analysis and visualization capabilities in an easy-to-use and integrated environment. In this paper, we briefly report on these tools and their major features and capabilities, along with examples of how they can facilitate molecular research in computational biophysics that might be otherwise prohibitively complex.

Factor VII (FVII) is a serine protease with a key role in initiating the coagulation cascade. It is part of a family of vitamin K-dependent clotting proteins, which require vitamin K for formation of their specialized membrane-binding domains (Gla domains). Membrane binding of the FVII Gla domain is critical to the activity of FVII, mediating the formation of its complex with other clotting factors. While Gla domains among coagulation factors are highly conserved in terms of amino acid sequence and structure, they demonstrate differential binding specificity toward anionic lipids. Although most Gla domain-containing clotting proteins display a strong preference for phosphatidylserine (PS), it has been demonstrated that FVII and protein C instead bind preferentially to phosphatidic acid (PA). We have developed the first model of the FVII Gla domain bound to PA lipids in membranes containing PA, the highly mobile membrane mimetic model, which accelerates slow diffusion of lipids in molecular dynamics simulations and therefore facilitates the membrane binding process and enhances sampling of lipid interactions. Simulations were performed using atomic level molecular dynamics, requiring a fixed charge to all atoms. The overall charge assigned to each PA lipid for this study was -1. We also developed an additional model of the FVII Gla domain bound to a 1:1 PS/PC membrane and compared the modes of binding of PS and PA lipids to FVII, allowing us to identify potential PA-specific binding sites.

py-MCMD, an open-source Python software, provides a robust workflow layer that manages communication of relevant system information between the simulation engines NAMD and GOMC and generates coherent thermodynamic properties and trajectories for analysis. To validate the workflow and highlight its capabilities, hybrid Monte Carlo/molecular dynamics (MC/MD) simulations are performed for SPC/E water in the isobaric-isothermal ( NPT ) and grand canonical (GC) ensembles as well as with Gibbs ensemble Monte Carlo (GEMC). The hybrid MC/MD approach shows close agreement with reference MC simulations and has a computational efficiency that is 2 to 136 times greater than traditional Monte Carlo simulations. MC/MD simulations performed for water in a graphene slit pore illustrate significant gains in sampling efficiency when the coupled-decoupled configurational-bias MC (CD-CBMC) algorithm is used compared with simulations using a single unbiased random trial position. Simulations using CD-CBMC reach equilibrium with 25 times fewer cycles than simulations using a single unbiased random trial position, with a small increase in computational cost. In a more challenging application, hybrid grand canonical Monte Carlo/molecular dynamics (GCMC/MD) simulations are used to hydrate a buried binding pocket in bovine pancreatic trypsin inhibitor. Water occupancies produced by GCMC/MD simulations are in close agreement with crystallographically identified positions, and GCMC/MD simulations have a computational efficiency that is 5 times better than MD simulations. py-MCMD is available on GitHub at https://github.com/GOMC-WSU/py-MCMD.

Excitatory amino acid transporters (EAATs) are glutamate transporters that belong to the solute carrier 1A (SLC1A) family. They couple glutamate transport to the cotransport of three sodium (Na + ) ions and one proton (H + ) and the counter-transport of one potassium (K + ) ion. In addition to this coupled transport, binding of cotransported species to EAATs activates a thermodynamically uncoupled chloride (Cl - ) conductance. Structures of SLC1A family members have revealed that these transporters use a twisting elevator mechanism of transport, where a mobile transport domain carries substrate and coupled ions across the membrane, while a static scaffold domain anchors the transporter in the membrane. We recently demonstrated that the uncoupled Cl - conductance is activated by the formation of an aqueous pore at the domain interface during the transport cycle in archaeal Glt Ph . However, a pathway for the uncoupled Cl - conductance has not been reported for the EAATs, and it is unclear if such a pathway is conserved. Here, we employ all-atom molecular dynamics (MD) simulations combined with enhanced sampling, free-energy calculations, and experimental mutagenesis to approximate large-scale conformational changes during the transport process and identified a Cl - -conducting conformation in human EAAT1 (hEAAT1). Sampling the large-scale structural transitions in hEAAT1 allowed us to capture an intermediate conformation formed during the transport cycle with a continuous aqueous pore at the domain interface. The free-energy calculations performed for the conduction of Cl - and Na + ions through the captured conformation highlight the presence of two hydrophobic gates that control low-barrier movement of Cl - through the aqueous pathway. Overall, our findings provide insights into the mechanism by which a human neurotransmitter transporter supports functional duality of active transport and passive Cl - permeation and confirm the commonality of this mechanism in different members of the SLC1A family.

The string method with swarms of trajectories (SMwST) is an algorithm that identifies a physically meaningful transition pathway─a one-dimensional curve, embedded within a high-dimensional space of selected collective variables. The SMwST algorithm leans on a series of short, unbiased molecular dynamics simulations spawned at different locations of the discretized path, from whence an average dynamic drift is determined to evolve the string toward an optimal pathway. However conceptually simple in both its theoretical formulation and practical implementation, the SMwST algorithm is computationally intensive and requires a careful choice of parameters for optimal cost-effectiveness in applications to challenging problems in chemistry and biology. In this contribution, the SMwST algorithm is presented in a self-contained manner, discussing with a critical eye its theoretical underpinnings, applicability, inherent limitations, and use in the context of path-following free-energy calculations and their possible extension to kinetics modeling. Through multiple simulations of a prototypical polypeptide, combining the search of the transition pathway and the computation of the potential of mean force along it, several practical aspects of the methodology are examined with the objective of optimizing the computational effort, yet without sacrificing accuracy. In light of the results reported here, we propose some general guidelines aimed at improving the efficiency and reliability of the computed pathways and free-energy profiles underlying the conformational transitions at hand.

Cytochrome c 2 (cyt. c 2 ) is a major element in electron transfer between redox proteins in bioenergetic membranes. While the interaction between cyt. c 2 and anionic lipids abundant in bioenergetic membranes has been reported, their effect on the shuttling activity of cyt. c 2 remains elusive. Here, the effect of anionic lipids on the interaction and binding of cyt. c 2 to the cytochrome bc 1 complex ( bc 1 ) is investigated using a combination of molecular dynamics (MD) and Brownian dynamics (BD) simulations. MD is used to generate thermally accessible conformations of cyt. c 2 and membrane-embedded bc 1 , which were subsequently used in multireplica BD simulations of diffusion of cyt. c 2 from solution to bc 1 , in the presence of various lipids. We show that, counterintuitively, anionic lipids facilitate association of cyt. c 2 with bc 1 by localizing its diffusion to the membrane surface. The observed lipid-mediated bc 1 association is further enhanced by the oxidized state of cyt. c 2 , in line with its physiological function. This lipid-mediated enhancement is salinity-dependent, and anionic lipids can disrupt cyt. c 2 - bc 1 interaction at nonphysiological salt levels. Our data highlight the importance of the redox state of cyt. c 2 , the lipid composition of the chromatophore membrane, and the salinity of the chromatophore in regulating the efficiency of the electron shuttling process mediated by cyt. c 2 . The conclusions can be extrapolated to mitochondrial systems and processes, or any bioenergetic membrane, given the structural similarity between cyt. c 2 and bc 1 and their mitochondrial counterparts.

Molecular dynamics (MD) simulations of biological membranes have achieved such levels of sophistication that are commonly used to predict unresolved structures and various properties of lipids and to substantiate experimental data. While achieving sufficient sampling of lipid dynamics remains a major challenge, a commonly used method to improve lipid sampling, e.g., in terms of specific interactions with membrane-associated proteins, is to randomize the initial arrangement of lipid constituents in multiple replicas of simulations, without changing the overall lipid composition of the membrane of interest. Here, we introduce a method that can rapidly generate multiple replicas of lipid bilayers with different spatial and conformational configurations for any given lipid composition. The underlying algorithm, which allows one to shuffle lipids at any desired level, relies on the application of an external potential, here referred to as the "carving potential", that removes clashes/entanglements before lipid positions are exchanged (shuffled), thereby minimizing the energy penalty due to abrupt lipid repositioning. The method is implemented as "Membrane Mixer Plugin (MMP) 1.0" in VMD, with a convenient graphical user interface that guides the user in setting various options and parameters. The plugin is fully automated and generates new membrane replicas more rapidly and conveniently than other analogous tools. The plugin and its capabilities introduced here can be extended to include additional features in future versions.

The march toward exascale computing will enable routine molecular simulation of larger and more complex systems, for example, simulation of entire viral particles, on the scale of approximately billions of atoms─a simulation size commensurate with a small bacterial cell. Anticipating the future hardware capabilities that will enable this type of research and paralleling advances in experimental structural biology, efforts are currently underway to develop software tools, procedures, and workflows for constructing cell-scale structures. Herein, we describe our efforts in developing and implementing an efficient and robust workflow for construction of cell-scale membrane envelopes and embedding membrane proteins into them. A new approach for construction of massive membrane structures that are stable during the simulations is built on implementing a subtractive assembly technique coupled with the development of a structure concatenation tool (fastmerge), which eliminates overlapping elements based on volumetric criteria rather than adding successive molecules to the simulation system. Using this approach, we have constructed two "protocells" consisting of MARTINI coarse-grained beads to represent cellular membranes, one the size of a cellular organelle and another the size of a small bacterial cell. The membrane envelopes constructed here remain whole during the molecular dynamics simulations performed and exhibit water flux only through specific proteins, demonstrating the success of our methodology in creating tight cell-like membrane compartments. Extended simulations of these cell-scale structures highlight the propensity for nonspecific interactions between adjacent membrane proteins leading to the formation of protein microclusters on the cell surface, an insight uniquely enabled by the scale of the simulations. We anticipate that the experiences and best practices presented here will form the basis for the next generation of cell-scale models, which will begin to address the addition of soluble proteins, nucleic acids, and small molecules essential to the function of a cell.

Cytochrome P450 2D6 (CYP2D6) is primarily expressed in the liver and in the central nervous system. It is known to be highly polymorphic in nature. It metabolizes several endogenous substrates such as anandamide (AEA). Concomitantly, it is involved in phase 1 metabolism of several antidepressants, antipsychotics, and other drugs. Research in the field of phytocannabinoids (pCBs) has recently accelerated owing to their legalization and increasing medicinal use for pain and inflammation. The primary component of cannabis is THC, which is well-known for its psychotropic effects. Since CYP2D6 is an important brain and liver P450 and is known to be inhibited by CBD, we investigated the interactions of four important highly prevalent CYP2D6 polymorphisms with selected phytocannabinoids (CBD, THC, CBDV, THCV, CBN, CBG, CBC, β-carophyllene) that are rapidly gaining popularity. We show that there is differential binding of CYP2D6*17 to pCBs as compared to WT CYP2D6. We also perform a more detailed comparison of WT and *17 CYP2D6, which reveals the possible regulation of AEA metabolism by CBD. Furthermore, we use molecular dynamics to delineate the mechanism of this binding, inhibition, and regulation. Taken together, we have found that the interactions of CYP2D6 with pCBs vary by polymorphism and by specific pCB class.

Alzheimer's Disease (AD) is the most common neurodegenerative disease, and efficient therapeutic and early diagnostic agents for AD are still lacking. Herein, we report the development of a novel amphiphilic compound, LS-4, generated by linking a hydrophobic amyloid-binding distyrylbenzene fragment with a hydrophilic triazamacrocycle, which dramatically increases the binding affinity toward various amyloid β (Aβ) peptide aggregates, especially for soluble Aβ oligomers. Moreover, upon the administration of LS-4 to 5xFAD mice, fluorescence imaging of LS-4-treated brain sections reveals that LS-4 can penetrate the blood-brain barrier and bind to the Aβ oligomers in vivo . In addition, the treatment of 5xFAD mice with LS-4 reduces the amount of both amyloid plaques and associated phosphorylated tau aggregates vs the vehicle-treated 5xFAD mice, while microglia activation is also reduced. Molecular dynamics simulations corroborate the observation that introducing a hydrophilic moiety into the molecular structure of LS-4 can enhance the electrostatic interactions with the polar residues of the Aβ species. Finally, exploiting the Cu 2+ -chelating property of the triazamacrocycle, we performed a series of imaging and biodistribution studies that show the 64 Cu-LS-4 complex binds to the amyloid plaques and can accumulate to a significantly larger extent in the 5xFAD mouse brains vs the wild-type controls. Overall, these results illustrate that the novel strategy, to employ an amphiphilic molecule containing a hydrophilic moiety attached to a hydrophobic amyloid-binding fragment, can increase the binding affinity for both soluble and insoluble Aβ aggregates and can thus be used to detect and regulate various Aβ species in AD.

Harnessing the power of graphics processing units (GPUs) to accelerate molecular dynamics (MD) simulations in the context of free-energy calculations has been a longstanding effort toward the development of versatile, high-performance MD engines. We report a new GPU-based implementation in NAMD of free-energy perturbation (FEP), one of the oldest, most popular importance-sampling approaches for the determination of free-energy differences that underlie alchemical transformations. Compared to the CPU implementation available since 2001 in NAMD, our benchmarks indicate that the new implementation of FEP in traditional GPU code is about four times faster, without any noticeable loss of accuracy, thereby paving the way toward more affordable free-energy calculations on large biological objects. Moreover, we have extended this new FEP implementation to a code path highly optimized for a single-GPU node, which proves to be up to nearly 30 times faster than the CPU implementation. Through optimized GPU performance, the present developments provide the community with a cost-effective solution for conducting FEP calculations. The new FEP-enabled code has been released with NAMD 3.0.

ApoB lipoproteins (apo B-Lp) are produced in hepatocytes, and their secretion requires the cargo receptor sortilin. We examined the secretion of apo B-Lp-containing very low-density lipoprotein (VLDL), an LDL progenitor. Sortilin also regulates the trafficking of the subtilase PCSK9, which when secreted binds the LDL receptor (LDLR), resulting in its endocytosis and destruction at the lysosome. We show that the site 2 binding compound (cpd984) has multiple effects in hepatocytes, including (1) enhanced Apo-Lp secretion, (2) increased cellular PCSK9 retention, and (3) augmented levels of LDLR at the plasma membrane. We postulate that cpd984 enhances apo B-Lp secretion in part through binding the lipid phosphatidylinositol 3,4,5-trisphosphate (PIP 3 ), which is present at higher levels on circulating VLDL form fed rats relative to after fasting. We attribute the enhanced VLDL secretion to its increased binding affinity for sortilin site 1 induced by cpd984 binding site 2. This hinders PCSK9 binding and secretion, which would subsequently prevent its binding to LDLR leading to its degradation. This suggests that site 2 is an allosteric regulator of site 1 binding. This effect is not limited to VLDL, as cpd984 augments binding of the neuropeptide neurotensin (NT) to sortilin site 1. Molecular dynamics simulations demonstrate that the C-terminus of NT (Ct-NT) stably binds site 1 through an electrostatic interaction. This was bolstered by the ability of Ct-NT to disrupt lower-affinity interactions between sortilin and the site 1 ligand PIP 3 . Together, these data show that binding cargo at sortilin site 1 is allosterically regulated through site 2 binding, with important ramifications for cellular lipid homeostasis involving proteins such as PCSK9 and LDLR.

We utilize various computational methodologies to study menthol's interaction with multiple organic phases, a lipid bilayer, and the human α4β2 nicotinic acetylcholine receptor (nAChR), the most abundant nAChR in the brain. First, force field parameters developed for menthol are validated in alchemical free energy perturbation simulations to calculate solvation free energies of menthol in water, dodecane, and octanol and compare the results against experimental data. Next, umbrella sampling is used to construct the free energy profile of menthol permeation across a 1-palmitoyl-2-oleoyl- sn -glycero-3-phosphocholine (POPC) bilayer. The results from a flooding simulation designed to study the water-membrane partitioning of menthol in a POPC lipid bilayer are used to determine the penetration depth and the preferred orientation of menthol in the bilayer. Finally, employing both docking and flooding simulations, menthol is shown to bind to different sites on the human α4β2 nAChR. The most likely binding mode of menthol to a desensitized membrane-embedded α4β2 nAChR is identified to be via a membrane-mediated pathway in which menthol binds to the sites at the lipid-protein interface after partitioning in the membrane. A rare but distinct binding mode in which menthol binds to the extracellular opening of receptor's ion permeation pore is also reported.

The assembly of proteins into amyloid fibrils has become linked not only with the progression of myriad human diseases, but also important biological functions. Understanding and controlling the formation, structure, and stability of amyloid fibrils is therefore a major scientific goal. Here we utilize electron microscopy-based approaches combined with quantitative statistical analysis to show how recently developed kind of amyloid modulators-multivalent polymer-peptide conjugates (mPPCs)-can be applied to control the structure and stability of amyloid fibrils. In doing so, we demonstrate that mPPCs are able to convert 40-residue amyloid beta fibrils into ordered nanostructures through a combination of fragmentation and bundling. Fragmentation is shown to be consistent with a model where the rate constant of fibril breakage is independent of the fibril length, suggesting a local and specific interaction between fibrils and mPPCs. Subsequent bundling, which was previously not observed, leads to the formation of sheet-like nanostructures which are surprisingly much more uniform than the starting fibrils. These nanostructures have dimensions independent of the molecular weight of the mPPC and retain the molecular-level ordering of the starting amyloid fibrils. Collectively, we reveal quantitative and nanoscopic understanding of how mPPCs can be applied to control amyloid structure and stability, and demonstrate approaches to elucidate nanoscale amyloid phase behavior in the presence of functional macromolecules and other modulators.

Protein aggregation is implicated in multiple deposition diseases including Alzheimer's Disease, which features the formation of toxic aggregates of amyloid beta (Aβ) peptides. Many inhibitors have been developed to impede or reverse Aβ aggregation. Multivalent inhibitors, however, have been largely overlooked despite the promise of high inhibition efficiency endowed by the multivalent nature of Aβ aggregates. In this work, we report the success of multivalent polymer-peptide conjugates (mPPCs) as a general class of inhibitors of the aggregation of Aβ 40 . Significantly delayed onset of fibril formation was realized using mPPCs prepared from three peptide/peptoid ligands covering a range of polymer molecular weights (MWs) and ligand loadings. Dose dependence studies showed that the nature of the ligands is a key factor in mPPC inhibition potency. The negatively charged ligand LPFFD (LD) leads to more efficient mPPCs compared to the neutral ligands, and is most effective at 7% ligand loading across different MWs. Molecular dynamics simulations along with dynamic light scattering experiments suggest that mPPCs form globular structures in solution due to ligand-ligand interactions. Such interactions are key to the spatial proximity of ligands and thus to the multivalency effect of mPPC inhibition. Excess ligand-ligand interactions, however, reduce the accessibility of mPPC ligands to Aβ peptides, and impair the overall inhibition potency.

Transcription factor STAT3 has been shown to regulate genes that are involved in stem cell self-renewal and thus represents a novel therapeutic target of great biological significance. However, many small-molecule agents with potential effects through STAT3 modulation in cancer therapy lack aqueous solubility and high off-target toxicity, hence impeding efficient bioavailability and activity. This work, for the first time, reports a prodrug-based strategy for selective and safer delivery of STAT3 inhibitors designed toward metastatic and drug-resistant breast cancer. We have synthesized a novel lipase-labile SN-2 phospholipid prodrug from a clinically investigated STAT3 inhibitor, nifuroxazide (Pro-nifuroxazide), which can be regioselectively cleaved by the membrane-abundant enzymes in cancer cells. Pro-nifuroxazide self-assembled to sub 20 nm nanoparticles (NPs), and the cytotoxic ability was screened in ER(+)-MCF-7 and ER(-)-MD-MB231 cells at 48-72 h using 3-(4,5-dimethylthiazol-2-yl)-2,5-diphenyltetra-zolium bromide proliferation assay. Results indicated that Pro-nifuroxazide NPs are multifold more effective toward inhibiting cancer cells in a time-dependent manner compared to parent nifuroxazide. A remarkable improvement in the local concentration of drugs to as high as ∼240 fold when assembled into NPs is presumably the reason for this functional improvement. We also introduced molecular dynamics simulations to generate Pro-nifuroxazide nano-assembly, as a model assembly from triggerable anti-cancer drugs, to provide molecular insights correlating physicochemical and anti-cancer properties. In silico properties of Pro-nifuroxazide including size, chemistry of NPs and membrane interactions with individual molecules could be validated by in vitro functional activities in cells of breast cancer origin. The in vivo anti-cancer efficiencies of Pro-nifuroxazide NPs in nude mice xenografts with MCF-7 revealed remarkable growth inhibition of as high as 400%. Histopathological analysis corroborated these findings to show significantly high nuclear fragmentation and retracted cytoplasm. Immunostaining on tumor section demonstrated a significantly lower level of pSTAT-3 by Pro-nifuroxazide NP treatment, establishing the inhibition of STAT-3 phosphorylation. Our strategy for the first time proposes a translatable prodrug agent self-assembled into NPs and demonstrates remarkable enhancement in IC 50 , induced apoptosis, and reduced cancer cell population through STAT-3 inhibition via reduced phosphorylation.

The cellular membrane constitutes one of the most fundamental compartments of a living cell, where key processes such as selective transport of material and exchange of information between the cell and its environment are mediated by proteins that are closely associated with the membrane. The heterogeneity of lipid composition of biological membranes and the effect of lipid molecules on the structure, dynamics, and function of membrane proteins are now widely recognized. Characterization of these functionally important lipid-protein interactions with experimental techniques is however still prohibitively challenging. Molecular dynamics (MD) simulations offer a powerful complementary approach with sufficient temporal and spatial resolutions to gain atomic-level structural information and energetics on lipid-protein interactions. In this review, we aim to provide a broad survey of MD simulations focusing on exploring lipid-protein interactions and characterizing lipid-modulated protein structure and dynamics that have been successful in providing novel insight into the mechanism of membrane protein function.

The role of Phe213 in the allosteric mechanism of human cytochrome P450 CYP3A4 was studied using a combination of progesterone (PGS) and carbamazepine (CBZ) as probe substrates. We expressed, purified, and incorporated into POPC Nanodiscs three mutants, F213A, F213S, and F213Y, and compared them with wild-type (WT) CYP3A4 by monitoring spectral titration, the rate of NADPH oxidation, and steady-state product turnover rates with pure substrates and substrate mixtures. All mutants demonstrated higher activity with CBZ, lower activity with PGS, and a reduced level of activation of CBZ epoxidation by PGS, which was most pronounced in the F213A mutant. Using all-atom molecular dynamics simulations, we compared the dynamics of WT CYP3A4 and the F213A mutant incorporated into the lipid bilayer and the effect of the presence of the PGS molecule at the allosteric peripheral site and evaluated the critical role of Phe213 in mediating the heterotropic allosteric interactions in CYP3A4.

The plasma membrane of the cell is a complex, tightly regulated, heterogeneous environment shaped by proteins, lipids, and small molecules. Ca 2+ ions are important cellular messengers, spatially separated from anionic lipids. After cell injury, disease, or apoptotic events, anionic lipids are externalized to the outer leaflet of the plasma membrane and encounter Ca 2+ , resulting in dramatic changes in the plasma membrane structure and initiation of signaling cascades. Despite the high chemical and biological significance, the structures of lipid-Ca 2+ nanoclusters are still not known. Previously, we demonstrated by solid-state nuclear magnetic resonance (NMR) spectroscopy that upon binding to Ca 2+ , individual phosphatidylserine lipids populate two distinct yet-to-be-characterized structural environments. Here, we concurrently employ extensive all-atom molecular dynamics (MD) simulations with our accelerated membrane mimetic and detailed NMR measurements to identify lipid-Ca 2+ nanocluster conformations. We find that major structural characteristics of these nanoclusters, including interlipid pair distances and chemical shifts, agree with observable NMR parameters. Simulations reveal that lipid-ion nanoclusters are shaped by two characteristic, long-lived lipid structures induced by divalent Ca 2+ . Using ab initio quantum mechanical calculations of chemical shifts on MD-captured lipid-ion complexes, we show that computationally observed conformations are validated by experimental NMR data. Both NMR measurements of diluted specifically labeled lipids and MD simulations reveal that the basic structural unit that reshapes the membrane is a Ca 2+ -coordinated phosphatidylserine tetramer. Our combined computational and experimental approach presented here can be applied to other complex systems in which charged membrane-active molecular agents leave structural signatures on lipids.

The human body contains endogenous cannabinoids (endocannabinoids) that elicit effects similar to those of Δ 9 -tetrahydrocanabinol, the principal bioactive component of cannabis. The endocannabinoid virodhamine (O-AEA) is the constitutional isomer of the well-characterized cardioprotective and anti-inflammatory endocannabinoid anandamide (AEA). The chemical structures of O-AEA and AEA contain arachidonic acid (AA) and ethanolamine; however, AA in O-AEA is connected to ethanolamine via an ester linkage, whereas AA in AEA is connected through an amide linkage. O-AEA is involved in regulating blood pressure and cardiovascular function. We show that O-AEA is found at levels 9.6-fold higher than that of AEA in porcine left ventricle. On a separate note, the cytochrome P450 (CYP) epoxygenase CYP2J2 is the most abundant CYP in the heart where it catalyzes the metabolism of AA and AA-derived eCBs to bioactive epoxides that are involved in diverse cardiovascular functions. Herein, using competitive binding studies, kinetic metabolism measurements, molecular dynamics, and wound healing assays, we have shown that O-AEA is an endogenous inhibitor of CYP2J2 epoxygenase. As a result, the role of O-AEA as an endogenous eCB inhibitor of CYP2J2 may provide a new mode of regulation to control the activity of cardiovascular CYP2J2 in vivo and suggests a potential cross-talk between the cardiovascular endocannabinoids and the cytochrome P450 system.

Carboxysomes are closed polyhedral cellular microcompartments that increase the efficiency of carbon fixation in autotrophic bacteria. Carboxysome shells consist of small proteins that form hexameric units with semipermeable central pores containing binding sites for anions. This feature is thought to selectively allow access to RuBisCO enzymes inside the carboxysome by HCO 3 - (the dominant form of CO 2 in the aqueous solution at pH 7.4) but not O 2 , which leads to a nonproductive reaction. To test this hypothesis, here we use molecular dynamics simulations to characterize the energetics and permeability of CO 2 , O 2 , and HCO 3 - through the central pores of two different shell proteins, namely, CsoS1A of α-carboxysome and CcmK4 of β-carboxysome shells. We find that the central pores are in fact selectively permeable to anions such as HCO 3 - , as predicted by the model.

Cytochrome aa 3 is the terminal respiratory enzyme of all eukaryotes and many bacteria and archaea, reducing O 2 to water and harnessing the free energy from the reaction to generate the transmembrane electrochemical potential. The diffusion of O 2 to the heme-copper catalytic site, which is buried deep inside the enzyme, is the initiation step of the reaction chemistry. Our previous molecular dynamics (MD) study with cytochrome ba 3 , a homologous enzyme of cytochrome aa 3 in Thermus thermophilus, demonstrated that O 2 diffuses from the lipid bilayer to its reduction site through a 25 Å long tunnel inferred by Xe binding sites detected by X-ray crystallography [Mahinthichaichan, P., Gennis, R., and Tajkhorshid, E. (2016) Biochemistry 55, 1265-1278]. Although a similar tunnel is observed in cytochrome aa 3 , this putative pathway appears partially occluded between the entrances and the reduction site. Also, the experimentally determined second-order rate constant for O 2 delivery in cytochrome aa 3 (∼10 8 M -1 s -1 ) is 10 times slower than that in cytochrome ba 3 (∼10 9 M -1 s -1 ). A question to be addressed is whether cytochrome aa 3 utilizes this X-ray-inferred tunnel as the primary pathway for O 2 delivery. Using complementary computational methods, including multiple independent flooding MD simulations and implicit ligand sampling calculations, we probe the O 2 delivery pathways in cytochrome aa 3 of Rhodobacter sphaeroides. All of the O 2 molecules that arrived in the reduction site during the simulations were found to diffuse through the X-ray-observed tunnel, despite its apparent constriction, supporting its role as the main O 2 delivery pathway in cytochrome aa 3 . The rate constant for O 2 delivery in cytochrome aa 3 , approximated using the simulation results, is 10 times slower than in cytochrome ba 3 , in agreement with the experimentally determined rate constants.

Heterotropic interactions between atorvastatin (ARVS) and dronedarone (DND) have been deciphered using global analysis of the results of binding and turnover experiments for pure drugs and their mixtures. The in vivo presence of atorvastatin lactone (ARVL) was explicitly taken into account by using pure ARVL in analogous experiments. Both ARVL and ARVS inhibit DND binding and metabolism, while a significantly higher affinity of CYP3A4 for ARVL makes the latter the main modulator of activity (effector) in this system. Molecular dynamics simulations reveal significantly different modes of interactions of DND and ARVL with the substrate binding pocket and with a peripheral allosteric site. Interactions of both substrates with residues F213 and F219 at the allosteric site play a critical role in the communication of conformational changes induced by effector binding to productive binding of the substrate at the catalytic site.

Doxorubicin (DOX) is a chemotherapeutic that is used in the treatment of a wide variety of cancers. However, it causes cardiotoxicity partly because of the formation of reactive oxygen species. CYP2J2 is a human cytochrome P450 that is strongly expressed in cardiomyocytes. It converts arachidonic acid (AA) into four different regioisomers of epoxyeicosatrienoic acids (EETs). Using kinetic analyses, we show that AA metabolism by CYP2J2 is modulated by DOX. We show that cytochrome P450 reductase, the redox partner of CYP2J2, metabolizes DOX to 7-deoxydoxorubicin aglycone (7-de-aDOX). This metabolite then binds to CYP2J2 and inhibits and alters the preferred site of metabolism of AA, leading to a change in the ratio of the EET regioisomers. Furthermore, molecular dynamics simulations indicate that 7-de-aDOX and AA can concurrently bind to the CYP2J2 active site to produce these changes in the site of AA metabolism. To determine if these observations are unique to DOX/7-de-aDOX, we use noncardiotoxic DOX analogues, zorubicin (ZRN) and 5-iminodaunorubicin (5-IDN). ZRN and 5-IDN inhibit CYP2J2-mediated AA metabolism but do not change the ratio of EET regioisomers. Altogether, we demonstrate that DOX and 7-de-aDOX inhibit CYP2J2-mediated AA metabolism and 7-de-aDOX binds close to the active site to alter the ratio of cardioprotective EETs. These mechanistic studies of CYP2J2 can aid in the design of new alternative DOX derivatives.

The mechanics of the protein-lipid interactions of transmembrane proteins are difficult to capture with conventional atomic molecular dynamics, due to the slow lateral diffusion of lipids restricting sampling to states near the initial membrane configuration. The highly mobile membrane mimetic (HMMM) model accelerates lipid dynamics by modeling the acyl tails nearest to the membrane center as a fluid organic solvent while maintaining an atomic description of the lipid headgroups and short acyl tails. The HMMM has been applied to many peripheral protein systems; however, the organic solvent used to date caused deformations in transmembrane proteins by intercalating into the protein and disrupting interactions between individual side chains. We ameliorate the effect of the solvent on transmembrane protein structure through the development of two new in silico Lennard-Jones solvents. The parameters for the new solvents were determined through an extensive parameter search in order to match the bulk properties of alkanes in a highly simplified model. Using these new solvents, we substantially improve the insertion free energy profiles of 10 protein side chain analogues across the entire bilayer. In addition, we reduce the intercalation of solvent into transmembrane systems, resulting in native-like transmembrane protein structures from five different topological classes within a HMMM bilayer. The parametrization of the solvents, in addition to their computed physical properties, is discussed. By combining high lipid lateral diffusion with intact transmembrane proteins, we foresee the developed solvents being useful to efficiently identify membrane composition inhomogeneities and lipid binding caused by the presence of membrane proteins.

Synaptotagmin (Syt) is a membrane-associated protein involved in vesicle fusion through the SNARE complex that is found throughout the human body in 17 different isoforms. These isoforms have two membrane-binding C2 domains, which sense Ca 2+ and thereby promote anionic membrane binding and lead to vesicle fusion. Through molecular dynamics simulations using the highly mobile membrane mimetic acclerated bilayer model, we have investigated how small protein sequence changes in the Ca 2+ -binding loops of the C2 domains may give rise to the experimentally determined difference in binding kinetics between Syt-1 and Syt-7 isoforms. Syt-7 C2 domains are found to form more close contacts with anionic phospholipid headgroups, particularly in loop 1, where an additional positive charge in Syt-7 draws the loop closer to the membrane and causes the anchoring residue F167 to insert deeper into the bilayer than the corresponding methionine in Syt-1 (M173). By performing additional replica exchange umbrella sampling calculations, we demonstrate that these additional contacts increase the energetic cost of unbinding the Syt-7 C2 domains from the bilayer, causing them to unbind more slowly than their counterparts in Syt-1.

Cytochrome P450 (CYP) 2J2 is the primary epoxygenase in the heart and is responsible for the epoxidation of arachidonic acid (AA), an ω-6 polyunsaturated fatty acid (PUFA), into anti-inflammatory epoxide metabolites. It also epoxidizes other PUFAs such as docosahexaenoic acid (DHA), linoleic acid (LA), and eicosapentaenoic acid (EPA). Herein, we have performed detailed thermodynamic and kinetic analyses to determine how DHA, LA, and EPA modulate the metabolism of AA by CYP2J2. We use the Nanodisc system to stabilize CYP2J2 and its redox partner, CYP reductase (CPR). We observe that DHA strongly inhibits CYP2J2-mediated AA metabolism, LA only moderately inhibits AA metabolism, and EPA exhibits insignificant inhibition. We also characterized the binding of these molecules using ebastine competitive binding assays and show that DHA binds significantly tighter to CYP2J2 than AA, EPA, or LA. Furthermore, we utilize a combined approach of molecular dynamics (MD) simulations and docking to predict key residues mediating the tight binding of DHA. We show that although all the tested fatty acids form similar contacts to the active site residues, the affinity of DHA for CYP2J2 is tighter because of the interaction of DHA with residues Arg-321, Thr-318, and Ser-493. To demonstrate the importance of these residues in binding, we mutated these residues to make two mutant variants, CYP2J2-T318A and CYP2J2-T318V/S493A. Both mutant variants showed weaker binding than the wild type (WT) to DHA and AA; DHA inhibition of AA was also mitigated in the mutants compared to the WT. Therefore, using a combined experimental and MD simulation approach, we establish that CYP2J2 inhibition of AA metabolism by DHA, EPA, and LA is asymmetric because of tighter binding of DHA to select residues in the active site.

The respiratory cytochrome bo 3 ubiquinol oxidase from Escherichia coli has a high-affinity ubiquinone binding site that stabilizes the one-electron reduced ubisemiquinone (SQ H ), which is a transient intermediate during the electron-mediated reduction of O 2 to water. It is known that SQ H is stabilized by two strong hydrogen bonds from R71 and D75 to ubiquinone carbonyl oxygen O1 and weak hydrogen bonds from H98 and Q101 to O4. In this work, SQ H was investigated with orientation-selective Q-band (∼34 GHz) pulsed 1 H electron-nuclear double resonance (ENDOR) spectroscopy on fully deuterated cytochrome (cyt) bo 3 in a H 2 O solvent so that only exchangeable protons contribute to the observed ENDOR spectra. Simulations of the experimental ENDOR spectra provided the principal values and directions of the hyperfine (hfi) tensors for the two strongly coupled H-bond protons (H1 and H2). For H1, the largest principal component of the proton anisotropic hfi tensor T z' = 11.8 MHz, whereas for H2, T z' = 8.6 MHz. Remarkably, the data show that the direction of the H1 H-bond is nearly perpendicular to the quinone plane (∼70° out of plane). The orientation of the second strong hydrogen bond, H2, is out of plane by ∼25°. Equilibrium molecular dynamics simulations on a membrane-embedded model of the cyt bo 3 Q H site show that these H-bond orientations are plausible but do not distinguish which H-bond, from R71 or D75, is nearly perpendicular to the quinone ring. Density functional theory calculations support the idea that the distances and geometries of the H-bonds to the ubiquinone carbonyl oxygens, along with the measured proton anisotropic hfi couplings, are most compatible with an anionic (deprotonated) ubisemiquinone.

Small diffusible redox proteins facilitate electron transfer in respiration and photosynthesis by alternately binding to their redox partners and integral membrane proteins and exchanging electrons. Diffusive search, recognition, binding, and unbinding of these proteins often amount to kinetic bottlenecks in cellular energy conversion, but despite the availability of structures and intense study, the physical mechanisms controlling redox partner interactions remain largely unknown. The present molecular dynamics study provides an all-atom description of the cytochrome c2-docked bc1 complex in Rhodobacter sphaeroides in terms of an ensemble of favorable docking conformations and reveals an intricate series of conformational changes that allow cytochrome c2 to recognize the bc1 complex and bind or unbind in a redox state-dependent manner. In particular, the role of electron transfer in triggering a molecular switch and in altering water-mediated interface mobility, thereby strengthening and weakening complex formation, is described. The results resolve long-standing discrepancies between structural and functional data.

Molecular dynamics (MD) simulation engines use a variety of different approaches for modeling molecular systems with force fields that govern their dynamics and describe their topology. These different approaches introduce incompatibilities between engines, and previously published software bridges the gaps between many popular MD packages, such as between CHARMM and AMBER or GROMACS and LAMMPS. While there are many structure building tools available that generate topologies and structures in CHARMM format, only recently have mechanisms been developed to convert their results into GROMACS input. We present an approach to convert CHARMM-formatted topology and parameters into a format suitable for simulation with GROMACS by expanding the functionality of TopoTools, a plugin integrated within the widely used molecular visualization and analysis software VMD. The conversion process was diligently tested on a comprehensive set of biological molecules in vacuo. The resulting comparison between energy terms shows that the translation performed was lossless as the energies were unchanged for identical starting configurations. By applying the conversion process to conventional benchmark systems that mimic typical modestly sized MD systems, we explore the effect of the implementation choices made in CHARMM, NAMD, and GROMACS. The newly available automatic conversion capability breaks down barriers between simulation tools and user communities and allows users to easily compare simulation programs and leverage their unique features without the tedium of constructing a topology twice.

Human islet amyloid polypeptide (hIAPP) is a 37-residue peptide hormone, which upon misfolding changes from the physiologically active monomer into pathological amyloid fibril aggregates in the pancreas of type 2 diabetes mellitus patients. During this process, the insulin-producing pancreatic β-cells are damaged; however, the underlying mechanism of this mode of cytotoxicity remains elusive. It is known that anionic lipids accelerate amyloid fibril formation, implicating the importance of the cellular membrane in the process, and that a pH close to the level in the β-cell secretory granules (pH 5.5) inhibits amyloid fibril formation. Using all-atom molecular dynamics simulations, we have investigated the membrane-associated monomer state of α-helical hIAPP, analyzed specific interactions of hIAPP with a mixed anionic-zwitterionic lipid membrane and examined the influence of pH on the structure and dynamics of hIAPP and its interaction with the membrane. We find that hIAPP primarily interacts with the membrane by forming favorable interactions between anionic lipids and the positively charged residues in the N-terminal part of the peptide. Rationalizing experimental findings, the simulations show that the N-terminal part of the peptide interacts with the membrane in the lipid headgroup region. At neutral pH, the C-terminal part of the peptide, which contains the residues that initiate fibril formation, displays a highly dynamic, unfolded state, which interacts with the membrane significantly less than the N-terminal part. Such an unfolded form can be proposed to contribute to the acceleration of fibril formation. At low pH, protonation of His18 mediates a stronger interaction of the C-terminal part with the membrane, resulting in the immobilization of the C-terminal part on the membrane surface that might constitute a mechanism by which low pH inhibits fibril formation.

Noncovalent second-shell interactions are important in controlling metal-binding affinity and activity in metalloenzymes, but fine-tuning these interactions in designed metalloenzymes has not been fully explored. As a result, most designed metalloenzymes have low metal-binding affinity and activity. Here we identified three mutations in the second coordination shell of an engineered Mn(II)-binding site in cytochrome c peroxidase (called MnCcP.1, containing Glu45, Glu37, and Glu181 ligands) that mimics the native manganese peroxidase (MnP), and explored their effects on both Mn(II)-binding affinity and MnP activity. First, removing a hydrogen bond to Glu45 through Tyr36Phe mutation enhanced Mn(II)-binding affinity, as evidenced by a 2.8-fold decrease in the KM of Mn(II) oxidation. Second, introducing a salt bridge through Lys179Arg mutation improved Glu35 and Glu181 coordination to Mn(II), decreasing KM 2.6-fold. Third, eliminating a steric clash that prevented Glu37 from orienting toward Mn(II) resulted in an 8.6-fold increase in kcat/KM, arising primarily from a 3.6-fold decrease in KM, with a KM value comparable to that of the native enzyme (0.28 mM vs 0.19 mM for Pleurotus eryngii MnP PS3). We further demonstrated that while the effects of Tyr36Phe and Lys179Arg mutations are additive, because involved in secondary-shell interactions to different ligands, other combinations of mutations were antagonistic because they act on different aspects of the Mn(II) coordination at the same residues. Finally, we showed that these MnCcP variants are functional models of MnP that mimic its activity in both Mn(II) oxidation and degradation of a phenolic lignin model compound and kraft lignin. In addition to achieving KM in a designed protein that is similar to the that of native enzyme, our results offer molecular insight into the role of noncovalent interactions around metal-binding sites for improving metal binding and overall activity; such insight can be applied to rationally enhance these properties in other metalloenzymes and their models.

Cl–/H+ transporters of the CLC superfamily form a ubiquitous class of membrane proteins that catalyze stoichiometrically coupled exchange of Cl– and H+ across biological membranes. CLC transporters exchange H+ for halides and certain polyatomic anions, but exclude cations, F–, and larger physiological anions, such as PO43– and SO42–. Despite comparable transport rates of different anions, the H+ coupling in CLC transporters varies significantly depending on the chemical nature of the transported anion. Although the molecular mechanism of exchange remains unknown, studies on bacterial ClC-ec1 transporter revealed that Cl– binding to the central anion-binding site (Scen) is crucial for the anion-coupled H+ transport. Here, we show that Cl–, F–, NO3–, and SCN– display distinct binding coordinations at the Scen site and are hydrated in different manners. Consistent with the observation of differential bindings, ClC-ec1 exhibits markedly variable ability to support the formation of the transient water wires, which are necessary to support the connection of the two H+ transfer sites (Gluin and Gluex), in the presence of different anions. While continuous water wires are frequently observed in the presence of physiologically transported Cl–, binding of F– or NO3– leads to the formation of pseudo-water-wires that are substantially different from the wires formed with Cl–. Binding of SCN–, however, eliminates the water wires altogether. These findings provide structural details of anion binding in ClC-ec1 and reveal a putative atomic-level mechanism for the decoupling of H+ transport to the transport of anions other than Cl–.

Cytochrome ba3 is a proton-pumping heme-copper oxygen reductase from the extreme thermophile Thermus thermophilus. Despite the fact that the enzyme's active site is buried deep within the protein, the apparent second order rate constant for the initial binding of O2 to the active-site heme has been experimentally found to be 10(9) M(-1) s(-1) at 298 K, at or near the diffusion limit, and 2 orders of magnitude faster than for O2 binding to myoglobin. To provide quantitative and microscopic descriptions of the O2 delivery pathway and mechanism in cytochrome ba3, extensive molecular dynamics simulations of the enzyme in its membrane-embedded form have been performed, including different protocols of explicit ligand sampling (flooding) simulations with O2, implicit ligand sampling analysis, and in silico mutagenesis. The results show that O2 diffuses to the active site exclusively via a Y-shaped hydrophobic tunnel with two 25-Å long membrane-accessible branches that coincide with the pathway previously suggested by the crystallographically identified xenon binding sites. The two entrances of the bifurcated tunnel of cytochrome ba3 are located within the lipid bilayer, where O2 is preferentially partitioned from the aqueous phase. The largest barrier to O2 migration within the tunnel is estimated to be only 1.5 kcal/mol, allowing O2 to reach the enzyme active site virtually impeded by one-dimensional diffusion once it reaches a tunnel entrance at the protein surface. Unlike other O2-utilizing proteins, the tunnel is "open" with no transient barriers observed due to protein dynamics. This unique low-barrier passage through the protein ensures that O2 transit through the protein is never rate-limiting.

A recently developed membrane-mimetic model was applied to study membrane interaction and binding of the two anchoring C2-like discoidin domains of human coagulation factor VIIIa (FVIIIa), the C1 and C2 domains. Both individual domains, FVIII C1 and FVIII C2, were observed to bind the phospholipid membrane by partial or full insertion of their extruding loops (the spikes). However, the two domains adopted different molecular orientations in their membrane-bound states; FVIII C2 roughly was positioned normal to the membrane plane, while FVIII C1 displayed a multitude of tilted orientations. The results indicate that FVIII C1 may be important in modulating the orientation of the FVIIIa molecule to optimize the interaction with FIXa, which is anchored to the membrane via its γ-carboxyglutamic acid-rich (Gla) domain. Additionally, a structural change was observed in FVIII C1 in the coiled main chain leading the first spike. A tight interaction with one lipid per domain, similar to what has been suggested for the homologous FVa C2, is characterized. Finally, we rationalize known FVIII antibody epitopes and the scarcity of documented hemophilic missense mutations related to improper membrane binding of FVIIIa, based on the prevalent nonspecificity of ionic interactions in the simulated membrane-bound states of FVIII C1 and FVIII C2.

Hemagglutinin (HA) is a protein located on the surface of the influenza virus that mediates viral fusion to the host cellular membrane. During the fusion process the HA fusion peptide (HAfp), formed by the first 23 N-terminal residues of HA and structurally characterized by two alpha helices (Helix A and Helix B) tightly packed in a hairpin-like arrangement, is the only part of the virus in direct contact with the host membrane. After encountering the host cell, HAfp is believed to insert into the membrane, thereby initiating the fusion of the viral and host membranes. Detailed characterization of the interactions between the HAfp and cellular membrane is therefore of high relevance to the mechanism of viral entry into the host cell. Employing HMMM membrane representation with enhanced lipid mobility, we have performed a large set of independent simulations of unbiased membrane binding of HAfp. We have been able to capture spontaneous binding and insertion of HAfp consistently in nearly all the simulations. A reproducible membrane-bound configuration emerges from these simulations, despite employing a diverse set of initial configurations. Extension of several of the simulations into full membrane systems confirms the stability of the membrane-bound form obtained from HMMM binding simulations. The resulting model allows for the characterization of important interactions between the peptide and the membrane and the details of the binding process of the peptide for the first time. Upon membrane binding, Helix A inserts much deeper into the membrane than Helix B, suggesting that the former is responsible for hydrophobic anchoring of the peptide into the membrane. Helix B, in contrast, is found to establish major amphipathic interactions at the interfacial region thereby contributing to binding strength of HAfp.

Using Nanodiscs, we quantitate the heterotropic interaction between two different drugs mediated by monomeric CYP3A4 incorporated into a nativelike membrane environment. The mechanism of this interaction is deciphered by global analysis of multiple-turnover experiments performed under identical conditions using the pure substrates progesterone (PGS) and carbamazepine (CBZ) and their mixtures. Activation of CBZ epoxidation and simultaneous inhibition of PGS hydroxylation are measured and quantitated through differences in their respective affinities for both a remote allosteric site and the productive catalytic site near the heme iron. Preferred binding of PGS at the allosteric site and a stronger preference for CBZ binding at the productive site give rise to a nontrivial drug-drug interaction. Molecular dynamics simulations indicate functionally important conformational changes caused by PGS binding at the allosteric site and by two CBZ molecules positioned inside the substrate binding pocket. Structural changes involving Phe-213, Phe-219, and Phe-241 are thought to be responsible for the observed synergetic effects and positive allosteric interactions between these two substrates. Such a mechanism is likely of general relevance to the mutual heterotropic effects caused by biologically active compounds that exhibit different patterns of interaction with the distinct allosteric and productive sites of CYP3A4, as well as other xenobiotic metabolizing cytochromes P450 that are also involved in drug-drug interactions. Importantly, this work demonstrates that a monomeric CYP3A4 can display the full spectrum of activation and cooperative effects that are observed in hepatic membranes.

Ubiquinone forms an integral part of the electron transport chain in cellular respiration and photosynthesis across a vast number of organisms. Prior experimental results have shown that the photosynthetic reaction center (RC) from Rhodobacter sphaeroides is only fully functional with a limited set of methoxy-bearing quinones, suggesting that specific interactions with this substituent are required to drive electron transport and the formation of quinol. The nature of these interactions has yet to be determined. Through parameterization of a CHARMM-compatible quinone force field and subsequent molecular dynamics simulations of the quinone-bound RC, we have investigated and characterized the interactions of the protein with the quinones in the Q(A) and Q(B) sites using both equilibrium simulation and thermodynamic integration. In particular, we identify a specific interaction between the 2-methoxy group of ubiquinone in the Q(B) site and the amide nitrogen of GlyL225 that we implicate in locking the orientation of the 2-methoxy group, thereby tuning the redox potential difference between the quinones occupying the Q(A) and Q(B) sites. Disruption of this interaction leads to weaker binding in a ubiquinone analogue that lacks a 2-methoxy group, a finding supported by reverse electron transfer electron paramagnetic resonance experiments of the Q(A)⁻Q(B)⁻ biradical and competitive binding assays.

While small molecules have been used to induce anesthesia in a clinical setting for well over a century, a detailed understanding of the molecular mechanism remains elusive. In this study, we utilize ab initio calculations to develop a novel set of CHARMM-compatible parameters for the ubiquitous modern anesthetics desflurane, isoflurane, sevoflurane, and propofol for use in molecular dynamics (MD) simulations. The parameters generated were rigorously tested against known experimental physicochemical properties including dipole moment, density, enthalpy of vaporization, and free energy of solvation. In all cases, the anesthetic parameters were able to reproduce experimental measurements, signifying the robustness and accuracy of the atomistic models developed. The models were then used to study the interaction of anesthetics with the membrane. Calculation of the potential of mean force for inserting the molecules into a POPC bilayer revealed a distinct energetic minimum of 4-5 kcal/mol relative to aqueous solution at the level of the glycerol backbone in the membrane. The location of this minimum within the membrane suggests that anesthetics partition to the membrane prior to binding their ion channel targets, giving context to the Meyer-Overton correlation. Moreover, MD simulations of these drugs in the membrane give rise to computed membrane structural parameters, including atomic distribution, deuterium order parameters, dipole potential, and lateral stress profile, that indicate partitioning of anesthetics into the membrane at the concentration range studied here, which does not appear to perturb the structural integrity of the lipid bilayer. These results signify that an indirect, membrane-mediated mechanism of channel modulation is unlikely.

Characterizing large-scale structural transitions in biomolecular systems poses major technical challenges to both experimental and computational approaches. On the computational side, efficient sampling of the configuration space along the transition pathway remains the most daunting challenge. Recognizing this issue, we introduce a knowledge-based computational approach toward describing large-scale conformational transitions using (i) nonequilibrium, driven simulations combined with work measurements and (ii) free energy calculations using empirically optimized biasing protocols. The first part is based on designing mechanistically relevant, system-specific reaction coordinates whose usefulness and applicability in inducing the transition of interest are examined using knowledge-based, qualitative assessments along with nonequilirbrium work measurements which provide an empirical framework for optimizing the biasing protocol. The second part employs the optimized biasing protocol resulting from the first part to initiate free energy calculations and characterize the transition quantitatively. Using a biasing protocol fine-tuned to a particular transition not only improves the accuracy of the resulting free energies but also speeds up the convergence. The efficiency of the sampling will be assessed by employing dimensionality reduction techniques to help detect possible flaws and provide potential improvements in the design of the biasing protocol. Structural transition of a membrane transporter will be used as an example to illustrate the workings of the proposed approach.

The potassium channel KcsA offers a unique opportunity to explicitly study the dynamics of the moving parts of ion channels, yet our understanding of the extent and dynamic behavior of the physiologically relevant structural changes at the inner gate in KcsA remains incomplete. Here, we use electron paramagnetic resonance, nuclear magnetic resonance, and molecular dynamics simulations to characterize the extent of pH-dependent conformational changes of the inner gate in lipid bilayers or detergent micelles. Our results show that under physiological conditions the inner gate experiences a maximal diagonal opening of ∼24 Å with the largest degree of dynamics near the pKa of activation (pH ∼3.9). These results extend the observation that the C-terminus is necessary to limit the extent of opening and imply that the inner gate regulates the extent of conformational change at the zone of allosteric coupling and at the selectivity filter.