Membrane binding and lipid-protein interaction of the C2 domain from coagulation factor V

Anchoring of coagulation factors to anionic regions of the membrane involves the C2 domain as a key player. The rate of enzymatic reactions of the coagulation factors is increased by several orders of magnitude upon membrane binding. However, the precise mechanisms behind the rate acceleration remain unclear, primarily because of a lack of understanding of the conformational dynamics of the C2-containing factors and corresponding complexes. We elucidate the membrane-bound form of the C2 domain from human coagulation factor V (FV-C2) by characterizing its membrane binding the specific lipid-protein interactions. Employing all-atom molecular dynamics simulations and leveraging the highly mobile membrane-mimetic (HMMM) model, we observed spontaneous binding of FV-C2 to a phosphatidylserine (PS)-containing membrane within 2-25 ns across twelve independent simulations. FV-C2 interacted with the membrane through three loops (spikes 1-3), achieving a converged, stable orientation. Multiple HMMM trajectories of the spontaneous membrane binding provided extensive sampling and ample data to examine the membrane-induced effects on the conformational dynamics of C2 as well as specific lipid-protein interactions. Despite existing crystal structures representing presumed "open" and "closed" states of FV-C2, our results revealed a continuous distribution of structures between these states, with the most populated structures differing from both "open" and "closed" states observed in crystal environments. Lastly, we characterized a putative PS-specific binding site formed by K23, Q48, and S78 located in the groove enclosed by spikes 1-3 (PS-specificity pocket), suggesting a different orientation of a bound headgroup moiety compared to previous proposals based upon analysis of static crystal structures.

Staphylococcus aureus FtsZ and PBP4 bind to the conformationally dynamic N-terminal domain of GpsB

In the Firmicutes phylum, GpsB is a membrane associated protein that coordinates peptidoglycan synthesis with cell growth and division. Although GpsB has been studied in several bacteria, the structure, function, and interactome of Staphylococcus aureus GpsB is largely uncharacterized. To address this knowledge gap, we solved the crystal structure of the N-terminal domain of S. aureus GpsB, which adopts an atypical, asymmetric dimer, and demonstrates major conformational flexibility that can be mapped to a hinge region formed by a three-residue insertion exclusive to Staphylococci. When this three-residue insertion is excised, its thermal stability increases, and the mutant no longer produces a previously reported lethal phenotype when overexpressed in Bacillus subtilis. In S. aureus, we show that these hinge mutants are less functional and speculate that the conformational flexibility imparted by the hinge region may serve as a dynamic switch to fine-tune the function of the GpsB complex and/or to promote interaction with its various partners. Furthermore, we provide the first biochemical, biophysical, and crystallographic evidence that the N-terminal domain of GpsB binds not only PBP4, but also FtsZ, through a conserved recognition motif located on their C-termini, thus coupling peptidoglycan synthesis to cell division. Taken together, the unique structure of S. aureus GpsB and its direct interaction with FtsZ/PBP4 provide deeper insight into the central role of GpsB in S. aureus cell division.

Distal Protein-Protein Interactions Contribute to SARS-CoV-2 Main Protease Substrate Binding and Nirmatrelvir Resistance

SARS-CoV-2 main protease, M pro , is responsible for the processing of the viral polyproteins into individual proteins, including the protease itself. M pro is a key target of anti-COVID-19 therapeutics such as nirmatrelvir (the active component of Paxlovid). Resistance mutants identified clinically and in viral passage assays contain a combination of active site mutations (e.g. E166V, E166A, L167F), which reduce inhibitor binding and enzymatic activity, and non-active site mutations (e.g. P252L, T21I, L50F), which restore the fitness of viral replication. Although the mechanism of resistance for the active site mutations is apparent, the role of the non-active site mutations in fitness rescue remains elusive. In this study, we use the model system of a M pro triple mutant (L50F/E166A/L167F) that confers not only nirmatrelvir drug resistance but also a similar fitness of replication compared to the wild-type both in vitro and in vivo. By comparing peptide and full-length M pro protein as substrates, we demonstrate that the binding of M pro substrate involves more than residues in the active site. In particular, L50F and other non-active site mutations can enhance the M pro dimer-dimer interactions and help place the nsp5-6 substrate at the enzyme catalytic center. The structural and enzymatic activity data of M pro L50F, L50F/E166A/L167F, and others underscore the importance of considering the whole substrate protein in studying M pro and substrate interactions, and offers important insights into M pro function, resistance development, and inhibitor design.

Overlay databank unlocks data-driven analyses of biomolecules for all

Tools based on artificial intelligence (AI) are currently revolutionising many fields, yet their applications are often limited by the lack of suitable training data in programmatically accessible format. Here we propose an effective solution to make data scattered in various locations and formats accessible for data-driven and machine learning applications using the overlay databank format. To demonstrate the practical relevance of such approach, we present the NMRlipids Databank-a community-driven, open-for-all database featuring programmatic access to quality-evaluated atom-resolution molecular dynamics simulations of cellular membranes. Cellular membrane lipid composition is implicated in diseases and controls major biological functions, but membranes are difficult to study experimentally due to their intrinsic disorder and complex phase behaviour. While MD simulations have been useful in understanding membrane systems, they require significant computational resources and often suffer from inaccuracies in model parameters. Here, we demonstrate how programmable interface for flexible implementation of data-driven and machine learning applications, and rapid access to simulation data through a graphical user interface, unlock possibilities beyond current MD simulation and experimental studies to understand cellular membranes. The proposed overlay databank concept can be further applied to other biomolecules, as well as in other fields where similar barriers hinder the AI revolution.

Dynamic Nature of Staphylococcus aureus Type I Signal Peptidases

Molecular dynamics simulations are used to interrogate the dynamic nature of Staphylococcus aureus Type I signal peptidases, SpsA and SpsB, including the impact of the P29S mutation of SpsB. Fluctuations and plasticity- rigidity characteristics vary among the proteins, particularly in the extracellular domain. Intriguingly, the P29S mutation, which influences susceptibility to arylomycin antibiotics, affect the mechanically coupled motions in SpsB. The integrity of the active site is crucial for catalytic competency, and variations in sampled structural conformations among the proteins are consistent with diverse peptidase capabilities. We also explored the intricate interactions between the proteins and the model S. aureus membrane. It was observed that certain membrane-inserted residues in the loop around residue 50 (50s) and C-terminal loops, beyond the transmembrane domain, give rise to direct interactions with lipids in the bilayer membrane. Our findings are discussed in the context of functional knowledge about these signal peptidases, offering additional understanding of dynamic aspects relevant to some cellular processes with potential implications for drug targeting strategies.

A 19F-qNMR-Guided Mathematical Model for G Protein-Coupled Receptor Signaling

G protein-coupled receptors (GPCRs) exhibit a wide range of pharmacological efficacies, yet the molecular mechanisms responsible for the differential efficacies in response to various ligands remain poorly understood. This lack of understanding has hindered the development of a solid foundation for establishing a mathematical model for signaling efficacy. However, recent progress has been made in delineating and quantifying receptor conformational states and associating function with these conformations. This progress has allowed us to construct a mathematical model for GPCR signaling efficacy that goes beyond the traditional ON/OFF binary switch model. In this study, we present a quantitative conformation-based mathematical model for GPCR signaling efficacy using the adenosine A2A receptor (A2AR) as a model system, under the guide of 19F quantitative nuclear magnetic resonance experiments. This model encompasses two signaling states, a fully activated state and a partially activated state, defined as being able to regulate the cognate Gα s nucleotide exchange with respective G protein recognition capacity. By quantifying the population distribution of each state, we can now in turn examine GPCR signaling efficacy. This advance provides a foundation for assessing GPCR signaling efficacy using a conformation-based mathematical model in response to ligand binding. SIGNIFICANCE STATEMENT: Mathematical models to describe signaling efficacy of GPCRs mostly suffer from considering only two states (ON/OFF). However, research indicates that a GPCR possesses multiple active-(like) states that can interact with Gαβγ independently, regulating varied nucleotide exchanges. With the guide of 19F-qNMR, the transitions among these states are quantified as a function of ligand and Gαβγ, serving as a foundation for a novel conformation-based mathematical signaling model.

Deciphering specificity and cross-reactivity in tachykinin NK1 and NK2 receptors

The tachykinin receptors neurokinin 1 (NK1R) and neurokinin 2 (NK2R) are G protein-coupled receptors that bind preferentially to the natural peptide ligands substance P and neurokinin A, respectively, and have been targets for drug development. Despite sharing a common C-terminal sequence of Phe-X-Gly-Leu-Met-NH2 that helps direct biological function, the peptide ligands exhibit some degree of cross-reactivity toward each other's non-natural receptor. Here, we investigate the detailed structure-activity relationships of the ligand-bound receptor complexes that underlie both potent activation by the natural ligand and cross-reactivity. We find that the specificity and cross-reactivity of the peptide ligands can be explained by the interactions between the amino acids preceding the FxGLM consensus motif of the bound peptide ligand and two regions of the receptor: the β-hairpin of the extracellular loop 2 (ECL2) and a N-terminal segment leading into transmembrane helix 1. Positively charged sidechains of the ECL2 (R177 of NK1R and K180 of NK2R) are seen to play a vital role in the interaction. The N-terminal positions 1 to 3 of the peptide ligand are entirely dispensable. Mutated and chimeric receptor and ligand constructs neatly swap around ligand specificity as expected, validating the structure-activity hypotheses presented. These findings will help in developing improved agonists or antagonists for NK1R and NK2R.

Quantitative Comparison against Experiments Reveals Imperfections in Force Fields' Descriptions of POPC-Cholesterol Interactions

Cholesterol is a central building block in biomembranes, where it induces orientational order, slows diffusion, renders the membrane stiffer, and drives domain formation. Molecular dynamics (MD) simulations have played a crucial role in resolving these effects at the molecular level; yet, it has recently become evident that different MD force fields predict quantitatively different behavior. Although easily neglected, identifying such limitations is increasingly important as the field rapidly progresses toward simulations of complex membranes mimicking the in vivo conditions: pertinent multicomponent simulations must capture accurately the interactions between their fundamental building blocks, such as phospholipids and cholesterol. Here, we define quantitative quality measures for simulations of binary lipid mixtures in membranes against the C-H bond order parameters and lateral diffusion coefficients from NMR spectroscopy as well as the form factors from X-ray scattering. Based on these measures, we perform a systematic evaluation of the ability of commonly used force fields to describe the structure and dynamics of binary mixtures of palmitoyloleoylphosphatidylcholine (POPC) and cholesterol. None of the tested force fields clearly outperforms the others across the tested properties and conditions. Still, the Slipids parameters provide the best overall performance in our tests, especially when dynamic properties are included in the evaluation. The quality evaluation metrics introduced in this work will, particularly, foster future force field development and refinement for multicomponent membranes using automated approaches.

Cholesterol and M2 Rendezvous in Budding and Scission of Influenza A Virus

The cholesterol of the host cell plasma membrane and viral M2 protein plays a crucial role in multiple stages of infection and replication of the influenza A virus. Cholesterol is required for the formation of heterogeneous membrane microdomains (or rafts) in the budozone of the host cell that serves as assembly sites for the viral components. The raft microstructures act as scaffolds for several proteins. Cholesterol may further contribute to the mechanical forces necessary for membrane scission in the last stage of budding and help to maintain the stability of the virus envelope. The M2 protein has been shown to cause membrane scission in model systems by promoting the formation of curved lipid bilayer structures that, in turn, can lead to membrane vesicles budding off or scission intermediates. Membrane remodeling by M2 is intimately linked with cholesterol as it affects local lipid composition, fluidity, and stability of the membrane. Thus, both cholesterol and M2 protein contribute to the efficient and proper release of newly formed influenza viruses from the virus-infected cells.

Mechanistic basis of GPCR activation explored by ensemble refinement of crystallographic structures

G protein‐coupled receptors (GPCRs) are important drug targets characterized by a canonical seven transmembrane (TM) helix architecture. Recent advances in X‐ray crystallography and cryo‐EM have resulted in a wealth of GPCR structures that have been used in drug design and formed the basis for mechanistic activation hypotheses. Here, ensemble refinement (ER) of crystallographic structures is applied to explore the impact of binding of agonists and antagonist/inverse agonists to selected structures of cannabinoid receptor 1 (CB1R), β₂ adrenergic receptor (β₂AR), and A_2A adenosine receptor (A_2AAR). To assess the conformational flexibility and its role in GPCR activation, hydrogen bond (H‐bond) networks are analyzed by calculating and comparing H‐bond propensities. Mapping pairwise propensity differences between agonist‐ and inverse agonist/antagonist‐bound structures for CB1R and β₂AR shows that agonist binding destabilizes H‐bonds in the intracellular parts of TM 5–7, forming the G protein binding cavity, while H‐bonds of the extracellular segment of TMs surrounding the orthosteric site are conversely stabilized. Certain class A GPCRs, for example, A_2AAR, bind an allosteric sodium ion that negatively modulates agonist binding. The impact of sodium‐excluding mutants (D52^2.50N, S91^3.39A) of A_2AAR on agonist binding is examined by applying ER analysis to structures of wildtype and the two mutants in complex with a full agonist. While S91^3.39A exhibits normal activity, D52^2.50N quenches the downstream signaling. The mainchain H‐bond pattern of the latter is stabilized in the intracellular part of TM 7 containing the NPxxY motif, indicating that an induced rigidity of the mutation prevents conformational selection of G proteins resulting in receptor inactivation.

19F NMR: A promising tool for dynamic conformational studies of G protein-coupled receptors

Advances in X-ray crystallography and cryoelectron microscopy enabled unprecedented insights into the activation processes of G protein-coupled receptors (GPCRs). However, these static receptor structures provide limited information about dynamics and conformational transitions that play pivotal roles in mediating signaling diversity through the multifaceted interactions between ligands, receptors, and transducers. Developing NMR approaches to probe the dynamics of conformational transitions will push the frontier of receptor science toward a more comprehensive understanding of these signaling processes. Although much progress has been made during the last decades, it remains challenging to delineate receptor conformational states and interrogate the functions of the individual states at a quantitative level. Here we cover the progress of ¹⁹F NMR applications in GPCR conformational and dynamic studies during the past 20 years. Current challenges and limitations of ¹⁹F NMR for studying GPCR dynamics are also discussed, along with experimental strategies that will drive this field forward.

An in-membrane NMR spectroscopic approach probing native ligand-GPCR interaction

Conventional approaches to study ligand-receptor interactions using solution-state NMR often involve laborious sample preparation, isotopic labeling, and receptor reconstitution. Each of these steps remains challenging for membrane proteins such as G protein-coupled receptors (GPCRs). Here we introduce a combinational approach integrating NMR and homogenized membrane nano-discs preparation to characterize the ligand-GPCR interactions. The approach will have a great potential for drug screening as it benefits from minimal receptor preparation, minimizing non-specific binding. In addition, the approach maintains receptor structural heterogeneity essential for functional diversity, making it feasible for probing a more reliable ligand-GPCR interaction that is vital for faithful ligand discovery.

A systematic approach for evaluating the role of surface-exposed loops in trypsin-like serine proteases applied to the 170 loop in coagulation factor VIIa

Proteases play a major role in many vital physiological processes. Trypsin-like serine proteases (TLPs), in particular, are paramount in proteolytic cascade systems such as blood coagulation and complement activation. The structural topology of TLPs is highly conserved, with the trypsin fold comprising two β-barrels connected by a number of variable surface-exposed loops that provide a surprising capacity for functional diversity and substrate specificity. To expand our understanding of the roles these loops play in substrate and co-factor interactions, we employ a systematic methodology akin to the natural truncations and insertions observed through evolution of TLPs. The approach explores a larger deletion space than classical random or directed mutagenesis. Using FVIIa as a model system, deletions of 1–7 amino acids through the surface exposed 170 loop, a vital allosteric regulator, was introduced. All variants were extensively evaluated by established functional assays and computational loop modelling with Rosetta. The approach revealed detailed structural and functional insights recapitulation and expanding on the main findings in relation to 170 loop functions elucidated over several decades using more cumbersome crystallization and single deletion/mutation methodologies. The larger deletion space was key in capturing the most active variant, which unexpectedly had a six-amino acid truncation. This variant would have remained undiscovered if only 2–3 deletions were considered, supporting the usefulness of the methodology in general protease engineering approaches. Our findings shed further light on the complex role that surface-exposed loops play in TLP function and supports the important role of loop length in the regulation and fine-tunning of enzymatic function throughout evolution.

Inverse Conformational Selection in Lipid-Protein Binding

Interest in lipid interactions with proteins and other biomolecules is emerging not only in fundamental biochemistry but also in the field of nanobiotechnology where lipids are commonly used, for example, in carriers of mRNA vaccines. The outward-facing components of cellular membranes and lipid nanoparticles, the lipid headgroups, regulate membrane interactions with approaching substances, such as proteins, drugs, RNA, or viruses. Because lipid headgroup conformational ensembles have not been experimentally determined in physiologically relevant conditions, an essential question about their interactions with other biomolecules remains unanswered: Do headgroups exchange between a few rigid structures, or fluctuate freely across a practically continuous spectrum of conformations? Here, we combine solid-state NMR experiments and molecular dynamics simulations from the NMRlipids Project to resolve the conformational ensembles of headgroups of four key lipid types in various biologically relevant conditions. We find that lipid headgroups sample a wide range of overlapping conformations in both neutral and charged cellular membranes, and that differences in the headgroup chemistry manifest only in probability distributions of conformations. Furthermore, the analysis of 894 protein-bound lipid structures from the Protein Data Bank suggests that lipids can bind to proteins in a wide range of conformations, which are not limited by the headgroup chemistry. We propose that lipids can select a suitable headgroup conformation from the wide range available to them to fit the various binding sites in proteins. The proposed inverse conformational selection model will extend also to lipid binding to targets other than proteins, such as drugs, RNA, and viruses.

Uncovering Membrane-Bound Models of Coagulation Factors by Combined Experimental and Computational Approaches

In the life sciences, including hemostasis and thrombosis, methods of structural biology have become indispensable tools for shedding light on underlying mechanisms that govern complex biological processes. Advancements of the relatively young field of computational biology have matured to a point where it is increasingly recognized as trustworthy and useful, in part due to their high space–time resolution that is unparalleled by most experimental techniques to date. In concert with biochemical and biophysical approaches, computational studies have therefore proven time and again in recent years to be key assets in building or suggesting structural models for membrane-bound forms of coagulation factors and their supramolecular complexes on membrane surfaces where they are activated. Such endeavors and the proposed models arising from them are of fundamental importance in describing and understanding the molecular basis of hemostasis under both health and disease conditions. We summarize the body of work done in this important area of research to drive forward both experimental and computational studies toward new discoveries and potential future therapeutic strategies.

Conformational Plasticity-Rigidity Axis of the Coagulation Factor VII Zymogen Elucidated by Atomistic Simulations of the N-Terminally Truncated Factor VIIa Protease Domain

The vast majority of coagulation factor VII (FVII), a trypsin-like protease, circulates as the inactive zymogen. Activated FVII (FVIIa) is formed upon proteolytic activation of FVII, where it remains in a zymogen-like state and it is fully activated only when bound to tissue factor (TF). The catalytic domains of trypsin-like proteases adopt strikingly similar structures in their fully active forms. However, the dynamics and structures of the available corresponding zymogens reveal remarkable conformational plasticity of the protease domain prior to activation in many cases. Exactly how ligands and cofactors modulate the conformational dynamics and function of these proteases is not entirely understood. Here, we employ atomistic simulations of FVIIa (and variants hereof, including a TF-independent variant and N-terminally truncated variants) to provide fundamental insights with atomistic resolution into the plasticity-rigidity interplay of the protease domain conformations that appears to govern the functional response to proteolytic and allosteric activation. We argue that these findings are relevant to the FVII zymogen, whose structure has remained elusive despite substantial efforts. Our results shed light on the nature of FVII and demonstrate how conformational dynamics has played a crucial role in the evolutionary adaptation of regulatory mechanisms that were not present in the ancestral trypsin. Exploiting this knowledge could lead to engineering of protease variants for use as next-generation hemostatic therapeutics.

Trifluorinated Keto-Enol Tautomeric Switch in Probing Domain Rotation of a G Protein-Coupled Receptor

Conformational dynamics and transitions of biologically active molecules are pivotal for understanding the physiological responses they elicit. In the case of receptor activation, there are major implications elucidating disease mechanisms and drug discovery innovation. Yet, incorporation of these factors into drug screening systems remains challenging in part due to the lack of suitable approaches to include them. Here, we present a novel strategy to probe the GPCR domain rotation by utilizing the ¹⁹fluorine signal variability of a trifluorinated keto-enol (TFKE) chemical equilibrium. The method takes advantage of the high sensitivity of the TFKE tautomerism toward microenvironmental changes resulting from receptor conformational transitions upon ligand binding. We validated the method using the adenosine A_2AR receptor as a model system in which the TFKE was attached to two sites exhibiting opposing motions upon ligand binding, namely, V229C^6.31 on transmembrane domain VI (TM6) and A289C^7.54 on TM7. Our results demonstrated that the TFKE switch was an excellent reporter for the domain rotation and could be used to study the conformational transition and dynamics of relative domain motions. Although further studies are needed in order to establish a quantitative relationship between the rotational angle and the population distribution of different components in a particular system, the research presented here provides a foundation for its application in studying receptor domain rotation and dynamics, which could be useful in drug screening efforts.

Cholesterol Alters the Orientation and Activity of the Influenza Virus M2 Amphipathic Helix in the Membrane

The influenza virus M2 amphipathic helix (M2AH) alters membrane curvature in a cholesterol-dependent manner, mediating viral membrane scission during influenza virus budding. Here, we have investigated the biophysical effects of cholesterol on the ability of an M2AH peptide to manipulate membrane properties. We see that the ability of the M2AH to interact with membranes and form an α-helix is independent of membrane cholesterol concentration; however, cholesterol affects the angle of the M2AH peptide within the membrane. This change in membrane orientation affects the ability of the M2AH to alter lipid order. In low-cholesterol membranes, the M2AH is inserted near the level of the lipid head groups, increasing lipid order, which may contribute to generation of the membrane curvature. As the cholesterol content increases, the M2AH insertion becomes flatter and slightly deeper in the membrane below the lipid headgroups, where the polar face can continue to interact with the headgroups while the hydrophobic face binds cholesterol. This changed orientation minimizes lipid packing defects and lipid order changes, likely reducing the generation of membrane curvature. Thus, cholesterol regulates M2 membrane scission by precisely modulating M2AH positioning within the membrane. This has implications for the understanding of many of amphipathic-helix-driven cellular budding processes that occur in specific lipid environments.

Water-Intake and Water-Molecule Paths to the Active Site of Secretory Phospholipase A2 Studied Using MD Simulations and the Tracking Tool AQUA-DUCT

Secretory phospholipases A₂ (sPLA₂s) are a subclass of enzymes that catalyze the hydrolysis at the sn-2 position of glycerophospholipids, producing free fatty acids and lysophospholipids. In this study, different phospholipids with structural modifications close to the scissile sn-2 ester bond were studied to determine the effect of the structural changes on the formation of the Michaelis-Menten complex and the water entry/exit pathways using molecular dynamics simulations and the computational tracking tool AQUA-DUCT. Structural modifications include methylation, dehydrogenation, and polarization close to the sn-2 scissile bond. We found that all water molecules reaching the active site of sPLA₂-IIA pass by the aromatic residues Phe⁵ and Tyr⁵¹ and enter the active site through an active-site cleft. The relative amount of water available for the enzymatic reaction of the different phospholipid-sPLA₂ complexes was determined together with the distance between key atoms in the catalytic machinery. The results showed that (Z)-unsaturated phospholipid is a good substrate for sPLA₂-IIA. The computational results are in good agreement with previously reported experimental data on the ability of sPLA₂-IIA to hydrolyze liposomes made from the different phospholipids, and the results provide new insights into the necessary active-site solvation of the Michaelis-Menten complex and can pave the road for rational design in engineering applications.

Allostery in Coagulation Factor VIIa Revealed by Ensemble Refinement of Crystallographic Structures

A critical step in injury-induced initiation of blood coagulation is the formation of the complex between the trypsin-like protease coagulation factor VIIa (FVIIa) and its cofactor tissue factor (TF), which converts FVIIa from an intrinsically poor enzyme to an active protease capable of activating zymogens of downstream coagulation proteases. Unlike its constitutively active ancestor trypsin, FVIIa is allosterically activated (by TF). Here, ensemble refinement of crystallographic structures, which uses multiple copies of the entire structure as a means of representing structural flexibility, is applied to explore the impacts of inhibitor binding to trypsin and FVIIa, as well as cofactor binding to FVIIa. To assess the conformational flexibility and its role in allosteric pathways in these proteases, main-chain hydrogen bond networks are analyzed by calculating the hydrogen-bond propensity. Mapping pairwise propensity differences between relevant structures shows that binding of the inhibitor benzamidine to trypsin has a minor influence on the protease flexibility. For FVIIa, in contrast, the protease domain is "locked" into the catalytically competent trypsin-like configuration upon benzamidine binding as indicated by the stabilization of key structural features: the nonprime binding cleft and the oxyanion hole are stabilized, and the effect propagates from the active site region to the calcium-binding site and to the vicinity of the disulphide bridge connecting with the light chain. TF binding to FVIIa furthermore results in stabilization of the 170 loop, which in turn propagates an allosteric signal from the TF-binding region to the active site. Analyses of disulphide bridge energy and flexibility reflect the striking stability difference between the unregulated enzyme and the allosterically activated form after inhibitor or cofactor binding. The ensemble refinement analyses show directly, for the first time to our knowledge, whole-domain structural footprints of TF-induced allosteric networks present in x-ray crystallographic structures of FVIIa, which previously only have been hypothesized or indirectly inferred.

Systematic Coarse-Grained Lipid Force Fields with Semiexplicit Solvation via Virtual Sites

Despite the central role of lipids in many biophysical functions, the molecular mechanisms that dictate macroscopic lipid behavior remain elusive to both experimental and computational approaches. As such, there has been much interest in the development of low-resolution, implicit-solvent coarse-grained (CG) models to dynamically simulate biologically relevant spatiotemporal scales with molecular fidelity. However, in the absence of solvent, a key challenge for CG models is to faithfully emulate solvent-mediated forces, which include both hydrophilic and hydrophobic interactions that drive lipid aggregation and self-assembly. In this work, we provide a new methodological framework to incorporate semiexplicit solvent effects through the use of virtual CG particles, which represent structural features of the solvent-lipid interface. To do so, we leverage two systematic coarse-graining approaches, multiscale coarse-graining (MS-CG) and relative entropy minimization (REM), in a hybrid fashion to construct our virtual-site CG (VCG) models. As a proof-of-concept, we focus our efforts on two lipid species, 1,2-dioleoyl-sn-glycero-3-phosphocholine (DOPC) and 1,2-dipalmitoyl-sn-glycero-3-phosphocholine (DPPC), which adopt a liquid-disordered and gel phase, respectively, at room temperature. Through our analysis, we also present, to our knowledge, the first direct comparison between the MS-CG and REM methods for a complex biomolecule and highlight each of their strengths and weaknesses. We further demonstrate that VCG models recapitulate the rich biophysics of lipids, which enable self-assembly, morphological diversity, and multiple phases. Our findings suggest that the VCG framework is a powerful approach for investigation into macromolecular biophysics.

Acid activation mechanism of the influenza A M2 proton channel

The homotetrameric influenza A M2 channel (AM2) is an acid-activated proton channel responsible for the acidification of the influenza virus interior, an important step in the viral lifecycle. Four histidine residues (His37) in the center of the channel act as a pH sensor and proton selectivity filter. Despite intense study, the pH-dependent activation mechanism of the AM2 channel has to date not been completely understood at a molecular level. Herein we have used multiscale computer simulations to characterize (with explicit proton transport free energy profiles and their associated calculated conductances) the activation mechanism of AM2. All proton transfer steps involved in proton diffusion through the channel, including the protonation/deprotonation of His37, are explicitly considered using classical, quantum, and reactive molecular dynamics methods. The asymmetry of the proton transport free energy profile under high-pH conditions qualitatively explains the rectification behavior of AM2 (i.e., why the inward proton flux is allowed when the pH is low in viral exterior and high in viral interior, but outward proton flux is prohibited when the pH gradient is reversed). Also, in agreement with electrophysiological results, our simulations indicate that the C-terminal amphipathic helix does not significantly change the proton conduction mechanism in the AM2 transmembrane domain; the four transmembrane helices flanking the channel lumen alone seem to determine the proton conduction mechanism.

Molecular Basis of Enhanced Activity in Factor VIIa-Trypsin Variants Conveys Insights into Tissue Factor-mediated Allosteric Regulation of Factor VIIa Activity

The complex of coagulation factor VIIa (FVIIa), a trypsin-like serine protease, and membrane-bound tissue factor (TF) initiates blood coagulation upon vascular injury. Binding of TF to FVIIa promotes allosteric conformational changes in the FVIIa protease domain and improves its catalytic properties. Extensive studies have revealed two putative pathways for this allosteric communication. Here we provide further details of this allosteric communication by investigating FVIIa loop swap variants containing the 170 loop of trypsin that display TF-independent enhanced activity. Using x-ray crystallography, we show that the introduced 170 loop from trypsin directly interacts with the FVIIa active site, stabilizing segment 215-217 and activation loop 3, leading to enhanced activity. Molecular dynamics simulations and novel fluorescence quenching studies support that segment 215-217 conformation is pivotal to the enhanced activity of the FVIIa variants. We speculate that the allosteric regulation of FVIIa activity by TF binding follows a similar path in conjunction with protease domain N terminus insertion, suggesting a more complete molecular basis of TF-mediated allosteric enhancement of FVIIa activity.

Membrane Interaction of the Factor VIIIa Discoidin Domains in Atomistic Detail

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.

Factor VIII Interacts with the Endocytic Receptor Low-density Lipoprotein Receptor-related Protein 1 via an Extended Surface Comprising "Hot-Spot" Lysine Residues

Lysine residues are implicated in driving the ligand binding to the LDL receptor family. However, it has remained unclear how specificity is regulated. Using coagulation factor VIII as a model ligand, we now study the contribution of individual lysine residues in the interaction with the largest member of the LDL receptor family, low-density lipoprotein receptor-related protein (LRP1). Using hydrogen-deuterium exchange mass spectrometry (HDX-MS) and SPR interaction analysis on a library of lysine replacement variants as two independent approaches, we demonstrate that the interaction between factor VIII (FVIII) and LRP1 occurs over an extended surface containing multiple lysine residues. None of the individual lysine residues account completely for LRP1 binding, suggesting an additive binding model. Together with structural docking studies, our data suggest that FVIII interacts with LRP1 via an extended surface of multiple lysine residues that starts at the bottom of the C1 domain and winds around the FVIII molecule.

Tissue factor activates allosteric networks in factor VIIa through structural and dynamic changes

BACKGROUND

Tissue factor (TF) promotes colocalization of enzyme (factor VIIa) and substrate (FX or FIX), and stabilizes the active conformation of FVIIa. Details on how TF induces structural and dynamic changes in the catalytic domain of FVIIa to enhance its efficiency remain elusive.

OBJECTIVE

To elucidate the activation of allosteric networks in the catalytic domain of the FVIIa protease it is when bound to TF.

METHODS

Long-timescale molecular dynamics simulations of FVIIa, free and in complex with TF, were executed and analyzed by dynamic network analysis.

RESULTS

Allosteric paths of correlated motion from the TF contact point, Met306, in FVIIa to the active site triad can be described and quantified. In particular, the shortest paths from Met306 to Ser344 and His193 are 16% and 8% longer in free FVIIa than in TF-FVIIa, and they encompass previously undiscovered residue-residue interactions that are not likely to be inferred from mutagenesis studies. Furthermore, paths from Met306 to Ile153 (N-terminus) and Trp364, both representing hallmark residues of allostery, are 7% and 37% longer, respectively, in free FVIIa. Thus, there is significantly weaker coupling between the TF contact point and key residues in the catalytic domain of FVIIa, causing the active site triad to disintegrate in the simulation when TF is not present.

CONCLUSIONS

These findings complement our current understanding of how the protease FVIIa is stimulated by TF. We demonstrate allosteric networks in the catalytic domain that are activated by TF and help to make FVIIa an efficient catalyst of FIX and FX activation.

The length of the linker between the epidermal growth factor-like domains in factor VIIa is critical for a productive interaction with tissue factor

Formation of the factor VIIa (FVIIa)-tissue factor (TF) complex triggers the blood coagulation cascade. Using a structure-based rationale, we investigated how the length of the linker region between the two epidermal growth factor (EGF)-like domains in FVIIa influences TF binding and the allosteric activity enhancement, as well as the interplay between the γ-carboxyglutamic acid (Gla)-containing and protease domains. Removal of two residues from the native linker was compatible with normal cofactor binding and accompanying stimulation of the enzymatic activity, as was extension by two (Gly-Ser) residues. In sharp contrast, truncation by three or four residues abolished the TF-mediated stabilization of the active conformation of FVIIa and abrogated TF-induced activity enhancement. In addition, FVIIa variants with short linkers associated 80-fold slower with soluble TF (sTF) as compared with wild-type FVIIa, resulting in a corresponding increase in the equilibrium dissociation constant. Molecular modeling suggested that the shortest FVIIa variants would have to be forced into a tense and energetically unfavorable conformation in order to be able to interact productively with TF, explaining our experimental observations. We also found a correlation between linker length and the residual intrinsic enzymatic activity of Ca(2+)-free FVIIa; stepwise truncation resulting in gradually higher activity with des(83-86)-FVIIa reaching the level of Gla-domainless FVIIa. The linker appears to determine the average distance between the negatively charged Gla domain and a structural element in the protease domain, presumably of opposite charge, and proximity has a negative impact on apo-FVIIa activity.

Interaction of neurotransmitters with a phospholipid bilayer: a molecular dynamics study

We have performed a series of molecular dynamics simulations to study the interactions between the neurotransmitters (NTs) γ-aminobutyrate (GABA), glycine (GLY), acetylcholine (ACH) and glutamate (GLU) as well as the amidated/acetylated γ-aminobutyrate (GABA(neu)) and the osmolyte molecule glycerol (GOL) with a dipalmitoylphosphatidylcholine (DPPC) bilayer. In agreement with previously published experimental data, we found the lowest membrane affinity for the charged molecules and a moderate affinity for zwitterionic and polar molecules. The affinity can be ranked as follows: ACH-GLU<<GABA<GLY<<GABA(neu)<<GOL. The latter three penetrated the bilayer at most with the deepest location being close to the glycerol backbone of the phospholipids. Even at that position, these solutes were noticeably hydrated and carried ∼30-80% of the bulk water along. The mobility of hydration water at the solute is also affected by the penetration into the bilayer. Two time scales of exchanging water molecules could be determined. In the bulk phase, the hydration layer contains ∼20% slow exchanging water molecules which increases 2-3 times as the solutes entered the bilayer. Our results indicate that there is no intermediate exchange of slow moving water molecules from the solutes to the lipid atoms and vice versa. Instead, the exchange relies on the reservoir of unbounded ("free") water molecules in the interfacial bilayer region. Results from the equilibrium simulations are in good agreement with the results from umbrella sampling simulations, which were conducted for the four naturally occurring NTs. Free energy profiles for ACH and GLU show a minimum of ∼2-3 kJ/mol close to the bilayer interface, while for GABA and GLY, a minimum of respectively ∼2 kJ/mol and ∼5 kJ/mol is observed when these NTs are located in the vicinity of the lipid glycerol backbone. The most important interaction of NTs with the bilayer is the charged amino group of NTs with the lipid phosphate group.

Binding of serotonin to lipid membranes

Serotonin (5-hydroxytryptamine, 5-HT) is a prevalent neurotransmitter throughout the animal kingdom. It exerts its effect through the specific binding to the serotonin receptor, but recent research has suggested that neural transmission may also be affected by its nonspecific interactions with the lipid matrix of the synaptic membrane. However, membrane-5-HT interactions remain controversial and superficially investigated. Fundamental knowledge of this interaction appears vital in discussions of putative roles of 5-HT, and we have addressed this by thermodynamic measurements and molecular dynamics (MD) simulations. 5-HT was found to interact strongly with lipid bilayers (partitioning coefficient ~1200 in mole fraction units), and this is highly unusual for a hydrophilic solute like 5-HT which has a bulk, oil-water partitioning coefficient well below unity. It follows that membrane affinity must rely on specific interactions, and the MD simulations identified the salt-bridge between the primary amine of 5-HT and the lipid phosphate group as the most important interaction. This interaction anchored cationic 5-HT in the membrane interface with the aromatic ring system pointing inward and a prevailing residence between the phosphate and the carbonyl groups of the lipid. The unprotonated form of 5-HT shows the opposite orientation, with the primary amine pointing toward the membrane core. Partitioning of 5-HT was found to decrease lipid chain order. These distinctive interactions of 5-HT and model membranes could be related to nonspecific effects of this neurotransmitter.

Secretory phospholipase A2 activity toward diverse substrates

We have studied secretory phospholipase A(2)-IIA (sPLA(2)) activity toward different phospholipid analogues by performing biophysical characterizations and molecular dynamics simulations. The phospholipids were natural substrates, triple alkyl phospholipids, a prodrug anticancer etherlipid, and an inverted ester. The latter were included to study head group-enzyme interactions. Our simulation results show that the lipids are optimally placed into the binding cleft and that water molecules can freely reach the active site through a well-defined pathway; both are indicative that these substrates are efficiently hydrolyzed, which is in good agreement with our experimental data. The phospholipid analogue with three alkyl side chains forms aggregates of different shapes with no well-defined sizes due to its cone-shape structure. Phosphatidylglycerol and phosphatidylcholine head groups interact with specific charged residues, but relatively large fluctuations are observed, suggesting that these interactions are not necessarily important for stabilizing substrate binding to the enzyme.

Transmission of Salmonella between wildlife and meat-production animals in Denmark

AIMS

To investigate the transmission of Salmonella spp. between production animals (pigs and cattle) and wildlife on production animal farms in Denmark.

METHODS AND RESULTS

In the winter and summer of 2001 and 2002, 3622 samples were collected from Salmonella-infected and noninfected herds of pigs and cattle and surrounding wildlife. Salmonella was detected in wildlife on farms carrying Salmonella-positive production animals and only during the periods when Salmonella was detected in the production animals. The presence of Salmonella Typhimurium in wild birds significantly correlated to their migration pattern and food preference.

CONCLUSIONS

Salmonella was transmitted from infected herds of production animals (cattle and pigs) to wildlife that lived amongst or in close proximity to them.

SIGNIFICANCE AND IMPACT OF THE STUDY

Salmonella in animal food products is associated with the occurrence of Salmonella in primary animal production. Strategies to control the introduction and spread of infection should include wildlife management, as the nearby wildlife may act as reservoirs for Salmonella spp. and/or may be passive carriers of the bacteria.

Verocytotoxin-producing Escherichia coli in wild birds and rodents in close proximity to farms

Wild animals living close to cattle and pig farms (four each) were examined for verocytotoxin-producing Escherichia coli (VTEC; also known as Shiga toxin-producing E. coli). The prevalence of VTEC among the 260 samples from wild animals was generally low. However, VTEC isolates from a starling (Sturnus vulgaris) and a Norway rat (Rattus norvegicus) were identical to cattle isolates from the corresponding farms with respect to serotype, virulence profile, and pulsed-field gel electrophoresis type. This study shows that wild birds and rodents may become infected from farm animals or vice versa, suggesting a possible role in VTEC transmission.