Impact of AlphaFold on structure prediction of protein complexes: The CASP15-CAPRI experiment

We present the results for CAPRI Round 54, the 5th joint CASP-CAPRI protein assembly prediction challenge. The Round offered 37 targets, including 14 homodimers, 3 homo-trimers, 13 heterodimers including 3 antibody-antigen complexes, and 7 large assemblies. On average ~70 CASP and CAPRI predictor groups, including more than 20 automatics servers, submitted models for each target. A total of 21 941 models submitted by these groups and by 15 CAPRI scorer groups were evaluated using the CAPRI model quality measures and the DockQ score consolidating these measures. The prediction performance was quantified by a weighted score based on the number of models of acceptable quality or higher submitted by each group among their five best models. Results show substantial progress achieved across a significant fraction of the 60+ participating groups. High-quality models were produced for about 40% of the targets compared to 8% two years earlier. This remarkable improvement is due to the wide use of the AlphaFold2 and AlphaFold2-Multimer software and the confidence metrics they provide. Notably, expanded sampling of candidate solutions by manipulating these deep learning inference engines, enriching multiple sequence alignments, or integration of advanced modeling tools, enabled top performing groups to exceed the performance of a standard AlphaFold2-Multimer version used as a yard stick. This notwithstanding, performance remained poor for complexes with antibodies and nanobodies, where evolutionary relationships between the binding partners are lacking, and for complexes featuring conformational flexibility, clearly indicating that the prediction of protein complexes remains a challenging problem.

Prediction of homoprotein and heteroprotein complexes by protein docking and template-based modeling: A CASP-CAPRI experiment

We present the results for CAPRI Round 30, the first joint CASP-CAPRI experiment, which brought together experts from the protein structure prediction and protein-protein docking communities. The Round comprised 25 targets from amongst those submitted for the CASP11 prediction experiment of 2014. The targets included mostly homodimers, a few homotetramers, and two heterodimers, and comprised protein chains that could readily be modeled using templates from the Protein Data Bank. On average 24 CAPRI groups and 7 CASP groups submitted docking predictions for each target, and 12 CAPRI groups per target participated in the CAPRI scoring experiment. In total more than 9500 models were assessed against the 3D structures of the corresponding target complexes. Results show that the prediction of homodimer assemblies by homology modeling techniques and docking calculations is quite successful for targets featuring large enough subunit interfaces to represent stable associations. Targets with ambiguous or inaccurate oligomeric state assignments, often featuring crystal contact-sized interfaces, represented a confounding factor. For those, a much poorer prediction performance was achieved, while nonetheless often providing helpful clues on the correct oligomeric state of the protein. The prediction performance was very poor for genuine tetrameric targets, where the inaccuracy of the homology-built subunit models and the smaller pair-wise interfaces severely limited the ability to derive the correct assembly mode. Our analysis also shows that docking procedures tend to perform better than standard homology modeling techniques and that highly accurate models of the protein components are not always required to identify their association modes with acceptable accuracy. Proteins 2016; 84(Suppl 1):323-348. © 2016 Wiley Periodicals, Inc.

The effect of R249S carcinogenic and H168R-R249S suppressor mutations on p53-DNA interaction, a multi scale computational study

In this study we have undertaken the theoretical analysis of the effect of R249S carcinogenic and H168R-R249S suppressor mutation at core domain of the tumor suppressor protein p53, on its natural interaction with DNA using a newly developed method. The results show that the carcinogenic mutation R249S affects the flexibility of L3 loop region in p53, inducing the loss of important hydrogen bonds observed at interaction in the wild-type with DNA, on the other hand the suppressor mutation H168R on the R249S assists in maintaining the wild-type like flexibility of the L3 region in p53 and thus recover the interaction terms lost in the carcinogenic mutation alone. The present study sets a new direction in the development of new drugs that may restore the interactions that lost as a consequence of the carcinogenic mutations in p53.

Tribochemical reaction dynamics of molybdenum dithiocarbamate on nascent iron surface: a hybrid quantum chemical/classical molecular dynamics study

Using a hybrid quantum chemical/classical molecular dynamics method, we have studied the tribochemical reaction dynamics of molybdenum dithiocarbamate (MoDTC), a commonly used friction modifier in automobile engine oils. MoDTC molecule adsorbed on rubbing nascent iron surface was situated. We firstly investigated the dynamic behavior of MoDTC molecule on the rubbing Fe(001) surface. During the friction simulation, the elongation of Mo-O bonds was observed, forming the Mo2S4 and thiocarbamic acid molecules. To unveil the detailed mechanism of this bond elongation, the electronic states of the MoDTC molecule and Fe(001) surface were computed, and the catalytic effects of Fe(001) surface to the molecule was found. We also found that extreme friction would influence the complete Mo-O bond dissociation. By using the hybrid quantum chemical/classical molecular dynamics method, we successfully simulated the tribochemical reaction dynamics of MoDTC as a friction modifier and obtained the influences of nascent iron surface and friction on its chemical reaction.

Quantum chemistry study on absorption spectra, electronic and electrical properties of organic dye on anatase(001)

As the most reactive surface, the stoichiometric O-bridge terminated anatase(001) surface attracted considerable attentions in many application fields. The interfacial electron transfer in dye-sensitized anatase(001) plays a principal role in a variety of photoinduced reactions. In the present work, the UV-vis absorption spectrum of TiO2 bulk and different surface models were calculated by means of tight-binding quantum chemical molecular dynamics program "Colors-excite" for the first time. The thickness dependence on electronic and electrical properties of anatase(001) surface was achieved. The anatase(001) surface with a thickness of 1.0 nm shows excellent electronic and electrical properties. Moreover, the most suitable binding mode (dissociative adsorption) and absorption spectra of perylene with acrylic acid (PAA) on the optimum anatase(001) were investigated. A significant red-shift was observed from the UV-vis absorption spectrum of PAA/anatase(001) system. The red-shift occurring when PAA adsorbed on anatase(001) surface suggests that PAA/anatase(001) may be potential candidate for dye-sensitized solar cell. This study also proposed an effective computational tool "Colors-excite" to study of the electronic excitation properties for both molecular and periodic systems.

First-principles study on proton dissociation properties of fluorocarbon- and hydrocarbon-based membranes in low humidity conditions

We present a theoretical study on the proton dissociation properties of the membranes for polymer electrolyte fuel cells. A density functional theory method is used to study the influence of fluorocarbon and hydrocarbon backbones on proton dissociation, the interaction of water molecules with the sulfonic acid group, and the energy barriers for proton dissociation. Better proton dissociation properties of CH(3)SO(3)H compared to CF(3)SO(3)H are observed from statistical analyses of the optimized structures for both systems. However, the calculated energy barriers for proton dissociation are lower for CF(3)SO(3)H than for the CH(3)SO(3)H system. At the same time, the interaction of water molecules is stronger for CH(3)SO(3)H than for CF(3)SO(3)H. Also, the analysis of the hydrogen-bonding network in both systems shows that the number of hydrogen bonds formed around the sulfonic acid group in CH(3)SO(3)H is larger than that in CF(3)SO(3)H. Therefore, the decrease of the energy barrier with increasing number of coordinating water molecules, pronounced in the case of CH(3)SO(3)H, may lower the barrier, which enhances good proton conductivity of a hydrocarbon-based polymer in low humidity conditions. Thus the hydration ability of a sulfonic acid group is an important factor for realizing better proton dissociation in low humidity conditions.

Oxidation mechanism in the metabolism of (S)-N-[1-(3-morpholin-4-ylphenyl)ethyl]-3-phenylacrylamide on oxyferryl active site in CYP3A4 Cytochrome: DFT modeling

The metabolism mechanism of (S)-N-[1-(3-morpholin-4ylphenyl)ethyl]-3-phenylacrylamide, mediated by CYP3A4 Cytochrome has been investigated by density functional QM calculations aided with molecular mechanics/molecular dynamics simulations. Two different orientations of phenyl ring for substrate approach toward oxyferryl center, imposing two subsequent rearrangement pathways have been investigated. Starting from sigma-complex in perpendicular orientation enzymatic mechanism involves consecutive proton shuttle intermediate, which further leads to the formation of alcohol and ketone. Parallel conformation leads solely to ketone product by 1,2 hydride shift. Although parallel and perpendicular sigma-complexes are energetically equivalent both for the gas phase or PCM solvent model, molecular dynamics studies in full CYP3A4 environment show that perpendicular conformation of the sigma-complex should be privileged, stabilized by hydrophobic interactions of phenylacrylamide chain. After assessing probability of the two conformations we postulate that the alcohol, accessible with the lowest energy barriers should be the major metabolite for studied substrate and CYP3A4 enzyme.

Combinatorial Computational Chemistry Approach for Materials Design:Applications in deNOx Catalysis, Fischer-Tropsch Synthesis, Lanthanoid Complex, and Lithium Ion Secondary Battery

Computational chemistry can provide fundamental knowledge regarding various aspects of materials. While its impact in scientific research is greatly increasing, its contributions to industrially important issues are far from satisfactory. In order to realize industrial innovation by computational chemistry, a new concept "combinatorial computational chemistry" has been proposed by introducing the concept of combinatorial chemistry to computational chemistry. This combinatorial computational chemistry approach enables theoretical high-throughput screening for materials design. In this manuscript, we review the successful applications of combinatorial computational chemistry to deNO(x) catalysts, Fischer-Tropsch catalysts, lanthanoid complex catalysts, and cathodes of the lithium ion secondary battery.

Identification of amino acid residues in the Ah receptor involved in ligand binding

The Ah receptor (AhR) is a ligand-activated transcription factor. Five amino acids as candidate amino acids necessary for ligand binding within or near the ligand-binding domain were selected based on their evolutional conservation and their aromatic nature that could interact with xenobiotic ligands. These amino acids were changed to Ala, and the mutated AhRs were subjected to a test of their transactivation activity in HeLa cells. Mutation of Phe318 completely lost its activity whereas other mutations only weakly impaired activity. The Leu-substituted mutant, AhR(Phe318Leu), activated the luciferase activity to the level comparable to wild type in the cells treated with 3-methylcholanthrene (MC) but not at all with beta-naphthoflavone (beta-NF). Ligand-binding activity of mutants was examined with [3H]MC in vitro. AhR(Phe318Ala) could not bind to [3H]MC. [3H]MC bound by AhR(Phe318Leu) was competed with unlabeled MC but not with beta-NF. A structural model of the ligand-binding domain was constructed.

Three-dimensional quantitative structure-activity relationship (3 D-QSAR) and docking studies on (benzothiazole-2-yl) acetonitrile derivatives as c-Jun N-terminal kinase-3 (JNK3) inhibitors

Three-dimensional quantitative structure-activity relationship (3D-QSAR) models were developed for 44 (benzothiazole-2-yl) acetonitrile derivatives, inhibiting c-Jun N-terminal kinase-3 (JNK3). It includes molecular field analysis (MFA) and receptor surface analysis (RSA). The QSAR model was developed using 34 compounds and its predictive ability was assessed using a test set of 10 compounds. The predictive 3D-QSAR models have conventional r2 values of 0.849 and 0.766 for MFA and RSA, respectively; while the cross-validated coefficient r(cv)2 values of 0.616 and 0.605 for MFA and RSA, respectively. The results of the QSAR model were further compared with a structure-based analysis using docking studies with crystal structure of JNK3. Ligands bind in the ATP pocket and the hydrogen bond with GLN155 was found to be crucial for selectivity among other kinases. The results of 3D-QSAR and docking studies validate each other and hence, the combination of both methodologies provides a powerful tool directed to the design of novel and selective JNK3 inhibitors.

Tribochemical Reaction Dynamics of Phosphoric Ester Lubricant Additive by Using a Hybrid Tight-Binding Quantum Chemical Molecular Dynamics Method

To study the atomistic behavior of the phosphoric ester molecule on the nascent Fe surface under boundary lubrication conditions, we adopted a hybrid tight-binding quantum chemical molecular dynamics method. First, we investigated chemical interactions between phosphoric ester and the nascent Fe surface. Phosphoric ester was shown to interact with the nascent Fe surface, forming both covalent and ionic bonds. Formation and dissociation dynamics of covalent bonds during tribochemical reaction was clearly observed during the simulation. The effect of friction condition on the tribochemical reaction dynamics was then studied, and it was indicated that friction would influence the formation and the dissociation of covalent bonds. By using a hybrid tight-binding quantum chemical molecular dynamics method, we obtained insights on initial tribochemical reaction processes for the formation of tribofilm from the phosphoric ester molecule on the nascent Fe surface.

Molecular dynamics study on the ligand recognition by tandem SH3 domains of p47phox, regulating NADPH oxidase activity

The phagocyte NADPH oxidase complex plays a crucial role in host defense against microbial infection through the production of superoxides. Chronic granulomatous disease (CGD) is an inherited immune deficiency caused by the absence of certain components of the NADPH oxidase. Key to the activation of the NADPH oxidase is the cytoplasmic subunit p47phox, which includes the tandem SH3 domains (N-SH3 and C-SH3). In active phagocytes, p47phox forms a stable complex with the cytoplasmic region of membrane subunit p22phox that forms a left-handed polyproline type-II (PPII) helix conformation. In this report, we have analyzed the conformational changes of p47phox-p22phox complexes of wild-type and three mutants, which have been detected in CGD patients, using molecular dynamics simulations. We have found that in the wild-type, two basal planes of PPII prism in cytoplasmic region of p22phox interacted with N-SH3 and C-SH3. In contrast, in the modeled mutants, the residue at the ape of PPII helix, which interacts simultaneously with both of the tandem SH3 domains in the wild-type, moved toward C-SH3. Furthermore, interaction energies of the cytoplasmic region of p22phox with C-SH3 tend to decrease in these mutants. All these findings led us to conclude that interactions between N-SH3 of p47phox and PPII helix, which is formed by cytoplasmic region of p22phox, may play a significant role in the activation of the NADPH oxidase.

Robotic path planning and protein complex modeling considering low frequency intra-molecular loop and domain motions

A novel algorithm is introduced to deal with intra-molecular motions of loops and domains that undergo proteins at interaction with other proteins. The methodology is based on complex energy landscape sampling and robotic motion planning. Mapping high flexibility regions on the protein underlies the proposed algorithm. This is the first time this type of research has been reported. Application of the methodology to several protein complexes where remarkable backbone rearrangement is observed shows that the new algorithm is able to deal with the problem of change of backbone conformation at protein interaction. We have implemented the module within the system MIAX (Macromolecular interaction assessment computer system) and together with our already reported soft and flexible docking algorithms we have developed a powerful tool for protein function analysis as part of wide genome function evaluation.

A combined bioinformatic approach oriented to the analysis and design of peptides with high affinity to MHC class I molecules

We report on a new method to compute the antigenic degree of peptides from available experimental data on peptide binding affinity to class I MHC molecules. The methodology is a combination of two strategies at different levels of information. The first, at the primary structure level, consists in expressing the peptides binding activity as a profile of amino acid contributions, amino acid similarity being accounted for by their characteristic physicochemical properties and their position within the sequence. The higher level of the strategy is based on a meticulous analysis of the contact interface of the peptides with the cleft constituting the receptor region of a particular class I MHC molecule. Interaction interfaces are inferred by docking the peptide onto the receptor groove of the MHC molecule; evaluation of the affinity of the peptide to the receptor is then performed by analysis of the electrostatic and hydrophobic energies on points of the interaction interface. The result is a robust system for analysis of peptide affinity to class I MHC molecules since while the first analysis dictates the composition of active sequences at the amino acid level, the second translates this information to the atomic level, where the molecular interaction can be analyzed in terms of the intrinsic interatomic forces and energies. Evaluation results for the methodology are encouraging since high affinity peptides are reflected by high scores at both levels of information, and are proportionally lower for peptides of medium and lower affinity for which interaction surfaces show relatively lower electrostatic complementarity and hydrophobic correlation than for the former.

Detection of Chlamydia trachomatis infection with DNA extracted from formalin-fixed paraffin-embedded tissues

It is a well-known fact that tubal stenosis and/or peritubal adhesion are often associated with Chlamydia trachomatis infection. Although tubal pregnancy may be attributed to this infection, there are only extremely rare cases in which the presence of C. trachomatis has been confirmed by immumo-histochemical examination on tissues isolated from patients with tubal pregnancy. We thus tried to confirm the presence of C. trachomatis infection by detecting DNA of the organism in tissues surgically isolated from patients with tubal pregnancy. Two detection methods, a ligase chain reaction (LCR) method and an immuno-histochemical staining which detects an antigen of C. trachomatis, were compared using paraffin-embedded tissue samples which were surgically isolated from patients with tubal pregnancy or hydrosalpinx. The LCR method was capable of detecting DNA of C. trachomatis in tissue samples in which the C. trachomatis-specific antigen could not be detected using immuno-histochemical staining. This was due to the fact that immuno-histochemical staining methods are applicable to the analysis of sequences the length of which range from 200 to 400 bp (base pairs), while the LCR method used here allows the analysis of sequences as small as 48 bp. This fact makes the LCR method a very convenient tool, as compared with immuno-histochemical methods, for analysis of the paraffin embedded tissue samples where by effects of formalin--used for fixation for pathologic diagnosis--the risk of fragmenting the DNA samples is relatively high. Presence of C. trachomatis DNA as detected by LCR method in surgically isolated samples from patients with tubal pregnancy supports the involvement of Chlamydia trachomatis infection in the occurrence of tubal pregnancy. Accordingly the LCR method is capable of detecting DNA of C. trachomatis in paraffin-embedded samples of tubal tissue in which presence of C. trachomatis could not be confirmed by an immuno-histochemical staining method.