Comparison of the pH- and thermally-induced fluctuations of a therapeutic antibody Fab fragment by molecular dynamics simulation.

Successful development of protein therapeutics depends critically on achieving stability under a range of conditions. A deeper understanding of the drivers of instability across different stress conditions, will enable the engineering of more robust protein scaffolds. We compared the impacts of low pH and high temperature stresses on the structure of a humanized antibody fragment (Fab) A33, using atomistic molecular dynamics simulations, using a recent 2.5 Å crystal structure. This revealed that low-pH induced the loss of native contacts in the domain CL. By contrast, thermal stress led to 5-7% loss of native contacts in all four domains, and simultaneous loss of >30% of native contacts in the VL-VH and CL-CH interfaces. This revealed divergent destabilising pathways under the two different stresses. The underlying cause of instability was probed using FoldX and Rosetta mutation analysis, and packing density calculations. These agreed that mutations in the CL domain, and CL-CH1 interface have the greatest potential for stabilisation of Fab A33. Several key salt bridge losses underpinned the conformational change in CL at low pH, whereas at high temperature, salt bridges became more dynamic, thus contributing to an overall destabilization. Lastly, the unfolding events at the two stress conditions exposed different predicted aggregation-prone regions (APR) to solvent, which would potentially lead to different aggregation mechanisms. Overall, our results identified the early stages of unfolding and stability-limiting regions of Fab A33, and the VH and CL domains as interesting future targets for engineering stability to both pH- and thermal-stresses simultaneously.

An Expanded Conformation of an Antibody Fab Region by X-Ray Scattering, Molecular Dynamics, and smFRET Identifies an Aggregation Mechanism.

Protein aggregation is the underlying cause of many diseases, and also limits the usefulness of many natural and engineered proteins in biotechnology. Better mechanistic understanding and characterization of aggregation-prone states is needed to guide protein engineering, formulation, and drug-targeting strategies that prevent aggregation. While several final aggregated states-notably amyloids-have been characterized structurally, very little is known about the native structural conformers that initiate aggregation. We used a novel combination of small-angle x-ray scattering (SAXS), atomistic molecular dynamic simulations, single-molecule Förster resonance energy transfer, and aggregation-prone region predictions, to characterize structural changes in a native humanized Fab A33 antibody fragment, that correlated with the experimental aggregation kinetics. SAXS revealed increases in the native state radius of gyration, Rg, of 2.2% to 4.1%, at pH 5.5 and below, concomitant with accelerated aggregation. In a cutting-edge approach, we fitted the SAXS data to full MD simulations from the same conditions and located the conformational changes in the native state to the constant domain of the light chain (CL). This CL displacement was independently confirmed using single-molecule Förster resonance energy transfer measurements with two dual-labeled Fabs. These conformational changes were also found to increase the solvent exposure of a predicted APR, suggesting a likely mechanism through which they promote aggregation. Our findings provide a means by which aggregation-prone conformational states can be readily determined experimentally, and thus potentially used to guide protein engineering, or ligand binding strategies, with the aim of stabilizing the protein against aggregation.

Probing the early stages of prion protein (PrP) aggregation with atomistic molecular dynamics simulations.

Prions are self-replicating infectious proteinaceous agents whose conformations are capable of forming amyloid-like aggregate fibrils. Here we present molecular dynamics simulations aimed at investigating the aggregation process of the ß-rich H2H3 domain of the ovine prion protein (H2H3-OvPrPSc), known to be the portion of prion protein carrying oligomerization activity.

Computational Design To Reduce Conformational Flexibility and Aggregation Rates of an Antibody Fab Fragment.

Computationally guided semirational design has significant potential for improving the aggregation kinetics of protein biopharmaceuticals. While improvement in the global conformational stability can stabilize proteins to aggregation under some conditions, previous studies suggest that such an approach is limited, because thermal transition temperatures ( Tm) and the fraction of protein unfolded ( fT) tend to only correlate with aggregation kinetics where the protein is incubated at temperatures approaching the Tm. This is because under these conditions, aggregation from globally unfolded protein becomes dominant. However, under native conditions, the aggregation kinetics are presumed to be dependent on local structural fluctuations or partial unfolding of the native state, which reveal regions of high propensity to form protein-protein interactions that lead to aggregation. In this work, we have targeted the design of stabilizing mutations to regions of the A33 Fab surface structure, which were predicted to be more flexible. This Fab already has high global stability, and global unfolding is not the main cause of aggregation under most conditions. Therefore, the aim was to reduce the conformational flexibility and entropy of the native protein at various locations and thus identify which of those regions has the greatest influence on the aggregation kinetics. Highly dynamic regions of structure were identified through both molecular dynamics simulation and B-factor analysis of related X-ray crystal structures. The most flexible residues were mutated into more stable variants, as predicted by Rosetta, which evaluates the ΔΔ GND for each potential point mutation. Additional destabilizing variants were prepared as controls to evaluate the prediction accuracy and also to assess the general influence of conformational stability on aggregation kinetics. The thermal conformational stability, and aggregation rates of 18 variants at 65 °C, were each examined at pH 4, 200 mM ionic strength, under which conditions the initial wild-type protein was <5% unfolded. Variants with decreased Tm values led to more rapid aggregation due to an increase in the fraction of protein unfolded under the conditions studied. As expected, no significant improvements were observed in the global conformational stability as measured by Tm. However, 6 of the 12 stable variants led to an increase in the cooperativity of unfolding, consistent with lower conformational flexibility and entropy in the native ensemble. Three of these had 5-11% lower aggregation rates, and their structural clustering indicated that the local dynamics of the C-terminus of the heavy chain had a role in influencing the aggregation rate.

Mapping the Aggregation Kinetics of a Therapeutic Antibody Fragment.

The analytical characterization of biopharmaceuticals is a fundamental step in the early stages of development and prediction of their behavior in bioprocesses. Protein aggregation in particular is a common issue as it affects all stages of product development. In the present work, we investigate the stability and the aggregation kinetics of A33Fab, a therapeutically relevant humanized antibody fragment at a wide range of pH, ionic strength, and temperature. We show that the propensity of A33Fab to aggregate under thermally accelerated conditions is pH and ionic-strength dependent with a stronger destabilizing effect of ionic strength at low pH. In the absence of added salts, A33Fab molecules appear to be protected from aggregation due to electrostatic colloidal repulsion at low pH. Analysis by transmission electron microscopy identified significantly different aggregate species formed at low and high pH. The correlations between apparent midpoints of thermal transitions (Tm,app values), or unfolded mole fractions, and aggregation rates are reported here to be significant only at the elevated incubation temperature of 65 °C, where aggregation from the unfolded state predominates. At all other conditions, particularly at 4-45 °C, aggregation of A33 Fab was predominantly from a native-like state, and the kinetics obeyed Arrhenius behavior. Despite this, the rank order of aggregation rates observed at 45 °C, 23 and 4 °C still did not correlate well to each other, indicating that forced degradation at elevated temperatures was not a good screen for predicting behavior at low temperature.

Decrypting Prion Protein Conversion into a ß-Rich Conformer by Molecular Dynamics.

Prion diseases are fatal neurodegenerative diseases characterized by the formation of ß-rich oligomers and the accumulation of amyloid fibrillar deposits in the central nervous system. Understanding the conversion of the cellular prion protein into its ß-rich polymeric conformers is fundamental to tackling the early stages of the development of prion diseases. In this paper, we have identified unfolding and refolding steps critical to the conversion into a ß-rich conformer for different constructs of the ovine prion protein by molecular dynamics simulations. By combining our results with in vitro experiments, we show that the folded C-terminus of the ovine prion protein is able to recurrently undergo a drastic conformational change by displacement of the H1 helix, uncovering of the H2H3 domain, and formation of persistent ß-sheets between H2 and H3 residues. The observed ß-sheets refold toward the C-terminus exposing what we call a "bending region" comprising residues 204-214. This is strikingly coincident with the region harboring mutations determining the fate of the prion oligomerization process. The ß-rich intermediate is used here for the construction of a putative model for the assembly into an oligomeric aggregate. The results presented here confirm the importance of the H2H3 domain for prion oligomer formation and therefore its potential use as molecular target in the design of novel prion inhibitors.

Protein-water interactions in MD simulations: POPS/POPSCOMP solvent accessibility analysis, solvation forces and hydration sites.

The effects of solvation on molecular recognition are investigated from different perspectives, ranging from methods to analyse explicit solvent dynamical behaviour at the protein surface to methods for the implicit treatment of solvent effects associated with the conformational behaviour of biomolecules. The here presented implicit solvation method is based on an analytical approximation of the Solvent Accessible Surface Area (SASA) of solute molecules, which is computationally efficient and easy to parametrise. The parametrised SASA solvation method is discussed in the light of protein design and ligand binding studies. The POPS program for the SASA computation on single molecules and complex interfaces is described in detail. Explicit solvent behaviour is described here in the form of solvent density maps at the protein surface. We highlight the usefulness of that approach in defining the organisation of specific water molecules at functional sites and in determining hydrophobicity scores for the identification of potential interaction patches.

The oligomerization properties of prion protein are restricted to the H2H3 domain.

The propensity of the prion protein (PrP) to adopt different structures is a clue to its pathological behavior. The determination of the region involved in the PrP(C) to PrP(Sc) conversion is fundamental for the understanding of the mechanisms underlying this process at the molecular level. In this paper, the polymerization of the helical H2H3 domain of ovine PrP (OvPrP) was compared to the full-length construct (using chromatography and light scattering). We show that the oligomerization patterns are identical, although the H2H3 domain has a higher polymerization rate. Furthermore, the depolymerization kinetics of purified H2H3 oligomers compared to those purified from the full-length PrP reveal that regions outside H2H3 do not significantly contribute to the oligomerization process. By combining rational mutagenesis and molecular dynamics to investigate the early stages of H2H3 oligomerization, we observe a conformationally stable beta-sheet structure that we propose as a possible nucleus for oligomerization; we also show that single point mutations in H2 and H3 present structural polymorphisms and oligomerization properties that could constitute the basis of species or strain variability.

Genome sequence of Vibrio splendidus: an abundant planctonic marine species with a large genotypic diversity.

Vibrio splendidus is a dominant Vibrio species in seawater presenting a remarkable genetic diversity; several strains have been linked to invertebrate's mortality. We report the complete genome sequence of V. splendidus LGP32, an oyster pathogen, and its comparison with partial genome sequences from related strains. As is typical for the genus, V. splendidus LGP32 contains two chromosomes (3.29 and 1.67 Mb) and most essential cellular processes are encoded by chromosome 1. Comparison with two other V. splendidus partial genome sequences (strains 12B01 and Med222) confirms the previously suggested high genotypic diversity within this species and led to the identification of numerous strain-specific regions that could frequently not be assigned to a specific mechanisms of recombination. Surprisingly, the chromosomal integron, the most variable genetic element in all other Vibrio species analysed to date, is absent from 12B01 and inactivated by a mobile element in Med222, while in LGP32 it only contains a limited number of cassettes. Finally, we found that the LGP32 integron contains a new dfrA cassette, related to those found in resistance integrons of gram-negative clinical isolates. Those results suggest that marine Vibrio can be a source of antibiotic resistance genes.

Regulatory role of UvrY in adaptation of Photorhabdus luminescens growth inside the insect.

We report global expression profiling of a uvrY-deficient mutant of Photorhabdus luminescens. We found that the regulator moiety of the two-component regulatory system BarA/UvrY regulated more than 500 target genes coding for functions involved in the synthesis of major compartments and metabolic pathways of the cell. This regulation appeared to be in part indirect as UvrY affected the expression of several regulators. Indeed, the flagellum biosynthesis transcription activator FlhC and the flagella regulon were induced in the absence of UvrY, leading to a hyperflagellated phenotype and an increase in motility and biofilm formation. Two major regulatory systems were also altered: the type 2 quorum-sensing inducer AI-2 was activated by UvrY, and the CsrA regulator function appeared to be repressed by the increase of the small-untranslated RNA csrB, the CsrA activity inhibitor TldD and the chaperonin GroESL. Both through and independently of these systems, UvrY regulated oxidative stress resistance; bioluminescence; iron, sugar and peptide transport; proteases; polyketide synthesis enzymes and nucleobases recycling, related to insect degradation and assimilation by bacteria. As a consequence, the uvrY-deficient strain exhibited a decreased killing of insect cells and a reduced growth on insect cells culture, suggesting a UvrY role in the adaptation of P. luminescens inside the insect.

Pleiotropic role of quorum-sensing autoinducer 2 in Photorhabdus luminescens.

Bacterial virulence is an integrative process that may involve quorum sensing. In this work, we compared by global expression profiling the wild-type entomopathogenic Photorhabdus luminescens subsp. laumondii TT01 to a luxS-deficient mutant unable to synthesize the type 2 quorum-sensing inducer AI-2. AI-2 was shown to regulate more than 300 targets involved in most compartments and metabolic pathways of the cell. AI-2 is located high in the hierarchy, as it controls the expression of several transcriptional regulators. The regulatory effect of AI-2 appeared to be dose dependent. The luxS-deficient strain exhibited decreased biofilm formation and increased type IV/V pilus-dependent twitching motility. AI-2 activated its own synthesis and transport. It also modulated bioluminescence by regulating the synthesis of spermidine. AI-2 was further shown to increase oxidative stress resistance, which is necessary to overcome part of the innate immune response of the host insect involving reactive oxygen species. Finally, we showed that the luxS-deficient strain had attenuated virulence against the lepidopteran Spodoptera littoralis. We concluded that AI-2 is involved mainly in early steps of insect invasion in P. luminescens.