Insights into scFv:drug binding using the molecular dynamics simulation and free energy calculation

Rheumatoid arthritis autoantibodies benefit from inflammation temperatures

BACKGROUND

Symptoms of rheumatoid arthritis are associated with local inflammation and may include low-grade fever.

OBJECTIVE

To establish the role of fever/inflammation temperatures (38℃-39℃) on the activity of autoantibodies and therapeutic antibodies relevant for rheumatoid arthritis.

METHODS

Through the use of molecular dynamics and free energy calculations, we investigated the role of temperature on the formation of pertinent immune complexes.

RESULTS

Rheumatoid arthritis autoantibodies bind with higher affinity at febrile/inflammation temperatures.

CONCLUSIONS

Fever may modulate binding affinity of autoantibodies relevant for rheumatoid arthritis.

Designing antibody against highly conserved region of dengue envelope protein by in silico screening of scFv mutant library

Dengue being one of the deadliest diseases of tropical regions, enforces to put continuous efforts for the development of vaccine and effective therapeutics. Most of the antibodies generated during dengue infection are non-neutralizing and cause antibody dependent enhancement. Hence, making a potent neutralizing antibody against all four dengue serotypes could be very effective for the treatment. However, designing a single antibody for all serotypes is difficult due to variation in protein sequences. Therefore, the objective is to identify conserved region of dengue envelope protein and then develop an antibody against that conserved region. Before advancing to the development of such an antibody, it is desirable to validate the interactions between antibody and dengue envelope protein. In silico analysis of such interactions provides a good platform to find out a suitable region to design and construct an antibody against it by analyzing antigen-antibody interaction before synthesizing the antibody. In this study, two highly conserved regions of dengue envelope protein were identified and an scFv was constructed against it. Both scFv and FuBc proteins were expressed in bacterial expression system and binding efficiency was analyzed by SPR analysis with KD value 2.3 μM. In order to improve binding efficiency, an in silico scFv mutant library was created which was virtually screened for higher binding efficiency. Six mutants with high binding efficiency were selected for further analysis. The binding ability of these mutants were predicted using simulation analysis which shows these mutations were stabilizing scFv-FuBc complex.

Thinking outside the Laboratory: Analyses of Antibody Structure and Dynamics within Different Solvent Environments in Molecular Dynamics (MD) Simulations

Monoclonal antibodies (mAbs) have revolutionized the biomedical field, directly influencing therapeutics and diagnostics in the biopharmaceutical industry, while continuing advances in computational efficiency have enabled molecular dynamics (MD) simulations to provide atomistic insight into the structure and function of mAbs. Despite the success of MD tools, further optimizations are still required to enhance the computational efficiency of complex mAb simulations. This issue can be tackled by changing the way the solvent system is modelled to reduce the number of atoms to be tracked but must be done without compromising the accuracy of the simulations. In this work, the structure of the IgG_2a antibody was analyzed in three solvent systems: explicit water and ions, implicit water and ions, and implicit water and explicit ions. Root-mean-square distance (RMSD), root-mean-square fluctuations (RMSF), and interchain angles were used to quantify structural changes. The explicit system provides the most atomistic detail but is ~6 times slower in its exploration of configurational space and required ~4 times more computational time on our supercomputer than the implicit simulations. Overall, the behavior of the implicit and explicit simulations is quantifiably similar, with the inclusion of explicit ions in the implicit simulation stabilizing the antibody to reproduce well the statistical fluctuations of the fully explicit system. Therefore, this approach holds promise to maximize the use of computational resources to explore antibody behavior.

The conformational feasibility for the formation of reaching dimer in ASV and HIV integrase: a molecular dynamics study

Retroviral integrases are reported to form alternate dimer assemblies like the core-core dimer and reaching dimer. The core-core dimer is stabilized predominantly by an extensive interface between two catalytic core domains. The reaching dimer is stabilized by N-terminal domains that reach to form intermolecular interfaces with the other subunit's core and C-terminal domains (CTD), as well as CTD-CTD interactions. In this study, molecular dynamics (MD), Brownian dynamics (BD) simulations, and free energy analyses, were performed to elucidate determinants for the stability of the reaching dimer forms of full-length Avian Sarcoma Virus (ASV) and Human Immunodeficiency Virus (HIV) IN, and to examine the role of the C-tails (the last ~16-18 residues at the C-termini) in their structural dynamics. The dynamics of an HIV reaching dimer derived from small angle X-ray scattering and protein crosslinking data, was compared with the dynamics of a core-core dimer model derived from combining the crystal structures of two-domain fragments. The results showed that the core domains in the ASV reaching dimer express free dynamics, whereas those in the HIV reaching dimer are highly stable. BD simulations suggest a higher rate of association for the HIV core-core dimer than the reaching dimer. The predicted stability of these dimers was therefore ranked in the following order: ASV reaching dimer < HIV reaching dimer < composite core-core dimer. Analyses of MD trajectories have suggested residues that are critical for intermolecular contacts in each reaching dimer. Tests of these predictions and insights gained from these analyses could reveal a potential pathway for the association and dissociation of full-length IN multimers.

Production of a soluble single-chain variable fragment antibody against okadaic acid and exploration of its specific binding

Okadaic acid is a lipophilic marine algal toxin commonly responsible for diarrhetic shellfish poisoning (DSP). Outbreaks of DSP have been increasing and are of worldwide public health concern; therefore, there is a growing demand for more rapid, reliable, and economical analytical methods for the detection of this toxin. In this study, anti-okadaic acid single-chain variable fragment (scFv) genes were prepared by cloning heavy and light chain genes from hybridoma cells, followed by fusion of the chains via a linker peptide. An scFv-pLIP6/GN recombinant plasmid was constructed and transformed into Escherichia coli for expression, and the target scFv was identified with IC-CLEIA (chemiluminescent enzyme immunoassay). The IC15 was 0.012 ± 0.02 μg/L, and the IC50 was 0.25 ± 0.03 μg/L. The three-dimensional structure of the scFv was simulated with computer modeling, and okadaic acid was docked to the scFv model to obtain a putative structure of the binding complex. Two predicted critical amino acids, Ser32 and Thr187, were then mutated to verify this theoretical model. Both mutants exhibited significant loss of binding activity. These results help us to understand this specific scFv-antigen binding mechanism and provide guidance for affinity maturation of the antibody in vitro. The high-affinity scFv developed here also has potential for okadaic acid toxin detection.

Insight into conformational change for 14-3-3σ protein by molecular dynamics simulation

14-3-3σ is a member of a highly conserved family of 14-3-3 proteins that has a double-edged sword role in human cancers. Former reports have indicated that the 14-3-3 protein may be in an open or closed state. In this work, we found that the apo-14-3-3σ is in an open state compared with the phosphopeptide bound 14-3-3σ complex which is in a more closed state based on our 80 ns molecular dynamics (MD) simulations. The interaction between the two monomers of 14-3-3σ in the open state is the same as that in the closed state. In both open and closed states, helices A to D, which are involved in dimerization, are stable. However, large differences are found in helices E and F. The hydrophobic contacts and hydrogen bonds between helices E and G in apo-14-3-3σ are different from those in the bound 14-3-3σ complex. The restrained and the mutated (Arg56 or Arg129 to alanine) MD simulations indicate that the conformation of four residues (Lys49, Arg56, Arg129 and Tyr130) may play an important role to keep the 14-3-3σ protein in an open or closed state. These results would be useful to evaluate the 14-3-3σ protein structure-function relationship.

Molecular dynamics simulation of the crystallizable fragment of IgG1-insights for the design of Fcabs

An interesting format in the development of therapeutic monoclonal antibodies uses the crystallizable fragment of IgG1 as starting scaffold. Engineering of its structural loops allows generation of an antigen binding site. However, this might impair the molecule’s conformational stability, which can be overcome by introducing stabilizing point mutations in the CH3 domains. These point mutations often affect the stability and unfolding behavior of both the CH2 and CH3 domains. In order to understand this cross-talk, molecular dynamics simulations of the domains of the Fc fragment of human IgG1 are reported. The structure of human IgG1-Fc obtained from X-ray crystallography is used as a starting point for simulations of the wild-type protein at two different pH values. The stabilizing effect of a single point mutation in the CH3 domain as well as the impact of the hinge region and the glycan tree structure connected to the CH2 domains is investigated. Regions of high local flexibility were identified as potential sites for engineering antigen binding sites. Obtained data are discussed with respect to the available X-ray structure of IgG1-Fc, directed evolution approaches that screen for stability and use of the scaffold IgG1-Fc in the design of antigen binding Fc proteins.

Structural characterization of a therapeutic anti-methamphetamine antibody fragment: oligomerization and binding of active metabolites

Vaccines and monoclonal antibodies (mAb) for treatment of (+)-methamphetamine (METH) abuse are in late stage preclinical and early clinical trial phases, respectively. These immunotherapies work as pharmacokinetic antagonists, sequestering METH and its metabolites away from sites of action in the brain and reduce the rewarding and toxic effects of the drug. A key aspect of these immunotherapy strategies is the understanding of the subtle molecular interactions important for generating antibodies with high affinity and specificity for METH. We previously determined crystal structures of a high affinity anti-METH therapeutic single chain antibody fragment (scFv6H4, K_D = 10 nM) in complex with METH and the (+) stereoisomer of 3,4-methylenedioxymethamphetamine (MDMA, or “ecstasy”). Here we report the crystal structure of scFv6H4 in homo-trimeric unbound (apo) form (2.60Å), as well as monomeric forms in complex with two active metabolites; (+)-amphetamine (AMP, 2.38Å) and (+)-4-hydroxy methamphetamine (p-OH-METH, 2.33Å). The apo structure forms a trimer in the crystal lattice and it results in the formation of an intermolecular composite beta-sheet with a three-fold symmetry. We were also able to structurally characterize the coordination of the His-tags with Ni²⁺. Two of the histidine residues of each C-terminal His-tag interact with Ni²⁺ in an octahedral geometry. In the apo state the CDR loops of scFv6H4 form an open conformation of the binding pocket. Upon ligand binding, the CDR loops adopt a closed formation, encasing the drug almost completely. The structural information reported here elucidates key molecular interactions important in anti-methamphetamine abuse immunotherapy.

Ligand selectivity of estrogen receptors by a molecular dynamics study

Estrogen receptors α (ERα) and β (ERβ) have different physiological functions and expression levels in different tissues. ERα and ERβ are highly homologous and have only two residue substitutions in the binding pocket. This high similarity at the active site stimulates the interests for discovering subtype selective ligands. In this study, molecular dynamics (MD) simulations combined with molecular mechanics generalized Born surface area (MM-GBSA) method have been carried out to analyze the basis of selectivity of three ligands (659, 818 and 041). The calculated binding free energies show that all the ligands bind more tightly to ERβ than to ERα. The dominant free energy components of selectivity for 659 are similar to that for 041, but different from that for 818. The decompositions of free energy contributions and structural analysis imply that there are eight residues primarily contributing to the selectivity for 659, five residues for 041, as well as two residues for 818. The structural analysis implies that two residue substitutions in binding packet cause the position of 659 in ERβ-659 complex to shift relative to that in ERα-659 complex and also cause the conformational changes of other residues in the binding pocket. The higher selectivity for 041 is mainly caused by three residues, Ile373 (Met421), His475 (His524) and Leu476 (Leu525).

Pharmacophore based 3D-QSAR modeling and free energy analysis of VEGFR-2 inhibitors

VEGFR-2, a transmembrane tyrosine kinase receptor is responsible for angiogenesis and has been an attractive target in treating cancers. The inhibition mechanism of structurally diverse urea derivatives, reported as VEGFR-2 inhibitors, was explored by pharmacophore modeling, QSAR, and molecular dynamics based free energy analysis.The pharmacophore hypothesis AADRR, resulted in a highly significant atom based 3D-QSAR model (r(2) = 0.94 and q(2) = 0.84). Binding free energy analysis of the docked complexes of highly active and inactive compounds, after 7 ns MD simulation, revealed the importance of van der Waals interaction in VEGFR-2 inhibition. The decomposition of binding free energy on a per residue basis disclosed that the residues in hinge region and hydrophobic pocket play a role in discriminating the active and inactive inhibitors. Thus, the present study proposes a pharmacophore hypothesis representing the identified interactions pattern and its further application as a template in screening databases to identify novel VEGFR-2 inhibitor scaffolds.