Integration of Molecular Dynamics Based Predictions into the Optimization of De Novo Protein Designs: Limitations and Benefits.

Recent advances in de novo protein design have gained considerable insight from the intrinsic dynamics of proteins, based on the integration of molecular dynamics simulations protocols on the state-of-the-art de novo protein design protocols used nowadays. With this protocol we illustrate how to set up and run a molecular dynamics simulation followed by a functional protein dynamics analysis. New users will be introduced to some useful open-source computational tools, including the GROMACS molecular dynamics simulation software package and ProDy for protein structural dynamics analysis.

Mild and cost-effective green fluorescent protein purification employing small synthetic ligands.

The green fluorescent protein (GFP) is a useful indicator in a broad range of applications including cell biology, gene expression and biosensing. However, its full potential is hampered by the lack of a selective, mild and low-cost purification scheme. In order to address this demand, a novel adsorbent was developed as a generic platform for the purification of GFP or GFP fusion proteins, giving GFP a dual function as reporter and purification tag. After screening a solid-phase combinatorial library of small synthetic ligands based on the Ugi-reaction, the lead ligand (A4C7) selectively recovered GFP with 94% yield and 94% purity under mild conditions and directly from Escherichia coli extracts. Adsorbents containing the ligand A4C7 maintained the selectivity to recover other proteins fused to GFP. The performance of A4C7 adsorbents was compared with two commercially available methods (immunoprecipitation and hydrophobic interaction chromatography), confirming the new adsorbent as a low-cost viable alternative for GFP purification.

Comparison of the Internal Dynamics of Metalloproteases Provides New Insights on Their Function and Evolution.

Metalloproteases have evolved in a vast number of biological systems, being one of the most diverse types of proteases and presenting a wide range of folds and catalytic metal ions. Given the increasing understanding of protein internal dynamics and its role in enzyme function, we are interested in assessing how the structural heterogeneity of metalloproteases translates into their dynamics. Therefore, the dynamical profile of the clan MA type protein thermolysin, derived from an Elastic Network Model of protein structure, was evaluated against those obtained from a set of experimental structures and molecular dynamics simulation trajectories. A close correspondence was obtained between modes derived from the coarse-grained model and the subspace of functionally-relevant motions observed experimentally, the later being shown to be encoded in the internal dynamics of the protein. This prompted the use of dynamics-based comparison methods that employ such coarse-grained models in a representative set of clan members, allowing for its quantitative description in terms of structural and dynamical variability. Although members show structural similarity, they nonetheless present distinct dynamical profiles, with no apparent correlation between structural and dynamical relatedness. However, previously unnoticed dynamical similarity was found between the relevant members Carboxypeptidase Pfu, Leishmanolysin, and Botulinum Neurotoxin Type A, despite sharing no structural similarity. Inspection of the respective alignments shows that dynamical similarity has a functional basis, namely the need for maintaining proper intermolecular interactions with the respective substrates. These results suggest that distinct selective pressure mechanisms act on metalloproteases at structural and dynamical levels through the course of their evolution. This work shows how new insights on metalloprotease function and evolution can be assessed with comparison schemes that incorporate information on protein dynamics. The integration of these newly developed tools, if applied to other protein families, can lead to more accurate and descriptive protein classification systems.

A theoretical and experimental approach toward the development of affinity adsorbents for GFP and GFP-fusion proteins purification.

The green fluorescent protein (GFP) is widely employed to report on a variety of molecular phenomena, but its selective recovery is hampered by the lack of a low-cost and robust purification alternative. This work reports an integrated approach combining rational design and experimental validation toward the optimization of a small fully-synthetic ligand for GFP purification. A total of 56 affinity ligands based on a first-generation lead structure were rationally designed through molecular modeling protocols. The library of ligands was further synthesized by solid-phase combinatorial methods based on the Ugi reaction and screened against Escherichia coli extracts containing GFP. Ligands A4C2, A5C5 and A5C6 emerged as the new lead structures based on the high estimated theoretical affinity constants and the high GFP binding percentages and enrichment factors. The elution of GFP from these adsorbents was further characterized, where the best compromise between mild elution conditions, yield and purity was found for ligands A5C5 and A5C6. These were tested for purifying a model GFP-fusion protein, where ligand A5C5 yielded higher protein recovery and purity. The molecular interactions between the lead ligands and GFP were further assessed by molecular dynamics simulations, showing a wide range of potential hydrophobic and hydrogen-bond interactions.

A tailor-made "tag-receptor" affinity pair for the purification of fusion proteins.

A novel affinity "tag-receptor" pair was developed as a generic platform for the purification of fusion proteins. The hexapeptide RKRKRK was selected as the affinity tag and fused to green fluorescent protein (GFP). The DNA fragments were designed, cloned in Pet-21c expression vector and expressed in E. coli host as soluble protein. A solid-phase combinatorial library based on the Ugi reaction was synthesized: 64 affinity ligands displaying complementary functionalities towards the designed tag. The library was screened by affinity chromatography in a 96-well format for binding to the RKRKRK-tagged GFP protein. Lead ligand A7C1 was selected for the purification of RKRKRK fusion proteins. The affinity pair RKRKRK-tagged GFP with A7C1 emerged as a promising solution (Ka of 2.45×10(5) M(-1) ). The specificity of the ligand towards the tag was observed experimentally and theoretically through automated docking and molecular dynamics simulations.

Structural evaluation of an alternative Protein A biomimetic ligand for antibody purification.

Affinity chromatography is one of the most common techniques employed at the industrial-scale for antibody purification. In particular, the purification of human immunoglobulin G (hIgG) has gained relevance with the immobilization of its natural binding counterpart-Staphylococcus aureus Protein A (SpA) or with the recent development of biomimetic affinity ligands, namely triazine-based ligands. These ligands have been developed in order to overcome economic and leaching issues associated to SpA. The most recent triazine-based ligand-TPN-BM, came up as an analogue of 2-(3-amino-phenol)-6-(4-amino-1-naphthol)-4-chloro-sym-triazine ligand also known as ligand 22/8 with improved physico-chemical properties and a greener synthetic route. This work intends to evaluate the potential of TPN-BM as an alternative affinity ligand towards antibody recognition and binding, namely IgG, at an atomic level, since it has already been tested, after immobilization onto chitosan-based monoliths and demonstrated interesting affinity behaviour for this purpose. Herein, combining automated molecular docking and molecular dynamics simulations it was predicted that TPN-BM has high propensity to bind IgG through the same binding site found in the crystallographic structure of SpA_IgG complex, as well as theoretically predicted for ligand 22/8_IgG complex. Furthermore, it was found that TPN-BM established preferential interactions with aromatic residues at the Fab domain (Trp 50, Tyr 53, Tyr 98 and Trp 100), while in the Fc domain the main interactions are based on hydrogen bonds with pH sensitive residues at operational regime for binding and elution like histidines (His 460, His 464, His 466). Moreover, the pH dependence of TPN-BM_IgG complex formation was more evident for the Fc domain, where at pH 3 the protonation state and consequently the charge alteration of histidine residues located at the IgG binding site induced ligand detachment which explains the optimal elution condition at this pH observed experimentally.

Understanding the molecular recognition between antibody fragments and protein A biomimetic ligand.

Affinity chromatography with protein A from Staphylococcus aureus (SpA) is the most widespread and accepted methodology for antibody capture during the downstream process of antibody manufacturing. A triazine based ligand (ligand 22/8) was previously developed as an inexpensive and robust alternative to SpA chromatography (Li et al. and Teng et al.). Despite the experimental success, there is no structural information on the binding modes of ligand 22/8 to antibodies, namely to Immunoglobulin G (IgG) molecules and fragments. In this work, we addressed this issue by a molecular docking approach allied to molecular dynamics simulations. Theoretical results confirmed the preference of the synthetic ligand to bind IgG through the binding site found in the crystallographic structure of the natural complex between SpA and the Fc fragment of IgG. Our studies also suggested other unknown "hot-spots" for specific binding of the affinity ligand at the hinge between V(H) and C(H)1 domains of Fab fragment. The best docking poses were further analysed by molecular dynamics studies at three different protonation states (pH 3, 7 and 11). The main interactions between ligand 22/8 and the IgG fragments found at pH 7 were weaker at pH 3 and pH 11 and in these conditions the ligand start losing tight contact with the binding site, corroborating the experimental evidence for protein elution from the chromatographic adsorbents at these pH conditions.

Molecular mechanism of the hydration of Candida antarctica lipase B in the gas phase: Water adsorption isotherms and molecular dynamics simulations.

Hydration is a major determinant of activity and selectivity of enzymes in organic solvents or in gas phase. The molecular mechanism of the hydration of Candida antarctica lipase B (CALB) and its dependence on the thermodynamic activity of water (a(w)) was studied by molecular dynamics simulations and compared to experimentally determined water sorption isotherms. Hydration occurred in two phases. At low water activity, single water molecules bound to specific water binding sites at the protein surface. As the water activity increased, water networks gradually developed. The number of protein-bound water molecules increased linearly with a(w), until at a(w)=0.5 a spanning water network was formed consisting of 311 water molecules, which covered the hydrophilic surface of CALB, with the exception of the hydrophobic substrate-binding site. At higher water activity, the thickness of the hydration shell increased up to 10 A close to a(w)=1. Above a limit of 1600 protein-bound water molecules the hydration shell becomes unstable and the formation of pure water droplets occurs in these oversaturated simulation conditions. While the structure and the overall flexibility of CALB was independent of the hydration state, the flexibility of individual loops was sensitive to hydration: some loops, such as those part of the substrate-binding site, became more flexible, while other parts of the protein became more rigid upon hydration. However, the molecular mechanism of how flexibility is related to activity and selectivity is still elusive.

Anchoring effects in a wide binding pocket: the molecular basis of regioselectivity in engineered cytochrome P450 monooxygenase from B. megaterium.

The molecular basis of regioselectivity of cytochrome P450 monooxygenases from Bacillus megaterium (CYP102A1) with its flexible and widely opened active site is still not well understood. In the present work (-)-alpha-pinene bound complexes with two triple mutants were modeled to elucidate the contribution of the three major factors that mediate selectivity: active site shape, protein flexibility, and chemical reactivity of the substrate. For the triple mutant A74G F87V L188Q (GVQ), one stable, productive conformation of the substrate (conformation I) was identified by multiple molecular dynamics simulations. The model predicts pinene epoxide as a major product (42% pinene oxide, 23% verbenol) which is in agreement with the experimental product profile (70% pinene oxide, 20% verbenol). In contrast, for the triple mutant A74G F87G L188Q (GGQ) two stable productive substrate conformations were identified (conformations IIa and IIb), and verbenol was predicted as major product (81% verbenol, 16% myrtenol), which is in agreement with experimental results (77% verbenol, 10% myrtenol). The effect of chemical reactivity of the substrate was demonstrated by comparison of (-)-alpha-pinene to its regioisomer (-)-beta-pinene, where the product profile is shifted from 68% pinocarveol and 32% myrtanal in mutant GVQ, to 40% pinocarveol and 60% myrtanal in mutant GGQ. Our results strongly suggest a major role of residue 87 in anchoring (-)-alpha-pinene during substrate binding which provides a simple and elegant rationalization of the dynamic structure of this enzyme-substrate complex.

Molecular dynamics simulations of the enzyme Cu, Zn superoxide dismutase.

The enzyme Cu, Zn superoxide dismutase (Cu,Zn-SOD) is a ubiquitous oxireductase, which is responsible for the cellular defense against oxidative stress caused by the high toxicity of the superoxide radical, and has been also linked to some cases of familiar amyotrophic lateral sclerosis. In the present study a set of molecular mechanics parameters for the active site of Cu,Zn-SOD has been derived. Afterward, an extensive molecular dynamics simulation has been carried out in an aqueous environment. The obtained results shed a further light on the structural flexibility of the backbone, where the active site is nested, and the solvation shell occupancy. The relatively small backbone deviation, shown by a root-mean-square deviation below 1.0 A, confirms the accuracy of the parameters. The solvent shell analysis has shown that the first solvation shell is located at about 5 A from the copper ion, generating an empty cavity with enough space to accommodate the superoxide radical. The low residence time means that a high permutation rate of water molecules in both solvation shells is consistent with the efficiency of this catalytic mechanism. Hybrid studies using ONIOM methodologies can now be done to evaluate the mechanistic implications of the explicit inclusion of the whole system.