Molecular Dynamics Simulations Reveal Isoform Specific Contact Dynamics between the Plexin Rho GTPase Binding Domain (RBD) and Small Rho GTPases Rac1 and Rnd1.

The Plexin family of transmembrane receptors are unique in that their intracellular region interacts directly with small GTPases of the Rho family. The Rho GTPase binding domain of plexin (RBD)-which is responsible for these interactions-can bind with Rac1 as well as Rnd1 GTPases. GTPase complexes have been crystallized with the RBDs of plexinA1, -A2, and -B1. The protein association is thought to elicit different functional responses in a GTPase and plexin isoform specific manner, but the origin of this is unknown. In this project, we investigated complexes between several RBD and Rac1/Rnd1 GTPases using multimicrosecond length all atom molecular dynamics simulations, also with reference to the free forms of the RBDs and GTPases. In accord with the crystallographic data, the RBDs experience more structural changes than Rho-GTPases upon complex formation. Changes in protein dynamics and networks of correlated motions are revealed by analyzing dihedral angle fluctuations in the proteins. The extent of these changes differs between the different RBDs and also between the Rac1 and Rnd1 GTPases. While the RBDs in the free and bound states have similar-if not decreased-correlations, correlations within the GTPases are increased upon binding. Mapping highly correlated residues to the structures, it is found that the plexinA1, -B1, and -A2 RBDs all have similar communication pathways within the ubiquitin fold, but that different residues are involved. Dynamic network analyses indicate that plexinA1 and -B1 RBDs interact with small GTPases in a similar manner, whereas complexes with the plexinA2 RBD display different features. Importantly complexes with Rnd1 have a considerable number of dynamic correlations and network connections between the proteins, whereas such features are missing in the RBD-Rac1 complexes. Overall, the simulations suggest mechanisms that are consistent with the experimental data on plexinB1 and indicate RBD and GTPase isoform specific changes in protein dynamics upon complex formation.